在MATLAB中求解一個非線性方程組的錯誤

Question

這段代碼有什麼問題？

function F = Final_Project_2(x) 
F = [(1/x(1)) + (1/x(2))- (2/(7*x(3))); 
     (x(3)+2*x(4))*(15*x(2))/((x(1)+x(2))*x(3)-0.7/x(3))-14; 
     (x(3)*((0.576*x(2)/(x(1)+x(2))*x(3)) - 0.27/x(3)))/(x(3)+(0.576*x(2)/(x(1)+x(2))*x(3)) - 0.27/x(3)) -20; 
     ((0.576*x(2)*x(3)/(x(1)+x(2))*x(3)) - 0.27*x(4)/x(3))-100]; 

這裏是工作區代碼：

x0 = [1; 1; 1; 1];   % Make a starting guess at the solution 
options=optimset('Display','iter'); % Option to display output 
[x,fval] = fsolve(@Final_Project_2,x0,options) % Call optimizer 
??? Error using ==> vertcat 
CAT arguments dimensions are not consistent. 

Error in ==> Final_Project_2 at 5 
F = [(1/x(1)) + (1/x(2))- (2/(7*x(3))); 

Error in ==> fsolve at 254 
      fuser = feval(funfcn{3},x,varargin{:}); 

Caused by: 
    Failure in initial user-supplied objective function evaluation. FSOLVE cannot 
    continue. 

Answer 1

Matlab的未連接以下的數組：

F = [(1/x(1)) + (1/x(2))- (2/(7*x(3))); 
     (x(3)+2*x(4))*(15*x(2))/((x(1)+x(2))*x(3)-0.7/x(3))-14; 
     (x(3)*((0.576*x(2)/(x(1)+x(2))*x(3)) - 0.27/x(3)))/(x(3)+(0.576*x(2)/(x(1)+x(2))*x(3)) - 0.27/x(3)) -20; 
     ((0.576*x(2)*x(3)/(x(1)+x(2))*x(3)) - 0.27*x(4)/x(3))-100]; 

第三行，其中最後-20被解釋爲出現該問題的矢量中的單個值，而不是您計算的術語的一部分。爲了證明這一點，你可以括號這個陣列的每一行，並看看結果：

x = [1 2 3 4]; % just for demonstration purpuse 

[(1/x(1)) + (1/x(2))- (2/(7*x(3)));] 
[  (x(3)+2*x(4))*(15*x(2))/((x(1)+x(2))*x(3)-0.7/x(3))-14;] 
[  (x(3)*((0.576*x(2)/(x(1)+x(2))*x(3)) - 0.27/x(3)))/(x(3)+(0.576*x(2)/(x(1)+x(2))*x(3)) - 0.27/x(3)) - 20;] 
[  ((0.576*x(2)*x(3)/(x(1)+x(2))*x(3)) - 0.27*x(4)/x(3))-100] 

這導致

ans = 
    1.4048 

ans = 
    23.6426 

ans = 
    0.7843 -20.0000 

ans = 
    -96.9040 

請注意，第三行中的1x2的載體。

爲了解決您的問題，或者更重視在計算中放置空格（即省略所有空格，或者在您的示例中寫入- 20而不是-20），或將所有屬於一個詞的項放在括號內。