第二章 MATLAB.docVIP

2 MATLAB——一个强有力的科学计算与 控制系统CAD软件 2.1 概述 MATLAB的名字由MATrix和LABoratory两词的前三个字母组合而成，始创者是时任美国新墨西哥大学计算机科学系主任的Cleve Moler教授。于1984年由MathWorks公司推出（DOS版）。1993年推出MATLAB 4.0（Windows版）。1997年MATLAB 5.0版问世。1999年春已到了5.3版。今天MATLAB已成为国际上最优秀的科技应用软件之一，其强大的科学计算与可视化功能、简单易用的开放式可推展环境以及多达30余个面向不同领域而扩展的工具箱（Toolbox）支持，使得MATLAB在许多学科领域在成为科学计算、计算机辅助设计与分析的基本工具和首选平台。 在国内外的许多大学中，MATLAB已成为自动控制，信号分析等课程和教科书的内容，成为许多学科攻读学位的大学生、硕士生、博士生的必须掌握的工具。在国际学术界，MATLAB已被确认作为一个准确、可靠的科学计算软件而得到承认。在许多国际一流的学术刊物上都可以看到MATLAB的应用，在研究单位和工业部门，MATLAB已成为一些分析统计软件的工作平台和主要支撑。许多硬件、仪器都接受MATLAB的支持。 面向控制工程应用一直是MATLAB的主要功能之一，早期的版本就提供了控制系统设计工具箱（Control Toolbox）。90年代初的3.5版推出Robust Toolbox， 4.0版推出基于模块图的控制系统仿真软件 Simulink。到目前为止，MATLAB 5.x版本中包含的控制工程类工具箱已超过10个。MATLAB所具备的强有力的计算功能和图形表现，以及各种工具箱提供的丰富的专用函数，为设计研究人员避免重复繁琐的计算和编程，更快、更好、更准确地进行控制系统分析和设计提供了极大的帮助。 本书的后叙章节中将结合控制系统理论与设计方法介绍几个控制类工具箱的应用。本章中，我们将通过一些简单例子表现MATLAB的编程和功能。 2.2 MATLAB基本功能及举例 本节中，我们通过一些简单例题来说明MATLAB的基本计算与绘图功能。 2.2.1 矩阵运算 既然以Matrix Laboratory的缩写来命名，MATLAB无疑将在矩阵的表示和运算方面独具特色。 矩阵的MATLAB表示： 对于矩阵，可以简单地表示为A=[1 2 3 ；4 5 6；7 8 9] 此外，a=[-1.3/7 sqrt(4.5) 2+3i]也表示一个矩阵。 在MATLAB中，可以方便地用C=[A；a]来构成一个新的矩阵。 也可以用 x=C（：，2）等命令来提取一个矩阵的某一部分。 MATLAB的矩阵运算： MATLAB在矩阵运算方面的优势则更是明显。如，可以直接用 C=A+B， D=A*B等计算式形式的命令来完成矩阵的四则运算。 如有 AX=B 则求解指令为 X=A\B XA=B 则求解指令为 X=B/A 矩阵转置： B=A’ 矩阵求逆： B=inv(A) 行列式： det（A） 此外，在矩阵求特征值，LU分解等诸多方面，MATLAB都提供了大量丰富的专用函数和指令，使得用户可以用计算公式的形式来解决Fortran、C等语言中须通过复杂编程来求解的问题。 读者可以通过运行本书提供的chp2_1.m程序来体会MATLAB中矩阵计算的一些基本方法。 %======= Program chp2_1 ========= %-------- Basioc Module -------- clc disp(Matrix expression) disp(A=[1 2 3;4 5 6;7 8 9] :) A=[1 2 3;4 5 6;7 8 9] disp(a=[-1.3/7 sqrt(4.5) 2+3i] :) a=[-1.3/7 sqrt(4.5) 2+3i] disp(Matrix Turn :) B=A pause clc A=[1 2 3;4 5 6;7 8 9] a=[-1.3/7 sqrt(4.5) 2+3i] disp( C=[A;a] :) C=[A;a] disp( x=C(:,2) :) x=C(:,2) pause clc A=[1 2 3;4 5 6;7 8 9] disp( C=A+A :) C=A+A disp( B=A*E(E is an 3x3 Identity matrix) :) B=A*eye(3) pause clc A=[1 2 3;3 2 1;1 2 5] B=[-2 9 -15;6 15 -21;-6 11 -2] disp(Ax=B ===== x=A\B :) x=inv(A)*B disp(the determinant