malab程序求解线性方程组

《malab程序求解线性方程组》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、题目：编写程序，求解线性方程组基于MATLAB软件平台编写了利用高斯消元法解线性方程组的程序，过程如下：一、变量说明：解线性方程组Ax=b，其中：A－线性方程组系数矩阵；b－列向量n－系数矩阵行数；x－解向量i,j－分别控制系数矩阵的行和列；m－消元过程的乘数二、基于MATLAB软件平台编写程序，如下：functionx=gauss(A,b)%定义Gauss消元函数[n,m]=size(A);%获得系数矩阵A的行、列数nb=length(b);%获得矩阵b的长度det=1;%存储行列式值x=zeros(n,1);%

2、预定义解向量x为n维1列的全0数组fork=1:n-1fori=k+1:nifA(k,k)==0returnend%以上判断系数矩阵A的主对角元素是否为0m=A(i,k)/A(k,k);%定义消元过程乘数forj=k+1:n%以下为消元过程，将系数矩阵转化为上三角矩阵A(i,j)=A(i,j)-mA(k,j);endb(i)=b(i)-mb(k);enddet=detA(k,k);enddet=detA(n,n);fork=n:-1:1%以下为逐步回代求解解向量xforj=k+1:nb(k)=b(k)-A(k,j)

3、x(j);endx(k)=b(k)/A(k,k);end三、程序应用例如：解线性方程组：打开MATLAB软件平台后，输入：>>A=[2,-1,5,1,1;4,2,4,8,-1;1,2,0,7,3;3,1,3,3,2;2,4,-3,0,4];>>b=[9;10;12;20;5];>>x=gauss(A,b)输出结果为：x=19.8840-21.9279-12.04703.33864.2006