MATLAB教程 基础准备及入门

1、MATLAB (MATrix LABoratory)具有用法简易、可灵活运用、程序结构强又兼具延展性。以下为其几个特色： · 功能强的数值运算 - 在MATLAB环境中，有超过500种数学、统计、科学及工程方面的函数可使用，函数的标示自然，使得问题和解答像数学式子一般简单明了，让使用者可全力发挥在解题方面，而非浪费在电脑操作上。 · 先进的资料视觉化功能 - MATLAB的物件导向图形架构让使用者可执行视觉数据分，并制作高品质的图形，完成科学性或工程 性图文并茂的文章。 · 高阶但简单的程式环境 - 做为一种直译式的程式语言，MATLAB容许使用者在短时间内写完程

2、式，所花的时间约为用 FORTRAN 或 C 的几分之一，而且不需要编译(compile)及联结 (link) 即能执行，同时包含了更多及更容易使用的内建 功能。 · 开放及可延伸的架构 - MATLAB容许使用者接触它大多数的数学原码，检视运算法，更改现存函数，甚至加入自己的函数使 MATLAB成为使用者所须要的环境。 · 丰富的程式工具箱 - MATLAB的程式工具箱融合了套装前软体的优点，与一个灵活的开放但容易操作之环境，这些工具箱提 供了使用者在特别应用领域所需之许多函数。现有工具箱有：符号运算（利用Maple V的计算核心执行 ）、影像处理、统计分析、讯号处理、神

3、经网路、模拟分析、控制系统、即时控制、系统确认、强建控 制、弧线分析、最佳化、模糊逻辑、mu分析及合成、化学计量分析。 第一章 基础准备及入门1.1 Desktop操作桌面的启动l MATLAB的启动l Desktop操作桌面简介1.2 Command Window运行入门l Command Window指令窗简介l 最简单的计算器使用法【例1.2-1】求的算术运算结果。>> (12+2*(7-4)/32 【例1.2-2】简单矩阵的输入步骤。A = 1,2,3; 4,5,6; 7,8,9 A = 1 2 3 4 5 6 7 8 9 【例1.2-3】矩阵的分行输入。A=1,2,34,

4、5,67,8,9 A = 1 2 3 4 5 6 7 8 9 【例1.2-4】指令的续行输入S=11/2+1/31/4+1/51/6+1/7 -1/8 S =0.6345 1.2.1 数值、变量和表达式 一 数值的记述3 -99 0.001 9.456 1.3e-3 4.5e32 二 变量命名规则变量的第一个字符必须是英文字母，最多可包括31个字符（英文、数字和下连符），区分大小写。如Mas_2，mas_2为两个不同的变量。 三 MATLAB默认的预定义变量help 在线上说明, 如 help，quitwho 列出所有定义过的变量名称ans 预设的计算结果的变量名pi=3.1415926 i=

5、j=inf 无穷大eps 机器中的最小数NaN 不定值nargin 所用函数的输入变量数目nargout 所用函数的输出变量数目realmin 最小可用正实数realmax 最大可用正实数 四 运算符和表达式数学表达式MATLAB运算符MATLAB表达式加a+b+a+b减a-b-a-b乘a×b*a*b除a÷b/或a/b或 ba幂abab说明：l 所有运算定义在复数域上。对于方根问题，运算只返回处于第一象限的解。l MATLAB用左斜杠或右斜杠分别表示“左除”或“右除”运算。对标量而言，两者的作用没有差别；但对矩阵来说，“左除”或“右除”运算的意义不同，结果也不同。 五 复数

6、和复数矩阵【例1.2-5】复数表达，及计算。（1）z1= 3 + 4i z1 = 3.0000 + 4.0000i （2）z2 = 1 + 2 * iz3=2*exp(i*pi/6)z=z1*z2/z3 z2 = 1.0000 + 2.0000iz3 = 1.7321 + 1.0000iz = 0.3349 + 5.5801i 【例1.2-6】复数矩阵的生成及运算A=1,3;2,4-5,8;6,9*iB=1+5i,2+6i;3+8*i,4+9*i C=A*B A = 1.0000 - 5.0000i 3.0000 - 8.0000i 2.0000 - 6.0000i 4.0000 - 9.00

7、00iB = 1.0000 + 5.0000i 2.0000 + 6.0000i 3.0000 + 8.0000i 4.0000 + 9.0000iC = 1.0e+002 * 0.9900 1.1600 - 0.0900i 1.1600 + 0.0900i 1.3700 【例1.2-7】求上例复数矩阵C的实部、虚部、模和相角。C_real=real(C)C_imag=imag(C)C_magnitude=abs(C)C_phase=angle(C)*180/pi C_real = 99 116 116 137C_imag = 0 -9 9 0C_magnitude = 99.0000 116

8、.3486 116.3486 137.0000C_phase = 0 -4.4365 4.4365 0 【例1.2-8】 用MATLAB计算能得到 2 吗？（1）a=-8;r=a(1/3) r = 1.0000 + 1.7321i （2）的全部方根计算如下m=0,1,2;R=abs(a)(1/3);Theta=(angle(a)+2*pi*m)/3;rrr=R*exp(i*Theta) rrr = 1.0000 + 1.7321i -2.0000 + 0.0000i 1.0000 - 1.7321i （3）图形表示t=0:pi/20:2*pi;x=R*sin(t);y=R*cos(t);plo

9、t(x,y,'b:'),gridhold onplot(rrr(1),'.','MarkerSize',30,'Color','r')plot(rrr(2,3),'o','MarkerSize',15,'Color','b')axis(-3,3,-3,3),axis squarehold off 图 1.2-2 （-8）的全部三次方根分布1.2.2 计算结果的图形表示【例1.2-9】画出衰减振荡曲线及其它的包络线。的取值范围是。（图1.2-3）t=0:

10、pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b') 图 1.2-3 衰减振荡曲线与包络【例1.2-10】画出所表示的三维曲面（图1.2-4）。的取值范围是。clear;x=-8:0.5:8;y=x'X=ones(size(y)*x;Y=y*ones(size(x);R=sqrt(X.2+Y.2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot) xlabel('x'),ylabe

11、l('y'),zlabel('z') 图 1.2-4 z = f (x,y) 的三维网线图1.3 Command Window操作要旨数值计算结果的显示格式format short（缺省形式）5 位定点表示format long 15 位定点表示format short e 5 位浮点表示format long e 1 5 位浮点表示指令行中的标点符号MATLAB常用标点的功能名称标点作 用空格用作输入量之间的分隔符；数组元素分隔符逗号，用作要显示计算结果的指令与其后指令之间的分隔符；用作输入量与输入量之间的分隔符，用作数组元素之间的分割符黑点.表示数值中的小数

12、点。分号；用作不显示技算结果指令的“结尾”标志，用作不显示计算结果指令与其后指令的分隔；用作数组行间分隔符。冒号：用来生成一维数组；用作单下标援引时，表示全部元素构成的长列；用作多下标援引时，表示那维上的全部。注释号%标注非执行行单引号字符串记述符圆括号（）数组援引时用方括号输入数组时用下连符_用作变量的连字符续行号.由4个以上黑点组成，把其下的物理行看作该行的“逻辑”行。 指令窗的常用控制指令常见的通用操作指令指令含义指令含义cd设置当前工作目录exit关闭/退出 MATLABclf清除图形窗quit关闭/退出 MATLABclc清除指令窗中显示的内容mkdir('c:',&

13、#39;my_dir0')创建子目录clear清除MATLAB工作空间中保存的变量more使其后的显示内容分页进行dir列出指定目录下的文件和子目录清单type显示指定M文件的内容edit打开M文件编辑器which指出其后文件所在的目录指令窗中指令行的编辑MATLAB指令窗中实施指令行编辑的常用操作键键名作用键名作用前寻式调回已输入过的指令行Home使光标移动到当前行的首端后寻式调回已输入过的指令行End使光标移动到当前行的未端在当前行中左移光标Delete删去光标右边的字符在当前行中右移光标Backspace删去光标左边的字符PageUp前寻式翻阅当前窗的内容Esc清除当前行的全部内

14、容PageDown后寻式翻阅当前窗的内容【例1.4-1】指令行操作过程示例。（1）若用户想计算的值，那末用户应依次键入以下字符y1=2*sin(0.3*pi)/(1+sqrt(5) y1 = 0.5000 （3）通过反复按键盘的箭头键，实现指令回调和编辑，进行新的计算。y2=2*cos(0.3*pi)/(1+sqrt(5) y2 =0.3633 1.4 Command History1.4.1 Command History历史指令窗简介1.4.2 历史指令行的再运行【例1.5-1】再运行图1.5-2所示历史指令窗中的三行指令。图1.5-2 再运行历史指令的演示1.5 Current Dire

15、ctory、路径设置器和文件管理1.5.1 Current Directory当前目录浏览器简介 一 用户目录和当前目录设置 二 借助当前目录浏览器获取M和MAT文件信息【例1.6-1】从图1.2-1所示MATLAB缺省桌面开始，叙述引出图1.6-1所示面貌的当前目录浏览器的操作过程。1.5.2 MATLAB的搜索路径和路径设置器 一 搜索路径及使用要点1.5.3 MATLAB搜索路径的扩展和修改 一 何时需要修改搜索路径 二 利用设置路径对话框修改搜索路径 三 利用指令path 设置路径1.6 Workspace Browser和Array Editor1.6.1 Workspace Bro

16、wser工作空间浏览器简介1.6.2 现场菜单用于内存变量的查阅和删除 一 内存变量查阅、删除的指令操作法【例1.7-1】在指令窗中运用who, whos查阅MATLAB内存变量。who Your variables are:A Bnumber D R XYZ Z y B C DD X Y x whos Name Size Bytes Class A 2x2 230 cell array B 1x1 264 struct array Bnumber 1x1 8 double array C 2x2 408 sym object D 1x2 4 char array DD 2x2 8 char

17、array R 33x33 8712 double array X 33x33 8712 double array XYZ 33x33x3 26136 double array Y 33x33 8712 double array Z 33x33 8712 double array x 1x33 264 double array y 33x1 264 double arrayGrand total is 7722 elements using 62434 bytes 【例1.7-2】在指令窗中运用clear指令可以删除内存中的变量。clear Bnumberwho Your variables

18、are:A B C D DD R X XYZ Y Z x y 二 内存变量查阅和删除的现场菜单操作法【例1.7-3】通过“工作空间浏览器”的运作，采用图形显示内存变量Z 。（注意：本例是在例1.3-10运行后进行的，因此内存中保存着由例1.3-10产生的全部变量。）图 1.7-3 数组Z 的三维网线图【例1.7-4】通过“工作空间浏览器”删除内存变量。（1）删除部分内存变量的方法（2）删除全部内存变量的方法1.6.3 Array Editor数组编辑器和大数组的输入1.6.4 数据文件的存取 一 存取数据文件的指令操作法 二 通过内存变量浏览器实现数据文件的存取【例1.7-5】数据的存取。（假

19、定内存中已经存在变量X,Y,Z）（1）建立用户目录，并使之成为当前目录，保存数据mkdir('c:','my_dir');cd c:my_dirsave saf x y zdir . . saf.mat （2）清空内存，从saf.mat 向内存装载变量Zclearload saf Zwho Your variables are:Z 1.7 Editor/Debugger和脚本编写初步1.7.1 Editor/Debugger M文件编辑器简介图 1.9-1 M文件编辑器示图1.7.2 M脚本文件编写初步【例1.9-1】编写解算例1.3-9题目的M脚本文件，并运行

20、之。选中M文件编辑器下拉菜单项【Debug：Run】，就可在图形窗中看到如图1.3-3所示的曲线。1.8 帮助系统1.8.1 帮助方式概述 一 “纯文本”帮助【例1.10-1】在指令窗中运行help 的示例。（1）help help HELP On-line help, display text at command line. HELP, by itself, lists all primary help topics. Each primary topic corresponds to a directory name on the MATLABPATH. . . . （2）help HE

21、LP topics:matlabgeneral - General purpose commands.matlabops - Operators and special characters.matlablang - Programming language constructs.matlabelmat - Elementary matrices and matrix manipulation.matlabelfun - Elementary math functions. . . . For more help on directory/topic, type "help topi

22、c". （3）help elmat Elementary matrices and matrix manipulation. Elementary matrices. zeros - Zeros array. ones - Ones array. eye - Identity matrix. . . . （4）help eye EYE Identity matrix. EYE(N) is the N-by-N identity matrix. EYE(M,N) or EYE(M,N) is an M-by-N matrix with 1's on the diagonal a

23、nd zeros elsewhere. EYE(SIZE(A) is the same size as A.See also ONES, ZEROS, RAND, RANDN. 【例1.10-2】在指令窗中，运用lookfor找H1行（M函数文件的第一注释行）包含关键词fourier的所有M函数文件名。（注意把本例与例1.10-4、例1.10-5比较。）lookfor fourier FFT Discrete Fourier transform.FFT2 Two-dimensional discrete Fourier Transform.FFTN N-dimensional discrete Fourier Transform.IFFT Inverse discrete Fourier transform.IFFT2 Two-dimensional inverse discrete Fourier transform.IFFTN N-dimens