[問題]關於程式控制的問題
不好意思最近在寫程式,想利用程式來解線性聯立方程式
可是不知道哪裡出了問題........怎麼跑都只會出現cannot handle這行...
以下是我所寫的....請教各位大大了...
function result=matrix(A,b)
[NR,NC]=size(A)
[nr,nc]=size(b)
if (nc ~= 1|NR~=nr)
disp('b的行數須為1且A的列數須等於b的列數,否則無法運算')
result='cannot handle';
else
if(NR == NC&NR == nr&nc == 1)
if (det(A)==0);
C=pinv(A);
D=C*b;
disp('第一個數為X的最接近範數解、第二個數為Y的最接近範數解...以此類推')
result=D
else
disp('第一個為X的解,第二個為Y的解......以此類推')
E=inv(A);
F=E*b;
result=F
end
elseif(NR < NC&NR == nr&nc == 1)
disp('此為其中可能的解')
result=A\b
elseif(NR > NC&NR == nr&nc == 1)
Am=[A,b];
disp('最後一行皆為解')
result=rref(Am)
else
disp('cannot handle')
end
end
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.115.203.129
推
03/12 15:43, , 1F
03/12 15:43, 1F
→
03/12 15:53, , 2F
03/12 15:53, 2F
推
03/12 17:00, , 3F
03/12 17:00, 3F
→
03/12 17:12, , 4F
03/12 17:12, 4F