马克,晚上有空发 。
方法有很多的,说说高斯列主元消去法解一般线性方程组的做法,以下是liezy.m文件,文件名不要修改就要用这个
function[RA,RB,n,x]=liezy(A,b);
B=[A b];n=length(b);
RA=rank(A);RB=rank(B);
zhicha=RA-RB;
if zhicha>0,
disp('因为RA~=RB,所以次方程无解')
return
end
if RA==RB
if RA==n
disp('方程有唯一一组解')
x=zeros(n,1); c=zeros(1,n+1);
for p=1:n-1
[Y,j]max(abs(B(p:n,p))),C=B(p,:);
B(p,:)=B(j+p-1,:);
然后再命令窗口输入:
系数矩阵A=[......,.......,.......,.......]
b=[0;0;0;0.........]
[RA,RB,n,x]=liezy(A,b)
然后就有结果了
A = [1,1,1; 2,-1,4; 1,4,-1];
r=rank(A)
y = null(A,'r')
结果
r =
2
y =
-1.6667
0.6667
1.0000