数值分析报告编程及运行结果(高斯顺序消元法)

数值分析报告编程及运行结果(高斯顺序消元法)

ID:29446540

大小:369.50 KB

页数:28页

时间:2018-12-19

数值分析报告编程及运行结果(高斯顺序消元法)_第1页
数值分析报告编程及运行结果(高斯顺序消元法)_第2页
数值分析报告编程及运行结果(高斯顺序消元法)_第3页
数值分析报告编程及运行结果(高斯顺序消元法)_第4页
数值分析报告编程及运行结果(高斯顺序消元法)_第5页
资源描述:

《数值分析报告编程及运行结果(高斯顺序消元法)》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库

1、实用标准文案高斯消元法1.程序:clearformatratA=input('输入增广矩阵A=')[m,n]=size(A);fori=1:(m-1)numb=int2str(i);disp(['第',numb,'次消元后的增广矩阵'])forj=(i+1):mA(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i);endAend%回代过程disp('回代求解')x(m)=A(m,n)/A(m,m);fori=(m-1):-1:1x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m)')/A(i,i);endx精彩文档实用标准文案2.运行结果:精彩

2、文档实用标准文案高斯选列主元消元法1.程序:clearformatratA=input('输入增广矩阵A=')[m,n]=size(A);fori=1:(m-1)numb=int2str(i);disp(['第',numb,'次选列主元后的增广矩阵'])temp=max(abs(A(i:m,i)));[a,b]=find(abs(A(i:m,i))==temp);tempo=A(a(1)+i-1,:);A(a(1)+i-1,:)=A(i,:);A(i,:)=tempodisp(['第',numb,'次消元后的增广矩阵'])forj=(i+1):mA(j,:)=A(j,:)

3、-A(i,:)*A(j,i)/A(i,i);endAend%回代过程disp('回代求解')精彩文档实用标准文案x(m)=A(m,n)/A(m,m);fori=(m-1):-1:1x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m)')/A(i,i);endx2.运行结果:精彩文档实用标准文案精彩文档实用标准文案追赶法1.程序:function[x,L,U]=zhuiganfa(a,b,c,f)a=input('输入矩阵-1对角元素a=');b=input('输入矩阵对角元素b=');c=input('输入矩阵+1对角元素c=');f=input('输入增广矩

4、阵最后一列元素f=');n=length(b);%对A进行分解u(1)=b(1);fori=2:nif(u(i-1)~=0)l(i-1)=a(i-1)/u(i-1);u(i)=b(i)-l(i-1)*c(i-1);elsebreak;endendL=eye(n)+diag(l,-1);U=diag(u)+diag(c,1);x=zeros(n,1);y=x;精彩文档实用标准文案%求解Ly=by(1)=f(1);fori=2:ny(i)=f(i)-l(i-1)*y(i-1);end%求解Ux=yif(u(n)~=0)x(n)=y(n)/u(n);endfori=n-1:-1

5、:1x(i)=(y(i)-c(i)*x(i+1))/u(i);end2.运行结果:精彩文档实用标准文案高斯-塞德尔迭代格式1.程序:functionx=Gauss_Seidel(a,b)a=input('输入系数矩阵a=')b=input('输入增广矩阵最后一列b=');e=0.5e-7;n=length(b);N=50;x=zeros(n,1);t=zeros(n,1);fork=1:Nsum=0;E=0;t(1:n)=x(1:n);fori=1:nx(i)=(b(i)-a(i,1:(i-1))*x(1:(i-1))-a(i,(i+1):n)*t((i+1):n))/a

6、(i,i);endifnorm(x-t)

7、))^2;endifsqrt(sum)

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

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

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