Question

我正在尝试在Matlab中实现以下代码：

v0=eps;
t0=0; tf=10;
[t,v]=ode45(@const_pow, [t0 tf], v0);

const_pow函数是：

function f=const_pow(t,v)
   f=(P/v-F0-c*v^2)/m; 
end

在运行程序时，出现以下错误：

 Error using ==> odearguments at 103

CONST_POW returns a vector of length 0, but the length of initial 
conditions vector is 1. The vector returned by CONST_POW and the   
initial conditions vector must have the same number of elements.

有关错误可能位置的任何建议吗？

Answer 1

从您的评论中，您将全局定义一些变量。要在const_pow内访问它们，您需要告诉该函数使用全局变量，例如：

function f=const_pow(t,v)
   global P F0 c m
   f=(P/v-F0-c*v^2)/m; 
end

如果不这样做，该函数将不会触及全局变量（这是为了阻止您意外更改其值）。例如：

function f=const_pow(t,v)
   global P F0 c
   m = 23.7; % this does not change the global value of m
   f=(P/v-F0-c*v^2)/m; 
end

在不使用全局变量的情况下传递大量变量的替代方法是使用所需数据定义结构：

data.P = % initial values
data.F0 = 
data.c = 
data.m =

如果有多组选项，这也会为您提供更多灵活性。或者，您可以直接传递变量。要使用ode45执行此操作，您需要使用匿名函数parameterize函数：

function f=const_pow(t,v,P,F0,c,m)
    % and so on
end

P,F0,c,m预先定义：

odefun = @(t,v) const_pow(t,v,P,F0,c,m); % anonymous function

然后：

[t,v]=ode45(odefun, [t0 tf], v0);

运行

1 个答案: