课件第二matlab的数值计算.ppt

a ^ p —— a 自乘p次幂 方阵 1的整数 3. 矩阵乘方—— a^n,a^p,p^a 对于p的其它值,计算将涉及特征值 和特征向量，如果p是矩阵，a是标量 a^p使用特征值和特征向量自乘到p次 幂；如a,p都是矩阵，a^p则无意义。 a=[1,2,3;4,5,6;7,8,9];a^2 ans =30 36 42 66 81 96 102 126 150 ※当一个方阵有复数特征值或负实特征值时，非整数幂是复数阵。 a^0.5 ans = 0.4498 + 0.7623i 0.5526 + 0.2068i 0.6555 -0.3487i 1.0185 + 0.0842i 1.2515 + 0.0228i 1.4844 - 0.0385i 1.5873 - 0.5940i 1.9503 - 0.1611i 2.3134 + 0.2717i inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 diag —— 对角矩阵 ’ —— 矩阵转置 sqrt —— 矩阵开方 4. 矩阵的其它运算 5. 矩阵的范数 矩阵范数的函数为：(1) norm(V)或norm(V,2)：计算矩阵V的 2—范数。(2) norm(V,1)：计算矩阵V的1—范数。(3) norm(V,inf)：计算矩阵V的∞—范数。 6.矩阵的一些特殊操作 矩阵的变维 a=[1:12];b=reshape(a,3,4) c=zeros(3,4);c(:)=a(:) 矩阵的变向 rot90:旋转; fliplr:左右翻; flipud:上下翻 矩阵的抽取 diag:抽取主对角线;(对于非方阵的情况?) tril: 抽取主下三角; triu:抽取主上三角，然后其余补零元素 矩阵的扩展 关系运算 关系符号 意义 = = == ~= 小于 小于或等于 大于 大于或等于 等于 不等于 关系运算符的运算法则 当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为1，否则为0。 当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 (3) 当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 注意：其书写方法与数学中的不等式符号不尽相同。 数组运算指元素对元素的算术运算， 与通常意义上的由符号表示的线性代数 矩阵运算不同。 数组加减(.+,.-) a.+b a.- b 7. 矩阵的数组运算 对应元素相加减（与矩阵加减等效） 2. 数组乘除(??，./，.\） a??b —— a，b两数组必须有相同的行 和列两数组相应元素相乘。 a=[1 2 3;4 5 6;7 8 9]; b=[2 4 6;1 3 5;7 9 10]; a.*b ans = 2 8 18 4 15 30 49 72 90 a=[1 2 3;4 5 6;7 8 9]; b=[2 4 6;1 3 5;7 9 10]; a*b ans = 25 37 46 55 85 109 85 133 172 a./b=b.\a a.\b=b./a a./b=b.\a — 都是a的元素被b的对应元 素除， “/”是斜杠 a.\b=b./a — 都是b的元素被a的对应元 素除， “\”是反斜杠 例: a=[1 2 3];b=[4 5 6]; c1=a.\b; c2=b./a c1 = 4.0000 2.5000 2.0000 c2 = 4.0000 2.5000 2.0000 —— 给出a,b对应元素