数值分析1_列主元&全主元消去法.pdf

作业二：分别编写用列主元消去法和全主元消去法求解线性方程组的

标准程序，求下列方程的解。

0.8320.4480.19311

[][]=[]

0.7840.4210.207x20

0.784−0.4210.293x30

高斯全主元消去法：

流程图

开

判断有没有解，并输出

寻找绝对值最大的项

并记录位置

换主行，换主列并记录

求解，并换回交换的列位子

结束

解题思路

首先寻找绝对值最大的项并记录位置，然后换主行，换主列并记录，最后再

求解，并换回之前交换的矩阵的列位。

程序

function[x,qa]=gaussq(a,b)%输出想x为解，qa为全主元变换后的a矩

阵

a=[0.8320.4480.193;0.7840.4210.207;0.784-0.4210.293];

b=[1,0,0]';

d=[ab];

RA=rank(a);RD=rank(d);L=length(b);n=size(a);pos=1:n(1);

ifRA~=RD

fprintf('无解')

els

ifRA~=L

fprintf('有无数多个解')

els

fprintf('有唯一解')

forq=1:n

big=max(max(abs(a(q:n,q:n))));

forr=q:n

fort=q:n

ifbig==abs(a(r,t))

zhuh=r;

zhul=t;

end

end

end

p=a(q,:);%换主

a(q,:)=a(zhuh,:);

a(zhuh,:)=p;

bb=b(q);

b(q)=b(zhuh);

b(zhuh)=bb;

p=a(:,q);%换主列

a(:,q)=a(:,zhul);

a(:,zhul)=p;

p=pos(q);%记录由于换主列而造成的解的位置的变化

pos(q)=pos(zhul);

pos(zhul)=p;

end

c=[ab];

forj=1:L-1%化为上三角阵

fori=(j+1):L

m=c(i,j)/c(j,j);

c(i,:)=c(i,:)-c(j,:)*m;

end

end

x(L,1)=c(L,L+1)/c(L,L);

fork=L-1:-1:1%求解x

x(k,1)=(c(k,L+1)-c(k,k+1:L)*x(k+1:L))/c(k,k);

end

y=[1:n(1)]';%交换被列调换时打乱的解的位置

forw=1:

forv=1:

if(pos(v)==w)

y(w)=x(v);

end

end

end

x=y;

end

end

qa=a

程序结果

高斯列主元消去法：

流程图

开始

判断有没有解，并输出

寻找当前列中绝对值

最大的项

换主行

求解

结束

程序

%function[x,qa]=gaussq(a,b)

a=[0.8320.4480.193;0.7840.4210.207;0.784-0.4210.293];

b=[1,0,0]';

d=[ab];

L=length(b);n=size(a);

forq=1:n

big=max(max(abs(a(q:n,q:n))));

forr=q:n

ifbig==abs(