function [x,XA]= GaussXQAllMain (A,b)
N = size(A);
n = N(1);
index_l = 0;
index_r = 0;
order = 1:n; %记录未知数顺序的向量
for i=1:(n-1)
me = max(max(abs(A(i:n,i:n)))); %选取全主元
for k=i:n
for r=i:n
if(abs(A(k,r))==me)
index_l = k;
index_r = r; %保存主元所在的行和列
k=n;
break;
end
end
end
temp = A(i,1:n);
A(i,1:n) = A(index_l,1:n);
A(index_l,1:n) = temp;
bb = b(index_l);
b(index_l)=b(i);
b(i) = bb; %交换主行
temp = A(1:n,i);
A(1:n,i) = A(1:n,index_r);
A(1:n,index_r) = temp; %交换主列
pos = order(i);
order(i) = order(index_r);
order(index_r) = pos; %主列的交换会造成未知数顺序的变化
for j=(i+1):n
if(A(i,i)==0)
disp('对角元素为0!');
return;
end
l = A(j,i);
m = A(i,i);
A(j,1:n)=A(j,1:n)-l*A(i,1:n)/m;
b(j)=b(j)-l*b(i)/m;
end
end
x=SolveUpTriangle(A,b);
y=zeros(n,1);
for i=1:n
for j=1:n
if(order(j)==i)
y(i)=x(j);
end
end
end %恢复未知数原来的顺序
x=y;
XA = A;
matlab解线性方程组的直接法.rar
版权申诉
90 浏览量
2023-06-06
16:59:33
上传
评论 1
收藏 6KB RAR 举报
Matlab科研辅导帮
- 粉丝: 1w+
- 资源: 7491
最新资源
- IMG_0694.GIF
- 基于图像的三维模型重建C++源代码+文档说明(高分课程设计)
- 基于聚焦法的工件立体测量方案,根据数据进行三维重建 使用HALCON处理图像,MATLAB拟合数据+源代码+数据集+效果图
- 锄战三国村 修改:货币使用不减 v1.10(2) 原创 (中文).apk
- 基于python实现的单目双目视觉三维重建+源代码+图像图片(高分课程设计)
- 基于C+++OPENCV的全景图像拼接源码(课程设计)
- 基于Python+OpenCV对多张图片进行全景图像拼接,消除鬼影,消除裂缝+源代码+文档说明+界面截图(高分课程设计)
- 基于C++实现的全景图像拼接源码(课程设计)
- 基于SIFT特征点提取和RASIC算法实现全景图像拼接python源码+文档说明+界面截图+详细注释(95分以上课程大作业)
- 基于matlab实现眼部判别的疲劳检测系统+源代码+全部数据+文档说明+详细注释+使用说明+截图(高分课程设计)
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈