MATLAB与系统仿真实验手册2010

实验一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 语句绘曲线。