[其它]Matlab课件.ppt

例: a =[1, 2, 3;4, 5, 6;7,8,9] g=[1,2,0;4,0,0;9,8,7] a =[1, 2, 3;4, 5, 6;7,8,9] h =[2,1,2;3,6,8;6,8,8] 这些操作符与逻辑操作符配合使用，可使程序设计更加灵活。 例如： 又如： 输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母。 fs=input(请输入分数：) if fs60 disp(不及格) elseif fs=60 fs75 disp(一般) elseif fs=75 fs90 disp(良好) elseif fs=90 fs=100 disp(优秀) else disp(分数错误) end \ 左除 A\B=inv（A）*B 其中inv（A）=A-1 左除对于求解线性方程组非常有用:对于AX=B, 若已知A-1,X= A-1B eg: 求解线性方程组的解 6X1+3X2+4X3=3 -2X1+5X2+7X3=-4 8X1-4X2-3X3=-7 ● ?A^B完成矩阵幂，它有两种形式： X^p(指数为标量)和x^P(底数为标量)； X^p是矩阵X的p次幂(p为标量)，当p为整数时，可通过连续相乘来完成；当p为负整数时，则先对X求逆；当p为其它值时，幂函数计算要涉及到特征值和特征向量，这样当[V，D]=eig(X)时，有X^p=V*D.^p/V。 x^P是标量x的矩阵指数函数，这也要用到矩阵P的特征值和特征向量。当指数和底数均为矩阵时，无法求解。 a =[1, 2, 3;4, 5, 6;7,8,9] h =[2,1,2;3,6,8;6,8,8] 例: a =[1, 2, 3;4, 5, 6;7,8,9] g=[1,2,0;4,0,0;9,8,7] d=zeros(3，2，2)； ％产生3×2×2的零矩阵(多维阵列) d(:,:,1) = 0 0 0 0 0 0 d(:,:,2) = 0 0 0 0 0 0 d=rand(3,2,2) d(:,:,1) = 0.9218 0.4057 0.7382 0.9355 0.1763 0.9169 d(:,:,2) = 0.4103 0.3529 0.8936 0.8132 0.0579 0.0099 fliplr(diag(d)) ％非零元素位于反主对角线 ans = 0 0 2 0 ?10 0 8 0 0 利用下标修改矩阵中的个别元素是很方便的，例如： (1) 连接操作符[ ]： 像分块矩阵构造大矩阵一样，通过连接操作符[ ]将小矩阵连接成大矩阵。例如： a=fix(10*rand(2,4)) ％产生[0，10]之间均匀分布的随机矩阵(fix为取整函数) 　　首次使用cputime时，表示启动了CPU时间计数器，因此可利用cputime计算出某一程序段的执行时间。例如： t1=cputime; x=rand(2048,1); fft(x); t=cputime-t1 结果： t = 0.0310 12．calendar 功能：日历。 格式： c=calendar c=calendar(y,m) c=calendar(d) calendar(…) 　　说明： 　　c=calendar可得到一个6×7的矩阵，它包含当月的日历；c=calendar(d)可得到由d指定的日期所在月的日历，其中d可以是日期数值或日期字符串；c=calendar(y,m)可得到由y，m指定年份、月份的日历。 　　　calendar(…)可直接在屏幕上显示出日历。 2.5.4 矩阵操作 　　1．diag 功能：对角矩阵和矩阵的对角化。 格式： X=diag(V,k) %当V为n元向量时，X=diag(V,k)可得到n+abs(k)阶的方阵X，其V的元素处于第k条对角线上，k=0表示主对角线，k0表示在主对角线之上，k0表示在主对角线之下。 V=diag(X,k) %当X为矩阵时，V=diag(X,k)可得到列向量V，它取自于X的第k个对角线上的元素； X=diag(V) %产生V的元素处于主对角线的对角方阵。 V=diag(X) %V=diag(X)相当于