计算机仿真实验报告.docVIP

PAGE PAGE 9 目 录 TOC \o 1-3 \h \z \u 实验一 Matlab语言编程 1 一. 实验目的 1 二. 具体实验内容、步骤、要求： 1 实验二 数值积分算法及函数调用练习 3 一． 实验目的 3 二． 实验实例介绍： 3 实验三 控制工具箱与SIMULINK软件应用 9 一． 实验目的 9 二． 实验预习要求： 9 三． 学会调出、运行已由SIMULINK建立的仿真模型。 9 四． 实验设计题目与要求： 10 实验一 Matlab语言编程 实验目的 熟悉Matlab语言及其编程环境，掌握编程方法 要求认真听取实验指导老师讲解与演示 具体实验内容、步骤、要求： 运行交互式学习软件，学习Matlab语言 在Matlab的命令窗口下输入如下命令：INTRO,然后根据显示出来的幻灯片右面按钮进行操作，可按START—NEXT—NEXT按钮，一步步运行，观察。 自编程序并完成上机编辑、调试、运行，存盘。 (1). 用Matlab命令完成矩阵的各种运算，例如： 求出下列运算结果，并上机验证。 A(:,1)，A(2,:)，A(1:2,2:3)，A(2:3,2:3)，A(:,1:2)，A(2:3), A(:),A(:,:),ones(2,2), eye(2)  (2). 绘制数学函数图形 t=0:0.1:8; y=1-2*t.*sin(t); plot(t,y) 理解命令文件和函数文件的区别，并自编函数文件并调用。 学会通过Help,熟悉Matlab中为用户提供的功能各异的命令和函数。 实验二 数值积分算法及函数调用练习 实验目的 理解数值积分算法，熟练掌握Matlab的函数调用。 实验实例介绍： 1. 用Euler法求初值问题的数值解。 设方程如下： （2-1式） t=[0,1] 取步长h=0.1，名为FZSYZ1.M 上机用如下程序名FZSYZ1.M可求出数值解。 t0=0; t0=0; tf=1; x0=1; h1=0.1; t=[t0:h1:tf]; n=length(t); u=x0; uu(1)=u; for i=2:n du=u-2*t(i-1)/u; u=du*h1+u; uu(i)=u; end plot(t,uu,b-) 2. 在Matlab中提供了现成的数值积分的函数，如ode1,ode23,ode45求解微分方程组，下面介绍ode23它为二阶/三阶的RKF法在MATLAB的ToolBox文件夹的MATLAB/funfun下的M文件，在此介绍其调用方法与应用例子如下： 其中，系统函数名为描述系统状态方程的M函数的名称，该函数名在调用时应该用引号括起来（文件名字串） T0,tf为起始和终止时间 X0系统初始状态变量的值（列向量） T01:控制解的精度。（缺省值为在ode23中） Trace:输出形式控制变量，非零则程序运行的每步都显示出来。 T:输出参数返回积分的时间离散值（列向量） X0:输出参量，返回每个时间点值的解的列向量 注意：系统函数的编写格式为固定的。 Function xdot=系统函数名（t,x） 例如：若MATLAB求解初值问题的解，其方程如下： （2-2式）；x(0)=1;t=[0 3] 第一步：编写如下程序： function xdot=fun22(t,x) xdot=x-2*t/x 并以fun21.m存盘 第二步：编写如下程序并以fzsy22.m存盘 t0=0;tf=1;tol=1e-6 x0=1;trace=1; [t,x]=ode23(fun22,t0,tf,x0,tol,trace) plot(t,x,g-.) 第三步：在命令窗口中运行fzsy22即可求出x的解，并画出曲线，如下： 图2-1 Euler方程和ode23算法求解的解曲线 3．实验具体内容、步骤、要求： 运行交互式软件中函数调用，学习程序； 试将（2-2）方程改为用Euler方程求解是比较用ode23求解结果。 Euler法程序： t0=0; t0=0; tf=3; x0=1; h1=0.1; t=[t0:h1:tf]; n=length(t); u=x0; uu(1)=u; for i=2:n du=u-t(i-1).^2; u=du*h1+u; uu(i)=u; end plot(t,uu,b-) 图形如2-1所示 试将（2-1）方程改为用ode23算法调用函数求解，并试比较结果。 F Fun22.m文件 function xdot=fun22(t,x) xdot=x-2*t/x Fun22_.m文件 t0=0;tf=1;tol=1e-6 x0=1;trace=1; [t,x]=ode23(fun