第三章 矩阵及其运算.pptVIP

3.3 特殊矩阵 3.4 矩阵基本运算_1 3.5 矩阵高级运算_1 3.6 求解线性方程组_1 3.6 求解线性方程组_2 3.6 求解线性方程组_2 第3章 矩阵线性代数算法实现 矩阵的生成 矩阵的修改 特殊矩阵 矩阵基本运算 矩阵高级运算 求解线性方程组 加/减/乘 矩阵有数值/符号/特殊矩阵 生成实数数值矩阵方法: (1) 由命令窗口直接输入 同一行用，或空格分隔(个数不限) 不同行用；分隔或分行输入； 所有元素置于一[ ]内。 %例3-1 x = [1 2 3 4;2 3 4 5;3 4 5 6] (2) 由m文件生成 调用时 run ‘…’ (3) 由文本文件生成 txt文件不含变量名称 文件名为矩阵变量名 每行数值个数必须相等 调用：load ‘d:\...’ 部分扩充 D=[A;B C] 部分删除 A(:,n)=[];A(m,:)=[] 部分修改 A(m,:)=[a b …]; A(:,n)=[ a b …] 结构改变 左右翻转 fliplr(A) 上下翻转 flipud(A) 逆时针旋转rot90(A,k); 按指定维数翻转矩阵 flipdim(A,dim) 平铺矩阵 B=repmat(A,m,n) 矩阵的变维 B(:)=A(:); B = reshape(A,m,n) 矩阵数据变换 取整数 floor; ceil; round; fix 有理数[n,d]=rat(A) 余数 B=rem (A, x) 常用特殊矩阵函数 特殊矩阵的生成方法 1 单位阵 eye 2 1矩阵 ones 3 零矩阵 zeros 4 随机阵 randn 5 魔方阵 magic 6 对角阵 diag 7 三角阵 triu 8 Hilbert阵 hilb 9 托普利茲阵 toeplitz 除/乘方 逆/范数/条件数/秩 除法运算 1 左除(\),右除(/),点除 B./A x=A\b是方程A*x =b的解 x=b/A是方程x*A=b的解 A=[1 0 3;4 13 6;7 4 9]; b=[4;7;1];C=A\b 乘方运算 1 矩阵乘方^ 2 矩阵的数量乘方 .^ 转置/方阵/矩阵函数 矩阵函数 1 方阵的指数 expm(A) 2 矩阵的对数 B=logm(A)； 3 方阵的函数 F = funm(A,fun) 4 矩阵的方根 X = sqrtm(A) 5 矩阵A的多项式 polyvalm(P, A) 矩阵转置 ′ 方阵的运算 1 方阵行列式 d=det(A) 2 方阵的迹 trace 矩阵的逆与伪逆 1 方阵的逆矩阵 inv(A) 2 方阵的伪逆矩阵 pinv(A) 矩阵和向量的范数 1 向量的范数 norm(X) 2 矩阵的范数 norm(A) 矩阵的条件数 cond(A) 矩阵的秩 rank (A) 矩阵元素个数 numel(A) 矩阵分解 1 Cholesky分解 chol(X) 2 LU分解 lu(X) 3 QR分解 qr(A) 4 schur分解 schur(A) 5 实Schur分解转化成复Schur [U,T]=rsf2csf(u,t) 6 特征值分解 eig(A) 7 奇异值分解 svd (X) 8 特征值问题的QZ分解 qz(A,B) 9 海森伯格形式的分解 hess(A) 线性方程组一般求解可分为两类： 求方程组唯一解(特解) 求方程组无穷解(通解) 通过系数矩阵的秩r(rank)判断: (n为未知变量个数)） r =n，有唯一解； rn，可有无穷解(只给1个解。 rn，给出LSM意义上的解。 线性方程组无穷解 = 齐次方程组通解+ 非齐次方程组1个特解 一般方法 特殊方法 唯一解解法 1 矩阵除法解法 AX=b =X=A\b (方法1) 2 矩阵LU、QR和cholesky分解解法 齐次线性方程组通解解法 非齐次线性方程组通解的解法 1,一般: AX=b通解: AX=0通解+ AX=b特解