如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
1.矩阵的的生成MATLAB在线性代数计算中的常用命令例1a=[1,2,3;4,5,6;7,8,10]a=1234567810(2)向量的冒号定义:a:d:b形式的语句生成一个行向量,范围在a与b之间,a为第一个元素,d为间隔,d的取值不能为0。(3)语句定义zeros(m,n)产生m行n列的元素全为0的矩阵;ones(m,n)产生m行n列的元素全为1的矩阵;eye(n)产生n阶单位矩阵;diag(u)产生一个对角矩阵,其对角线元素与向量u的元素一致.2.矩阵的元素操作MATLAB利用下标访问矩阵的元素。运行以上语句得到的结果为:a(:,4)=[3,2,1]%将矩阵a添上第4列b5=reshape(a,2,6)%将a重写为2行6列的矩阵c=find(b3==8)%求b3中等于8的元素的位置%将a的所有列按照从左到右的次序排列,b4求第7个元素。3.矩阵的行列式与转置a2=[1,2,3;4,5,6;7,8,9];a20=det(a2)在MATLAB中矩阵的转置阵的调用函数是A'。4.矩阵的加、减、乘与乘方运算(2)矩阵A与B相乘C=A*B,其结果与代数中矩阵相乘也是一致的,要求A的列数等于B的行数。例9已知矩阵a=[1,2,3;4,5,6;7,8,0],b=[1,2,2;3,4,1;5,6,2],试求a.*b。例10已知矩阵a=[1,2,3;4,5,6;7,8,0],试求a3。a=[123;456;780];a3=a*a*a;a3在MATLAB中,对矩阵还有另一种乘方。若A为一个一般的矩阵或向量,A.^n表示A的每个元素求n次方。5.矩阵的初等变换与矩阵的秩clearA=[1,1,1;1,2,-5;2,3,4];rref(A)rank(A)6.求向量组的极大线性无关组[Rjb]=rref(A)A(:,jb)7.矩阵的求逆运算A=[-16-4-6;15-39;1809];V=inv(A)它警告说:“此矩阵接近奇异,数据尺度很差,结果可能不准确。”矩阵方程组AX=B以及XA=B的解可以分别用A\B与A/B来表示。A=[303-6;5-11-5;-314-9;1-34-4];运行各语句得到的结果都为:a=[1-1;11];b=[5;1];inv(a)*bA=[1-1;11];b=[5,1];b*inv(A)9.求解一般线性方程组rref(A)clearA=[1-111;1-11-2;1-1-21];formatratn=4;RA=rank(A)symskX=k*bclearA=[-1-24;211;11-1];rank(A)(2)非齐次线性方程组的求解10.求过渡矩阵clearA=[111;100;1-11];B=[0-11;112;101];K=inv(A)*B,X=[1;2;-1],Y=inv(K)*X11.向量的正交规范化在线性代数中,矩阵的特征值与特征向量的计算是一个比较复杂繁琐的过程。clearA=[5-20;-262;027];[P,R]=eig(A)P=2/32/3-1/32/3-1/32/3-1/32/32/313.实对称矩阵的对角化B=inv(P)*A*P14.计算标准化二次型的转换矩阵x=P*y15.判别二次型的正定性