数学建模中常用的线性代数基础知识在MATLAB中的实现

MATLAB中向量和矩阵的基本运算

设A,B两个矩阵

1.A+B，A-B
2.k*A
3.A*B
4.A\B 左除 A-1B，A必须为方阵
5.A/B 右除 AB-1，B必须为方阵
6.det(A) 求|A|，A必须为方阵
7.inv(A)或A-1
8.A^n
9.A’或transpose(A)
10.rank(A)
rref(A) 矩阵行变化化简，求矩阵A阶梯形的行最简形式
例题:

矩阵的变换与分解及其在MATLAB中的实现

矩阵的对角元素

函数diag将一个矩阵的对角元素提取出来

diag(A) 由矩阵A的对角线元素得到一个列向量

用该函数来产生第k阶对角线上的元素

diag(A,k)
其中,k=0表示主对角线;
k>0表示在主对角线以上;
k<0表示在主对角线一下

函数blkdiag可以根据输入创建一个分块对角矩阵

S=blkdiag(a,b,c,…)
功能：根据输入的a,b,c等参数构造一个分块对角矩阵

矩阵的分解

(1) 矩阵的奇异值分解 s=svd(A)
[U,S,V]=svd(A)
[U,S,V]=svd(A,0)
(2) 矩阵的LU分解 [L,U]=lu(A)
(2) 矩阵的QR分解 [Q,R]=qr(A)
(2) 矩阵的Cholesky分解 C=chol(A)

MATLAB中矩阵特征值和特征向量的求解方法

矩阵的特征值和特征向量

d=eig(A) 返回方阵A的全部特征值组成的特征向量d
[V,D]=eig(A) 返回矩阵A的特征值D与特征向量矩阵V,满足AV=VD
poly(A) 求矩阵A的特征多项式

注：求矩阵A的特征根 D=roots(poly(A)) 。

矩阵的相似对角化

线性方程组的直接求解法在MATLAB中的实现

齐次线性方程组的直接求解法

A=[1 1 1 1 1;3 2 1 1 -3;0 1 2 2 6;5 4 3 3 -1]; B=null(A,‘r’);

syms k1 k2 k3; x=k1B(:,1)+k2B(:,2)+k3*B(:,3)

非齐次线性方程组的直接求解法

求逆法

对AX=B，若|A|不等于0，则解由inv(A)*B获得

左除法和右除法

A=[2 3 1 4;1 -2 4 -5;3 8 -2 13; 4 -1 9 -6];
rref(A) %转化为矩阵简最简阶梯式

符号方程组的求解方法

线性方程组AX=B的符号解

X=linesolve(A,B) %只给出特解

syms a b c;
A=[a 0 0;0 b 0]; B=[1;c]
X=linesolve(A,B)

非线性方程组的符号姐

[x1,x2,x3,…]=solve(e1,e2,e3,…)

e1=sym(‘a+b+x=y’);
e2=sym(‘2ax-by=-1’);
e3=sym('2(a+b)=x+y’);
e4=sym(‘ay+bx=4’);
[a,b,x,y]=solve(e1,e2,e3,e4)