æ¹æ³æå¾å¤çï¼è¯´è¯´é«æ¯å主å
æ¶å»æ³è§£ä¸è¬çº¿æ§æ¹ç¨ç»çåæ³ï¼ä»¥ä¸æ¯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)
ç¶åå°±æç»æäº
温馨提示:答案为网友推荐,仅供参考