控制新版系统数字仿真试验参考指导书.docVIP

淮 阴 工 学 院 控制系统数字仿真试验指导书 编者：王建国 适用学院：电子和电气工程学院 电子和电气工程学院 6 月 试验一 面向一阶微分方程组连续系统数字仿真试验 一、试验目标 1．掌握以一阶微分方程组形式描述连续系统数字仿真方法和步骤。 2. 了解怎样编程实现对一阶微分方程组模型仿真研究。 二、试验内容 单位反馈控制系统结构图A.1所表示。 图A.1 系统结构 1．依据系统闭环传输函数，求出对应状态空间模型，并利用MATLAB中ODE解函数编程，对系统进行仿真研究，并绘制对应单位阶跃响应曲线。 2．考虑微分方程模型 试编程对其进行仿真（0≤t≤100），并用不一样色彩、线型绘制仿真结果曲线。 三、预习要求 1. 复习ODE解函数使用方法和调用格式。 2．复习一阶微分方程组仿真程序编制方法。 四、试验汇报要求 1. 整理多种试验条件下数据和曲线。 2. 写出试验汇报。 五、教案 1．利用MATLAB编程求出图A.1所表示系统状态空间模型程序（文件名为Appl1_1_1.m）以下： % 这是试验一中第一个系统求取状态空间模型程序 clear % 输入三个步骤系统传输函数 sys1=tf([1],[0.5,1]); sys2=tf([1],[0.2,1]); sys3=tf([1],[1,0]); % 求系统闭环传输函数 sys=sys1*sys2*sys3; sys_cl_tf=feedback(sys, -1); % 求系统状态空间模型 sys_cl_ss=ss(sys_cl_tf) 运行Appl1_1_1.m后，得到系统状态空间模型为： 据此，编制仿真程序（文件名为Appl1_1_1_fang.m）以下： % 这是试验一中第一个系统仿真程序 clear x0=[0,0,0]; tspan=[0,10]; [t,x]=ode45(Appl_1_1_func, tspan, x0); plot(t, 1.25*x(:,3)); xlabel(t); ylabel(y(t)); grid; title(试验一 第一个系统单位阶跃响应); 其中，Appl_1_1_func.m程序以下： % 这是试验一中第一个系统对应数学模型 function xdot=Appl_1_1_func(t, x) xdot1(1)=-7*x(1)-2.5*x(2)-1.25*x(3)+1; xdot1(2)=4*x(1); xdot1(3)=2*x(2); xdot=xdot1; 运行Appl1_1_1_fang.m后，得到系统单位阶跃响应曲线图A.1_1所表示。 图A.1_1 试验一第一个系统单位阶跃响应 2. 利用MATLAB编程编制求解系统 仿真程序（文件名为Appl1_1_2_fang.m）以下： % 这是试验一中第二个系统仿真程序 clear x0=[10,10]; tspan=[0,100]; [t,x]=ode45(Appl_1_2_func, tspan, x0); plot(t, x(:,1), k, LineWidth,2); hold plot(t, x(:,2), r-., LineWidth,2); legend(x(1), x(2)); xlabel(t); ylabel(状态向量); grid; title(试验一 第二个系统仿真输出); 其中，Appl_1_2_func.m程序以下： % 这是试验一中第二个系统对应数学模型 function xdot=Appl_1_2_func(t, x) xdot1(1)=- x(1) *(x(1)*x(1)+x(2)*x(2))+x(2); xdot1(2)=-x(1)-x(2) * (x(1)*x(1)+x(2)*x(2)); xdot=xdot1; 运行Appl1_1_2_fang.m后，得到系统仿真曲线图A.1_2所表示。 图A.1_2 试验一第二个系统仿真曲线 试验2 面向系统结构图连续系统数字仿真试验 一、试验目标 （1）掌握以系统结构图形式描述连续系统数字仿真方法和步骤。 （2）初步了解怎样用仿真方法来分析系统动态性能。 （3）了解不一样数值积分算法和仿真计算精度之间关系。 （4）学会一个初步寻求合理仿真步长方法。 （5）了解RK4法计算稳定性和步长关系。 二、试验内容 单位反馈控制系统结构如附录图1所表示，其开环传输函数和闭环传输函数分别为 在零初始条件下闭环系统单位阶跃响应标准解为 按试验目标、要求和已知条件，建立系统Simulink模型。 2. 按下面经验公式（1）初选仿真步长 （为此，需要学习MATLAB中求取开环剪切频率 函数margin使用方法）。 （1） 为系统中反应最快小闭环开环剪切频率。