MATLAB数值运算讲义.ppt

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对应元素间的商. 3. 数组乘方(.^) — 元素对元素的幂 例: a=[1 2 3];b=[4 5 6]; z=a.^2 z = 1.00 4.00 9.00 z=a.^b z = 1.00 32.00 729.00 (1 .^4 2 .^5 3 .^6) matlab语言把多项式表达成一个行向量， 该向量中的元素是按多项式降幂排列的。 f(x)=a0xn+a1xn-1+……an-1x+an 可用行向量 p=[a0 a1 …… an-1 an ]表示 poly —— 产生特征多项式系数向量 特征多项式一定是n+1维的 特征多项式第一个元素一定是1 四、 多项式运算 例:a=[1 2 3;4 5 6;7 8 0]; p=poly(a) p =1.00 -6.00 -72.00 -27.00 p是多项式p(x)=x3-6x2-72x-27的系数 matlab描述方法，我们可用： p1=poly2str(p,‘x’) — 函数文件，显示 数学多项式的形式 p1 =x^3 - 6 x^2 - 72 x – 27 注意：多项式中缺少的幂次用‘0’补齐。