Matlab数值代数.pptVIP

数学软件 Matlab 数值代数简介 Some numerical software 线性方程组求解模块 特征值求解模块 矩阵市场： James Demmel 数值代数常用网站 * 特殊矩阵 Matlab 提供了一些函数，用于生成一些常见的特殊矩阵 例： A=magic(3) B=ones(4) 特殊矩阵的生成 A=magic(5) a1=diag(A) a2=diag(A,1) a3=diag(A,-2) a=[1,2,3,4] A1=diag(a) A2=diag(a,1) A3=diag(a,-2) 注意 diag 的用法 * 常见矩阵生成函数 zeros(m,n) zeros(n) 生成一个 m 行 n 列的零矩阵 m=n 时可简写为 zeros(n) ones(m,n) ones(n) 生成一个 m 行 n 列的元素全为 1 的矩阵 m=n 时可简写为 ones(n) eye(m,n) eye(n) 生成一个主对角线全为 1 的 m 行 n 列矩阵 m=n 时可简写为 eye(n)，即为 n 维单位矩阵 diag(X) diag(X,k) 若 X 是矩阵，则 diag(X) 为 X 的主对角线向量 若 X 是向量，diag(X) 产生以 X 为主对角线的对角矩阵 tril(A) 提取一个矩阵的下三角部分 triu(A) 提取一个矩阵的上三角部分 rand(m,n) rand(n) 产生 0～1 间均匀分布的随机矩阵 m=n 时简写为 rand(n) randn(m,n) randn(n) 产生均值为0，方差为1的标准正态分布随机矩阵 m=n 时简写为 randn(n) 其它特殊矩阵生成函数：magic、hilb、pascal 等 toeplitz，hankel T = toeplitz(r): Hermitian Toeplitz matrix formed from vector r. T = toeplitz(c,r):Toeplitz matrix T having c as its first column and r as its first row. (注意c和r维数一样，且第一个元素相同（不同的话，以列为主）) 类似地：H = hankel(c)；H = hankel(c,r) 例如： c = [1 2 3 4 5]; r = [1.5 2.5 3.5 4.5 5.5]; toeplitz(c,r) H = hankel(c) * 翻转与旋转 矩阵的翻转与旋转 例： 注意矩阵旋转与转置的区别！ fliplr(A) 左右翻转 flipud(A) 上下翻转 rot90(A) rot90(A,k) 逆时针旋转 90 度 逆时针旋转 k×90 度 A = [1 2 3;4 5 6] B = fliplr(A) C = flipud(A) D = rot90(A) E = rot90(A,-1) * 查看矩阵的大小 查看矩阵的大小：size、length size(A) 返回矩阵 A 的行数和列数 size(A,1) 返回矩阵 A 的行数 size(A,2) 返回矩阵 A 的列数 length(x) length(X) 若 x 是向量，则返回 x 的长度 若 X 是矩阵，则返回行数和列数中大的一个 numel(A) 返回 A 的元素的个数 A=[1,2,3; 4,5,6] size(A) size(A,1) size(A,2) length(A), numel(A) x=1:5; length(x) 例： * 改变矩阵的形状 改变矩阵的形状：reshape(A,m,n) 将矩阵元素按 列方向 进行重新排列成一个 m×n 的新矩阵 新矩阵的元素个数必须与原矩阵元素个数相等！ A=[1,2,3; 4,5,6; 7,8,9; 10,11,12] B=reshape(A,4,3) 例： C=reshape(A,2,6) C=reshape(A,5,2) × * 矩阵的超越函数 Matlab 提供的矩阵函数：expm、sqrtm、logm 详情参见联机帮助（help expm / sqrtm / logm ） 更一般的矩阵函数： funm funm(A, @fun) fun 可以是 exp, log，cos，sin，... ... eA =1+A+A^2/2!+...+A^n/n!+... 矩阵函数 将矩阵作为一个整体参与运算 怎样计算 eA ? lu分解 [L,U,P] = lu(A)即LU=PA L单位下三角矩阵，U为上三角矩阵，P为置换矩阵 [L,U] = lu(A)：L is either a unit lower tr