如何求解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的六个微分方程？

回答

相关问题