文档介绍:14-15-2 《几何与代数》数学实验报告学号: __________ 姓名: _________ 得分: ___________ 实验 1 求解线性方程组实验内容: 用 MATLAB 求解如下线性方程组 Ax=b, 其中: 1 1 0 0 0 0 0 0 1 2 1 0 0 0 0 0 0 1 3 1 0 0 0 0 0 0 1 4 1 0 0 0 0 0 0 1 5 1 0 0 0 0 0 0 1 6 1 0 0 0 0 0 0 1 7 1 0 0 0 0 0 0 1 8 A ? ?? ?? ?? ?? ?? ??? ?? ?? ?? ?? ?? ?? ? b=[02114102] T. 实验目的: 1. 了解 MATLAB 软件, 学会 MATLAB 软件的一些基本操作. 2. 熟悉 MATLAB 软件的一些数值计算功能. 3. 练****编写简单的 MATLAB 程序. 实验原理: 1. 对于齐次线性方程组 Ax =b, 根据 gramer 法则, 其解为 X i =D i /D; 2. 当A 可逆,即|A| ≠0 时, 方程组 Ax =b 的解为 X=inv(A)*B; 3. 当A 可逆,即|A| ≠0 时, 对增广矩阵[A,b] 进行初等行变换, 把它化为最简形, 则原矩阵 b 对应的地方即为 A -1 *b , 取最简形的最后一列赋给 X 即得线性方程组的解。实验方案一:利用 Cramer 法则求解: 1. 在 MATLAB 命令窗口中输入如下命令: >> a=[1;1;0;0;0;0;0;0];b=[1;2;1;0;0;0;0;0];c=[0;1;3;1;0;0;0;0];d=[0;0;1;4;1;0;0;0];e=[0;0;0;1;5;1;0;0]; f=[0;0;0;0;1;6;1;0];g=[0;0;0;0;0;1;7;1];h=[0;0;0;0;0;0;1;8];i=[0;2;1;1;4;1;0;2]; % 输入矩阵 A 与i >> D=det([a,b,c,d,e,f,g,h]), % 把行列式的值求出来赋给 D D= 10241 >> D1=det([i,b,c,d,e,f,g,h]),D2=det([a,i,c,d,e,f,g,h]),D3=det([a,b,i,d,e,f,g,h]),D4=det([a,b,c,i,e,f,g,h]), D5=det([a,b,c,d,i,f,g,h]),D6=det([a,b,c,d,e,i,g,h]),D7=