高斯消元法matlab实现

高斯消元法matlab实现

ID:14290739

大小:69.00 KB

页数:5页

时间:2018-07-27

高斯消元法matlab实现_第1页
高斯消元法matlab实现_第2页
高斯消元法matlab实现_第3页
高斯消元法matlab实现_第4页
高斯消元法matlab实现_第5页
资源描述:

《高斯消元法matlab实现》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、《数值分析》实验报告实验五班级:姓名:学号:成绩:一、实验目的与要求1.掌握高斯消去法的基本思路和迭代步骤;2.培养编程与上机调试能力。二、实验内容1.编写用高斯消元法解线性方程组的MATLAB程序,并求解下面的线性方程组,然后用逆矩阵解方程组的方法验证.(1)(2)2.编写用列主元高斯消元法解线性方程组的MATLAB程序,并求解下面的线性方程组,然后用逆矩阵解方程组的方法验证.(1)(2)三.MATLAB计算源程序1.用高斯消元法解线性方程组的MATLAB程序输入的量:系数矩阵和常系数向量;输出的量:系数矩阵和增广矩阵的秩RA,

2、RB,方程组中未知量的个数n和有关方程组解及其解的信息.function[RA,RB,n,X]=gaus(A,b)B=[Ab];n=length(b);RA=rank(A);RB=rank(B);zhica=RB-RA;ifzhica>0,disp('请注意:因为RA~=RB,所以此方程组无解.')returnendifRA==RBifRA==ndisp('请注意:因为RA=RB=n,所以此方程组有唯一解.')X=zeros(n,1);C=zeros(1,n+1);forp=1:n-1fork=p+1:nm=B(k,p)/B(p,

3、p);B(k,p:n+1)=B(k,p:n+1)-mB(p,p:n+1);endendb=B(1:n,n+1);A=B(1:n,1:n);X(n)=b(n)/A(n,n);forq=n-1:-1:1X(q)=(b(q)-sum(A(q,q+1:n)X(q+1:n)))/A(q,q);endelsedisp('请注意:因为RA=RB

4、程组中未知量的个数n和有关方程组解及其解的信息.function[RA,RB,n,X]=liezhu(A,b)B=[Ab];n=length(b);RA=rank(A);RB=rank(B);zhica=RB-RA;ifzhica>0,disp('请注意:因为RA~=RB,所以此方程组无解.')returnendifRA==RBifRA==ndisp('请注意:因为RA=RB=n,所以此方程组有唯一解.')X=zeros(n,1);C=zeros(1,n+1);forp=1:n-1[Y,j]=max(abs(B(p:n,p)));

5、C=B(p,:);B(p,:)=B(j+p-1,:);B(j+p-1,:)=C;fork=p+1:nm=B(k,p)/B(p,p);B(k,p:n+1)=B(k,p:n+1)-mB(p,p:n+1);endendb=B(1:n,n+1);A=B(1:n,1:n);X(n)=b(n)/A(n,n);forq=n-1:-1:1X(q)=(b(q)-sum(A(q,q+1:n)X(q+1:n)))/A(q,q);endelsedisp('请注意:因为RA=RB

6、元法的MATLAB文件如下:clear;A=[0.1012.3043.555;-1.3473.7124.623;-2.8351.0725.643];b=[1.183;2.137;3.035];[RA,RB,n,X]=gaus(A,b)运行结果为:请注意:因为RA=RB=n,所以此方程组有唯一解.RA=3RB=3n=3X=-0.39820.01380.3351(2)编写高斯消元法MATLAB文件如下:clear;A=[521;28-3;1-3-6];b=[8;21;1;];[RA,RB,n,X]=gaus(A,b)运行结果为:请注意

7、:因为RA=RB=n,所以此方程组有唯一解.RA=3RB=3n=3X=12-1在MATLAB中利用逆矩阵法检验结果:(1)在commandwindows中直接运行命令:A=[0.1012.3043.555;-1.3473.7124.623;-2.8351.0725.643];b=[1.183;2.137;3.035];X=Ab运行结果为:X=-0.39820.01380.3351(2)在commandwindows中直接运行命令:A=[521;28-3;1-3-6];b=[8;21;1;];X=Ab运行结果为:X=12-1两小

8、题所得结果相同,检验通过2(1)编写列组高斯消元法MATLAB文件如下:clear;A=[0.1012.3043.555;-1.3473.7124.623;-2.8351.0725.643];b=[1.183;2.137;3.035];[RA,RB

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。