计算机仿真复习

计算机仿真复习（1）给定被控对象2321()6116ssGssss,（a）试利用Matlab建立其状态空间表达式；MATLAB程序为：num=[111];den=[16116];[A,B,C,D]=tf2ss(num,den)运行结果为：A=-6-11-6100010B=100C=111D=0（b）(b)在采样周期0.1时进行离散化。MATLAB程序为：K=1,Z=[-1/2+1/2*i*3^(1/2);-1/2-1/2*i*3^(1/2)]P=[-1;-2;-3];,T=0.1;[A,B,C,D]=zp2ss(Z,P,K)运行结果为：A=-1.0000001.0000-5.0000-2.449502.44950B=100C=1.0000-4.0000-2.0412D=0（2）已知两子系统的传递函数分别为11(),(1)(2)Gsss21()(3)Gsss试利用Matlab求两个子系统串联、并联时的传递函数。并用Simulink实现，两种方法进行对比分析。1.用matlab求法：Matlab程序为：num1=1;den1=[131];num2=1;den2=[13];[num3,den3]=series(num1,den1,num2,den2)[num4,den4]=parallel(num1,den1,num2,den2)运行结果为：num3=0001den3=16103num4=0144den4=161032.用Simulink仿真方法Simulink编辑（3）已知系统的方框图如下图所示，求系统的传递函数。Matlab程序为：num1=10;den1=[11];num2=1;den2=[20.5];[num5,den5]=series(num1,den1,num2,den2);num4=0.1;den4=1;[num6,den6]=feedback(num5,den5,num4,den4);num3=540;den3=1;[num7,den7]=series(num3,den3,num6,den6);num8=1;den8=1;[num,den]=feedback(num7,den7,num8,den8);printsys(num,den)运行结果为：num/den=5400----------------------2s^2+2.5s+5401.5（4）试用三种方法分别判断下述系统的稳定性。（a）11221123xxxx(b)代码1：A=[-11;2-3];p=poly(A);r=roots(P);ii=find(real(r)0);n=length(ii);if(n0)disp('SystemisUnstable');elsedisp('SystemisStable');432543222()521ssssGssssssend显示结果:r=-5.50640.4371+0.7643i0.4371-0.7643i-0.6839+0.4849i-0.6839-0.4849iSystemisStable代码2：num=[12112];den=[151121];r=roots(den),pzmap(num,den)显示结果：r=-4.81960.4133+0.6826i0.4133-0.6826i-0.5035+0.2689i-0.5035-0.2689i有上图和运行结果可以看出在右半s平面有两个极点，所以系统不稳定。程序段2：num=[12112];den=[151121];[num,den]=cloop(num,den);r=roots(den),zplane(num,den)r=-5.50640.4371+0.7643i0.4371-0.7643i-0.6839+0.4849i-0.6839-0.4849i由结果可知，系统在z平面的单位圆外有一个极点，故系统不稳定。(5)对于典型二阶系统222()2nnnGsss试绘制出无阻尼自然振荡频率3n，阻尼比ζ分别为0.2,0.6,1.0,1.4,1.8,2.2,2.6,3.0时系统的单位阶跃响应曲线。wn=3;zata=[0.2:0.4:3];figure(1)holdonforA=zatanum=wn.^2;den=[12*A*wnnum];step(num,den)endtitle('单位阶跃响应曲线')holdoff运行结果：02468101200.20.40.60.811.21.41.6单位阶跃响应曲线Time(sec)Amplitude仿真技术是一门利用物理模型或数学模型模拟实际环境进行科学实验的技术。仿真的优点是：经济、可靠、实用、安全、灵活和多次重复使用。它所遵循的基本原则是相似性原理。仿真分类：依据相似原理，分为：物理仿真，数学仿真和混合仿真依据被研究系统的特征分为：连续系统仿真和离散事件系统仿真控制系统辅助设计变换法（频域法）单输入单输出系统状态空间法（时域法）多变量控制系统矩阵运算1、矩阵的元素表示2、特殊矩阵表示3、矩阵基本运算绘图1、subplot2、plot系统模型的转换模型之间的连接（串、并、反馈）系统离散化simulink的启动（三种方法）基本模块给出模块结构图，写出表达式稳定性判断时域分析（单位阶跃响应step、脉冲响应impluse、零输入响应initial、任意输入响应lsim）根轨迹rlocus伯德图bode、奈奎斯特图nyquist能控性ctrb、能观测性obsvgram极点配置的状态反馈acker\place