MATLAB实验报告2 第九讲+线性方程组.pdf

统计与数学模型分析实验中心

实验报告

课程名称：经济数学实验与建模

课程号：F108003

实验班级：

指导教师：

学期：2010－2011(1)

实验室：

2011年月日

安徽财经大学统计与数学模型分析实验中心

《经济数学实验与建模》实验报告

班级：091投资1班学号：2009200025姓名：檀静实验时间5月21日实验地点

实验楼505教室

实验名称：

实学习矩阵的运算、线性方程组的解法以及线性代数在经济管理方面的应用。

验（1）掌握矩阵的加法、数乘、乘法、乘方、行列式的运算。

目（2）掌握矩阵的初等行（或列）变换,求矩阵的秩和逆矩阵。

的（3）理解矩阵的LU、QR、SVD分解。

1．求线性齐次方程组的通解

例4-6求下面方程组的通解。

xx2xxx0,

12345

3xxx4x3x0,

12345

x5x9x8xx0.



12345

解：在命令窗口中输入以下命令。

A=[11-2-11;3-1143;15-9-81];%系数矩阵

formatrat%指定以有理式格式输出

C=null(A,r)%求解空间的有理基向量

运行后显示结果如下。

实

C=

验1/4-3/4-1

7/47/40

内

100

容010

001

从而可得原方程组的基础解系为：

(1/4,7/4,1,0,0)T,(3/4,7/4,0,1,0)T,(1,0,0,0,1)T

123

写出方程通解，命令如下。

symsk1k2k3%k1，k2，k3是方程组通解中的任意常数

X=k1*C(:,1)+k2*C(:,2)+k3*C(:,3)%写出方程组的通解

X=

[1/4*k1-3/4*k2-k3]

[7/4*k1+7/4*k2]

[k1]

[k2]

[k3]

2求非齐次线性方程组的通解

xx2xxx1,

12345

例4-7求解方程组3xxx4x3x