资源简介
数值分析中的高斯消元程序,在Matlab中编写
代码片段和文件信息
function [RARBnX]=gaus(Ab)
B=[A b];
n=length(b);
RA=rank(A);
RB=rank(B);
zhica=RB-RA;
if zhica>0
disp(‘因为RA不等于RB,所以此方程组无解.‘)
return
end
if RA==RB
if RA==n
disp(‘请注意:因为RA=RB=n,所以此方程组有唯一解.‘)
X=zeros(n1);
for p= 1:n-1
for k=p+1:n
m= B(kp)/ B(pp); B(kp:n+1)= B(kp:n+1)-m* B(pp:n+1);
end
end
b=B(1:nn+1);A=B(1:n1:n); X(n)=b(n)/A(nn);
for q=n-1:-1:1
X(q)=(b(q)-sum(A(qq+1:n)*X(q+1:n)))/A(qq);
end
else
disp(‘请注意:因为RA=RB end
end
end
属性 大小 日期 时间 名称
----------- --------- ---------- ----- ----
文件 632 2014-10-14 21:35 Gaus\gaus\gaus.asv
文件 618 2014-10-14 22:00 Gaus\gaus\gaus.m
文件 189 2014-10-14 21:43 Gaus\gaus\mainn3.m
文件 30984 2014-10-14 21:58 Gaus\Gaus.docx
目录 0 2014-11-13 20:37 Gaus\gaus
目录 0 2014-11-13 20:38 Gaus
----------- --------- ---------- ----- ----
32423 6
评论
共有 条评论