如何用matlab解方程组

可以使用Matlab内置的函数来解决线性方程组，其中包括LU分解法、高斯消元法、追赶法等。其中，LU分解法是一种常用的方法，可以通过以下步骤来实现：
1. 创建方程组系数矩阵A和常数向量b；
2. 判断方程组是否有解，可以通过求解方程组的秩来判断；
3. 创建增广矩阵B=[A,b]，并求解其秩s；
4. 若r=s=n(未知数)，则该齐次线性方程组有唯一解，可以利用矩阵的逆求解；
5. 若方程组有无穷多解，则可以利用LU分解法求解，具体步骤为：先将系数矩阵A进行LU分解，得到LU=PA，然后解Ly=Pb，最后再解Ux=y得到原方程组的解。

下面是一个使用LU分解法解决线性方程组的Matlab程序范例：
function x=solvebyLU(A,b)
%该函数利用LU分解法求解线性方程组Ax=b的解
flag=isexist(A,b); %调用自定义函数isexist（）判断方程组解的情况
if flag==0
disp('该方程组无解！');
x=[];
return;
else
r=rank(A);
[m,n]=size(A);
[L,U,P]=lu(A);
b=P*b;
%解Ly=b
y(1)=b(1);
if m>1
for i=2:m
y(i)=b(i)-L(i,1:i-1)*y(1:i-1)';
end
end
y=y';
%解Ux=y得原方程组得一个特解
x0(r)=y(r)/U(r,r);
if r>1
for i=r-1:-1:1
x0(i)=(y(i)-U(i,i+1:r)*x0(i+1:r)')/U(i,i);
end
end
x0=x0';
if flag==1 %若方程组有唯一解
x=x0;
return;
else %若方程组有无穷多解
format rat;
Z=null(A,'r'); %求出对应齐次方程组的基础解系
[mZ,nZ]=size(Z);
x0(r+1:n)=0;
for i=1:nZ
t=sym(char([107 48+i])); %取k=[k1,k2...,];
k(i)=t;
end
x=x0;
for i=1:nZ
x=x+k(i)*Z(:,i); %将方程组的通解表示为特解加对应齐次通解形式
end
end
end