ÄÚÈÝ·¢²¼¸üÐÂʱ¼ä : 2025/5/2 14:56:44ÐÇÆÚÒ» ÏÂÃæÊÇÎÄÕµÄÈ«²¿ÄÚÈÝÇëÈÏÕæÔĶÁ¡£
1£©rlocfind():¼ÆËã¸ø¶¨¸ùµÄ¸ù¹ì¼£ÔöÒæ 2£©rlocus() º¯Êý£º¹¦ÄÜΪÇóϵͳ¸ù¹ì¼£¡£ 4¡¢¿ØÖÆÏµÍ³ÆµÓò·ÖÎö»ù±¾·½·¨ 1£©Bodeͼ£ºbode()º¯Êý
2£©Nyquistͼ£ºnyquist()º¯Êý 3£©Îȶ¨Ô£¶È¼ÆË㣺margin( )º¯Êý 5¡¢ÏßÐÔϵͳʱ¼äÏìÓ¦·ÖÎö
1£©step( )º¯Êý---Çóϵͳ½×Ô¾ÏìÓ¦
2£©impulse( )º¯Êý£ºÇóȡϵͳµÄÂö³åÏìÓ¦ 3£©lsim( )º¯Êý£ºÇóϵͳµÄÈÎÒâÊäÈëϵķÂÕæ
Èý¡¢ÊµÑé³ÌÐò¼°½á¹û
1¡¢±íʾÏÂÁд«µÝº¯ÊýÄ£ÐÍ,²¢×ª»¯³ÉÆäËûµÄÊýѧģÐÍ
1.(1)´«µÝº¯ÊýÒÔ¼°´«µÝº¯Êýת»¯ÎªÁ㼫µãÐÎʽ
´«µÝº¯Êýת»¯Îª×´¿Õ±í´ïʽ
£¨2£©´«µÝº¯ÊýÒÔ¼°´«µÝº¯Êýת»¯ÎªÁ㼫µãÐÎʽ ´«µÝº¯Êýת»¯Îª×´¿Õ±í´ïʽ
£¨3£©´«µÝº¯ÊýÒÔ¼°´«µÝº¯Êýת»¯ÎªÁ㼫µãÐÎʽ
´«µÝº¯Êýת»¯Îª×´¿Õ±í´ïʽ
(4)×´¿Õ±í´ïʽÒÔ¼°×´¿Õ±í´ïʽת»¯Îª´«µÝº¯Êý
×´¿Õ±í´ïʽת»¯ÎªÁ㼫µãÐÎʽ
2¡¢Ò»¸öµ¥Î»¸º·´À¡¿ª»·´«µÝº¯ÊýΪ
ÊÔ»æ³öϵͳ±Õ»·µÄ¸ù¹ì¼£Í¼£»²¢ÔÚ¸ú¹ì¼£Í¼ÉÏÈÎѡһµã£¬ÊÔ¼ÆËã¸ÃµãµÄÔöÒæK¼°ÆäËùÓм«µãµÄλÖá£
3¡¢ÇóÏÂÃæÏµÍ³ÔÚ½×Ô¾ÐźÅΪ(t)ʱϵͳµÄÏìÓ¦¡£
²¢ÇóϵͳÐÔÄÜÖ¸±ê£ºÎÈֵ̬¡¢ÉÏÉýʱ¼ä¡¢µ÷½Úʱ¼ä¡¢³¬µ÷Á¿¡£
ÎÈֵ̬ ¡¢ÉÏÉýʱ¼ä ¡¢µ÷½Úʱ¼ä 10 ¡¢³¬µ÷Á¿
ËÄ ÊµÑé×ܽá
ͨ¹ý±¾´ÎʵÑéÕÆÎÕ¿ØÖÆÏµÍ³ÊýѧģÐ͵Ļù±¾ÃèÊö·½·¨ºÍÏ໥ת»¯£¬ÁË½â¿ØÖÆÏµÍ³µÄÎȶ¨ÐÔ·ÖÎö·½·¨£¬ÕÆÎÕ¿ØÖÆÏµÍ³ÆµÓòÓëʱÓò·ÖÎö»ù±¾·½·¨¡£
ʵÑéÆß ¿ØÖÆÏµÍ³PIDУÕýÆ÷Éè¼Æ·¨
Ò»¡¢ÊµÑéÄ¿µÄ
1¡¢ÊìϤ³£¹æPID¿ØÖÆÆ÷µÄÉè¼Æ·½·¨
2¡¢ÕÆÎÕPID²ÎÊýµÄµ÷½Ú¹æÂÉ
3¡¢Ñ§Ï°±àд³ÌÐòÇóϵͳµÄ¶¯Ì¬ÐÔÄÜÖ¸±ê
¶þ¡¢ÊµÑé»ù±¾ÔÀí
1£®Ä£ÄâPID¿ØÖÆÆ÷
µäÐ͵ÄPID¿ØÖƽṹÈçͼ1Ëùʾ¡£ PID¿ØÖÆÆ÷ ±È Àý e(t) r(t) u(t) ` ¶ÔÏóÄ£ÐÍ »ý ·Ö
΢ ·Ö
ͼ1 µäÐÍPID¿ØÖƽṹ
PIDµ÷½ÚÆ÷µÄÊýѧÃèÊöΪ
y(t) u(t)?Kp[e(t)?1Ti?t0e(?)d??Tdde(t)] dt2 Êý×ÖPID¿ØÖÆÆ÷
ÔÚ¼ÆËã»úPID¿ØÖÆÖУ¬Á¬ÐøPID¿ØÖÆËã·¨²»ÄÜÖ±½ÓʹÓã¬ÐèÒª²ÉÓÃÀëÉ¢»¯·½·¨£¬Í¨³£Ê¹ÓÃÊý×ÖPID¿ØÖÆÆ÷¡£ÒÔһϵÁвÉÑùʱ¿ÌµãkT£¨TΪ²ÉÑùÖÜÆÚ£©´ú±íÁ¬ÐøÊ±¼ät£¬ÒÔ¾ØÐη¨ÊýÖµ»ý·Ö½üËÆ´úÌæ»ý·Ö£¬ÒÔÒ»½×ºóÏò²î·Ö½üËÆ´úÌæÎ¢·Ö£¬¼´£º
??t?kT?kkt? e(jT)?T?e(j)??0e(?)d??T?j?0j?0??de(t)e(kT)?e((k?1)T)e(k)?e(k?1)???TT?dtÀëÉ¢PID±í´ïʽ£º
1u(k)?Kp[e(k)?Ti?e(j)T?Tdj?0ke(k)?e(k?1)]
TÈý¡¢ÊµÑé³ÌÐò¼°½á¹û
1¡¢ÔÚSIMULINK´°¿Ú½¨Á¢ÈçÏÂÄ£ÐÍ
2¡¢Éè¼ÆPID¿ØÖÆÆ÷£¬´«µÝº¯ÊýÄ£ÐÍÈçÏÂ
¼ÆËãµÃ³ö M= p= tp= tr= ts= yss= Âú×ãÌõ¼þ¡£
´ËʱPID²ÎÊýÈçÏÂ±í£º
²ÎÊý ÊýÖµ 5 Kp 5 Ti 2 Td ³ÌÐòÈçÏ£º
%system analyse %M maxmun %p over
%tp maxmun-time %ts transition-time %tr rise_time %yss stable
function[M,p,tp,tr,ts,yss]=sys_ana(t,y); tm=length(t); yss=y(tm);
[M tpk]=max(y); tp=t(tpk); p=(M-yss)/yss;
while (y(tm)>*yss)&(y(tm) ts=t(tm); k=1; while y(k)<=yss* k=k+1; end tr=t(k); 6¡¢½«ÊäÈëÐźÅÉèΪr(k)?0.50sin(2?t)£¬µ÷½ÚPID¿ØÖÆÆ÷²ÎÊý£¬»æÖÆÕýÏÒ¸ú×ÙÇúÏß¡£ ʵÑé°Ë ÏßÐÔ·½³Ì×éÇó½â¼°º¯ÊýÇó¼«Öµ Ò»¡¢ÊµÑéÄ¿µÄ 1¡¢Àí½âÏßÐÔ·½³Ì×éÇó½â·½·¨ 2¡¢Àí½âº¯ÊýÇó¼«Öµ·½·¨ 3¡¢Í¨¹ýÁ·Ï°ÒÔÏÂÄÚÈÝÊìϤÇó½âÏßÐÔ·½³Ì×éµÄ·½·¨ºÍº¯ÊýÇó¼«ÖµµÄ·½·¨ ¶þ¡¢ÊµÑé˵Ã÷ 1£®×ÔÖ÷±àд³ÌÐò£¬±ØÒªÊ±²Î¿¼Ïà¹Ø×ÊÁÏ 2£®ÊµÑéǰӦд³ö³ÌÐò´óÖ¿ò¼Ü»òÍêÕûµÄ³ÌÐò´úÂë 5£®ÊµÑéѧʱ£º2ѧʱ Èý¡¢ÊµÑéÄÚÈݼ°ÊµÑé³ÌÐò 1¡¢ÇóÏÂÁз½³Ì×éµÄ½â ?6x1?5x2?2x3?5x4??4?2x?3y?5z?10?9x?x?4x?x?13??12341£© ?3x?7y?4z?3 ? 3x?4x?2x?2x?1234?x?7y?z?5?1???3x1?9x2?2x4?11 2¡¢ÇóÏÂÁк¯ÊýÔÚÖ¸¶¨Çø¼äµÄ×î´óÖµ 1?x2,x?(0,2) 1£©f(x)?41?x>> f='-(1+x^2)/(1+x^4)',x=fminsearch(f,0),ymax=(1+x^2)/(1+x^4) f = -(1+x^2)/(1+x^4) x = ymax = 2£©f(x)?sinx?cosx,x?(0,?) >> f='-sin(x)-cos(x^2)',x=fminsearch(f,0),ymax=sin(x)+cos(x^2) f = -sin(x)-cos(x^2) x = 2