《控制系统数字仿真与CAD习题.docx》由会员分享，可在线阅读，更多相关《控制系统数字仿真与CAD习题.docx（34页珍藏版）》请在金锄头文库上搜索。
1、控制系统计算机仿真及辅助设计课程上机作业1月5日上机作业一：2-1 思考题：（1）数学模型的微分方程，状态方程，传递函数，零极点增益和部分分式五种形式，各有什么特点？（2）数学模型各种形式之间为什么要互相转换？（3）控制系统建模的基本方法有哪些？他们的区别和特点是什么？（4）控制系统计算机仿真中的“实现问题”是什么含意？（5）数值积分法的选用应遵循哪几条原则？（1）微分方程是直接描述系统输入和输出量之间的制约关系，是连续控制系统其他数学模型表达式的基础。状态方程能够反映系统内部各状态之间的相互关系，适用于多输入多输出系统。传递函数是零极点形式和部分分式形式的基础。零极点增益形式可用于分析系统的
2、稳定性和快速性。利用部分分式形式可直接分析系统的动态过程。（2）不同的控制系统的分析和设计方法，只适用于特定的数学模型形式。（3）控制系统的建模方法大体有三种：机理模型法，统计模型法和混合模型法。机理模型法就是对已知结构，参数的物理系统运用相应的物理定律或定理，经过合理的分析简化建立起来的各物理量间的关系。该方法需要对系统的内部结构和特性完全的了解，精度高。统计模型法是采用归纳的方法，根据系统实测的数据，运用统计规律和系统辨识等理论建立的系统模型。该方法建立的数学模型受数据量不充分，数据精度不一致，数据处理方法的不完善，很难在精度上达到更高的要求。混合法是上述两种方法的结合。（4）“实现问题”
3、就是根据建立的数学模型和精度，采用某种数值计算方法，将模型方程转换为适合在计算机上运行的公式和方程，通过计算来使之正确的反映系统各变量动态性能，得到可靠的仿真结果。（5）数值积分法应该遵循的原则是在满足系统精度的前提下，提高数值运算的速度和并保证计算结果的稳定。2-2.用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式：(1) G（s）= (2) =y=0 2 0 2 X解：（1）程序代码如下：clear;clc;num=1 7 24 24;den=1 10 35 50 24;Z,P,K=tf2zp(num,den);A,B
4、,C,D=zp2ss(Z,P,K);R,P1,H=residue(num,den);G1=tf(num,den)G2=zpk(Z,P,K)G3=ss(A,B,C,D)程序的输出结果如下： Transfer function: s3 + 7 s2 + 24 s + 24-s4 + 10 s3 + 35 s2 + 50 s + 24 Zero/pole/gain:(s+1.539) (s2 + 5.461s + 15.6)-(s+4) (s+3) (s+2) (s+1)a = x1 x2 x3 x4 x1 -3 -1.414 0 0 x2 1.414 0 0 0 x3 1 1.088 -7 -3.
5、464 x4 0 0 3.464 0 b = u1 x1 1 x2 0 x3 0 x4 0 c = x1 x2 x3 x4 y1 1 1.088 -1.539 1.038 d = u1 y1 0 Continuous-time model.（2）程序代码如下：clear;clc;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=4 2 2 0;C=0 2 0 2;D=0;num,den=ss2tf(A,B,C,D);Z,P,K=ss2zp(A,B,C,D);R,P1,H
6、=residue(num,den);G1=tf(num,den)G2=zpk(Z,P,K)G3=ss(A,B,C,D)程序输出结果如下：Transfer function: 4 s3 + 14 s2 + 22 s + 15-s4 + 4 s3 + 6.25 s2 + 5.25 s + 2.25 Zero/pole/gain:4 (s+1.5) (s2 + 2s + 2.5)- (s+1.5)2 (s2 + s + 1) a = x1 x2 x3 x4 x1 2.25 -5 -1.25 -0.5 x2 2.25 -4.25 -1.25 -0.25 x3 0.25 -0.5 -1.25 -1 x4
7、 1.25 -1.75 -0.25 -0.75 b = u1 x1 4 x2 2 x3 2 x4 0 c = x1 x2 x3 x4 y1 0 2 0 2 d = u1 y1 0 Continuous-time model.2-3.用欧拉法求下面系统的输出响应y(t)在0t1上，h=0.1时的数值。 要求保留4位小数，并将结果与真解比较。解：程序代码如下：clear;clc;h=0.1;y=1;disp(欧拉法求得值为：);for i=0:h:1 m=y; disp(y); y=m-m*h;enddisp(真解为：);for i=0:h:1 y=exp(-i); disp(y);end输出结果
8、如下：欧拉法值10.90000.81000.72900.65610.59050.53140.47830.43050.38740.3487真值10.90480.81870.74080.67030.60650.54880.49660.44930.40660.36792-4用二阶龙格库塔法求解2-3的数值解，并于欧拉法求得的结果比较。解：程序代码如下：clear;clc;h=0.1;y=1;disp(二阶龙格库塔法求得值为：);for i=0:h:1 disp(y); k1=-y; k2=-(y+k1*h); y=y+(k1+k2)*h/2;enddisp(欧拉法求得值为：); y=1;for i=
9、0:h:1 m=y; disp(y); y=m-m*h;end程序输出结果如下：欧拉法值10.90000.81000.72900.65610.59050.53140.47830.43050.38740.3487龙格库塔法值10.90500.81900.74120.67080.60710.54940.49720.45000.40720.36852-5用四阶龙格-库塔法求解题2-3数值解，并与前两题结果相比较。解：程序代码如下：clear;clc;h=0.1;y=1;disp(四阶龙格库塔法求得值为：);for i=0:h:1 disp(y); k1=-y; k2=-(y+k1*h/2); k3=
10、-(y+k2*h/2); k4=-(y+k3*h); y=y+(k1+2*k2+2*k3+k4)*h/6;end程序结果如下：欧拉法值10.90000.81000.72900.65610.59050.53140.47830.43050.38740.3487二阶龙格库塔法值10.90500.81900.74120.67080.60710.54940.49720.45000.40720.3685四阶龙格库塔法值10.90480.81870.74080.67030.60650.54880.49660.44930.40660.36792-6已知二阶系统状态方程为写出取计算步长为h时，该系统状态变量X=
11、的四阶龙格-库塔法递推关系式。2-7单位反馈系统的开环传递函数已知如下用matlab语句 、函数求取系统闭环零极点，并求取系统闭环状态方程的可控标准型实现。解：程序代码如下：clear;clc;num=5 100;den1=1 0;den2=1 4.6;den3=1 3.4 16.35;den4=0 0 0 5 100;den=conv(den1,den2);den=conv(den,den3);den=den+den4Z,P,K=tf2zp(num,den);P1=roots(den)Z1=roots(num)ZP输出结果如下：den = 1.0000 8.0000 31.9900 80.2
12、100 100.0000P1 = -0.9987 + 3.0091i -0.9987 - 3.0091i -3.0013 + 0.9697i -3.0013 - 0.9697iZ1 = -20Z = -20P = -0.9987 + 3.0091i -0.9987 - 3.0091i -3.0013 + 0.9697i -3.0013 - 0.9697i2-8用matlab语言编制单变量系统三阶龙格-库塔法求解程序，程序入口要求能接收状态方程各系数阵（A,B,C,D）,和输入阶跃函数r(t)=R*1(t);程序出口应给出输出量y（t）的动态响应数值解序列。解：函数程序代码如下：function y = hs( A,B,C,D,R
