如何求解ODE 45的六個微分方程？

我想正確地在Matlab中輸入微分方程，所以我可以ode45解決它們，但我不能讓我的代碼工作。這是我進入代碼：如何求解ODE 45的六個微分方程？

function ydot=Untitledrt(z,y) 
ydot = zeros(6,1); 
%y(1)=A 
%y(2)=B 
%y(3)=C 
%y(4)=D 
%y(5)=P 
%y(6)=T 
m1 = 6; 
m2 = 9; 
m3 = 5; 
k1 = 6; 
k2 = 7; 
k3 = 4; 
k4 = 1; 
c1 = 1; 
c2 = 0.2; 
c3 = 0.1; 
c4 = 2; 
F1 = 3; 
F2 = 9; 
F3 = 12; 
ydot(1)=y(2) 
ydot(2)=((-((k1+k2)./m1)).*y(1))-(((c1+c2)./m1).*y(2))+((k2./m1).*y(3))+ ((c2./m1).*y(4))+(F1./m1) 
ydot(3)=y(4) 
ydot(4)=((k2./m2).*y(1))+((c2./m2).*y(2))-(((k2+k3)./m2).*y(3))-(((c2+c3).*m2).*y(4))+((k3/m2).*y(5))+((c3./m2).*y(6))+(F2./m2) 
ydot(5)=y(6) 
ydot(6)=((k2./m3).*y(3))+((c3./m3).*y(4))-(((k3+k4)./m3).*y(5))-(((c3+c4)./m3).*y(6))+(F3./m3)

MATLAB保持返回以下錯誤信息：

Not enough input arguments. 

Error in Untitledrt (line 24) 
ydot(1)=y(2)

，我知道這種情況正在發生，因爲我還沒有定義的任何Y項的。然而，我見過的其他所有代碼也沒有定義這些條款。有沒有我缺少的東西，或者讓代碼在沒有定義這些術語的情況下工作的方法？請幫幫我！先謝謝你。

來源

2016-04-26 J.Slaugh

你是如何調用ODE45？你的函數'y'必須是一個向量 –

另外請注意，這裏的'z'將是ODE的't'，你沒有使用（意味着你的系統是自治的）。在這種情況下，你應該可以用〜代替它，即'function ydot = Untitledrt（〜，y）'，在許多情況下，MATLAB將識別並通過在進行調用時不實例化該變量來加速代碼。這可能不是一個問題，但是是一個很好的開發習慣。 –

首先創建一個M文件，對任何給定的t, y1, y2, y3, y4, y5, y6評估系統f(t,y)的右側，並將其命名爲Untitledrt.m。

function ydot=Untitledrt(t,y) 
ydot = zeros(6,1); 
m1 = 6; 
m2 = 9; 
m3 = 5; 
k1 = 6; 
k2 = 7; 
k3 = 4; 
k4 = 1; 
c1 = 1; 
c2 = 0.2; 
c3 = 0.1; 
c4 = 2; 
F1 = 3; 
F2 = 9; 
F3 = 12; 
ydot(1)=y(2) 
ydot(2)=((-((k1+k2)./m1)).*y(1))-(((c1+c2)./m1).*y(2))+((k2./m1).*y(3))+ ((c2./m1).*y(4))+(F1./m1) 
ydot(3)=y(4) 
ydot(4)=((k2./m2).*y(1))+((c2./m2).*y(2))-(((k2+k3)./m2).*y(3))-(((c2+c3).*m2).*y(4))+((k3/m2).*y(5))+((c3./m2).*y(6))+(F2./m2) 
ydot(5)=y(6) 
ydot(6)=((k2./m3).*y(3))+((c3./m3).*y(4))-(((k3+k4)./m3).*y(5))-(((c3+c4)./m3).*y(6))+(F3./m3)

現在，通過線或鍵入MatLab中窗線以下命令將其粘貼到新的腳本文件

[t,Y]=ode45('Untitledrt',[0 10],[1;-1;0;0;0;0]); 

plot(t,Y(:,1),'+',t,Y(:,2),'x',t,Y(:,3),'o',t,Y(:,4),'o',t,Y(:,5),'o',t,Y(:,6),'o')

Plot of the six dependent variables

來源

2016-12-02 05:04:36 zhk

如何求解ODE 45的六個微分方程？

回答

相關問題