资源描述：
《《控制系统数字仿真与cad》张晓华版课后答案》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、2・2.用MATLAB语言求下列系统的状态方程(1)状态方程模型参数:编写matlab程序如下»num=[l72424];»den=[l10355024];»[ABCD]=tf2ss(num,den)(2)零极点增益:编写程序»num=[l72424];»den=[l10355024];»fZPK]=tf2zp(num,den)得到结果Z=-2.7306+2.8531,-2.7306-2.853li,-1.5388P=-4,-3,・2,-lK=1⑶部分分式形式:编写程序»num=[l72424];»den=[l10355024];»[RPH]=residue
2、(num,den)得到结果R=4.0000<6.0000,2.0000,1.0000P=-4.0000,-3.0000,-2.0000,-1.0000H=[](2)解：(1)传递函数模型参数：编写程序»A=[2.25-5-1.25-0.5;2.25-4.25-1.25-0.25;0.25-0.5-1.25・1;1.25-1.75-0.25-0.75]；»B=[4220];»C=[0202J;»D=[0];»fnumden]=ss2tf(A,B,C,D)得到结果num=04.000014.000022.000015.0000den=1.00004.00006.
3、25005.25002.2500(2)零极点增益模型参数：A=[2.25-5-1.25-0.5:2.25・4.25-1.25・0.25:0.25-0.5-1.251.25・1.75・0.25-0.75]；B=[4220];C=[0202];D=IOJ;[Z,P,K]=ss2zp(A,B,C,D)(3)部分分式形式的模型参数：A=[2.25-5-1.25-0.52.25・4.25-1.25・0.250.25・0.5-1.25・11.25-1.75-0.25-0.75]；B=f4220]：C=[0202);D二[0];[numdenJ=ss2tf(A,B,C,D
4、)[R,P,H]=residue(num,den)2-6(1)解：m文件程序为h=0.1;dispCa数的数值解为');％显示''中间的文字％dispCy=');%同上％y二1;fort=0:h:lm=y;disp(y);%显示y的当前值％y=m-m*h;end保存文件da2.m在matalb命令行屮键入>>da得到结果函数的数值解为y=10.90000.81000.72900.65610.59050.53140.47830.43050.38740.3487(2)另建一个m文件求解y二e在点[0,1]的数值程序为h=0.1;disp(屈数的离散时刻解为)；d
5、isp('yM);fort=0:h:ly=exp(-t);disp(y);end保存文件da3.m在matalb命令行中键入»da3函数的离散时刻解为y=10.90480.81870.74080.67030.60650.54880.49660.44930.40660.36792-5.解：⑴编辑m文件程序h=0.1;dispC四阶龙格■库塔方法求解两数数值解为dispO*);y=i;fort=0:h:ldisp(y);kl=-y;k2=-(y+kl*h/2);k3=-(y+k2*h/2);k4=-(y+k3*h);y=y+(k1+2*k2+2*k3+k4)*h
6、/6;end保存文件q5.m在matlab命令行里键入>>q5得到结果四阶龙格■库塔方法求解函数数值解为y=10.90480.81870.74080.67030.60650.54880.49660.44930.40660.36792-7»a=ll0];»b二[14.6];»c=[l3.416.35];»d=conv(a,b)；»e=conv(d,c)e=1.00008.000031.990075.21000»f=[oo05100];»g=e+fg=1.00008.000031.990080.2100100.0000%以上是计算闭环传递函数的特征多项式以上是计
7、算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式以上是计算闭环传递函数的特征多项式％»p=roots(g)%计算特征多项式的根计算特征多项式的根计算特征多项式的根计算特征多项式的根，，，就是闭环，就是闭环就是闭环就是闭环传递函数的极点传递函数的极点传递函数的极点传递函数的极点％P=-0.9987+3.009li-0.9987-3.009li-3.0013+0.9697i-3.0013-0.9697i»m=[5100];»z=roots(m)z=-20%计算零点计算零点计算零点计算零点％3-2.进行下列计算，给出不使用for和while等循环语句的
8、计算方法。解:根据等比数列求和方法，再利用matla
