《MATLAB语言与应用》实验测验1 专业学号姓名成绩 自动化200909150151 高成 要求：随机生成一个4×4矩阵A，生成一个4×4的魔术矩阵B，计算C=A.*B，D=A*B，E=A/B，F=A\B，G=A.^2，H=A^2。 矩阵C、D、E、F、G和H的实验结果分别为： >> A=rand(4) A = 0.2769 0.6948 0.4387 0.1869 0.0462 0.3171 0.3816 0.4898 0.0971 0.9502 0.7655 0.4456 0.8235 0.0344 0.7952 0.6463 >> B=magic(4) B = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 1 5 1 >> C=A.*B C = 4.4308 1.3897 1.3162 2.4293 0.2309 3.4881 3.8156 3.9181 0.8742 6.6516 4.5931 5.3470 3.2938 0.4822 11.9280 0.6463 >> D=A*B D = 12.6011 13.8844 13.2146 14.6104 7.7173 13.1080 12.9453 8.2055 14.9772 22.2435 21.0705 18.4963 23.0896 16.6406 17.2807 21.1692
>> E=A/B Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.306145e-017. E = 1.0e+014 * 0.8053 2.4159 -2.4159 -0.8053 -0.5976 -1.7929 1.7929 0.5976 0.1930 0.5789 -0.5789 -0.1930 -1.9751 -5.9254 5.9254 1.9751 >> F=A\B F = 32.4732 0.3697 5.9411 15.7588 34.5598 -1.5535 0.4005 28.6980 -48.3016 -4.7808 -8.8210 -36.1811 22.4018 27.1553 26.4707 24.4556 >> G=A.^2 G = 0.0767 0.4828 0.1925 0.0349 0.0021 0.1006 0.1456 0.2399 0.0094 0.9029 0.5860 0.1985 0.6781 0.0012 0.6323 0.4177 >> H=A^2 H = 0.3053 0.8361 0.8711 0.7083 0.4678 0.5121 0.8228 0.6505 0.5120 1.1116 1.3455 1.1126 0.8391 1.3610 1.4971 0.9428
图5. 利用 SIMULINK仿
4. 建立如图11-54所示的仿真模型，其中PID控 制器采用Simulink子系统封装形式，其内部 结构如图11-31（a)所示。试设置正弦波信号 幅值为5、偏差为0、频率为10πHz\始终相位 为0，PID控制器的参数为Kp=10.75、 Ki=1.2、Kd=5，采用变步长的ode23t算法、 仿真时间为2s，对模型进行仿真。 （6)观察仿真结果。系统放着结束后，双击仿真模型中的示波器模块，得到仿真结果。单击示波器窗口工具栏上的Autoscale按钮，可以自动调整坐标来 使波形刚好完整显示，这时的波形如图所示。 图3 2. 题操作步骤如下： (1） 打开一个模型编辑窗口。 (2） 将所需模块添加到模型中。在模块库浏览器中单击Sources,将 Clock（时钟）拖到模型编辑窗口。同样，在User-Defined Functions（用户定义模块库）中把Fcn(函数模块）拖到模型编辑窗口，在Continuous（连续系统模块库）中把 Integrator（积分模块）拖到模型编辑窗口，在Sinks中把Display模块编辑窗口。 (3） 设置模块参数并连接各个模块组成仿真模型。双击Fcn模块，打开Function Block operations中把Add模块拖到模型编辑窗口，在Sinks中把Scope模块拖到模型编辑窗口。 （3） 设置模块参数并连接各个模块组成仿真模型。先双击各个正弦源，打开其Block Parameters对话框，分别设置Frequency（频率）为2*pi、 6*pi、10*pi、 14*pi、18*pi，设置Amplitude（幅值）为1、1/3、1/5、1/7和1/9，其余参数不改变。对于求和模块，將符号列表List of signs设置为 +++++。 （4） 设置系统仿真参数。单击模型
系统仿真 § 4.1控制系统的数学模型 1、传递函数模型（tranfer function） 2、零极点增益模型（zero-pole-gain） 3、状态空间模型（state-space） 4、动态结构图（Simulink结构图） 一、传递函数模型(transfer fcn-----tf) 1、传递函数模型的形式 传函定义：在零初始条件下，系统输出量的拉氏变换C（S）与输入量的拉氏变换R（S）之比。 C(S) b1S m+b2S m-1+…+b m G（S）=----------- =- -------------------------------- R(S) a1S n + a2S n-1 +…+ a n num(S) = ------------ den(S) 2、在MATLAB命令中的输入形式 在MATLAB环境中，可直接用分子分母多项式系数构成的两个向量num、den表示系统： num = [b1, b2, ..., b m]; den = [a1, a2, ..., a n]; 注：1）将系统的分子分母多项式的系数按降幂的方式以向量的形式输入两个变量，中间缺项的用0补齐，不能遗漏。 2）num、den是任意两个变量名，用户可以用其他任意的变量名来输入系数向量。 3）当系统种含有几个传函时，输入MATLAB命令状态下可用n1,d1;n2,d2…….。 4）给变量num，den赋值时用的是方括号；方括号内每个系数分隔开用空格或逗号；num,den方括号间用的是分号。 3、函数命令tf( ) 在MATLAB中，用函数命令tf( )来建立控制系统的传函模型，或者将零极点增益模型、状态空间模型转换为传函模型。 tf( )函数命令的调用格式为： 圆括号中的逗号不能用空格来代替 sys = tf ( num, den ) [G= tf ( num, den )]
课程设计报告 题目PID控制器应用 课程名称控制系统仿真院部名称龙蟠学院 专业自动化 班级M10自动化 学生姓名 学号 课程设计地点 C208 课程设计学时一周 指导教师应明峰 金陵科技学院教务处制成绩
一、课程设计应达到的目的 应用所学的自动控制基本知识与工程设计方法，结合生产实际，确定系统的性能指标与实现方案，进行控制系统的初步设计。 应用计算机仿真技术，通过在MATLAB软件上建立控制系统的数学模型，对控制系统进行性能仿真研究，掌握系统参数对系统性能的影响。 二、课程设计题目及要求 1．单回路控制系统的设计及仿真。 2．串级控制系统的设计及仿真。 3．反馈前馈控制系统的设计及仿真。 4．采用Smith 补偿器克服纯滞后的控制系统的设计及仿真。 三、课程设计的内容与步骤 (1)．单回路控制系统的设计及仿真。 （a）已知被控对象传函W(s) = 1 / (s2 +20s + 1)。 （b）画出单回路控制系统的方框图。 （c）用MatLab的Simulink画出该系统。
（d）选PID调节器的参数使系统的控制性能较好，并画出相应的单位阶约响应曲线。注明所用PID调节器公式。PID调节器公式Wc（s）=50(5s+1)/(3s+1) 给定值为单位阶跃响应幅值为3。 有积分作用单回路控制系统PID控制器取参数分别为：50 2 5 有积分作用单回路控制系统PID控制器取参数分别为：50 0 5
大比例作用单回路控制系统PID控制器取参数分别为：50 0 0 （e）修改调节器的参数，观察系统的稳定性或单位阶约响应曲线，理解控制器参数对系统的稳定性及控制性能的影响？ 答：由上图分别可以看出无积分作用和大比例积分作用下的系数响应曲线，这两个PID调节的响应曲线均不如前面的理想。增大比例系数将加快系统的响应，但是过大的比例系数会使系统有比较大的超调，并产生振荡，使稳定性变坏；
控制系统的MATLAB仿真与设计课后答案
>>z=-4*sqrt(2)*sin(t); >>plot3(x,y,z,'p'); >>title('Line in 3-D Space'); >>text(0,0,0,'origin'); >>xlabel('X'),ylable('Y'),zlable('Z');grid; 4>>theta=0:0.01:2*pi; >>rho=sin(2*theta).*cos(2*theta); >>polar(theta,rho,'k'); 5>>[x,y,z]=sphere(20); >>z1=z; >>z1(:,1:4)=NaN; >>c1=ones(size(z1)); >>surf(3*x,3*y,3*z1,c1); >>hold on >>z2=z; >>c2=2*ones(size(z2)); >>c2(:,1:4)=3*ones(size(c2(:,1:4))); >>surf(1.5*x,1.5*y,1.5*z2,c2); >>colormap([0,1,0;0.5,0,0;1,0,0]); >>grid on >>hold off 第四章 1>>for m=100:999 m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2M文件：function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; 主程序： [s,p]=fcircle(10) 3>>y=0;n=100; for i=1:n y=y+1/i/i; end >>y
第 4 章 MatLab 的程序设计 MatLab 是一个工具、开发平台，同时它也是一门编程语言。与在命令窗口用交互的方式工作相比，通过程序运行来解决实际问题，其效率更高，因此，凡是复杂的、大型的应用都是以程序的方式执行。相对其它高级语言， MatLab 更简单、编程的效率更高、调试过程也更容易。 MatLab 中的程序文件是以 m 为后缀，所以通常将 MatLab 的程序文件称为 m 文件。MatLab提供了两种形式的m文件，即：脚本（Script）式m文件（就简称m文件）、函数型 m 文件。在 MatLab 中已经嵌入了一个功能强大的集成开发环境—— m 文件编辑器，用它来进行程序的编辑、修改、调试、运行等，完成应用开发工作。 4.1 MatLab 程序设计基础 通过前面内容的学习，大家对 MatLab 已经有了一个初步的认识和印象，到目前为止，我们都是在“命令”窗口中，以交互的方式运行，完成我们的工作。实际上简单的m 文件，就是一个批处理程序，它是若干条命令的集合。 例： 4.1.1 M 文件规则和属性 函数 M 文件必须遵循一些特定的规则。除此之外，它们有许多的重要属性，这其中包括： 1. 函数名和文件名必须相同。例如，函数 fliplr 存储在名为 fliplr.m 文件中。 2. MATLAB 头一次执行一函数个 M 文件时，它打开相应的文本文件并将命令编辑成存储器的内部表示，以加速执行以后所有的调用。如果函数包含了对其它函 数 M 文件的引用，它们也同样被编译到存储器。普通的脚本 M 文件不被编译，即使它们是从函数 M 文件内调用；打开脚本 M 文件，调用一次就逐行进行注释。 3. 在函数 M 文件中，到第一个非注释行为止的注释行是帮助文本。当需要帮助时，返回该文本。例如， ? help fliplr 返回上述前八行注释。 4. 第一行帮助行，名为 H1 行，是由 lookfor 命令搜索的行。 5. 函数可以有零个或更多个输入参量。函数可以有零个或更多个输出参量。
一、 控制系统的模型与转换 1． 请将下面的传递函数模型输入到matlab 环境。 ]52)1)[(2(24)(322 33++++++=s s s s s s s G ) 99.02.0)(1(568 .0)(22+--+=z z z z z H ，T=0.1s >> s=tf('s'); G=(s^3+4*s+2)/(s^3*(s^2+2)*((s^2+1)^3+2*s+5)); G Transfer function: s^3 + 4 s + 2 ------------------------------------------------------ s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3 >> num=[1 0 0.56]; den=conv([1 -1],[1 -0.2 0.99]); H=tf(num,den,'Ts',0.1) Transfer function: z^2 + 0.56 ----------------------------- z^3 - 1.2 z^2 + 1.19 z - 0.99 2． 请将下面的零极点模型输入到matlab 环境。请求出上述模型的零极点，并绘制其位置。 )1)(6)(5()1)(1(8)(22 +++-+++=s s s s j s j s s G ) 2.8() 6.2)(2.3()(1 511-++=----z z z z z H ，T=0.05s >>z=[-1-j -1+j]; p=[0 0 -5 -6 -j j]; G=zpk(z,p,8) Zero/pole/gain: 8 (s^2 + 2s + 2) -------------------------- s^2 (s+5) (s+6) (s^2 + 1) >>pzmap(G)
2013年春季学期 《MATLAB语言及应用》课程试卷 姓名： 学号： 学院： 专业： 必答题 1.常用的matlab界面由哪些窗口组成，各有什么主要作用？(4分) (1)菜单和工具栏功能：【File】菜单主要用于对文件的处理。【Edit】菜单主 要用于复制、粘贴等操作，与一般Windows程序的类似，在此不作详细介绍。【Debug】菜单用于调试程序。【Desktop】菜单用于设置主窗口中需要打开的窗口。【Window】菜单列出当前所有打开的窗口。【Help】菜单用于选择打开不同的帮助系统。 (2)命令窗口功能：用于输入命令并显示除图形以外的所有执行结果 (3)历史命令窗口功能：主要用于记录所有执行过的命令 (4)当前工作目录窗口功能：对文件和目录进行操作 (5)工作空间窗口功能：查看、载入和保存变量 2.如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么 区别？(2分) 方法一：在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/|【SetPath】菜单，进入“设置搜索路径”对话框，通过该对话框编辑搜索路径。 方法二：在命令窗口执行“path(path,…D:\Study ?)”，然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。 方法三：在命令窗口执行“addpath D:\Study- end”，将新的目录加到整个搜索路径的末尾。如果将end改为begin，可以将新的目录加到整个搜索路径的开始。 区别：当前文件目录是正在运行的文件的目录，显示文件及文件夹的详细信息，且只有将文件设置为当前目录才能直接调用。搜索路径中的文件可以来自多个 不同目录，在调用时不用将其都设置为当前目录，为同时调用多个文件提供 方便。 3.有几种建立矩阵的方法？各有什么优点？(4分) 1.在命令窗口中直接输入优点：适合输入不规则和较小的矩阵 2.通过语句和函数生成矩阵优点：适合输入规则且较大的矩阵 3.通过M文件来建立矩阵，从外部数据文件中导入矩阵 优点：方便创建和导入大型矩阵 4.说明break语句、continue语句和return语句的用法。(3分)
《MATLAB与控制系统仿真》 实验报告 班级： 学号： 姓名： 时间：2013 年 6 月
目录实验一 MATLAB环境的熟悉与基本运算（一）实验二 MATLAB环境的熟悉与基本运算（二）实验三 MATLAB语言的程序设计 实验四 MATLAB的图形绘制 实验五基于SIMULINK的系统仿真 实验六控制系统的频域与时域分析 实验七控制系统PID校正器设计法 实验八线性方程组求解及函数求极值
实验一 MATLAB环境的熟悉与基本运算（一） 一、实验目的 1．熟悉MATLAB开发环境 2．掌握矩阵、变量、表达式的各种基本运算 二、实验基本原理 1.熟悉MATLAB环境: MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。 2.掌握MA TLAB常用命令 表1 MA TLAB常用命令 3.MATLAB变量与运算符 3．1变量命名规则 3．2 MATLAB的各种常用运算符 表3 MATLAB关系运算符 表4 MATLAB逻辑运算符
| Or 逻辑或 ~ Not 逻辑非 Xor 逻辑异或 符号功能说明示例符号功能说明示例 ：1:1:4;1:2:11 . ；分隔行.. ，分隔列… （）% 注释 [] 构成向量、矩阵！调用操作系统命令 {} 构成单元数组= 用于赋值 4.MATLAB的一维、二维数组的寻访 表6 子数组访问与赋值常用的相关指令格式 三、主要仪器设备及耗材 计算机 四．实验程序及结果 1、新建一个文件夹（自己的名字命名，在机器的最后一个盘符） 2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。 3、学习使用help命令。
MATLAB就是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模与仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。 MATLAB与Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数与数据、实现算法、创建用户界面、连 matlab开发工作界面 接其她编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位就是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
在新的版本中也加入了对C,FORTRAN,C++ ,JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用。 一种语言之所以能如此迅速地普及,显示出如此旺盛的生命力,就是由于它有着不同于其她语言的特点。正如同FORTRAN与C等高级语言使人们摆脱了需要直接对计算机硬件资源进行操作一样,被称作为第四代计算机语言的MATLAB,利用其丰富的函数资源,使编程人员从繁琐的程序代码中解放出来。MATLAB的最突出的特点就就是简洁。MATLAB用更直观的、符合人们思维习惯的代码,代替了C与FORTRAN语言的冗长代码。MATLAB给用户带来的就是最直观、最简洁的程序开发环境。以下简单介绍一下MATLAB的主要特点。 ①语言简洁紧凑,使用方便灵活,库函数极其丰富。MATLAB程序书写形式自由,利用其丰富的库函数避开繁杂的子程序编程任务,压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写,用户不必担心函数的可靠性。可以说,用MATLAB进行科技开发就是站在专家的肩膀上。 具有FORTRAN与C等高级计算机语言知识的读者可能已经注意到,如果用FORTRAN或C语言去编写程序,尤其当涉及矩阵运算与画图时,编程会很麻烦。例如,如果用户想求解一个线性代数方程,就得编写一个程序块读入数据,然后再使用一种求解线性方程的算法(例如追赶法)编写一个程序块来求解方程,最后再输出计算结果。在求解过程中,
MATLAB控制系统与仿真 课 程 设 计 报 告 院（系）：电气与控制工程学院 专业班级：测控技术与仪器1301班 姓名：吴凯 学号：1306070127
指导教师：杨洁昝宏洋 基于MATLAB的PID恒温控制器 本论文以温度控制系统为研究对象设计一个PID控制器。PID控制是迄今为止最通用的控制方法，大多数反馈回路用该方法或其较小的变形来控制。PID控制器（亦称调节器）及其改进型因此成为工业过程控制中最常见的控制器(至今在全世界过程控制中用的84%仍是纯PID调节器，若改进型包含在内则超过90%)。在PID控制器的设计中，参数整定是最为重要的,随着计算机技术的迅速发展，对PID参数的整定大多借助于一些先进的软件，例如目前得到广泛应用的MATLAB仿真系统。本设计就是借助此软件主要运用Relay-feedback法，线上综合法和系统辨识法来研究PID控制器的设计方法，设计一个温控系统的PID控制器，并通过MATLAB中的虚拟示波器观察系统完善后在阶跃信号下的输出波形。 关键词：PID参数整定；PID控制器；MATLAB仿真。 Design of PID Controller based on MATLAB Abstract This paper regards temperature control system as the research object to design a pid controller. Pid control is the most common control method up until now; the great majority feedback loop is controlled by this method or its small deformation. Pid controller (claim regulator also) and its second generation so become the most common controllers in the industry process control (so far, about 84% of the controller being used is the pure pid controller, it’ll exceed 90% if the second generation included). Pid parameter setting is most important in pid controller designing, and with the rapid development of the computer technology, it mostly recurs to some advanced software, for example, mat lab simulation software widely used now. this design is to apply that soft mainly use Relay feedback law and synthetic method on the line to study pid
《MATLAB语言及应用》课程教学大纲 课程编号：21311105 总学时数：32 总学分数：2 课程性质：专业必修课 适用专业：电气工程及其自动化 一、课程的任务和基本要求： 控制算法是集中现代控制系统CAD 技术的本质反映，对于控制算法系统深入地学习，在扩展计算机技术在控制理论中的应用和发展，同时培养学生运用计算机技术进行思维和开发的能力。控制系统理论、计算方法与计算机技术的结合是当代控制理论发展的标志，因此在以MATLAB 为代表的软件平台上，对控制系统进行分析、设计与仿真将成为控制工程领域工程师必须熟练掌握的重要知识和技能。因此，深入透彻地分析和理解控制算法的思想和构造就必须系统学习典型控制系统应用软件的基本原理和控制算法，将成为本课程的目的和任务。 通过本课程的学习，要求学生掌握对于控制系统的分析和综合设计的方法和基本技巧，而控制算法在控制系统CAD 技术中占有相当大的比重，本课程要求较熟练掌握控制算法的基本思想；MATLAB 是一种解释性编程语言，因此，要求熟练掌握MATLAB 的基本编程手段和模块化编程方法，消化和理解控制语言描述的图形界面的设计过程。 二、基本内容和要求： 1.自动控制系统与仿真基础知识 （1）自动控制系统基本概念 （2）自动控制系统分类 （3）控制系统仿真基本概念 （4）MA TLAB与控制系统仿真 （5）MA TLAB 7中控制相关的工具箱 要求：了解自动控制系统与仿真的基础知识，包括自动控制系统的基本概念、分类，以及控制系统仿真的基本概念和Matlab工具。 2.MA TLAB计算及仿真基础 （1）MA TLAB概述 （2）MA TLAB桌面操作环境 （3）MA TLAB数值计算 （4）关系运算和逻辑运算 （5）符号运算 （6）复数和复变函数运算 （7）MA TLAB常用绘图命令 （8）MA TLAB程序设计 要求：了解MA TLAB计算及仿真基础，包括MATLAB的安装、界面及其数值计算、函数运算、程序设计及其绘图命令。 3. Simulink仿真基础 （1）Simulink仿真概述
《自动控制原理实验》 目录 第一部分实验箱的使用 第二部分经典控制实验 第一章基本实验 实验一典型环节及其阶跃响应 实验二二阶系统阶跃响应 实验三控制系统的稳定性分析 实验四控制系统的频率特性 实验五连续控制系统的串联校正 实验六数字PID控制实验 第二章综合实验 第三部现代控制理论实验 第一章基本实验 第二章综合实验
实验一 典型环节及其阶跃响应 预习要求： 1、复习运算放大器的工作原理；了解采用A μ741运算放大器构成各种运算电路的方法； 2、了解比例控制、微分控制、积分控制的物理意义。 一、实验目的 1、学习自动控制系统典型环节的电模拟方法，了解电路参数对环节特性的影响。 2、学习典型环节阶跃响应的测量方法； 3、学会根据阶跃响应曲线计算确定典型环节的传递函数。 二、实验内容 1、比例环节 电路模拟： 图1-1 传递函数： 2211 ()()()U s R G s U s R ==- 2、惯性环节 电路模拟： 图1-2 传递函数： 22112()/()()11 U s R R K G s U s Ts R Cs = =-=- ++ 3、积分环节 电路模拟： A/D1 D/A1 A/D1
图1-3 传递函数： 21()11 ()()U s G s U s Ts RCs = =-=- 4、微分环节 电路模拟： 图1-4 传递函数： 211() ()() U s G s s RC s U s τ= =-=- 5、比例微分 电路模拟： 图1-5 传递函数： 222111 ()()(1)(1)()U s R G s K s R C s U s R τ= =-+=-+ 6、比例积分 电路模拟： 图1-6 A/D1 2 R D/A1 A/D1 A/D1 A/D1 C
一、MATLAB简介 MATLAB 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 MATLAB是矩阵实验室（Matrix Laboratory）的简称，和Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。 二、应用 MATLAB 产品族可以用来进行以下各种工作： （1）数值分析； （2）数值和符号计算； （3）工程与科学绘图； （4）控制系统的设计与仿真； （5）数字图像处理技术； （6）数字信号处理技术； （7）通讯系统设计与仿真； （8）财务与金融工程。 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特
第4章MatLab的程序设计 MatLab是一个工具、开发平台，同时它也是一门编程语言。与在命令窗口用交互的方式工作相比，通过程序运行来解决实际问题，其效率更高，因此，凡是复杂的、大型的应用都是以程序的方式执行。相对其它高级语言，MatLab更简单、编程的效率更高、调试过程也更容易。 MatLab中的程序文件是以m为后缀，所以通常将MatLab的程序文件称为m文件。MatLab提供了两种形式的m文件，即：脚本（Script）式m文件（就简称m文件）、函数型m文件。在MatLab中已经嵌入了一个功能强大的集成开发环境——m文件编辑器，用它来进行程序的编辑、修改、调试、运行等，完成应用开发工作。 4.1 MatLab程序设计基础 通过前面内容的学习，大家对MatLab已经有了一个初步的认识和印象，到目前为止，我们都是在“命令”窗口中，以交互的方式运行，完成我们的工作。实际上简单的m文件，就是一个批处理程序，它是若干条命令的集合。 例： 4.1.1 M文件规则和属性 函数M文件必须遵循一些特定的规则。除此之外，它们有许多的重要属性，这其中包括： 1. 函数名和文件名必须相同。例如，函数fliplr存储在名为fliplr.m文件中。 2. MATLAB头一次执行一函数个M文件时，它打开相应的文本文件并将命令编辑成存储器的内部表示，以加速执行以后所有的调用。如果函数包含了对其它函数M文件的引用，它们也同样被编译到存储器。普通的脚本M文件不被编译，即使它们是从函数M文件内调用；打开脚本M文件，调用一次就逐行进行注释。 页脚内容1
3. 在函数M文件中，到第一个非注释行为止的注释行是帮助文本。当需要帮助时，返回该文本。例如，? help fliplr返回上述前八行注释。 4. 第一行帮助行，名为H1 行，是由lookfor命令搜索的行。 5. 函数可以有零个或更多个输入参量。函数可以有零个或更多个输出参量。 6. 函数可以按少于函数M文件中所规定的输入和输出变量进行调用，但不能用多于函数M文件中所规定的输入和输出变量数目。如果输入和输出变量数目多于函数M文件中function语句一开始所规定的数目，则调用时自动返回一个错误。 相对于函数m文件，脚本式m文件就简单多了，它没有严格的格式要求，只要将有关的命令或函数一一敲入即可，但是还是有几个问题需要注意： 1. m文件的名称不得与MatLab的内部函数同名、第一个字符不得为数字（这点与变量的命名规则相同）； 2. 最好在文件的头部加上注释，对该m文件的作用、功能作一个简要说明，而在一些重要命令行后也加上注释行，以方便使用者阅读、查找； 3. 要特别注意m文件的保存路径或位置，如果不是保存在MatLab默认的路径下，可以使用addpath函数来设置、添加路径，否则，m文件不能运行。 脚本式m文件与函数m文件还有一个重要区别：脚本式m文件中的变量均为全局变量，而函数m文件中的变量则是局部变量。这可以从这两种程序文件运行后在Workspace中留下痕迹看出。当然，在函数m文件中也可以专门将某些变量定义为全局变量（关键字是：global）。不过，在使用全局变量（函数m文件中）时应特别注意： ①.全局变量需要函数体的变量赋值语句之前定义或说明； ②.全局变量名最好是大写，而且要尽量长，能反映它本身的含义； 页脚内容2
5.6 控制系统的时域响应MATLAB 仿真实训 5.6.1实训目的 1. 学会利用MATLAB 绘制系统的单位阶跃响应曲线，掌握读取系统动态性能指标的方法； 2. 学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法； 3. 掌握利用MATLAB 绘制系统的零输入响应曲线的方法； 4. 掌握利用MATLAB 绘制系统的一般输入响应曲线的方法； 5. 学会通过仿真曲线读取相关信息，并依据有关信息进行系统的时域分析。 5.6.2实训内容 1.编写程序求取下列各系统的单位阶跃响应，完成表5-5并记录相关曲线。 162.316)(21++= s s s G 16 4.216 )(22 ++=s s s G 166.116)(2 3++=s s s G 1616 )(24++=s s s G 解：>> n1=16; >> d1=[1,3.2,16]; >> sys1=tf(n1,d1); >> step(sys1) >> n2=16; >> d2=[1,2.4,16]; >> sys2=tf(n2,d2); >> step(sys2)
>> n3=16; >> d3=[1,1.6,16]; >> sys3=tf(n3,d3); >> step(sys3) >> n4=16; >> d4=[1,1,16]; >> sys4=tf(n4,d4); >> step(sys4)
序号ξnω m ax c p t s t(% 5 = ?) 计算值实验计算值实验计算值实验值 1 0.4 4 1.2538 1.25 0.8569 0.863 2.1875 2.1 2 0. 3 4 1.3723 1.37 0.8233 0.828 2.9167 2.81 3 0.2 4 1.5266 1.53 0.8016 0.8 4.3750 4.9 4 0.12 5 4 1.6731 1.67 0.791 6 0.803 7.0000 7.33 w=4; cmax1=1+exp(-z1*pi/sqrt(1-z1^2)); tp1=pi/(w*sqrt(1-z1^2)); ts1=3.5/(z1*w); [cmax1,tp1,ts1] ans = 1.2538 0.8569 2.1875 >> z2=0.3; w=4; cmax2=1+exp(-z2*pi/sqrt(1-z2^2)); tp2=pi/(w*sqrt(1-z2^2)); ts2=3.5/(z2*w); [cmax2,tp2,ts2]
摘要 自动控制原理理论性强，现实模型在实验室较难建立，因此利用SIMULINK进行仿真实验，可以加深我们学生对课程的理解，调动我们学习的积极性，同时大大提高了我们深入思考问题的能力和创新能力。本文针对自动控制系统的设计很大程度上还依赖于实际系统的反复实验、调整的普遍现象，结合具体的设计实例，介绍了利用较先进的MATLAB软件中的SIMULINK仿真工具来实现对自动控制系统建模、分析与设计、仿真的方法。它能够直观、快速地分析系统的动态性能、和稳态性能。并且能够灵活的改变系统的结构和参数，通过快速、直观的仿真达到系统的优化设计。关键词：MATLAB；自动控制；系统仿真
Abstract Strong theory of automatic control theory, the reality is more difficult to establish in the laboratory model, thus using the SIMULINK simulation experiment, students can deepen our understanding of the course, to mobilize the enthusiasm of our study, while greatly increasing our ability to think deeply and Innovationcapacity.In this paper, the design of automatic control system is still largely dependent on the actual system of repeated experiments, adjustment of the universal phenomenon, with specific design example, introduced the use of more advanced software in the MATLAB SIMULINK simulation tools to achieve the automatic control systemModeling, Analysis and design, simulation methods.It can intuitively and quickly analyze the dynamic performance, and steady-state performance. Keywords:MATLAB; Automatic control; System simulation
控制系统的时域响应MATLAB 仿真实训 实训目的 1. 学会利用MATLAB 绘制系统的单位阶跃响应曲线，掌握读取系统动态性能指标的方法； 2. 学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法； 3. 掌握利用MATLAB 绘制系统的零输入响应曲线的方法； 4. 掌握利用MATLAB 绘制系统的一般输入响应曲线的方法； 5. 学会通过仿真曲线读取相关信息，并依据有关信息进行系统的时域分析。 实训内容 1.编写程序求取下列各系统的单位阶跃响应，完成表5-5并记录相关曲线。 162.316)(21++= s s s G 16 4.216 )(22 ++=s s s G 166.116)(2 3++=s s s G 1616 )(24++=s s s G 解：>> n1=16; >> d1=[1,,16]; >> sys1=tf(n1,d1); >> step(sys1) >> n2=16; >> d2=[1,,16]; >> sys2=tf(n2,d2); >> step(sys2)
>> n3=16; >> d3=[1,,16]; >> sys3=tf(n3,d3); >> step(sys3) >> n4=16; >> d4=[1,1,16]; >> sys4=tf(n4,d4); >> step(sys4)
序号ξnω m ax c p t s t(% 5 = ?)计算值实验计算值实验计算值实验值 14 24 34 44 w=4; cmax1=1+exp(-z1*pi/sqrt(1-z1^2)); tp1=pi/(w*sqrt(1-z1^2)); ts1=(z1*w); [cmax1,tp1,ts1] ans = >> z2=; w=4; cmax2=1+exp(-z2*pi/sqrt(1-z2^2)); tp2=pi/(w*sqrt(1-z2^2)); ts2=(z2*w); [cmax2,tp2,ts2]
第三章1. (1)A=eye(3) (2)A=100+100*rand(5,6) (3)A=1+sqrt(0.2)*randn(10,50) (4)B=ones(size(A)) (5)A+30*eye(size(A)) (6)B=diag(diag(A)) 2. B=rot90(A) C=rot90(A,-1) 3. B=inv(A) ;A的逆矩阵 C=det(A) ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(A)*b x = -6.0000 26.6667 27.3333 5. (1) diag(A) ；主对角线元素 ans = 1 1 5 9 triu(A) ；上三角阵
ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ；下三角阵 ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ；秩 ans = 4 norm(A) ；范数 ans = 21.3005 cond(A) ；条件数 ans = 11.1739 trace(A) ；迹 ans = 16 （2）略 6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000
[V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365