下载文档原格式
机械与汽车工程学院《Matlab控制系统仿真》实验指导书学院班级姓名学号浙江科技学院机械与汽车工程学院制实验一 MATLAB语言基本命令1 实验目的1. 掌握科学计算的有关方法，熟悉MA TLAB语言及其在科学计算中的运用；2. 掌握MATLAB的命令运行方式和M文件运行方式；3. 掌握矩阵在MA TLAB中的运用。
2 实验器材计算机WinXP、Matlab7.0软件3 实验内容(1). 输入A=[7 1 5;2 5 6;3 1 5]，B=[1 1 1; 2 2 2;3 3 3]，在命令窗口中执行下列表达式，掌握其含义：A(2, 3) A(:,2) A(3,:) A(:,1:2:3)A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*BA^2 A.^2 B/A B./A(2).输入C=1:2:20，则C（i）表示什么？其中i=1,2,3, (10)(3).查找已创建变量的信息，删除无用的变量；(4). 试用help命令理解下面程序各指令的含义：cleart =0:0.001:2*pi;subplot(2,2,1);polar(t, 1+cos(t))subplot(2,2,2);plot(cos(t).^3,sin(t).^3)subplot(2,2,3);polar(t,abs(sin(t).*cos(t)))subplot(2,2,4);polar(t,(cos(2*t)).^0.5)4 实验步骤：打开MA TLAB程序，将实验内容中的题目依次输入MATLAB中，运行得到并记录结果，最后再对所得结果进行验证。
5 实验报告要求记录实验数据,理解其含义实验二 MATLAB语言程序设计1 实验目的（1）掌握Matlab程序的编制环境和运行环境。
（2）掌握Matlab程序的编写方法。
（3）能编写基本的数据处理Matlab程序。
（4）能编写基本的数据可视化Matlab程序。
2 实验器材计算机WinXP、Matlab7.0软件3 实验内容(1) Matlab脚本文件编写和执行(2) Matlab 函数文件的编写和调用(3) nargm和nargout函数使用方法(4) 局部变量与全局变量使用4 实验步骤1、Matlab命令文件编写(1) 建立自己工作目录，如/Mywork。
实验总要求1、封面必须注明实验名称、实验时间和实验地点，实验人员班级、学号（全号）和姓名等。
2、内容方面：注明实验所用设备、仪器及实验步骤方法；记录清楚实验所得的原始数据和图像，并按实验要求绘制相关图表、曲线或计算相关数据；认真分析所得实验结果，得出明确实验结论。
3、图形可以打印出来并剪贴上去，文字必须用标准试验纸手写。
实验一MATLAB绘图基础一、实验目的了解MATLAB常用命令和常见的内建函数使用。
熟悉矩阵基本运算以及点运算。
掌握MATLAB绘图的基本操作：向量初始化、向量基本运算、绘图命令plot,plot3,mesh,surf 使用、绘制多个图形的方法。
二、实验内容建立并执行M文件multi_plot.m，使之画出如图的曲线。
三、实验方法（参考程序）024681012Plot of y=sin(2x) and its derivative四、实验要求1. 分析给出的MA TLAB 参考程序，理解MA TLAB 程序设计的思维方法及其结构。
2. 添加或更改程序中的指令和参数，预想其效果并验证，并对各语句做出详细注释。
对不熟悉的指令可通过HELP 查看帮助文件了解其使用方法。
达到熟悉MA TLAB 画图操作的目的。
3. 总结MATLAB 中常用指令的作用及其调用格式。
五、实验思考1、实现同时画出多图还有其它方法，请思考怎样实现，并给出一种实现方法。
(参考程序如下)t=0:pi/100:4*pi;y1=sin(2*t);y2=2*cos(2*t);plot(t,y1,'-b');hold on; %保持原图plot(t,y2,'-g');grid onaxis([0 4*pi -2 2])title('Plot of y=sin(2x) and its derivative')Plot of y=sin(2x) and its deriv ativ e024681012024681012-2-1012xyPlot of y=sin(2x)024681012-2-1012xyPlot derivative of y=sin(2x);y=2cos(2x)t=0:pi/100:4*pi; y1=sin(2*t); y2=2*cos(2*t);024681012-2-1.5-1-0.500.511.52Plot of y=sin(2x) and its deriv ativ et=0:pi/100:4*pi; y1=sin(2*t); y2=2*cos(2*t); plot(t,y1,'r--'); hold on ;plot(t,y2,'-b'); grid onaxis([0 4*pi -2 2])title('Plot of y=sin(2x) and its derivative')2468101214Plot of y=sin(2x)xyPlot of y=sin(2x) and its deriv ativ exyt=0:pi/100:4*pi; y1=sin(2*t); y2=2*cos(2*t); plot(t,y1,'r--');title('Plot of y=sin(2x)'); xlabel('x'),ylabel('y'); figure(2) plot(t,y2,'-b');title('Plot of y=sin(2x) and its derivative') xlabel('x'),ylabel('y'); grid onaxis([0 4*pi -2 2])2、思考三维曲线（plot3）与曲面(mesh, surf)的用法，（1）绘制参数方程233,)3cos(,)3sin()(t z e t t y e t t t x t t ===--的三维曲线；t=0:pi/30:10*pi;plot3(t.^3.*sin(3.*t).*exp(-t),t.^3.*cos(3.*t).*exp(-t),t.^2);2（2）绘制二元函数xyy xe x x y xf z ----==22)2(),(2，在XOY 平面内选择一个区域(-3:0.1:3,-2:0.1:2)，然后绘制出其三维表面图形。
上机实验报告
实验项目名称Simulink熟悉及其应用所属课程名称系统仿真与MATLAB 实验日期
班级
学号
姓名
成绩
斜坡输入模型
斜坡响应曲线
Subsystem（PID Controller）内容：
输出指令：
plot(ScopeData(:,1),ScopeData(:,2),'-r') xlabel('Time')
ylabel('Output')
title('阶跃响应输出曲线')
grid on
的数值解并绘制函数
，其初始值为：
微分方程的数值解为：x=0
Workspace输出图像：
操作代码：
plot(ScopeData_exp5.time,ScopeData_exp5.signals.values,'-r') xlabel('Time')
ylabel('Output')
实验五阶跃响应输出曲线')
grid on
输出图像：
仿真输出结果：
y(t)的图像如下所示，每条曲线对应不同的信号（F(t)）来源。
说明：
1、将每一道题的程序、建立的模型放置在该题目下方；
2、小结部分为对本次实验的心得体会、思考和建议。
《自动控制原理》MATLAB分析与设计仿真实验任务书（2010）一．仿真实验内容及要求：1．MATLAB软件要求学生通过课余时间自学掌握MATLAB软件的基本数值运算、基本符号运算、基本程序设计方法及常用的图形命令操作；熟悉MATLAB 仿真集成环境Simulink的使用。
2．各章节实验内容及要求1）第三章 线性系统的时域分析法对教材P136.3-5系统进行动态性能仿真，并与忽略闭环零点的系统动态性能进行比较，分析仿真结果；对教材P136.3-9系统的动态性能及稳态性能通过的仿真进行分析，说明不同控制器的作用；在MATLAB环境下完成英文讲义P153.E3.3。
对英文讲义中的循序渐进实例“Disk Drive Read System”，在时，试采用微分反馈使系统的性能满足给定的设计指标。
2）第四章 线性系统的根轨迹法在MATLAB环境下完成英文讲义P157.E4.5；利用MATLAB绘制教材P181.4-5-(3）；在MATLAB环境下选择完成教材第四章习题4-10或4-18，并对结果进行分析。
3）第五章 线性系统的频域分析法利用MATLAB绘制本章作业中任意2个习题的频域特性曲线；4）第六章 线性系统的校正利用MATLAB选择设计本章作业中至少2个习题的控制器，并利用系统的单位阶跃响应说明所设计控制器的功能。
5）第七章 线性离散系统的分析与校正利用MATLAB完成教材P383.7-20的最小拍系统设计及验证。
利用MATLAB完成教材P385.7-25的控制器的设计及验证。
二．仿真实验时间安排及相关事宜1．依据课程教学大纲要求，仿真实验共6学时，教师可随课程进度安排上机时间，学生须在实验之前做好相应的准备，以确保在有限的机时内完成仿真实验要求的内容；2．实验完成后按规定完成相关的仿真实验报告；3．仿真实验报告请参照有关样本制作并打印装订；4．仿真实验报告必须在本学期第15学周结束之前上交授课教师。
MATLAB与控制系统仿真实验指导书河北大学电子信息工程学院20010年9月目录实验一MATLAB基本操作与基本运算 (1)实验二M文件及数值积分仿真方法设计 (3)实验三MATLAB 的图形绘制 (4)实验四函数文件设计和控制系统模型的描述 (6)实验五控制系统的分析与设计 (7)实验六连续系统离散化仿真方法设计 (8)实验七SIMULINK 仿真 (9)实验八SIMULINK 应用进阶 (10)附录MATLAB常用函数 (12)实验一MATLAB基本操作与基本运算一、实验目的及要求：1．熟悉MATLAB6.5的开发环境；2．掌握MATLAB6.5的一些常用命令；3．掌握矩阵、变量、表达式的输入方法及各种基本运算。
二、实验内容：1.熟悉MATLAB6的开发环境：①MATLAB的各种窗口：命令窗口、命令历史窗口、启动平台窗口、工作空间窗口、当前路径窗口。
图1 MA TLAB界面窗口②路径的设置：建立自己的文件夹，加入到MA TLAB路径中，并保存。
设置当前路径，以方便文件管理。
③改变命令窗口数据的显示格式>> format short>> format long然后键入特殊变量：pi （圆周率），比较显示结果。
2.掌握MATLAB 常用命令>> who %列出工作空间中变量>> whos %列出工作空间中变量，同时包括变量详细信息 >>save test %将工作空间中变量存储到test.mat 文件中 >>load test %从test.mat 文件中读取变量到工作空间中 >>clear %清除工作空间中变量>>help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 >>lookfor %查找具有某种功能的函数但却不知道该函数的准确名称如： lookfor Lyapunov 可列出与Lyapunov 有关的所有函数。
实验一MATLAB基本操作一、实验目的：①通过上机实验操作，使学生熟悉MATLAB实验环境，练习MATLAB命令、m文件，进行矩阵运算、图形绘制、数据处理。
②通过上机操作，使得学生掌握Matlab变量的定义和特殊变量的含义，理解矩阵运算和数组运算的定义和规则。
③通过上机操作，使得学生掌握数据和函数的可视化，以及二维曲线、三维曲线、三维曲面的各种绘图指令。
二、实验原理和说明Matlab是Matrix 和Laboratory两词的缩写，是美国Mathworks公司推出的用于科学计算和图形处理的可编程软件，经历了基于DOS版和Windows版两个发展阶段。
三、实验设备和仪器：PC电脑，Matlab7.0仿真软件四、实验内容、方法和步骤：数组运算和矩阵运算数组“除、乘方、转置”运算符前的“.”决不能省略，否则将按矩阵运算规则进行运算；执行数组和数组之间的运算时，参和运算的数组必须同维，运算所得的结果也和参和运算的数组同维。
A=[ 1 2 3; 4 5 6; 7 8 9];B=[-1 -2 -3；-4 -5 -6；-7 -8 -9];X=A.*BY=A*Bplot用于二维曲线绘图，若格式为plot（X，Y，’s’），其中X为列向量，Y是和X等行的矩阵时，以X 为横坐标，按Y的列数绘制多条曲线；若X为矩阵，Y是向量时，以Y为纵坐标按X的列数（或行数）绘制多条曲线。
参考程序如下：t=(0:pi/100:pi)'y1=sin(t)*[-1 1];y2=sin(t).*sin(9*t);plot(t,y1, 'r:', t, y2, 'b-.')axis([0 pi, -1, 1])title('Drawn by Dong-yuan GE')程序运行界面如下：plot3用于三维曲线绘制，其使用格式和plot十分相似。
参考程序如下：t=0:0.02:2*pi;x=sin(t);y=cos(t);z=cos(2*t);plot3(x,y,z,'b-', x,y,z,'o')程序运行界面如下：mesh和surf用于三维空间网线和曲面的绘制。
实验1 数据和函数的可视化实验要求 ⒈掌握二维、三维绘图指令；⒉掌握图形标识和精细指令； ⒊学习用编程器编程、程序调试、程序保存和运行程序进行仿真 实验内容⑴在例3—3的程序的基础上，加入坐标控制、分格线和图形标识等指令，并运行； ⑵在例3—5的程序的基础上，采用精细指令，并加入必要的文本，运行；⑶绘制如下函数的立体网状图和立体曲面图；)(22y x xez +-=参考程序 %dfx=linspace(-2,2,25) ;%在X 轴上取25个点 y=linspace(-2,2,25) ; %在Y 轴上取25个点[xx,yy]=meshgrid(x,y); %xx 和yy 都是25×25的矩阵 zz=xx.*exp(-xx.^2-yy.^2);%zz 也是25×25的矩阵 figure('name','立体网状图',... 'numbertitle','off');mesh(xx,yy,zz); %绘制立体网状图 title('立体网状图');xlabel('X');ylabel('Y');zlabel('Z'); figure('name','立体曲面图',... 'numbertitle','off'); surf(xx,yy,zz); % title('立体曲面图');xlabel('X');ylabel('Y');zlabel('Z');⑷绘制peaks 函数的各种效果图用编程器编辑如下程序，加入注释，分析运行结果； %三维空间立体绘图基本指令演示 figure(1)peaks; %直接绘制peaks 函数[x,y,z]=peaks; %对peaks 函数取点 figure(2)meshc(x,y,z); %同时绘出网状图和等高线axis([-inf,inf,-inf,inf,-inf,inf])title('同时绘出网状图和等高线')figure(3)contour3(peaks,25); %绘出曲面在三维空间中的等高线axis([-inf,inf,-inf,inf,-inf,inf])title('曲面在三维空间中的等高线')figure(4)contour(peaks,25); %绘制等高线在X—Y平面的投影axis([-inf,inf,-inf,inf])title('等高线在X—Y平面的投影')figure(5)meshz(x,y,z); %将曲面加上围裙axis([-inf,inf,-inf,inf,-inf,inf])title('曲面加上围裙')figure(6)waterfall(x,y,z); %在X方向加上水流效果axis([-inf,inf,-inf,inf,-inf,inf])title('在X方向加上水流效果')figure(7)waterfall(x',y',z'); %在Y方向加上水流效果axis([-inf,inf,-inf,inf,-inf,inf])title('在Y方向加上水流效果')实验报告要求⑴按自己的设想编制程序，给出每句程序的注释；⑵附上仿真的结果；实验2 动态系统的时域仿真实验要求⒈建立动态系统的数学模型，掌握数学模型的MATLAB语言表示和模型转换；⒉掌握时域仿真的方法和主要的仿真指令；⒊学习用编程器编程、程序调试、程序保存和运行程序进行仿真。
实验一MATLAB基本操作与基本运算 (1)实验二M文件及数值积分仿真方法设计 (3)实验三MATLAB 的图形绘制 (4)实验四函数文件设计和控制系统模型的描述 (6)实验五控制系统的分析与设计 (7)实验六连续系统离散化仿真方法设计 (8)实验七SIMULINK 仿真 (9)实验八SIMULINK 应用进阶 (10)附录MATLAB常用函数 (11)实验一MATLAB基本操作与基本运算一、实验目的及要求：1．熟悉MATLAB6.5的开发环境；2．掌握MATLAB6.5的一些常用命令；3．掌握矩阵、变量、表达式的输入方法及各种基本运算。
二、实验内容：1.熟悉MATLAB6的开发环境：①MATLAB的各种窗口：命令窗口、命令历史窗口、启动平台窗口、工作空间窗口、当前路径窗口。
图1 MA TLAB 界面窗口②路径的设置：建立自己的文件夹，加入到MA TLAB 路径中，并保存。
设置当前路径，以方便文件管理。
③改变命令窗口数据的显示格式>> format short >> format long然后键入特殊变量：pi （圆周率），比较显示结果。
2.掌握MATLAB 常用命令>> who %列出工作空间中变量>> whos %列出工作空间中变量，同时包括变量详细信息 >>save test %将工作空间中变量存储到test.mat 文件中 >>load test %从test.mat 文件中读取变量到工作空间中 >>clear %清除工作空间中变量>>help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 >>lookfor %查找具有某种功能的函数但却不知道该函数的准确名称如： lookfor Lyapunov 可列出与Lyapunov 有关的所有函数。
3. 在MATLAB 的命令窗口计算:1) )2sin(π2) 5.4)4.05589(÷⨯+4. 设计M 文件计算: 已知求出： 1) a+ba-b a+b*5 a-b+I （单位阵）2) a*b a.*b a/b 3) a^2 a.^2注意：点运算 . 的功能，比较结果。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=300120101b))(,(12h k x y h x f k m m ++=)(2)(211k k h x y y mm ++=+⎪⎩⎪⎨⎧))(,(1m m x y x f k =5. 设计M 文件计算:x=0:0.1:10当sum>1000时停止运算，并显示求和结果及计算次数。
三、实验报告要求：1）体会1、2的用法；2）对3、4、5写出程序及上机的结果。
实验二 M 文件及数值积分仿真方法设计一、实验目的及要求：1. 掌握MA TLAB 程序设计方法，会编写M 文件；2. 掌握用数值积分方法对连续系统建立仿真模型的方法。
二、实验内容：M 文件设计：仿真时间0～1秒，计算步长h=0.1，初始值y(0)=1，对连续系统y x dxdy+= 采用以下方法建立仿真模型： ① 前向Euler: n n n h f y y +=+1② 后向Euler: ③ 梯形法： ④ 改进Euler:⑤ 经典RK4：)2(1002i i i x x sum -∑==11+++=n n n hf y y ()112++++=n n n n f f hy y ))(,(1m m x y x f k =)2,2(12k hy h x f k m m ++=()4321122k k k k hy y m m ++++=+⎪⎪⎪⎪⎪⎪⎨⎧),(34h k y h x f k m m ++=)2,2(23k h y h x f k m m ++=⑥ Adams 法： ⑦ Milne 法：并与精确解x e x y x --=12)(的结果进行比较。
三、实验报告要求：写出程序及上机的结果（保存程序，以便下次实验使用）。
实验三 MATLAB 的图形绘制一、实验目的及要求：1.掌握MA TLAB 绘图的基本方法，熟悉各种绘图函数的使用；2.掌握图形的修饰方法和标注方法；3.了解MA TLAB 中图形窗口的操作。
二、实验内容：1、x=［－2π，2π］，y1=sinx 、y2=cosx 、y3=sin2x 、y4=cos 2x①用MA TLAB 语言分四个区域分别绘制的曲线，并且对图形标题及横纵坐标轴进行标注（如下图所示）。
()32115162312----+-+=n n n n n f f f hy y ()111143-+-++++=m m m m m f f f hyy图2 四分区绘制曲线②另建一个窗口，不分区，用不同颜色、线型绘出四条曲线，并标注图例注解。
图2 四分区绘制曲线图3 同一窗口绘制多条曲线③(选做)在图形窗口可利用Figure窗口菜单提供的功能进行操作，并保存成** .fig文件。
2、在实验二的基础上分别用分区域绘制和同一窗口绘制多条曲线的方法绘制7种算法和精确解的结果。
图4 仿真算法与精确解曲线分区域比较图5 仿真算法与精确解曲线同一窗口比较三、实验报告要求：写出1、①（对图形的修饰写出第一个子区域即可）和②的程序。
实验四 函数文件设计一、实验目的及要求：1. 掌握函数文件的设计方法；2. 掌握函数文件的调用方法。
二、实验内容：1、阶乘函数的设计： 设计一个函数文件实现一个阶乘运算 n n n n y ⨯-⨯⨯⨯⨯==)1(321!)( ， 并设计程序调用该函数。
为保证函数的通用性，当输入负数或小数时，显示出错提示：disp('Input parameter must be a positive integer!') 提示：fix(x) 对零方向取整数 ceil(x) 对+∞方向取整数round(x) 四舍五入取整数2、非线性限幅函数的设计：为，a a x a a x a a x x y ⎪⎩⎪⎨⎧-<->≤=正整数并要求：对输入的一个数列]2,2[a a x -∈，以x 为横坐标，y 为纵坐标绘制红色实线，图形的x 轴范围限定在[x 最小值-1 x 最大值+1]，y 轴范围限定在[-1.5*a 1.5*a]，如下左图：（右图为利用Figure 窗口的功能进行添加，选做）三、实验报告要求：写出程序及上机的结果。
实验五 控制系统的分析与设计一、实验目的及要求：1. 掌握控制系统数学模型的基本描述方法。
2. 掌握控制系统频域与时域分析基本方法。
3. 掌握现代控制理论的基本设计方法。
二、实验内容：1、已知两个传递函数分别为：ss x G s x G +=+=22132)(,131)(①在MATLAB 中分别用传递函数、零极点、和状态空间法表示； ②在MATLAB 中分别求出通过反馈、串联、并联后得到的系统模型；2．已知两个单位负反馈系统中前向通道的传递函数分别为：s s s s x G s s s x G +++=+++=23221312)(,1312)(，①分别绘制开环系统的bode 图和nyquist 曲线； ②分别绘制闭环系统的根轨迹；③分别绘制闭环系统的阶跃响应曲线，分析静态误差； ④分别绘制闭环系统的单位斜坡输入响应曲线，分析静态误差。
3．已知某被控对象的状态空间表达式为：()xy u x x101111101101121=⎥⎥⎦⎤⎢⎢⎣⎡+⎥⎥⎦⎤⎢⎢⎣⎡---=①判断能控性；②期望极点为p=(-1,-2,-3)； ③设计状态反馈K 进行极点配置；④分别绘制校正前后系统的阶跃响应曲线。
三、实验报告要求：写出程序及上机的结果，并对结果进行分析。
实验六 连续系统离散化仿真方法设计一、实验目的及要求：1. 掌握用替换法建立连续系统仿真模型的方法；2. 掌握用离散相似法建立连续系统仿真模型的方法。
二、实验内容：1、替换法建立连续系统仿真模型:①设连续系统传递函数为 351)()()(2++==S S s u s y s G ，用简单替换法求其Z 传递函数和差分方程（采样周期T=0.1s ）;② 若已知：输入作用为单位阶跃函数U(t)=1(t)，起始值y(0)=0、y(0.1)= 0.004251，仿真时间t=[0,10],在①的基础上设计程序求出响应结果Y(t)，并绘制响应曲线；并与控制系统step( )函数得出的响应结果比较。
2、离散相似法建立连续系统仿真模型:①若系统传递函数为52)()()(+==S s u s y s G ，当采用零阶保持器s e s G Tsh --=1)(，采样周期T=0.1S 时，用Z 域离散相似法求出系统的脉冲传递函数G(z)，并写出差分方程；② 若输入作用为斜坡函数U(t)=t ，仿真时间t=[0,10]，在①的基础上设计程序求出响应结果Y(t)，并绘制响应曲线；并与控制系统lsim( )函数得出的响应结果比较。
三、实验报告要求：写出用简单替换法和离散相似法建立差分方程的过程，并写出程序。
实验七 SIMULINK 仿真——单回路控制系统及PID 控制器参数整定一、实验目的及要求：1.熟悉SIMULINK 工作环境及特点；2.熟悉控制线性系统仿真常用基本模块的用法；3.掌握SIMULINK 的建模与仿真方法。
二、实验内容：用SIMULINK 建立被控对象的传递函数为11010)(23+++=s s s x G ，系统输入为单位阶跃，采用PID 控制器进行闭环调节。
①练习模块、连线的操作，并将仿真时间定为300 秒，其余用缺省值； ②试用稳定边界法（过程控制P5工程整定法之一）设置出合适的PID 参数，得出满意的响应曲线。
③设计M 文件在一个窗口中绘制出系统输入和输出的曲线，并加图解。
三、实验报告要求：①阐述用SIMULINK 进行控制系统仿真的一般过程； ②说明用工程整定法——稳定边界法整定PID 参数的过程。
实验八 SIMULINK 应用进阶 ——子系统封装及复杂系统仿真一、实验目的及要求：1. 掌握SIMULINK 子系统模块化和封装技术；2. 熟悉SIMULINK 对复杂控制系统仿真时常用的方法。
二、实验内容：1、用SIMULINK 建立系统：c bx ax y++=3，x 为输入，y 为输出，c b a ,,为常数。
对该系统进行封装，要求通过对话框能修改c b a ,,的值。
图6 系统封装后的效果图2、在第一题的基础上，①加入实验四建立的阶乘函数，实现ax^3+bx+c+n!的运算；②将模型中的a 、b 、c 、n 表示为字母形式的变量名，设计M 文件在文件中为模型参数赋值，然后运行仿真模型；③在②中的M 文件中实现用命令行运行仿真模型，并用plot 语句绘曲线。
