1.1 MATLAB 软件介绍 KLBU8%
; :4&nJ*qG
在欧美各高等学校,MATLAB 软件成为线性代数、自动控制理论、数字信号处理、时间序列分析、动态系统仿真、图像处理等诸多课程的基本教学工具,成为本科生、硕士生和博士生的必须掌握的基本技能。 ,.A@U*j
~:'tp28?
主要用途: .!e):&(8
(1)数值和符号计算; iyKAw
(2)绘图; Ye% e!
(3)一种语言体系,也可以方便地与 Fortran、C 等语言接口; Tty_P,
(4)工具箱(Toolbox),分为功能性和学科性两种类。同时,其工具箱的数据文件代码完全开发,用户也可以开发自己的工具箱。 DvB!-|ek
sC1Mwx
功能性工具箱主要用来扩充 MATLAB 软件的符号计算功能、图视建模仿真功能、文字处理功能以及硬件实时交互功能。这种功能性工具箱用于多种学科。而学科性工具箱是专业性比较强的,如控制工具箱(Control Toolbox)、信号处理工具箱(Signal Processing Toolbox)、通信工具箱(Communication Toolbox)等都属此类。 %UT5KYd!=N
bA!n;
MATLAB 语言的主要特点: zl^ %x1G
(1)语法规则简单。与其他编程语言相比更接近于常规数学表达,对于数组变量的使用,不需类型声明,也无需事先申请内存空间。 Yd'Fhvo8
(2)提供了数以千计的计算函数,极大的提高了用户的编程效率。 t~]n"zgovz
(3)是一种脚本式(scripted)的解释型语言,无论是命令、函数或变量,只要在命令窗口的提示符下键入并以回车键结束,则 MATLAB 都会予以解释执行。 #I8)|p?P
(4)可移植性,可跨平台运行。MATLAB 软件可以运行在很多不同的计算机系统平台上,包括大部分的 UNIX 和 Linux 系统,其编写的程序对应的数据文件是一致的,绘图功能也与平台无关。 LM\ H%=*L
Oi%\'biM
1.2 MATLAB 系统环境与数据操作 b+Vfi9<
l25_J.e
1.2.1 系统路径设置 P{fT5K|
(VV5SvdE
两个路径概念: )eIC5>#.
(1)当前目录窗口:指 MATLAB 运行时的工作目录。只有在当前目录和搜索路径下的文件、函数才可以被运行和调用;如果没有特殊指明,数据文件也将存放在当前目录下。 { RH&mu
(2)搜索路径:指 MATLAB 执行过程中对变量、函数和文件进行搜索的路径。 -FpZZ8=,M2
搜索路径是使用 MATLAB 时一个很重要的概念,MATLAB 具有大量的工具箱,用户在使用时也会设置大量的用户文件夹和文件,怎样让 MATLAB 能够识别出这些文件在哪里,这就是搜索路径需要确定的内容。 E6JfSH#
(注:MATLAB 启动后的默认目录是 C:\MATLAB\BIN 若不建子目录则MATLAB 环境产生的数据文件就登陆在这个缺省目录上。为了更加便于用户创建、修改 M 文件和其他文件,用户应建立自己的工作目录。) nsn
键入字符识别的搜寻次序为: 17P5Dr&
->>是否为变量 '5V}Z3zJ/
->>是否为函数 )Q= EmZbJz
->>是否为当前目录下的 M 文件 h K;9XJAf
->>是否为搜索路径下其他 M 文件 i<@"+~n~GK
A0X'|4I
系统路径设置方法 *U>"_h T0
(1)path 指令,可设置扩展(临时)搜索路径。例: jV{?.0/h|
path(path, ‘C:\tools\mydir’) PL} Wu=
(2)cd 指令,可设置当前工作目录,该目录必须已存在。例: C2}n &{T
cd E:\MATLAB xB-\yWDZe
(3)userpath 指令,可修改运行时默认路径。例: :j^IXZW
userpath(‘D:\MATLAB’) 4fauI%kc
(4)pathtool 指令,可直接更改、添加搜索路径。 I\e/
Bv^
idsBw!DB
1.2.2 变量 的基本表达 Znr6,[U+q
,Tp:. "
变量命名规则:以英文字母开头,可包括英文字母、数字和下划线,区分大小写,最多可包含 63 个字符(6.5 及以后版本)。 3bYPi^
变量书写特点: pPL)!=o!
(1)无需声明 m.F}9HI%hN
(2)末尾有分号,只执行不输出 jkAWRpOc)
(3)回车符代表这句结束 '3<T~t
(4)一行写不开,可用…续行 dCA| )
(5)注释以%开头 Lb]!TOl
d*$L$1S
变量赋值的两种常见方式: 5PY4PT=G
->>变量=表达式,例: 3JYhF)G
s = 1-1/2+3 CWBlDz
输出结果:s = 3.5000 Me 5Xd|
->>表达式(赋值给默认变量 ans) f$>KTb({B
2016/81 R7\T.;8+
输出结果:ans = 24.8889 A1Ru&fd!
*^y,Gg/
MATLAB 预定义的变量: B]2m(0Y>>v
(1)ans,最近表达式计算结果的变量名; <+y%k~("
(2)eps,预设正的极小值,其值为 2.2204e-16; n(uzqd
(3)pi,圆周率 π; 1;v,rs M
(4)inf,无穷大,∞值; T.(C`/VM
(5)i 或 j,虚数单元,sqrt(-1); N\l|3~
(6)NaN,非数,0/0,∞/∞。 G&Fe2&5!w
注:用户只能临时覆盖这些预定义变量的值,Clear 或重启 MATLAB 可恢复其值。 e"#QUc(
P2Jo^WS
变量查询方法: MO^Q 8v
(1)工作空间窗口显示; bG)EZ
(2)who 和 whos 指令,列出在 MATLAB 工作空间中已经驻留的变量名清单。其中,whos 指令在给出变量名的同时,还给出它们的维数及性质。 Mz7qC3Z
p40;@gUug
清除变量方法: 'VpzB
s#
->>clear 清除所有变量 g^l~AR
->>clear a 只清除变量 a FEH+ PKSc
NufRd/q
变量保存与调用: r01u3!
命令 quit 和 exit 都可退出 MATLAB,结束工作任务会删除工作空间中的变量,可采用如下命令来保存工作空间,以备再次调用这些变量。 ?B+]Ex(\B,
^HhV?Iqg
保存的指令格式: Sw.k,p*r
(1)save,将工作空间中的所有变量保存到 matlab.mat 文件中。 mzf~qV^T
(2)save [文件名] [变量名],将指定的变量保存在指定文件中。如: &w!(.uDO
save temp x y z _Z>ny&
把 x,y,z 这三个变量保存在文件 temp.mat 中。 +~K)
~
在下次加载 MATLAB 时,可以利用 load 指令将保存在文件中的变量恢复到工作空间中,其格式有: y(ceEV
(1)load,将保存在 matlab.mat 中的变量载入到 MATLAB 工作空间中。 'yo-`nNFD
(2)load [文件名] [变量名],从指定的文件中将指定的变量装入 MATLAB工作间。如: /IQ$[WR cx
load temp x 0fGt7 "Q
只将 temp.mat 文件的变量 x 载入到 MATLAB 工作空间中。 |.KB
=;"e Z
1.2.3 数据、操作与显示 2=
Y8$-
|2E:]wT}qg
数字的计述: hG~HV{6
MATLAB 的数值采用习惯的十进制表示,可以带小数点和负号,其缺省的数据类型为双精度浮点型(double)。例: _z=ytt9D
3 -99 0.0013 1.243e-6 J#IVu?B
表达式按常规相同的优先级从左至右执行运算,指数运算级别最高,乘除次之,加减最低。 XuoyB{U
加+ 减- 乘* 左除\ 右除/ 幂^ b_&:tE--]
注:对于标量而言,左除和右除的作用相同,但对矩阵运算产生不同的影响。 ?[]jJ
,|g&v/WlC%
在缺省的状态下 MATLAB 以短格式 short 格式显示计算结果。可以用MATLAB 命令窗口中 format 指令来改变数字的显示格式。由于 MATLAB 以双精度执行所有运算,显示格式的设置仅影响矩阵的显示,不影响矩阵的计算与存储。 g@'2 :'\
如果矩阵的所有元素都是整数,则矩阵以不带小数点的格式显示。如果有一个元素不是整数,则有几种输出格式。默认格式为 short格式,只显示 5 位有效数字,其他的显示格式可显示更多的有效数字,还可用科学表示法。 (A )f
r4
->>format short e 短格式科学表示; Nwj M=GG
->>format short g 在 short 和 short e 中自动选择最佳方式; G#Kw6
->>format long 长格式(15 位有效数字); cOf.z)kf6
->>format long e 长格式科学表示; wg+[T;0 S
->>format long g 在 long 和 long e 中自动选择最佳方式; qZdA%
->>format bank 银行格式; z }Lf]w?
->>format rat 近似有理数表示; m*wDJEKo
->>format hex 十六进制格式; KVevvy)W
->>format compact 紧凑格式。 Quwq_.DU
……
s)jNP\-
2#
字符串相关函数: =Fc}T%
class( )函数,可查看数据类型 #kk5{*`
double( )/abs( )函数,可将字符串转换成 ASCII 码 "/G]M&
str2num( )/num2str( )函数,实现字符串与数字之间的转换。 [!C!R$AMa
eval( )函数,执行字符串变量中的内容。 $O&N
#@' B\!<@=
1.3 MATLAB 帮助/命令查询 DUSQh+C
1pK(tm
MATLAB 的在线帮助系统相当完备,就查询系统的调用方式而言,可分为两种: O=eU38n:5u
(1)从 MATLAB 指令窗的 help 菜单选项中寻求帮助此与一般 windows的求助方法一样 ^6gEL~m|]
(2)在 MATLAB 指令窗中直接键入帮助指令。 6i;q=N$'
->>help 不带任何参数,显示出 MATLAB 的目录项产生清单信息; G6}!PEwM
->>help 目录名(lang/ matfun/ elfun),显示指定目录中的所有命令及其函数; ykRd+H-t
->>help 命令名/函数名/符号,显示指定的命令名/函数名/符号的详细信息。 EERCb%M8Z
lookfor 指令可以根据用户提供的完整或不完整的关键词,去搜索出一组与之有关的指令。 L#?mPF
exist 指令,检查指定名字的变量或函数文件的存在性; k]A=Q
what 指令,按扩展名分类列出在搜索路径中指定目录上的文件名; baBPf{<
which 指令,列出指定名字文件所在的目录。 F~0iJnF
6"[J[7up
注:help 的工作机理是把指定名字的那个 M 文件的第一段注释内容显示出来,以构成自己文件的再线求助。lookfor 指令的机制是对 MATLAB 中的每个 M文件注释区的第一行进行扫描,一旦发现包含要查询的字符串就显示出来。用户也可利用此机理建立自己文件的在线帮助。 m70AWG
!<HF764@`
1.4 MATLAB 的矩阵操作 \T[OF8yhW
,k1ns?i9KH
1.4.1 矩阵的建立 ~# \{'<
gL~3z'$
(1)直接输入法 \x<,Ma=D
->>在命令窗口直接输入矩阵元素; ^I9U<iNIL
->>按行输入矩阵元素,同一行元素使用空格或者逗号隔开,不同行元素使用分号或者回车作为间隔。 37biRXqLH
XTA:Y7"O
(2)M 文件建立矩阵 p(-EtxP
->>启动编辑器(输入 edit 命令,或者使用菜单栏按钮); f*:N*cC
->>输入待建立矩阵,输入格式要求和直接输入法类似。 Ga+Cb2$
/3.;sS]B
(3)zeros( )和 ones( )函数 CfO{KiM(2
z}iz~WZ
(4)冒号表达式建立向量 NiEz3ODSi
冒号是一个重要的运算符,语法: \vx'+}
e1:e2:e3 I]HLWF
e1 为初始值,e2 为步长(默认为 1),e3 为终止值,产生行向量; P8#_E{f
7L!}F;yT
(5)linspace 建立行向量 mhM;`dl
语法: wz@[rMf
linspace(a,b,n) &> _aY #
a 和 b 分别为行向量的第一个和最后一个元素,n 为总元素,省略 n 自动产生 100 个元素的行向量 FRL;fF
^Y!$WP
1.4.2 矩阵的简单操作 I %sw(uoE
]<ay_w;
(1)索引矩阵 3UF^Ff<wo
->>可以使用矩阵元素的行标和列标或者矩阵元素的位置索引来索引矩阵元素; Q-A_ 8
->>行标和列标共同使用的形式; m<LzgX
->>使用序号索引,涉及到 MATLAB 中矩阵的存储形式; R5mb4
矩阵元素的存放是按照列来排序 C8i4z
find(A==8)或 find(A,8)返回 A 矩阵内值为 8 的元素序号;find(A<8)返回所有小于 8 的序号; @Z7s3b
->>序号和下标的转换。 >l=jJTJ;q
sub2ind(size(A),m,n)返回序号 P8H2v_)X&
[m,n]=ind2sub(size(A),3)返回下标 *NM*
zlB[Eg^X
(2)重排矩阵 CKSs(-hkJ
->>A=[1,2,3;4,5,6;7,8,9] ~[kI![
->>E=reshape(A,9,1)重排为 9×1 列向量 M04u>|
,
->>I=A[:]同上效果 @\:@_}Z`_}
`Ba?4_>k
(3)矩阵的转置 0C3Y =F
->>L=K’ Z%b1B<u$
U6o]7j&6
(4)矩阵的拆分 e|>@ >F]K
->>K=A(1,:)获得第一行所有列的内容 +;)Xu}
->>K=A(1,1:end)效果同上
1D2RhM%
->>K=A(1,1:2) 获得第一行 1 到 2 列的内容 *v: .]_;
->>K=A(1,[1 2])效果同上 D(&Zq7]n
->>K=A([1 2],[1 2]) 0+)1KU)I
DjiI*HLNR
(5)删除矩阵元素 "uj@!SEs`?
->>赋空值,例如: oA:`=f%\
M(2)=[ ] EP90E^v^
A(1,:) =[ ] :T>OJ"p
&LG|YvMY6
(6)扩展矩阵元素 >]A#_p
->>Q=repmat(A,2,1) Bk@EQdn
将[1 2 3 4]扩展为[1 2 3 4;1 2 3 4] O2f-{jnTz,
B/mfm 7
(7)压缩矩阵元素 IL uQf-
->>R=[1 2 3 3 3 4 4]
hFan$W$
->>Q=unique(R)只保留不重复的元素 (=Oo=8\
sHV?njZd
1.4.3 矩阵的处理基础 eht>4)
rmFcSolt,f
MATLAB 运算以矩阵为单位,可以直接实现矩阵的加法、数乘等运算以及矩阵的行列式、矩阵的秩、矩阵的逆以及转置等运算。 % .ss
计算 A 的行列式; *(XGNp[0
->>使用 det(A)函数 g&Uu~;jq]
如果 A 的行列式不为 0,那么计算 A 的逆; bA'N2~.,
->>使用 inv(A)函数 Q ~n%c7
*.VNyay
向量可以看作是 1×n 矩阵或者 n×1 矩阵,因此向量的加法和数乘等运算和矩阵式一样的。 91nB?8ZE6,
求向量内积(a,b)=b H a; cXr_,>k
->>sum(conj(b).*a) dB QCr{7
->>a*b’,与上式一致 ccx0aC3@I
->>dot(b,a),与上式一致 S\GxLW@x
ka{!' ^
线性方程组求解 c+O:n:L
可以使用矩阵求逆的方法求解线性方程组,在一般线性代数中的解决方法类似。 wbk$(P'gN
->>x=inv(A)*b s;[=B
->>x=A\b,与上式一致 w'y,$gtX/
JW=q'ibR
矩阵的相似化简和分解 6\VZ6oS
求解 A 的 Jordan 标准型; ,% .)mf
->>使用 jordan(A)函数 T{N8 K K
求解 A 的特征值; t1$pl6&,
->>使用 eig(A)函数 9[
o$/x}
范数运算,包括 1 范数,2 范数和无穷范数以及 f 范数; ?iamo.0zN
->>使用 norm 函数,norm(A,1),norm(A,2),norm(A,inf),norm(A,’fro’) +CXq41g"c
nQg_1+
函数分析 M~y}0Ik
计算函数矩阵 A(x)的一阶和二阶导数。在 MATLAB 中函数的计算,比如求导数,默认对矩阵的每一个元素求导,这即为函数矩阵导数的计算方法。 G
c,
->>syms x 9Sa6v?sRor
->>A=[sin(x) exp(x) 1;cos(x) x^2+1 log(x)] <^942y-=
->>diff(A),一阶导数 znIS2{p/`
->> diff(A,2),二阶导数 n;:C{5
Ysw&J}6e
计算矩阵函数,e A ,sinA,cosA。在 MATLAB 中使用 funm 函数计算矩阵函数,funm 为通用矩阵函数计算函数,其使用方法如下: ta'wX
funm(A,@fun) ivt ~S
->>funm(A,@exp),求 e A i'1MZ%.
->>expm(A),同上 -3m!970
->>funm(A,@sin),求 sinA vTWm_ed+^
->>funm(A,@cos),求 cosA vfc,{F=Q
e)87
&
7
1.5 MATLAB 程序控制结构 gkca{BJ
41Ab,
1.5.1 M 文件 6!;D],,"#.
)M"xCO3a
分类: !-&;t7R
命令文件(脚本文件,Script File) 5{v uN)K3
函数文件(Function File) ^oClf(
P{: 5i%qC
区别: %&s4YD/{
命令文件没有输入,没有返回; &^Gp
命令文件可以对工作空间的变量操作,结果返回工作空间,而函数文件中的变量为局部变量,函数执行完毕,变量被清除; e+<9Sh7&
命令文件可以直接运行,函数文件见需要调用的形式运行(除了特殊的函数文件外)。 qG]0z_dPE~
.*J /F$
1.5.2 顺序结构 k/BlkjlNE
YlR9
1LX
按照代码顺序执行,即顺序结构,受控于代码前后位置; C9,Uwz<!]
使用于顺序结构中的简单输入和输出(命令交互): 1S yG
数据输入(input 函数) hZ"Sqm]
->>A=input(‘输入一个参数’) gyT0h?xDt
数据显示(disp 函数) FchO
6O
->>disp(A) %c8@
程序暂停(pause 函数) 8==_43
->>pause(3)暂停 3 秒 ["<'fq;PJ
.Sv/0&O
1.5.3 分支结构 "M^mJl&*b
IA`Lp3Z
分支结构,又称为选择结构,包括 if 分支,switch 分支和 try 分支。 'q@vTM'-
(1)if 分支语句: +/" \.wYv
if 条件 0F
2p4!@W
语句组 1 8`90a\t'Z
else I<+:Ho=6
语句组 2 p\.IP2+c
end *9EW&Ek
k^zU;
(2)switch 分支语句 1jaK N*
switch 表达式 [X>f;;h
case 表达式 1 0]^gT'
语句组 1 J@#rOOu
case 表达式 2 gP |>gy#e
语句组 2 i)$<j!L
… K6|R ;r5e{
case 表达式 n @Y}G,i
语句组 n jvo^I$|2h
otherwise r d)W+W9
语句组 m 432]yhQ
end 5Yww,s
8Dn~U:F/?
(3)try 分支语句 eo.B0NZsF
在经典的程序设计语言中未出现,这是一种试探性的分支语句,也就是如果这条语句执行有错误,则不执行,或者执行其他的语句,提高了程序的容错性能。 wyXQP+9G
try 'rA(+-.M;
语句组 1 xJA{Hws
catch t 6lwKK
语句组 2 Jb-.x_Bf
end (A "yE4rYK
\)ZCB7|
1.5.4 循环结构 77ztDQDtM
按照给定条件,重复执行某些语句,包括 for 循环和 while 循环。 MV07RjeS
(%"9LYv
(1)for 循环语句 k|U2Mp
for 循环变量=表达式 1:表达式 2:表达式 3(可以采用矩阵的形式) xV=Tmu6l
循环体 2.
v<pqn
end v-/vj/4>
注意循环变量自动增加,在循环体内对循环变量的赋值操作会带来不可预料的程序执行。 \1QY=}
;|:R*(2
(2)while 循环语句 BfD&