西电自动化系统仿真作业欧拉法仿真系统.doc

系统仿真作业四 问题：已知系统传递函数为 ,根据传递函数写出系统的状态方程，编写matlab程序，应用欧拉法计算系统在单位阶跃信号作用下的输出，并与准确输出做对比。（设定步长h1,选取不同的步长做多次计算） 解答：根据系统的传递函数，求得状态方程为： 输入U(t)为阶跃信号，故X1(0)=0，X2(0)=0, 由上式推导出相应的欧拉公式为： 根据上述欧拉公式，编写相应的.m文件，代码如下 x1=0; %x1初始条件 x2=0; %x2初始条件 h=0.01; %选取步长 for i=1:1000 %依照欧拉公式计算输出 x1(i+1)=(1-100*h)*x1(i)+100*h*x2(i); x2(i+1)=(1-h)*x2(i)+h; end t=0:0.01:10; %计算系统的准确输出并作图 subplot(2,2,1) step(100,[1,101,100],t) subplot(2,2,2) plot(t,x1,r) %用欧拉公式计算出来的点画图 xlabel(times) subplot(2,1,2) %将两幅图画在一起做对比 step(100,[1,101,100],t) hold on plot(t,x1,r) 当取h为0.1时，得到如下图像 观察图像得知，当步长取0.1时，利用欧拉法已无法计算出数值解，结算结果呈发散状态。仿真时间为50时，欧拉法计算出的最大值达到量级。 当取步长为0.05时 观察图像，当步长取0.05时，利用欧拉法已无法计算出数值解，结算结果呈发散状态。仿真时间为50时，欧拉法计算出的最大值达到量级。 当选取步长为0.02时，结果如下 将欧拉法计算的结果图像单独列出来 观察可得，当h取0.02时，欧拉法计算的数值总体变化趋势与准确值相仿，但可以看到，在具体某一点时，欧拉法计算出的数值已有一定的误差。 当取步长h为0.01时，结果如下 （左上为准确解图像，右上为欧拉公式数值解图像，下方为对比图像）观察图像，当步长为0.01时，欧拉法计算出来的解几乎与准确解相同。 选取h为0.001时，结果如下 观察图像，得知当h取0.001时，欧拉公式计算出的数值与准确解基本完全吻合，不存在误差。 结论：步长的选取影响欧拉公式计算结果的准确性，当步长h选取为0.01，欧拉公式计算结果与准确值基本吻合，无误差，当h大于0.01时，欧拉公式计算结果开始出现误差，随着步长的逐渐增大，欧拉公式计算结果的误差将越来越大，最终趋于发散。