MATLAB实验报告(二)矩阵代数.pdfVIP

实验名称：矩阵代数

姓名：学号：实验日期：

实验目的：

1、学习MATLAB的有关线性代数运算指令；

2、熟悉矩阵代数的MATLAB指令；

3、熟悉如何应用指令解决计算和建模实验问题。

实验项目：

1、学会应用矩阵代数的MATLAB指令解决问题；

2、通过实验解决矩阵除法、线性方程组的通解和矩相似对角化的问题；

3、通过建模实验研究投入产出分析、基因遗产等应用问题。

实验具体过程：

1、题目：用矩阵除法解如下线性方程组，并判断解的意义。

A=[],

4-33;32-6;1-53b=[-1;-2;1],x=[x1;x2;x3]，解Ax=b

解题思路：

（1）求出系数矩阵和增广矩阵的秩，判断解的个数。

（2）若秩（A）=秩(A，b)=n，存在唯一解；其解直接由x=A求得；若秩（A）

=秩(A，b)n，存在无穷多解，其解可表示为对应齐次方程Ax=0的一个

基础解系和Ax=b的一个特解x=A的叠加；若秩（A）≠秩(A，b)，则

无解，从而得到最小二乘近似解x=A。

对实验题目的解答：

A=[4-33;32-6;1-53];b=[-1;-2;1];

[rank(A),rank([A,b])]

ans=

33%秩相等且等于3，说明有唯一解

x=A

x=

-0.4706

-0.2941

0

2、题目：用矩阵除法解如下线性方程组，并判断解的意义。

A=[41;32;1-5];b=[1;1;1];x=x=[x1;x2],解Ax=b

对实验题目的解答：

A=[41;32;1-5];b=[1;1;1];

[rank(A),rank([A,b])]

ans=

23%秩不相等，即无解

x=A

x=%求得一最小二乘近似解

0.3311

-0.1219

3、题目：用矩阵除法解如下线性方程组，并判断解的意义。

A=[21-11;121-1;1121];b=[1;2;3];x=[x1;x2;x3;x4],解Ax=b

1

对实验题目的解答：

A=[21-11;121-1;1121];b=[1;2;3];

[rank(A),rank([A,b])]

ans=

33%秩相等且小于4，说明有无穷多解

x=A

x=%求得一特解

1.0000

0

1.0000

0.0000

x0=null(A)

x0=%得到Ax=0的基础解系

-0.6255

0.6255

-0.2085

0.4170

求得通解为kx0+x

改进或思考：

A=[21-11;121-1;1121];b=[1;2;3];

[rank(A),rank([A,b])]

ans=

33%秩相等且小于4，说明有无穷多解

rref([A,b])

ans=

1.0000001.50001.0000

01.00000-1.50000

001.00000.50001.0000

从而原方程等价于x1+1.5x4=1,x2-1.5x4=0,x3+0.5x4=