控制系统数字仿真与CAD第4章仿真题答案.doc

《控制系统数字仿真与CAD第4章仿真题答案.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、__________________________________________________《控制系统数字仿真与CAD》第四章上机实验报告院系：电信工程学院班级：11级控制工程基地1班姓名：张亚东学号：11220323时间：2014年5月27日例4-1求阶跃响应的数值解。解：先建立一个名为sp4-1.m的子程序如下：b=b/a(1);a=a/a(1);A=a(2:n+1);A=[rot90(rot90(eye(n-1,n)));-fliplr(A)];B=[zeros(1,n-1),1]';m1=length(b);C=[fliplr(b),z

2、eros(1,n-m1)];Ab=A-B*C*V;X=X0';y=0;t=T0;N=round(Tf-T0)/h;fori=1:NK1=Ab*X+B*R;K2=Ab*(X+h*K1/2)+B*R;K3=Ab*(X+h*K2/2)+B*R;K4=Ab*(X+h*K3)+B*R;X=X+h*(K1+2*K2+2*K3+K4)/6;y=[y,C*X];t=[t,t(i)+h];end[t',y']plot(t,y)_________________________________________________________________________

3、___________________________然后对上面的子程序进行调用，改变开环放大系数k，得到k取不同值时输出响应的曲线如下：>>k=2;>>a=conv([100],conv([0.251],[0.251]))a=0.06250.50001.000000>>b=[2*kk]b=42>>X0=[0000]X0=0000>>V=1;n=4;T0=0;Tf=10;h=0.25;R=1;>>sp4_1ans=000.25000.08850.50000.53930.75001.16501.00001.66901.25001.83321.50001.

4、61991.75001.17042.00000.71822.25000.47032.50000.51432.75000.79233.00001.14533.25001.39913.50001.44683.75001.29004.00001.0251____________________________________________________________________________________________________4.25000.78694.50000.68174.75000.74135.00000.91775.25001

5、.11435.50001.23735.75001.23826.00001.13106.25000.97826.50000.85586.75000.81657.00000.86837.25000.97607.50001.08287.75001.13918.00001.12508.25001.05548.50000.96918.75000.90819.00000.89759.25000.93609.50000.99999.75001.056310.00001.0801>>holdon>>k=1;>>a=conv([100],conv([0.251],[0.

6、251]))a=0.06250.50001.000000>>b=[2*kk]b=21>>X0=[0000]X0=0000____________________________________________________________________________________________________>>V=1;n=4;T0=0;Tf=10;h=0.25;R=1;>>sp4_1ans=000.25000.04430.50000.27790.75000.63921.00001.02371.25001.34101.50001.53611.

7、75001.59442.00001.53422.25001.39502.50001.22352.75001.06193.00000.93993.25000.87153.50000.85633.75000.88264.00000.93334.25000.99054.50001.03974.75001.07185.00001.08375.25001.07765.50001.05895.75001.03456.00001.01066.25000.99196.50000.98096.75000.97787.00000.98117.25000.98837.500

8、00.99687.75001.00438.00001.00958.25001.01188.50