采用GAUSS列主消元法求解线性方程组（MATLAB）

最新推荐文章于 2024-07-29 00:30:00 发布

Duke_Laurence

最新推荐文章于 2024-07-29 00:30:00 发布

阅读量9.6k

点赞数 7

文章标签： matlab gauss 数值理论计算

本文链接：https://blog.csdn.net/duke_laurence/article/details/52950974

版权

该博客介绍了一种使用MATLAB通过GAUSS列主消元法来求解线性方程组的方法。程序首先接收用户输入的方程组参数，然后通过矩阵增广和主元选择进行消元操作。如果出现小主元或矩阵不满秩的情况，程序会提示无解。最后，通过回代法计算出方程组的解并输出。

摘要由CSDN通过智能技术生成

程序清单：

%%求解任意线性方程组的解
clc;
clear all;
format long e
disp('线性方程组求解，请输入参数');
n=input('维数n=');
A=input('矩阵A=');
b=input('右端项b=');
eps=input('控制精度eps=');

b=b'; %%变为列向量
A=[A b]; %%矩阵增广

for k=1:n-1
B=A(k:n,k); %%先将第k列可能作为主元的元素取出方至矩阵B
P=max(abs(B)); %%选主元P
if(P<eps) %%控制小主元
disp('无解');
break;
else
u=find((abs(B))==P); %%计算主元所在行相对与k行的位置
if(u~=1)
A([k,u],:)=A([u,k],:); %%换行
end
m=A(k+1:n,k)/A(k,k); %%求出各行行乘数并放至矩阵m
for i=1:lengt