计算机仿真技术实验报告29633.docVIP

计算机仿真技术实验报告 学院： 专业： 班级： 姓名： 学号： 实验一MATLAB的ODE45算法分别求解下列二个方程。要求:1.编写出Matlab仿真程序；2.画出方程解的图形并对图形进行简要分析；3.分析下列二个方程的关系。 2． 程序代码 方程一： function fun=fun1(t,x) fun=-x^2; [t,x]=ode45(fun1,[0,1],1); plot(t,x,o,t,x,-);xlabel(time t0=0,tt=20);ylabel(x values x(0)=1);grid 执行结果 使用ode45 方程二： function fun=fun2(t,x) fun=x^2; [t,x]=ode45(fun2,[0,1],-1);plot(t,x,o,t,x,-);xlabel(time t0=0,tt=1);ylabel(x values x(0)=-1);grid 执行结果 使用ode45 题目二：下面方程组用在人口动力学中，可以表达为单一化的捕食者-被捕食者模式（例如，狐狸和兔子）。其中表示被捕食者， 表示捕食者。如果被捕食者有无限的食物，并且不会出现捕食者。于是有，则这个式子是以指数形式增长的。大量的被捕食者将会使捕食者的数量增长；同样，越来越少的捕食者会使被捕食者的数量增长。而且，人口数量也会增长。请分别调用ODE45、ODE23算法求解下面方程组。要求编写出Matlab仿真程序、画出方程组解的图形并对图形进行分析和比较。 程序代码 function fun=fun3(t,x) fun=[x(1)-0.1*x(1)*x(2)+0.01*t;-x(2)+0.02*x(1)*x(2)+0.04*t]; [t,x]=ode45(fun3,[0,20],[30,20]); plot(t,x);xlabel(time t0=0,tt=20);ylabel(x values x1(0)=30,x2(0)=20);grid 执行结果 [t,x]=ode23(fun3,[0,20],[30,20]); plot(t,x);xlabel(time t0=0,tt=20);ylabel(x values x1(0)=30,x2(0)=20);grid 2. Matlab中模型表示及模型转换指令的使用 题目三：若给定系统的的传递函数为 请用MATLAB编程求解其系统的极零点模型。 程序代码 num=[6 12 6 10]; den=[1 2 3 1 1]; sys=tf(num,den); H=zpk(sys) 执行结果： Zero/pole/gain: 6 (s+1.929) (s^2 + 0.07058s + 0.8638) ------------------------------------------------- (s^2 + 0.08663s + 0.413) (s^2 + 1.913s + 2.421) 题目四：习题2.4 系统传递函数为，求其对角标准型实现。 程序代码： num=[1 4 5]; den=[1 6 11 6]; sys=tf(num,den); canon(sys,modal) 执行结果： a = x1 x2 x3 x1 -3 0 0 x2 0 -2 0 x3 0 0 -1 b = u1 x1 -15.11 x2 -22.98 x3 8.958 c = x1 x2 x3 y1 -0.06619 0.04352 0.1116 d = u1 y1 0 Continuous-time model. 题目五：习题5.8 已知某离散系统，采样周期，先需要在下作数字仿真，求该系统的数字仿真模型。 分步执行： D=zpk(0,[0.1 0.2],1,0.1) Zero/pole/gain: z --------------- (z-0.1) (z-0.2) Sampling time: 0.1 D2=d2d(D,0.02) Zero/pole/gain: 0.278 (z-0.4926) -------------------- (z-0.631) (z-0.7248) Sampling time: 0.02 直接执行： D=zpk(0,[0.1 0.2],1,0.1); D2=d2d(D,0.02) Zero/pole/gain: 0.278 (