资源描述：
《控制工程基础实验报告20090907》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第一题：三、实验过程和结果：1：>>num=[1,3,5,4];den=[2,4,1,3,7];G1=tf(num,den)Transferfunction:s^3+3s^2+5s+4-----------------------------2s^4+4s^3+s^2+3s+7>>z=[1,4,5];p=[2,3,7,6];k=[2];G2=zpk(z,p,k)Zero/pole/gain:2(s-1)(s-4)(s-5)-----------------------(s-2)(s-3)(s-6)(s-7)>>G
2、3=G1*G2Zero/pole/gain:(s+1.453)(s-1)(s-4)(s-5)(s^2+1.547s+2.752)---------------------------------------------------------------------(s-2)(s-3)(s-6)(s-7)(s^2+3.101s+2.532)(s^2-1.101s+1.383)>>G4=G1+G2Zero/pole/gain:2.5(s+1.337)(s^2-9.187s+21.14)(s^2+1.713s+1.4
3、74)(s^2-3.263s+3.494)---------------------------------------------------------------------(s-2)(s-3)(s-6)(s-7)(s^2+3.101s+2.532)(s^2-1.101s+1.383)>>G5=feedback(G1,G2,-1)Zero/pole/gain:0.5(s+1.453)(s-2)(s-3)(s-6)(s-7)(s^2+1.547s+2.752)-------------------------
4、-----------------------------------------(s-6.544)(s-6.318)(s-2.72)(s-2.38)(s^2+3.059s+2.471)(s^2-1.096s+1.212)>>[z,p,p]=tf2zp([1,3],[3,2,4])z=-3p=0.3333p=0.3333>>[num]=zp2tf(z,p,k)num=262.>>simulink>>[n,d]=linmod('cdhs1')>Indlinmodat202Inlinmodat60n=0003.750
5、01.2500-0.0000d=1.00001.41670.7361-9.8333-4.5694-0.41>>[n,d]=linmod('cdhs2')>Indlinmodat202Inlinmodat60n=Columns1through600000.75000.2500Column7-0.0000d=Columns1through61.00001.61671.0194-9.6861-6.5361-1.3306Column7-0.0833第二题：三、实验过程和结果：>>G1=tf([1],[1,2]);G2=t
6、f([45],[22,3,1]);G3=zpk([1,2],[3,4,5],[1])Zero/pole/gain:(s-1)(s-2)-----------------(s-3)(s-4)(s-5)>>figure(1);subplot(2,2,1);step(G1);subplot(2,2,2);impulse(G1);subplot(2,2,3);G5=G1/tf([1,0],[1]);step(G5);subplot(2,2,4);t=0:0.001:15;u=1+sin(t);lsim(G1,u,t);>
7、>figure(2);subplot(3,2,1);step(G2);holdon;G4=tf([23],[16,2,5]);step(G4);subplot(3,3,3);impulse(G2);subplot(2,2,3);G6=G1/tf([1,0],[1]);step(G6);subplot(2,2,4);t=0:0.001:15;u=1+sin(t);lsim(G2,u,t);>>figure(3);subplot(3,2,1);step(G3);subplot(3,3,3);impulse(G3);s
8、ubplot(2,2,3);G7=G1/tf([1,0],[1]);step(G7);subplot(2,3,5);t=0:0.002:15;u=1+sin(t);lsim(G3,u,t);>>n=[45];d=[22,3,1];Gb=tf(n,d);sys1=1+Gb;sys2=tf(sys1.den,sys1.num);n3=[1,0];d3=[1];sys3=tf(
