控制系统仿真及CAD试题.docVIP

控制系统仿真及CAD试题（研2010） （20分）试论述系统仿真的目的、意义、分类及应用与发展概况。 解：系统仿真的目的：在分析系统各要素性质及其相互关系的基础上，建立能描述系统结构或行为过程的、且具有一定逻辑关系或数量关系的仿真模型，据此进行试验或定量分析，以获得正确决策所需的各种信息。的输出响应y(t)在0≤t≤1上，h=0.1时的数值。要求保留4位小数，并将结果与真解比较。 解：欧拉法（前向欧拉法,可以自启动）其几何意义：把f(t,y)在[]区间内的曲边面积用矩形面积近似代替。利用matlab提供的m文件编程，得到算法公式。如下所示 （1） m文件程序为 h=0.1; disp(函数的数值解为); %显示 ‘’中间的文字% disp(y=); %同上% y=1; for t=0:h:1 m=y; disp(y); %显示y的当前值% y=m-m*h; end 保存文件q2.m 在matalb命令行中键入 q2 得到结果 函数的数值解为 y= 1 0.9000 0.8100 0.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487 （2）另建一个m 文件求解在t[0,1]的数值 （ %是的真解%） 程序为h=0.1; disp(函数的离散时刻解为); disp(y=); for t=0:h:1 y=exp(-t); disp(y); end 保存文件q3.m 在matlab命令行中键入 q3 函数的离散时刻解为 y= 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 比较欧拉方法求解与真值的差别 欧拉 1 0.9000 0.8100 0.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.3874 0.3487 真值 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 误差 0 -0.0048 -0.007 –0.0118 –0.0142 –0.0160 –0.0174 –0.0183 –0.0188 -0.0192 -0.0192 显然误差与为同阶无穷小，欧拉法具有一阶计算精度，精度较低，但算法简单。 我们经常用到 预报-校正法 的二阶龙-格库塔法， 此方法可以自启动，具有二阶计算精度，几何意义：把f(t,y)在[]区间内的曲边面积用上下底为和、高为h的梯形面积近似代替。利用matlab提供的m文件编程，得到算法公式。如下所示 （1）m文件程序为 h=0.1; disp(函数的数值解为); disp(y=); y=1; for t=0:h:1 disp(y); k1=-y; k2=-(y+k1*h); y=y+(k1+k2)*h/2; end 保存文件q4.m 在matlab的命令行中键入 q4 显示结果为 函数的数值解为 y= 1 0.9050 0.8190 0.7412 0.6708 0.6071 0.5494 0.4972 0.4500 0.4072 0.3685 比较欧拉法与二阶龙格-库塔法求解.（误差为绝对值） 真值 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.4066 0.3679 龙库 1 0.9050 0.8190 0.7412 0.6708 0.6071 0.5494 0.4972 0.4500 0.4072 0.3685 误差 0 0.0002 0.0003 0.0004 0.0005 0.0006 0.0006 0.0006 0.0007 0.0006 0.0006 明显误差为的同阶无穷小，具有二阶计算精度，而欧拉法具有一阶计算精度，二阶龙格-库塔法比欧拉法计算精度高。 （20分）分别使用解微分方程方法、控制工具箱、simulink求解具有如下闭环传递函数的系统的阶跃响应。 解：（1）用解微分方程方法：将转化为状态方程，利用matlab语句 num=[10