【文章】**线性方程组求解在Excel和Matlab中的实现**

[spy1120]

一、用Excel求解线性方程组

　　 x1+x2+2x3+3x4=1
　　3x1-x2-x3-2x4=-4
　　2x1+3x2-x3-x4=-6
　　x1+2x2+3x3-x4=-4
　　可按如下的步骤来解这个方程组：
　　1.打开Excel。
　　2.由于在本方程组中未知数有4个，所以预留4个可变单元格的位置A1-A4。
　　3.将活动单元格移至B1处，从键盘键入：=A1＋A2＋2＊A3＋3＊A4：然后回车(此时B1显示0)。即在B1处输入方程组中第一个方程等号左边的表达式。
　　4.在B2处从键盘键入：=3＊A1－A2－A3－2＊A4；然后回车(此时B2显示0)。即在B2处输入方程组中第二个方程等号左边的表达式。
　　5.在B3处从键盘键入：=2＊A1＋3＊A2－A3－A4；然后回车(此时B3显示0)。即在B3处输入方程组中第三个方程等号左边的表达式。
　　6.在B4处从键盘键入：=A1＋2＊A2＋3＊A3－A4；然后回车(此时B4显示0)。即在B4处输入方程组中第四个方程等号左边的表达式。
　　7.点击工具?规划求解，出现规划求解参数对话框。
　　8.对话框中第一栏为：设置目标单元格，在相应的框中填入$B$1。
　　9.对话框中第二栏为：等于；后有三个选项，依次为最大值，最小值，值为。根据题意B1表示方程组中第一个方程等号左边的表达式，它的值应为1，因此点击值为前的圆圈，输入1。
　　10.对话框中第三栏为：可变单元格；我们预留的可变单元格为A1-A4，所以在可变单元格框内键入 A 1： A 4。
　　11.对话框中最后一栏为：约束；首先点击添加按钮，屏幕出现添加约束对话框。
　　12.在添加约束对话框的单元格引用位置键入：B2；在中间的下拉式菜单中选取=；在约束值处键入：－4；然后按添加按钮，屏幕出现空白的添加约束对话框。
　　13.在添加约束对话框的单元格引用位置键入：B3；在中间的下拉式菜单中选取=；在约束值处键入：－6；然后按添加按钮，屏幕出现空白的添加约束对话框。
　　14.在添加约束对话框的单元格引用位置键入：B4；在中间的下拉式菜单中选取=；在约束值处键入：－4；然后按确定键，返回规划求解参数对话框。特别注意在最后一个约束条件键入后，按确定键(而不是像前面一样按添加键)。
　　15.按求解键，出现求解结果对话框。此时在A1－A4的位置依次为：－1，－1，0，1；这就是说，原方程组的解为：X1=－1,X2=－1,X3=0,X4=1。这样我们就求出了方程组的解。

二、用Matlab求解线性方程组

MATLAB中线性方程组的求解用矩阵除法计算，只要输入系数矩阵A和常数矩阵B，方程组就成为矩阵方程AX=B，则该矩阵方程的解为X=A\B。
实例、求线性方程组 的解。
由线性方程组可知：
系统矩阵为 ，常数矩阵为 。操作如下：
>> A=[2 4 6;9 6 3;4 3 7]; %输入系数矩阵A
>> B=[7;9;6]; %输入常数矩阵B
>> X=A\B %求未知矩阵X
X =
0.0250
1.3250
0.2750
如果系数矩阵A是非奇异方阵，也可以先用inv命令求出A的逆阵NA，再用X=NA*B求出方程组的解。
>> NA=inv(A) %求矩阵A的逆阵
NA =
-0.2750 0.0833 0.2000
0.4250 0.0833 -0.4000
-0.0250 -0.0833 0.2000
>> X=NA*B %求未知矩阵X
X =
0.0250
1.3250
0.2750