计算方法-全主元消去法-matlab程序.pdf
版权申诉
51 浏览量
2022-11-16
20:32:45
上传
评论
收藏 17KB PDF 举报
%求四阶线性方程组的 MATLAB 程序
clear
Ab=[0.001 2 1 5 1;
3 -4 0.1 -2 2;
2 -1 2 0.01 3;
1.1 6 2.3 9 4];%增广矩阵
num=[1 2 3 4];%未知量 x 的对应序号
for i=1:3
A=abs(Ab(i:4,i:4));%系数矩阵取绝对值
[r,c]=find(A==max(A(:)));
r=r+i-1;%最大值对应行号
c=c+i-1;%最大值对应列号
q=Ab(r,:),Ab(r,:)=Ab(i,:),Ab(i,:)=q;%行变换
w=Ab(:,c),Ab(:,c)=Ab(:,i),Ab(:,i)=w;%列变换
n=num(i),num(i)=num(c),num(c)=n;%列变换引起未知量 x 次序变化
for j=i:3
Ab(j+1,:)=-Ab(j+1,i)*Ab(i,:)/Ab(i,i)+Ab(j+1,:);%消去过程
end
end %最后得到系数矩阵为上三角矩阵
%回代算法求解上三角形方程组
x(4)=Ab(4,5)/Ab(4,4);
x(3)=(Ab(3,5)-Ab(3,4)*x(4))/Ab(3,3);
x(2)=(Ab(2,5)-Ab(2,3)*x(3)-Ab(2,4)*x(4))/Ab(2,2);
x(1)=(Ab(1,5)-Ab(1,2)*x(2)-Ab(1,3)*x(3)-Ab(1,4)*x(4))/Ab(1,1);
for s=1:4
fprintf('未知量 x%g =%g',num(s),x(s))
end
%验证如下
%A=[0.001 2 1 5 1; 3 -4 0.1 -2 2;2 -1 2 0.01 3; 1.1 6 2.3 9 4];
%b=[1 2 3 4]';
%x=A;
%x1= 1.0308
%x2= 0.3144
%x3= 0.6267
%x4= -0.0513
资源评论
不吃鸳鸯锅
- 粉丝: 8314
- 资源: 2万+
最新资源
- 基于ARM的指令调度和延迟分支
- 基于Vue和TypeScript的极简聊天应用设计源码 - HasChat
- 基于Vue2全家桶和Zcool数据的图片收集网站设计源码 - cool-picture
- 基于C和C++的二维绘制工具设计源码 - DrawPro
- Object.defineProperty 的 IE 补丁object-defineproperty-ie-master.zip
- 整卷预览.mhtml
- MySQL是一种广泛使用的开源关系型数据库管理系统，它提供了丰富的SQL语句用于数据库的创建、查询、更新和管理 以下是一些常见的
- MySQL是一种广泛使用的开源关系型数据库管理系统，它提供了丰富的SQL语句用于数据库的创建、查询、更新和管理 以下是一些常见
- MySQL是一种广泛使用的开源关系型数据库管理系统，它提供了丰富的SQL语句用于数据库的创建、查询、更新和管理 以下是一些常见的
- 基于Javascript的结婚请帖设计源码 - Invitation
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
