在ode45中传递参数

Question

我想在我的ode45函数中传递一个简单的参数。我的职能如下：

function dxdt = state( t,x,vgth,vgval)
p=1;
k=10^0.7;
window1=1-((2*x)-1).^(2*p);
dxdt=k*(vgval-vgth+1.2)*window1;  
end

主要剧本是：

clear all
step=0.01;
t = 0:step:10;
f=2*0.157;
vg = 5*sin(2*f*t);
x0=0.01;
vgth=1.9;
[t,x] = ode45(@(t,x) state1 (t,x,vgth,vg(t)), t, x0);
plot(t,x)

在这里，当我在参数传递中排除参数vg(t)时，整个过程都很好。但在争论中，它不适用于vg(t)。

Answer 1

将评论汇集到一个答复中：

function dxdt = state( t,x,vgth,vgval)
    p=1;
    k=10^0.7;
    window1=1-((2*x)-1).^(2*p);
    dxdt=k*(vgval-vgth+1.2)*window1;  
end

function vgval = vg(t)
    f=2*0.157;
    vgval = 5*sin(2*f*t);  
end 

step=0.01;
t = 0:step:10;
x0=0.01;
vgth=1.9;
[t,x] = ode45(@(t,x) state(t,x,vgth,vg(t)), t, x0);

希望这能如愿以偿。