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,:);;;;;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 2 3;4 5 6;7 8 9;10 11 12],( b=[1 3 5 7]’) [Q,R]=qr(A),