控制系统数字仿真技术实验

实验一控制系统的模型转换、数据处理与曲线拟合（2学时）一．实验目的：掌握控制系统的微分方程、状态方程、传递函数、零极点增益、部分分式描述及转换；掌握常用数据拟合与插值方法。二．实验方法及预习内容：1．利用Matlab工具箱中常用的五种模型转换命令进行模型描述和转换；2．利用Matlab工具箱中的多项式拟合命令polyfit对实验数据进行拟合。三．实验内容：1．用Matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式：（1）G（s）=32432144848207010048sssssss（2）.X=2.25-5-1.25-0.542.25-4.25-1.25-0.2520.25-0.5-1.25-121.25-1.75-0.25-0.750Xuy=[0202]X2．已知元件的实验数据如下，拟合这一数据，并尝试给出其特性方程。X0.01001.01002.01003.01004.0100Y2.54377.88849.624211.607111.9727X5.01006.01007.01008.01009.0100y13.218914.267914.613415.404515.0805四．实验总体要求：1．每次实验前应做好预习和准备；2．实验后应及时提交仿真程序Word文档（包括：程序、实验结果和图示、实验分析与总结）；3．重点检查Matlab中M文件的运行；4．请认真撰写实验报告，实验结果可贴在报告相应位置。注：格式不符合要求、或字迹潦草的一律退回重写。五．本次实验要求：1．熟悉五种连续系统控制模型的Matlab转换命令，并得出相应数学模型表达式；2．熟悉常用数据拟合方法。（1）G（s）=32432144848207010048sssssss程序den=[1,20,70,100,48];num=[1,14,48,48];[A,B,C,D]=tf2ss(num,den)[z,p,k]=ss2zp(A,B,C,D)[R,P,H]=residue(num,den)结果A=-20-70-100-48100001000010B=1000C=1144848D=0z=-9.4641-2.5359-2.0000p=-16.0051-1.5269+0.9247i-1.5269-0.9247i-0.9412k=1.0000R=0.3892-0.0932-0.3754i-0.0932+0.3754i0.7973P=-16.0051-1.5269+0.9247i-1.5269-0.9247i-0.9412H=[]（2）.X=2.25-5-1.25-0.542.25-4.25-1.25-0.2520.25-0.5-1.25-121.25-1.75-0.25-0.750Xuy=[0202]X程序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,P,H]=residue(num,den)结果num=04.000014.000022.000015.0000den=1.00004.00006.25005.25002.2500z=-1.0000+1.2247i-1.0000-1.2247i-1.5000p=-0.5000+0.8660i-0.5000-0.8660i-1.5000-1.5000k=4.0000R=4.0000-0.00000.0000-2.3094i0.0000+2.3094iP=-1.5000-1.5000-0.5000+0.8660i-0.5000-0.8660iH=[]实验二基于Matlab的微分方程数值解法（2学时）一．实验目的：掌握欧拉法、四阶龙格库塔法的程序编制方法。二．实验方法及预习内容：对微分方程描述的控制系统，利用欧拉法、二阶龙格-库塔法、四阶龙格-库塔法分别编写M文件，进行数值计算和作图。三．实验内容：1．分别用欧拉法、二阶龙格-库塔法、四阶龙格-库塔法求下面系统的输出响应y(t)在0≤t≤1上，h=0.1时的数值解。'2,(0)1yyy要求保留4位小数，并将三种方法的结果与真解2()tyte进行比较。2.若为2',(0)1yyy如何编程计算？四．实验总体要求：1．每次实验应做好实验前的预习和准备；2．实验后应及时提交仿真程序（M文件与Word文档）、实验结果和图示、实验分析与总结；3．认真撰写实验报告。五．本次实验要求：熟悉微分方程数值解法（欧拉法、二阶与四阶龙格-库塔法）。1．分别用欧拉法、二阶龙格-库塔法、四阶龙格-库塔法求下面系统的输出响应y(t)在0≤t≤1上，h=0.1时的数值解。'2,(0)1yyy要求保留4位小数，并将三种方法的结果与真解2()tyte进行比较。程序欧拉法h=0.1;disp('欧拉法求函数的数值解为');disp('y=');y=1;fort=0:h:1m=y;disp(y);y=m-2*m*h;End结果欧拉法求函数的数值解为y=10.80000.64000.51200.40960.32770.26210.20970.16780.13420.1074二阶龙格—库塔法程序h=0.1;disp('二阶龙格—库塔法函数的数值解为');disp('y=');y=1;fort=0:h:1disp(y);k1=-2*y;k2=-2*(y+k1*h);y=y+(k1+k2)*h/2;End结果二阶龙格—库塔法函数的数值解为y=10.82000.67240.55140.45210.37070.30400.24930.20440.16760.1374四阶龙格—库塔法程序h=0.1;disp('四阶龙格—库塔法求解函数数值解为');disp('y=');y=1;fort=0:h:1disp(y);k1=-2*y;k2=-2*(y+k1*h/2);k3=-2*(y+k2*h/2);k4=-2*(y+k3*h);y=y+(k1+2*k2+2*k3+k4)*h/6;end结果四阶龙格—库塔法求解函数数值解为y=10.81870.67030.54880.44930.36790.30120.24660.20190.16530.1353真解程序h=0.1;disp('函数的离散时刻解为');disp('y=');y=1;fort=0:h:1y=exp(-2*t);disp(y);end结果函数的离散时刻解为y=10.81870.67030.54880.44930.36790.30120.24660.20190.16530.13532.若为2',(0)1yyy如何编程计算？欧拉法程序h=0.1;disp('欧拉法求函数的数值解为');disp('y=');y=1;fort=0:h:1m=y;disp(y);y=m+m*m*h;End结果欧拉法求函数的数值解为y=11.10001.22101.37011.55781.80052.12462.57603.23974.28926.1289二阶龙格—库塔法程序h=0.1;disp('二阶龙格—库塔法函数的数值解为');disp('y=');y=1;fort=0:h:1disp(y);k1=y*y;k2=y*(y+k1*h);y=y+(k1+k2)*h/2;end结果二阶龙格—库塔法函数的数值解为y=11.10501.23381.39551.60381.88162.26902.84233.76495.44929.2277四阶龙格—库塔法程序h=0.1;disp('四阶龙格—库塔法求解函数数值解为');disp('y=');y=1;fort=0:h:1disp(y);k1=y*y;k2=y*(y+k1*h/2);k3=y*(y+k2*h/2);k4=y*(y+k3*h);y=y+(k1+2*k2+2*k3+k4)*h/6;end结果四阶龙格—库塔法求解函数数值解为y=11.10521.23431.39651.60581.88552.27672.85873.80475.56599.7082实验三基于SIMULINK的控制系统时域分析（2学时）一．实验目的：掌握使用SIMULINK、控制工具箱求解系统的输入和输出响应的仿真方法。二．实验方法及预习内容：利用SIMULINK工具进行控制系统模型分析、系统设计与仿真的相关原理。三．实验内容：1．分别使用解微分方程方法、控制工具箱、Simulink求解具有如下闭环传递函数的系统的阶跃响应。43230()8364010sssss2．某小功率随动系统动态结构如图所示，已知：10.01T，20.05T，01K，1300K，21K，0.08cK。若系统输入分别为()1()srtt，srt，[1()1(1.5)]srt，试用Simulink分析系统的输出()sct分别如何？四．实验要求：1．熟悉Simulink法进行控制系统时域分析的基本步骤，并与微分方程法、控制工具箱法进行比较；2．熟悉Simulink法在输入不同信号（单位阶跃、单位斜坡、两阶跃叠加）下的控制系统输出响应实验方法、图形显示。控制工具箱法num=[30];den=[18364010];sys=tf(num,den);step(sys);Grid解微分方程法Num=[];den=[];[A,B,C,D]=tf2ss(num,den);建立wffc.m文件functiondx=wffc(t,x)u=1;dx=[-8*x(1)-36*x(2)-40*x(3)-10*x(4)+u;x(1);x(2);x(3)];主程序[t,x]=ode45('wffc',[0,20],[0,0,0,0]);y=30*x(:,4);plot(t,y);grid实验四面向系统结构图的仿真验证（2学时）一．实验目的：从控制系统常见的结构形式拓扑描述入手，掌握面向连续控制系统结构图的计算机仿真方法及其程序实现。二．实验原理及预习内容：1．原理：任何复杂连接结构的线性控制系统都是由一些简单的线性环节组合而成，按照它们之间相互连接的拓扑关系列出连接矩阵，可以得到能清晰地描述复杂连接系统的仿真模型。2．预习内容：利用连接矩阵进行复杂控制系统建模的方法和原理。三．实验步骤：1．对具有复杂连接闭环结构形式的系统，可用一阶环节作为典型环节，再运用拓扑描述中联接矩阵的表达方法得出此类系统结构的仿真模型；2．再通过数值积分法求取各环节的动态响应。注意：所确定的典型环节中，参数0iB，以保证系统仿真求解的基本条件。四．实验内容：习题3-2．设典型闭环结构控制系统如下图所示，当阶跃输入幅值R=20时，用面向系统结构图的数字仿真法sp3-1.m求取系统的输出响应。五．实验要求：1．列出复杂连接闭环系统的仿真程序框图；2．列出MATLAB程序实现的主要程序，包括输入数据、形成闭环各系统阵、数值积分求解、以及输出结果。解：sp3_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),zeros(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+