高斯消元法（matlab)

高斯部分主元消元法

高斯列主元消元法

高斯部分主元消去法：

原理:将线性方程组的系数即为矩阵A(n,n),对应的值即为 B(n,1),记增广矩阵C为（A，B)；

第一步:找出系数中绝对值最大的元素，将其交换到C（1，1）,通过线性运算，使得第一列C(1,1)下面的元素都消为0；

第二步：找出除第一行第一列元素，系数中绝对值最大的元素，将其交换到C（2，2），通过线性运算，使得第二列C（2，2）下面的元素全部消为0；并以此类推，得到一个上三角矩阵；

注意：在高斯全主元消元时，可能交换了列，再得到上三角后，需要将列交换回来；

例：

%高斯主元消去法
clc;
clear;
close all;
A=[1 2 3;5 4 10;3 -0.1 1];%输入矩阵系数
B=[1;0;2];
C=[A B];
[m,n]=size(C); %计算出增广矩阵的行数和列数
%找主元->消去—>找主元->消去
%高斯主元消元法会更改列数，需要记录下来，消元完毕后再改回来；
%生成一个行向量r来记录列数的更改；
r=1:m
H=m;L=n;
for q=1:m-1%找出主元max=0;
for i=q:Hfor j=q:L-1if(abs(C(i,j))>max)max=abs(C(i,j));a=i;b=j;%利用a,b来记录最大值的位置，方便后续的交换位置；endendend%将主元换至第q行第q列
C([q,a],:)=C([a,q],:);%交换第q行和最大值所在的第a行
C(:,[q,b])=C(:,[b,q]);%交换第q列和最大值所在的第b列
%r记录列的更换
r(:,[q,b])=r(:,[b,q]);%交换第q列和最大值所在的第b列
disp(C(q,q));
disp(C);
%消元过程,消去主元下面一列
for k=q+1:mtemp=C(q,q)/C(k,q);for d=1:nC(k,d)=temp*C(k,d)-C(q,d);end
end
disp(C);
end
disp(r(2));%反带回方程式
x=zeros(1,m)%生成1行m列的0矩阵
for s=m:-1:1;x(s)=C(s,n);for t=1:m-s;x(s)=x(s)-C(s,n-t)*x(m+1-t);%disp(C(s,n-t));endx(s)=x(s)/C(s,s);
end
%将更换的列更改回来,按照r给定的顺序重新排列
for l=1:mX(r(l))=x(l);
end
disp(X);

运行结果：

可以解得：x1=1.200 x2=2.00 x3=-1.40

高斯列主元消元法：

原理：将线性方程组的系数即为矩阵A(n,n),对应的值即为 B(n,1),记增广矩阵C为（A，B)；

第一步：找出第一列绝对值最大的元素，将其换到C（1，1),通过线性运算使得第一列C（1，1)下所有元素都消为0；

第二步：找出第二列绝对值最大的元素，将其换到C（2，2),通过线性运算使得第一列C（2，2）下所有元素都消为0；并以此类推：

例：

%高斯列主元消元法
clc;
clear;
close all;
A=[1/4 1/5 1/6;1/3 1/4 1/5;1/2 1 2];%输入矩阵系数
B=[9;8;8];C=[A B];
[m,n]=size(C); %计算出增广矩阵的行数和列数
H=m;L=n;%找出第一列中最大元素->消元—>找出第二列中最大元素—>消元...for i=1:n-2%第i列max=0;for j=i:m%第j行if(abs(C(j,i))>max)max=abs(C(j,i));a=i;%记录i列中最大在第a列b=j;%记录i列中最大在第b行endenddisp(max);%找到一列中最大元素后，交换位置C([i,b],:)=C([b,i],:);%交换第i行和最大值所在的第b行使得列最大元素在C（i，i）；
disp(C);
%消元
for k=i+1:mtemp=C(i,i)/C( k,i)for d=1:nC(k,d)=temp*C(k,d)-C(i,d);end
end
disp(C);end%反带回方程
x=zeros(1,m)%生成1行m列的0矩阵
x=zeros(1,m)%生成1行m列的0矩阵
for s=m:-1:1;x(s)=C(s,n);for t=1:m-s;x(s)=x(s)-C(s,n-t)*x(m+1-t);%disp(C(s,n-t));endx(s)=x(s)/C(s,s);
end
disp(x);

运行结果：

x1=-227.0769 x2=476.9231 x3=-177.6923

高斯消元法（matlab)相关推荐

高斯消元法 matlab程序
function[A,b,x]=jvzhen(A,b); [z,n]=size(A); %A为方阵,所以也可以为m=size(A,2) for i=1:z-1for t=i+1:z%从第i+1行始消第 ...
Jacobi迭代法分量形式matlab,mtalab中jacobi迭代法
一.实验目的及题目 1.1 实验目的: (1)学会用高斯列主元消去法,LU 分解法,Jacobi 迭代法和 Gauss-Seidel 迭代法解线性方程组. (2)学会用 Matlab 编写..... ...
matlab编jacobi方法,jacobi迭代法matlab
h Matlab 线性方程组的迭代解法(Jacobi 迭代法 Gauss-Seidel 迭代法)实验报告 2008年11月09日星期日 12:49 1.熟悉 Jacobi 迭代法,并编写 Matla ...
matlab 高斯迭代法求解,高斯迭代法matlab算例
Matlab 线性方程组的迭代解法 Gauss-Seidel 迭代法 Matlab 线性方程组的迭代解法 Gauss-Seidel 迭代法实验报告 1.熟悉 Gauss-Seidel 迭代法,并编写 ...
matlab用高斯消元法解线性方程组,用matlab利用高斯消元法求解线性方程组
用matlab利用高斯消元法求解线性方程组Tag内容描述: 1.用C语言编写软件完成以下任务:请用高斯列主元消元法解下列线性方程组:方法说明(以4阶为例):第1步消元在增广矩阵(A,b)第一列中找到绝 ...
高斯消元法原理与Matlab实现
直接法解线性方程组-高斯消元法 1.高斯消元法思想设有线性方程组如下所示: {a11x1+a12x2+⋯+a1nxn=b1,a21x1+a22x2+⋯+a2nxn=b2,⋮an1x1+an2x2+⋯ ...
MATLAB数值分析学习笔记：线性代数方程组的求解和高斯消元法
工程和科学计算的许多基本方程都是建立在守恒定律的基础之上的,比如质量守恒等,在数学上,可以建立起形如 [A]{x}={b} 的平衡方程.其中{x}表示各个分量在平衡时的取值,它们表示系统的状态或响应: ...
高斯消元法的MATLAB代码
高斯消元法的原理详见于线性代数课本,其MATLAB代码如下: function[x]=gauss_elim(A,b) % A为系数矩阵且必须为方阵 % b为常数项向量 n=size(A,1); x=z ...
利用Matlab进行高斯消元法
利用Matlab进行高斯消元法思路代码: 参考资料思路采取两种思路(1)不进行行行交换的高斯消元法(2)有行行交换(列主消元法) 算法如下: (1)没有行行交换的消元的算法图 (2)列主消元法 ...
MATLAB三角消元法,matlab中用的高斯消元法怎么使用！
matlab中用的高斯消元法怎么使用! mip版关注:201 答案:3 悬赏:30 解决时间 2021-01-25 19:12 已解决 2021-01-25 12:40 老师给我们一个程序让我 ...

高斯消元法（matlab)

高斯消元法（matlab)相关推荐

最新文章

热门文章