Matlab机电系统仿真-1.ppt

仍以前面的command_in_out 为例进行说明。仿真参数设置对话框内的设置与前面相同。使用四组不同的仿真时间对此系统进行仿真。创建m文件“command_in_out_m”,并输入代码如下： t=0:0.1:10;t=t; u=sin(t); sim_input=[t u]; [tout1,x1,yout1]=sim(command_in_out,5); %仿真时间范围0 ～5s ，输出时间向量tout1 由Simulimk的求解器步长决定 [tout2,x2,yout2]=sim(command_in_out,[1 8]); ％仿真时间范围1～8s，输出时间向量 tout2 由 Simulimk 的求解器步长决定 [tout3,x3,yout3]=sim(command_in_out,[1:8]); %仿真时间范围1 ～8s ，并且每隔1s 输出一次，即输出时间变量为[1 2 3 4 5 6 7 8] [tout4,x4,yout4]=sim(command_in_out,1:0.2:8); %仿真时间范围1 ～8s ，并且每隔0.2s 输出一次，即输出时间变量为[1 1.2 1.4 … 7.6 7.8 8 subplot(2,2,1); plot(t,u,tout1,yout1,*); subplot(2,2,2); plot(t,u,tout2,yout2,*); subplot(2,2,3); plot(t,u,tout3,yout3,*); subplot(2,2,4); plot(t,u,tout4,yout4,*); 外部输入变量设置 MATLAB 程序（程序名：command_in_out_m2.m） t=0:0.1:10; t=t; u=sin(t); sim_input=[t, u]; [tout1, x1, yout1]=sim (command_in_out,10); % 使用Simulink 仿真参数对话框中的workspace I/O 从MATLAB 工作空间中获得输入信号 u= cos(t); ut=[t, u]; %改变系统输入信号 [tout2, x2, yout2]=sim (command_in_out, 10, [], ut); %使用Sim 中的ut 参数获得系统输入信号，ut 的使用会覆盖由Workspace I/O 的系统输入设置，这一点可以在下面的系统仿真结果图形中反映出来 subplot(1, 2, 1); plot (tout1, yout1); grid subplot(1, 2, 2); plot (tout2, yout2); grid 图中左图表示系统输入为sin(t) 时的相应曲线，右图表示cos(t) 时的相应曲线。从图中可明显看出，当使用sim 命令的ut 参数时，Simulink 仿真参数设置对话框中的设置被覆盖。以前对话框中的外部输入是名为sim_input 的正弦信号，而采用ut 参数后执行的余弦输入信号。注意：这里指的“ 覆盖” ，并不是在Workspace I/O 对话框的Input 中，将sim_input 改变成了ut ，事实上并没有改变，只是不执行sim_input，而执行了命令行中的ut sim其他命令 simset 命令是用来创建和编辑options 结构的。 simget 命令用来获得系统模型的仿真参数设置。 simplot 命令 9、微分方程求解 求微分方程组的通解. y = 3*sin(5*x)*exp(-2*x) [y]=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x); y=simple(y) 创建m文件 解析解：dsolve 数值解：ode 求微分方程组的通解. [x,y,z]=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z, t); x=simple(x) ; % 化简 y=simple(y) ; z=simple(z) ; [x;y;z] ans = C8*exp(2*t) + C9*exp(-t) C8*exp(2*t) + C9*exp(-t) + C10*exp(-2*t) C8*exp(2*t) + C10*exp(-2*t) 创建m文件 求微分方程组的通解. 建立M文件vdp1000．m如下： function dy=vdp1000(t,y) dy=zeros