matlab - 用 ODE 45 求解器求解六个微分方程而不定义 y 项？MATLAB

Question

我正在尝试使用 ode45 求解器求解微分方程。但是，它没有按预期工作。

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 不断返回以下错误消息：

没有足够的输入参数。

Untitledrt 中的错误（第 24 行）

ydot(1)=y(2)

我了解该错误是由于未定义 y 术语。ODE45 求解器代码的其他一些示例也没有定义 y 项。有什么我想念的吗？或者有人可以在不定义 y 术语的情况下提出一些建议以使此代码正常工作。先感谢您。

Answer 1

f(t,y)首先创建一个 M 文件，用于评估任何给定的系统右侧，t, y1, y2, y3, y4, y5, y6并将其命名为 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')

六个因变量的绘图