(完整)高斯消元法MATLAB实现(2)_第1页
(完整)高斯消元法MATLAB实现(2)_第2页
(完整)高斯消元法MATLAB实现(2)_第3页
(完整)高斯消元法MATLAB实现(2)_第4页
免费预览已结束,剩余1页可下载查看

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、数值分析实验报告一、实验目的与要求1掌握高斯消去法的基本思路和迭代步骤;2培养编程与上机调试能力。二、实验内容1 编写用高斯消元法解线性方程组的 MATLA程序,并求解下面的线性方程组, 然后用逆矩阵解方程组的方法验证 .0.101x12.304x23.555 x31.1835x12x2x381) 1.347 x13.712x24.623x32.137(2) 2x18x23x3212.835x11.072 x25.643x33.035x13x26x31MATLA程序,并求解下面的线性方2编写用列主元高斯消元法解线性方程组的 程组,然后用逆矩阵解方程组的方法验证 .0.101x12.304x23

2、.555 x31.1835x12x2x381) 1.347 x13.712x24.623x32.137(2) 2x18x23x3212.835x11.072 x25.643x33.035x13x26x31三.MATLAB计算源程序1.用高斯消元法解线性方程组AX b的MATLAB程序输入的量:系数矩阵 A 和常系数向量 b ;输出的量:系数矩阵A和增广矩阵B的秩RA,RB方程组中未知量的个数n 和有关方程组解 X 及其解的信息 .function RA,RB,n,X= gaus(A,b)B=A b; n=length(b); RA=rank(A);RB=rank(B);zhica=RB-RA;

3、if zhica0, disp( 请注意:因为 RA=RB ,所以此方程组无解 .) returnendif RA=RBif RA=ndisp( 请注意:因为 RA=RB=n ,所以此方程组有唯一解 .) X=zeros(n,1); C=zeros(1,n+1);for p= 1:n-1for k=p+1:nm= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n);for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:

4、n)*X(q+1:n)/A(q,q);endelsedisp( 请注意:因为 RA=RB0,disp( 请注意:因为 RA=RB ,所以此方程组无解 .)returnendif RA=RB if RA=n disp( 请注意:因为 RA=RB=n ,所以此方程组有唯一解 .) X=zeros(n,1); C=zeros(1,n+1);for p= 1:n-1Y,j=max(abs(B(p:n,p); C=B(p,:); B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;for k=p+1:nm= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m*

5、 B(p,p:n+1); endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n);for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)/A(q,q);endelsedisp( 请注意:因为 RA=RBn ,所以此方程组有无穷多解 .) endend实验过程:1( 1)编写高斯消元法的 MATLAB 文件如下: clear;A=0.101 2.304 3.555;-1.347 3.712 4.623;-2.835 1.072 5.643; b=1.183;2.137;3.035;RA,RB,n,X =ga

6、us (A,b)运行结果为:请注意:因为 RA=RB=n ,所以此方程组有唯一解 . RA =3RB =3 n =3 X =-0.39820.01380.3351( 2)编写高斯消元法 MATLAB 文件如下: clear;A=5 2 1;2 8 -3;1 -3 -6; b=8;21;1;RA,RB,n,X =gaus (A,b)运行结果为:请注意:因为 RA=RB=n ,所以此方程组有唯一解 .RA =3RB =3 n =3X =12-1在 MATLAB 中利用逆矩阵法检验结果:(1) 在 command windows 中直接运行命令:A=0.101 2.304 3.555;-1.347

7、3.712 4.623;-2.835 1.072 5.643; b=1.183;2.137;3.035;X=Ab运行结果为:-0.3982X =0.01380.3351(2) 在 command windows 中直接运行命令:A=5 2 1;2 8 -3;1 -3 -6; b=8;21;1;X=Ab运行结果为:X =12-1两小题所得结果相同,检验通过2( 1)编写列组高斯消元法 MATLAB 文件如下:clear;A=0.101 2.304 3.555;-1.347 3.712 4.623;-2.835 1.072 5.643; b=1.183;2.137;3.035;RA,RB,n,X

8、=liezhu(A,b)运行结果:请注意:因为 RA=RB=n ,所以此方程组有唯一解 RA =3 RB =3 n =3X =-0.39820.01380.33512)编写列组高斯消元法的 MATLAB 文件如下: clear;A=5 2 1;2 8 -3;1 -3 -6;b=8;21;1;RA,RB,n,X =liezhu(A,b)运行结果为:请注意:因为 RA=RB=n ,所以此方程组有唯一解RA =3 RB =33n =X =12-1与题 1中逆矩阵计算所得结果相同,检验通过四 . 实验体会:通过实验我掌握了消元法解方程的一些基本算法以及用 matlab 实现矩阵的几种基本 计算。对 MATLAB

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论