主元消去法

最新推荐文章于 2021-05-20 01:21:06 发布

kill bert

最新推荐文章于 2021-05-20 01:21:06 发布

阅读量1k

点赞数

分类专栏：计算方法文章标签：主元消去法

本文链接：https://blog.csdn.net/Sweeping_Robot/article/details/106506843

版权

计算方法专栏收录该内容

31 篇文章 1 订阅

订阅专栏

在这里插入图片描述

matlab代码：

clear all
A=[1 0.8324 0.7675 0.9831;0.8324 0.6930 0.6400 0.8190;0.7675 0.6400 0.5911 0.7580;0.9831 0.8190 0.7580 0.0055];
b=[8.8997 7.4144 6.8428 4.9171]';
B= [A, b];
[m, n]=size (A) ;
X=zeros (m, 1);
for j=1:m-1
    Q=B(:,j);
    for k=1:j-1 Q(k)=0; end
    imax=find(abs(Q)==max(abs (Q)));
    i=j;
    P=B(i,:);
    B(i,:)=B(imax,:) ;
    B(imax,:)=P;
    for i=j:n-1
        B(i+1,: )=B(i+1,:)-(B(i+1,j)./B(j, j).*B(j,:));
    end
end
for i=m:-1:1
	sum=0;
    for j=m:-1:i+1
        sum= sum+X(j)*B(i,j);
    end
    X(i)= (B(i, m+1)-sum)/B(i,i);
end
X