Question

Matlab解多元一次方程组 调用m文件

如下是fce.m文件的内容

function [C,x]=fce(A,b)
[~,n]=size(A);
R=rank(A);
B=[A b];
Rr=rank(B);
if R==Rr&&R==n %n为未知数的个数，判断是否有唯一解
x=A\b;
elseif R==Rr&&R<n %判断是否有无穷解
x=A\b; %求特解?
C=null(A,'r');
%求AX=0的基础解系,所得C为n-R列矩阵，这n-R列即为对应的基础解系?
%这种情形方程组通解xx=k(p)*C(:,P)(p=1…n-R)
else x='No?solution!'; %判断是否无解
end
end

在command windows中进行调用，会出问题，但是不知道有什么问题。
例如：
>> a=[2 3;8 6]

a =

2 3
8 6

>> b=[8;20]

b =

8
20

>> [C,x]=fce(a,b)
Error in ==> fce at 2
[~,n]=size(A);

??? Output argument "C" (and maybe others) not assigned during call to
"C:\Users\203\Documents\MATLAB\fce.m>fce".

求大神指点啊

Answer 1

disp('输入第一个方程的系数')
a1=input('a1=');
b1=input('b1=');
c1=input('c1=');
d1=input('d1=');
disp('输入第二个方程的系数')
a2=input('a2=');
b2=input('b2=');
c2=input('c2=');
d2=input('d2=');
disp('输入第三个方程的系数')
a3=input('a3=');
b3=input('b3=');
c3=input('c3=');
d3=input('d3=');

A=[a1 b1 c1;a2 b2 c2;a3 b3 c3];
b=[d1 d2 d3]';
x=A\b