Matlab入门4.doc

传递函数形式转换传递函数模型可转化成零极点模型，也可转换成状态空间模型。例一、模型转换 % name:model_transform_01clear,clc % 清除此前运行程序的变量，清屏num1=5*1,-5,6 % 传函分子多项式den1=conv(conv(1,-1,1,-1),1,2) % 传函分母多项式tfGs=tf(num1,den1) % 传函的分式形式zpGs=zpk(tfGs) % 转换成传函的零极点形式z,p,k=zpkdata(zpGs,v) % 列出零极点和比例系数ssGs=ss(tfGs) % 转换成状态空间形式运行结果：num1 = 5 -25 30den1 = 1 0 -3 2 Transfer function:5 s2 - 25 s + 30- s3 - 3 s + 2 Zero/pole/gain:5 (s-3) (s-2)-(s+2) (s-1)2 z = 3.0000 2.0000p = -2.0000 1.0000 1.0000k = 5 a = x1 x2 x3 x1 0 0.75 -0.25 x2 4 0 0 x3 0 2 0 b = u1 x1 4 x2 0 x3 0 c = x1 x2 x3 y1 1.25 -1.563 0.9375 d = u1 y1 0 Continuous-time model.状态方程：输出方程：例二、传递函数模型的串联% series exampleclear,clcnum1=0,1;den1=1,1;num2=0,2;den2=1,2;num3=0,3;den3=1,3;Gs1=tf(num1,den1)Gs2=tf(num2,den2)Gs3=tf(num3,den3)Gs4=series(Gs1,Gs2) % 不允许三连串Gs4=series(Gs1,Gs2,Gs3)Gs5=series(Gs4,Gs3)Gs6=Gs1*Gs2*Gs3 % 允许三连乘运行结果：Transfer function： 1-s + 1 Transfer function: 2-s + 2 Transfer function: 3-s + 3 Transfer function: 2-s2 + 3 s + 2 Transfer function: 6-s3 + 6 s2 + 11 s + 6 Transfer function: 6-s3 + 6 s2 + 11 s + 6例三、传递函数模型的并联% parallel exampleclear,clcnum1=0,1;den1=1,1;num2=0,2;den2=1,2;num3=0,3;den3=1,3;Gs1=tf(num1,den1)Gs2=tf(num2,den2)Gs3=tf(num3,den3)Gs4= parallel (Gs1,Gs2) % 不允许三连并Gs4= parallel (Gs1,Gs2,Gs3)Gs5= parallel (Gs4,Gs3)Gs6=Gs1+Gs2+Gs3 % 允许三连加运行结果：Transfer function: 1-s + 1 Transfer function: 2-s + 2 Transfer function: 3-s + 3 Transfer function: 3 s + 4-s2 + 3 s + 2 Transfer function: 6 s2 + 22 s + 18-s3 + 6 s2 + 11 s + 6 Transfer function: 6 s2 + 22 s + 18-s3 + 6 s2 + 11 s + 6例四、传递函数模型的反馈% feedback_example clear,clcnum1=0,1;den1=1,1;num2=0,2;den2=1,2;num3=0,3;den3=1,3;Gs1=tf(num1,den1)Gs2=tf(num2,den2)Gs3=tf(num3,den3)Gs4= series (Gs1,Gs2) Gs5= feedback (Gs4,Gs3) % 建立反馈闭路的闭环传函Gs6= feedback (Gs1*Gs2,Gs3) % 建立反馈闭路的闭环传函 运行结果：Transfer function: 1-s + 1 Transfer function: 2-s + 2 Transfer function: 3-s + 3 Transfer function: 2-s2 + 3 s + 2 Transfer function: 2 s + 6-s3 + 6 s2 + 11 s + 12 Transfer function: 2 s + 6-s3 + 6 s2 + 11 s + 12例五、传递函数模型的单位反馈% unit_feedback_example clear,clcnum1=0,1;den1=1,1; num2=0,2;den2=1,2;Gs1=tf(num1,den1)Gs2=tf(num2,den2)num3,den3=series(num1,den1,num2,den2);num4,den4=cloop(num3,den3) % 求出单位反馈闭环传函的分子分母 Gs4=tf(num4,den4) % 求出单位反馈闭环传函运