matlab第一章MATLAB系统及语言简介.pptVIP

* * 1.5 矩阵运算与数组运算 1.5.1 矩阵的创建函数 创建向量对象 例如：x=1:5; y=0:pi/4:2*pi; z=6:-1:0; x=(0.0:0.2:3.0)’; y=exp(-x).*sin(x); [x y] 能创建向量的其它函数： logspace 函数，创建对数等距的向量 y=logspace(d1,d2) y=logspace(d1,d2,N) linspace 函数，创建指定长度的等距向量 k=linspace(-pi,pi,4) 创建特殊矩阵 MATLAB提供用于创建某些特殊矩阵的函数 zeros 元素全为零的矩阵 ones 元素全为 1 的矩阵 rand 元素为一致分布的随机矩阵 rand(4,3) randn 元素为正态分布的随机矩阵 eye 单位矩阵 compan 多项式的伴随矩阵 diag 向量的对角矩阵 gallery 检验矩阵 1.5.2 矩阵的角标 角标：矩阵元素在矩阵中的行列号，A(i,j) 矩阵角标可以是一个向量，如果 x 和 v都是长度为 n 的向量，则 x(v) 表示新的向量 [x(v(1)), x(v(2)), … , x(v(n))]。 若 A 是 10*10 的矩阵 语句 A(1:5,3) A(1:5, 7:10) A(:,3) A(1:5,:) A(:,[3 5 10])=B(:,1:3) A(:,n:-1:1) A(:) 矩阵A的元素按列排列为一个列向量 设A是3*2的矩阵， A(:)=11:16 矩阵的0-1角标 以0-1为元素的向量作为矩阵的角标向量 设A是一个m*n阶矩阵，L是一个长度为m的0-1向量，则 A(L,:) 表示由向量L的非零元素对应的矩阵A的行组成的子块 举例 x=[45,78,90;32,43,50;60,24,84;96,78,35] L=x(:,2)42; a=x(L,:) 设A是3*2的矩阵， A(:)=11:16； A矩阵为 11 14 12 15 13 16 median 矩阵各列的中值向量 std 标准偏差向量 sum 元素之和向量 prod 元素之积向量 cumsum 累加和向量 cumprod 累积向量 diff 差分向量 hist 直方图向量 cov 列向量卷积 sort 排序 corrcoef 矩阵相关系数 MATLAB5.0及以上版本中，上述函数（除corrcoef 之外）都可以应用于多维数组。 例如：sum(M,d) M是一个n维数组 M=zeros(3,3,4); for k=1:4 M(:,:,k)=k*ones(3,3); end sum(M,1); sum(M,2); sum(M,3) 1.5.4 线性代数 1. 矩阵分解和线性方程求解 （1）三角分解 A=L?U 是最基本的矩阵分解方法，是求解线性方程的基本算法，L: 下三角矩阵 U：上三角矩阵 lu(A) 在理论上L应当为对角线元素为1的下三角矩阵 调用命令：[L,U]=lu(A) L经过重新排列后才可以是一个下三角矩阵 调用命令：[L,U,P]=lu(A)， P ? A=L ? U P：置换阵 矩阵的除法以LU分解为基础（举例） （2）正交分解 QR分解 将一个方阵或长方阵表示为一个正交矩阵和个三角矩阵的乘积 A=Q?R 调用命令 [Q,R]=qr(A) 举例 A=[1