matlab4-线性代数问题求解.pptVIP

第四章 线性代数问题求解 矩阵 线性方程组的直接解法 线性方程组的迭代法 线性方程组的符号解法 稀疏矩阵技术 特征值与特征向量 4.1 矩阵4.1.1特殊矩阵的输入 数值矩阵的输入 零矩阵、幺矩阵及单位矩阵 生成n?n方阵：A=zeros(n), B=ones(n), C=eye(n) 生成m?n矩阵：A=zeros(m,n), B=ones(m,n), C=eye(m,n) 生成和矩阵B同样位数的矩阵：A=zeros(size(B)) 随机元素矩阵 若矩阵随机元素满足[0,1]区间上的均匀分布生成n?m阶标准均匀分布为随机数矩阵：A=rand(n,m)生成n?n阶标准均匀分布为随机数方阵：A=rand(n) 对角元素矩阵已知向量生成对角矩阵：A=diag(V) 已知矩阵提取对角元素列向量：V＝diag(A) 生成主对角线上第k条对角线为V的矩阵：A=diag(V,k) 例：diag( )函数的不同调用格式 C=[1 2 3]; V=diag(C) % 生成对角矩阵 V =100020003 V1=diag(V) % 将列向量通过转置变换成行向量 V1 =123 C=[1 2 3]; V=diag(C,2) % 主对角线上第 k条对角线为C的矩阵 V =0010000020000030000000000 生成三对角矩阵: V=diag([1 2 3 4])+diag([2 3 4],1)+diag([5 4 3],-1) V =1200523004340034 Hilbert矩阵及逆Hilbert矩阵生成n阶的Hilbert矩阵：A=hilb(n)求取逆Hilbert矩阵：B=invhilb(n) Hankel(汉克 ) 矩阵其中：第一列的各个元素定义为C向量，最后一行各个元素定义为R。H为对称阵。H1=hankel(C)由 Hankel 矩阵反对角线上元素相等得出一下三角阵均为零的Hankel 矩阵 Vandermonde(范德蒙)矩阵伴随矩阵其中：P(s)为首项系数为一的多向式。符号矩阵的输入数值矩阵A转换成符号矩阵：B=sym(A) 例： A=hilb(3) A =1.0000 0.5000 0.33330.5000 0.3333 0.25000.3333 0.2500 0.2000 B=sym(A) B = [ 1, 1/2, 1/3] [ 1/2, 1/3, 1/4] [ 1/3, 1/4, 1/5] 4.1.2 矩阵基本概念与性质 行列式格式 ：d=det(A) 例：求行列式 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; det(A) ans =0 例： tic, A=sym(hilb(20)); det(A), toc ans = 1/00000000000000000000000000000000000 elapsed_time =2.3140 高阶的Hilbert矩阵是接近奇异的矩阵。 矩阵的迹格式： t=trace(A) 矩阵的秩 格式：r=rank(A) ％用默认的精度求数值秩r=rank(A， ) ％给定精度下求数值秩矩阵的秩也表示该矩阵中行列式不等于0的子式的最大阶次。可证行秩和列秩（线性无关的）应相等。 例 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; rank(A) ans =3 该矩阵的秩为3，小于矩阵的阶次，故为非满秩矩阵。 例 H=hilb(20); rank(H) ％数值方法 ans = 13 H=sym(hilb(20)); rank(H) % 解析方法，原矩阵为非奇异矩阵 ans = 20 矩阵范数 矩阵的范数定义：格式：N=norm(A)％求解默认的2范数N=norm(A，选项)％选项可为1,2,inf等 例：求一向量、矩阵的范数 a=[16 2 3 13]; [norm(a), norm(a,2), norm(a,1), norm(a,Inf)] ans =2.0645635e+001 2.0645635e+001 3.0000000e+001 1.0000000e+001 A=[16 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1]; [norm(A), norm(A,2), norm(A,1), norm(A,Inf)] ans =34 34 34 34符号运算工具箱未提供norm( )函数，需先用double( )函数转换成双精