第二章MATLAB数值计算概要.ppt

例. 求矩阵A=[1,2,3;4,5,6;7,8,9]的特征值和特征向量。 所得结果为 解：键入并运行如下命令： A=[1,2,3;4,5,6;7,8,9]; E=eig(A),[V,D]=eig(A) E = 16.1168 -1.1168 -0.0000 V = -0.2320 -0.7858 0.4082 -0.5253 -0.0868 -0.8165 -0.8187 0.6123 0.4082 D = 16.1168 0 0 0 -1.1168 0 0 0 0 矩阵的一些特殊操作 1. 矩阵元素 MATLAB允许用户利用矩阵元素的序号对一个矩阵的元素进行提取和修改。 矩阵元素的序号标识方式为(双下标形式) 格式：A(p,q) 矩阵A的第p行第q列元素 矩阵元素的序号标识方式也可以用一个数字表示(单下标形式) ，这时将矩阵元素按列编号，先第一列，再第二列，依次类推。若A为m×n矩阵，元素A(p,q)的序号为(q-1)*m+p。 例. 运行如下MATLAB命令： A=[1 2 3;4 5 6;7 8 9] B1=A(3,2), B2=A(8) A(3,2)=0, A(8)=11 所得结果为 B1 = 8 B2 = 6 A = 1 2 3 4 5 6 7 0 9 A = 1 2 3 4 5 6 7 8 9 A = 1 2 3 4 5 11 7 0 9 2. 子矩阵 格式：A (u,v) u，v都是向量或冒号表达式，返回值为矩阵A的由u确定的行和由v确定的列所有元素构成的矩阵。 例如：u=[2,3,5]，v=[2,4]，则A(u,v)表示矩阵A的第2,3,5行中第2,4列构成的3×2矩阵。 矩阵的一些特殊操作 例. 运行如下MATLAB命令： A=[1:5;6:10;11:15;16:20;21:25;26:30] u=[2,3,5];v=[2,4]; B=A(u,v) 所得结果为 A = 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 B = 7 9 12 14 22 24 7 9 12 14 22 24 例：利用冒号表达式获得子矩阵 A(:,q)：矩阵A第q列全部元素构成的子矩阵 A(p,:)：矩阵A第p行的全部元素成的子矩阵； A(p:p+k,:)：矩阵A第p～p+k行的全部元素成的子矩阵； A(:,q:q+k)：矩阵A第q～q+k列的全部元素成的子矩阵； A(p:p+k,q:q+m)：矩阵A第p～p+k行内，并在第q～q+m列中的所有元素成的子矩阵。 2.2 多项式及其运算 多项式的表达方法 多项式运算 多项式的表达方法 (1) 多项式系数行向量的直接输入 多项式 a0xn+ a1xn-1+…+ an-1x+ an可用其系数行向量表示，即[a0 a1 … an-1 an] 注意：通过命令poly2str可以把多项式的向量表示转化成多项式的字符串表示 例：aa=[1,2,3,4]; p=poly2str(aa,x) 运行结果为 p = x^3 + 2 x^2 + 3 x + 4 (2) 利用方阵或向量生成多项式命令poly ① 利用方阵生成多项式 格式：P=poly(A) 其中A是一个n阶方阵，返回值P为n+1维行向量，P的元素是A的特征多项式（按降幂排列）的系数 例：A=[1,2;3,4]; Q=poly(A) 运行结果为： Q = 1 -5 -2 矩阵A的特征多项式为 多项式的表达方法 ② 利用向量生成多项式 例：v=[1,2,3,4]; S=poly(v) 运行结果为： S = 1 -10 35 -50 24 格式：P=poly(v) 其中v是一个n维向量，返回值P为n+1维行向量，P的元素是以向量v的元素为根且最高次项系数为1的多项式(按降