關於MATLAB中的inlineeval的錯誤

Question

作爲一個組項目的一部分，我們有一個2個非線性微分方程組，我們必須使用中點法繪製S = S（t），I = I（t）圖形。

而我得到試圖插入與相應的微分方程的矩陣時以下錯誤：在內嵌表達式==>矩陣（[[

「錯誤 - （I S）/ 1000]， [（I S）/ 1000 - （3 * I）/ 10]。]），用於類型 '雙重' 的輸入參數在內嵌/的subsref

誤差（第23行） 未定義的函數 '矩陣' INLINE_OUT_ = inlineeval（INLINE_INPUTS_，INLINE_OBJ_.inputExpr，INLINE_OBJ_.expr）;「

我做的代碼如下：

syms I S 
u=[S;I]; 
F=[-0.001*S*I;0.001*S*I-0.3*I]; 
F1=inline(char(F),'I','S'); 
h=100; %Valores aleatórios 
T=100000; 
ni=(T/h); 
u0=[799;1]; 
f=zeros(1,2); 
k=zeros(1,2); 
i=1; 

while i<=ni 
f(1)=F1(u0(1)); 
f(2)=F1(u0(2)); 
dx=h*f; 
k(1)=F1((u0(1)+h*(1/2)),(u0(2)+h*(1/2))); 
k(2)=F1((u0(1)+h*(1/2)),(u0(2)+h*(1/2))); 
u1=u0+h*k; 
disp('i:'),disp(i) 
disp('u= '),disp(u1) 
u0=u1; 
i=i+1; 
end 

我是新來這個所以它很可能是錯誤的算法，但如果有人可以幫助我的錯誤我會apreciate它。謝謝！

Answer 1

具體產生錯誤的問題是，您將兩個符號函數放入矩陣中，然後調用char（輸出matrix([[-(IS)/1000], [(IS)/1000 - (3*I)/10]])而不是很好地轉換爲字符串）。

第二個問題是您正試圖將兩個函數同時傳遞到inline。 inline從一個字符串創建一個函數（並且使用匿名函數而不是inline是首選）。你不能在其中放置多個功能。

這裏你不需要sym。事實上，如果你根本不需要操縱方程式，那就避免它（比它更值得的麻煩）。一種常見的方法是創建一個單元陣列：

F{1} = @(I,S) -0.001*S*I; 
F{2} = @(I,S) 0.001*S*I-0.3*I; 

然後，您可以通過在I和S像這樣：

F{1}(500,500) 

注意，您的功能包括I和S，所以他們總是必要。當只傳遞一個像這樣的變量時，請重新考慮您期待的內容：f(1)=F1(u0(1));，因爲這也會導致錯誤。