gramer法则matlab,东南大学几何与代数matlab实验报告

东南大学 几何与代数实验报告(利用matlab软件解决实际数学问题,几何与代数作业之一)

数学实验报告

学号: , 姓名: 吴雪松 , 得分:

实验1 求解线性方程组

实验内容: 用MATLAB求解如下线性方程组Ax = b, 其中

5 1 0 0A =

0 0 0 0

651000000651000000651000000651000000651000000651

0 0 0 0

, b = [09009320] T. 0 0 6 5

实验目的: 1. 了解MATLAB软件, 学会MATLAB软件的一些基本操作；

2. 熟悉MATLAB软件的一些数值计算功能，如在矩阵方面的运算； 3. 练习编写简单的MATLAB程序。

实验原理: 1. 对于齐次线性方程组Ax = b, 根据gramer法则, 其解为Xi=Di/D;

2. 当A可逆，即|A|≠0时, 方程组Ax = b的解为X=inv(A)*B;

3. 当……时, 对增广矩阵[A, b]进行初等行变换, 把它化为最简形, 则原矩阵b对应的地方即为A-1*b，取最简形的最后一列赋给X即得线性方程组的解。

实验方案: 1. 在MATLAB命令窗口中输入如下命令:

A=[5 6 0 0 0 0 0 0; 1 5 6 0 0 0 0 0; 0 1 5 6 0 0 0 0; 0 0 1 5 6 0 0 0; 0 0 0 1 5 6 0 0; 0 0 0 0 1 5 6 0; 0 0 0 0 0 1 5 6; 0 0 0 0 0 0 1 5;];

b=[0 9 0 0 9 3 2 0]'; %输入矩阵A与b X=[]; %声明一个名为X的空矩阵

D=det(A); %把A的行列式的值求出来赋给D for i=1:8 %开始循环语句 A=[5 6 0 0 0 0 0 0; 1 5 6 0 0 0 0 0; 0 1 5 6 0 0 0 0; 0 0 1 5 6 0 0 0;