实验十一.线性代数问题MATLAB实现

实验内容

1.对矩阵进行LU分解。（写出求解的MATLAB程序，及其求解结果） 2.对矩阵：进行Cholesky分解。（写出求解的MATLAB程序，及其求解结果） 3.对矩阵进行QR分解。（写出求解的MATLAB程序，及其求解结果） 4.对矩阵进行奇异值分解。（写出求解的MATLAB程序，及其求解结果） 5.已知矩阵，求A的逆，及其行列式值。（写出求解的MATLAB程序，及其求解结果） 6.已知，求矩阵A的行最简阶梯型。（写出求解的MATLAB程序，及其求解结果） 7.已知方程组，其中，，判断其是否有解？若有，求出其基础解系。（写出求解的MATLAB程序，及其求解结果） 8.已知矩阵，求其特征多项式、特征值与特征向量。（写出求解的MATLAB程序，及其求解结果）

1.

>> a=[1,2,3,-4;-3,-4,-12,13;2,10,0,-3;4,14,9,-13];
>> [L,U]=lu(a)
>> L*U

2.

>> a=[1,2,1;2,8,4;1,4,6];
>> R=chol(a)
>> R'*R
>> L=chol(a,'lower')
>> L*L'

3.

>> a=[1,-3,1;-3,1,-1;1,-1,5];
>> [Q,R]=qr(a)
>> Q*R

4.

>> a=[1,1,1,1;1,2,3,4;1,3,6,10;1,4,10,20];
>> [U,S,V]=svd(a)
>> U*S*V'

5.

>> a=[3,1,-1,2;-5,1,3,-4;2,0,1,-1;1,-5,3,-3];
>> rank(a) %求A的逆
>> det(a) %求A的行列式值

6.

>> a=[1,2,0,2,5;-2,-5,1,-1,-8;0,-3,3,4,1;3,6,0,-7,2];
>> format rat
>> rref(a)

7.

>> A=[1,5,-1,-1;1,-2,1,3;3,8,-1,1;1,-9,3,7];
>> b=[-1,3,1,7]';
>> B=[A,b];
>> r1=rank(A)
>> r2=rank(B)

8.

>> A=[-2,1,1;0,2,0;-4,1,3];
>> [v,l]=eig(A) %v的列表示特征向量，矩阵l对角线上元素为与之对应的特征值
>> poly(A) %矩阵A的特征多项式

答案供参考,如有错误请自行改正，本人只是随意分享之前的作业，代码均来源于网络+个人稍微修改，认为有用可参考哈～