2. Matlab入门.pptVIP

Matlab 入门 徐士河 shihexu03@ 1．指令行的编辑 1.1 启动Matlab后，就可以利用Matlab工作。 由于Matlab是一种交互式语言，随时输入指 令，即时给出运算结果是它的工作方式。 如：2*sin(0.3*pi)/(1+sqrt(5)) ans=0.5000 (ans是一个保留的Matlab字符 串，它表示上面一个式子的返回结果，用于 结果的缺省变量名) 2.简单的矩阵输入 ① 要直接输入矩阵时，矩阵元素用空格或逗号分隔；矩阵行用分号“；”隔离，整个矩阵放在方括号“[ ]”里。 A=[1,2,3;4,5,6;7,8,9] 说明：指令执行后，矩阵A被保存在Matlab的工作间（workspace）中，以备后用。如果用户不用clear指令清除它，或对它进行重新赋值，那末该矩阵会一直保存在工作间中，直到本 Matlab指令窗关闭为止。 2.简单的矩阵输入 ② 矩阵的分行输入 A=[1,2,3 4,5,6 7,8,9] 3.语句与变量 Matlab语句有两种最常见的形式 i) 表达式 ii)变量=表达式 [例1] 表达式的计算结果 1996/18 ans=110.8889 [例2] 计算结果的赋值 s=1-1/2+1/3-1/4+1/5-1/6+... 1/7-1/8; 说明：三个小黑点是“连行号”，分号“；”作用是：指令执行结果将不显示在屏幕上，但变量s将驻留在内存中。若用户想看s的值，可键入以下命令：s s= 0.6345 4.数据结构：向量与矩阵 ①向量的转置 z=[1+j,2+pi*i,-sqrt(-1)]’ z = 1.0000 - 1.0000i 2.0000 - 3.1416i 0 + 1.0000i 当对复数向量进行转置操作时,可以得到其共轭转置向量。 z. (非共轭转置向量) 4.数据结构：向量与矩阵 ②产生一个行向量 t=[0:0.1:10] %产生从0到10的行向量，元素之间间隔为0.1 t=linspace(n1,n2,n)%产生n1和n2之间线 性均匀分布的n个数 (缺省n时,产生100个数) t=logspace(n1,n2,n) (缺省n时,产生50个 点) 4.数据结构：向量与矩阵 ③矩阵的标号 A(m,n)表示矩阵A的第m行，第n列的元素； A(1:2,1:3)表示矩阵A的从第一行到第二行，从第一列到第三列的所有元素； A(:) 可以得到一个长向量,该向量的元素是按列一一叠加在一起的。例如 a=[1 2;3 4]; a(:) ans= 1 3 2 4 4.数据结构：向量与矩阵 ④ 特殊矩阵 i)单位矩阵 eye(m),eye(size(a)) 可以得到与矩阵a同样大小的单位矩 阵; eye(m,n)可得到一个可允许的最大单位矩阵而其余处补0。 ii)所有元素为1的矩阵 ones(n)，ones(size(a))，ones(m，n) iii）空矩阵是一个特殊矩阵,这在线性代数中是不存在的。 例如: q=[ ]矩阵q在工作空间之中,但它的大小为零。通过空 阵的办法可以删去行与列。例如 a=rand(5,5)； a(:,1:3)=[] 第一条指令执行后得到一个5′ 5的矩阵；第二条指令将矩阵 a的前三列删除。 4.数据结构：向量与矩阵 iv) 对角矩阵 当v是向量时，diag(v)得到以v的元素为对角线上元 素的对角矩阵； 当v是矩阵时，diag(v) 得到一个列向量，其 元素为矩阵v对角线上的元素， diag(v,1)得到矩阵v对角线上移一行的元素组 成的列向量，diag(v,-1)得到矩阵v对角线下 移一行的元素组成的列向量。 ⑤ 字符串要用单引号并用括号括在里面。例如： disp(text string) % disp 显示命令 7. 数学运算与函数 ① 基本代数运算操作 +，-，*，\，/，^, ② 矩阵运算函数：求行列式（det）,矩阵求逆（inv）,求秩（rank），求迹(trace)，求模（norm）,d=eig(A)求矩阵A的特征值[v,d]=eig(A)求矩阵A的特征向量和特征值，这里v的列向量是对应的特征向量。 ③矩阵基本运算： A\B, B/A, A.*B, A./B, A.\B, A.^B。 7. 数学运算与函数 ④ 多项式 任意多项式都可以用一个行向量来表示，即n +1维 的向量a表示多项式 y(x) = a(1)xn-1 + a(2)xn-2 +?+ a(n)x + a(n +1)，