ode45不传递参数(MATLAB)

Question

我在matlab中实现了一种滤波算法，在参数传递方面遇到了一些困难。

我的动力是由功能控制的：

function zdot = VortexDynamics(t, z, vorticity)
for ii=1:length(vorticity),
    zdot(ii,1) = 1i/(2*pi) * sum(  vorticity([1:ii-1,ii+1:end])  .* ( z(ii) - z([1:ii-1,ii+1:end])  ) ./ (abs(z(ii) - z([1:ii-1,ii+1:end])  ).^2 ) ) ;  
end
end

我最初用ode45计算轨迹，没有问题

 [t_true,z_true]=ode45(@VortexDynamics, tspan ,z0_true, OPTIONS, vorticity) ;

为了实现过滤器，我为每个时间步骤调用一个函数：

 [z_a , Gam, z_ens] = EnKF_nonlinear(num_tracers, Func_Observation, Observations(:,ii),vorticity,  [prev_t t], z_ens, Gam, W,R, VortexDynamics, H, N) ;

EnKF_nonlinear的定义如下：

 function [z_a, Gam_a, z_ens] = EnKF_nonlinear(num_tracers, Func_Observation, Observations,vorticity,  tspan, z0_est, Gam0_est, W,R, dynFunc, H,N )

在EnKF_nonlinear中，我需要再次将动态与ode45集成：

 for ii = 1:N

  [t_out, z_out] = ode45(dynFunc, tspan, z0_est(:,ii), vorticity);
  z_ens(:,ii) = z_out(end,:) ; 
 end

但我知道错误是：

 Error using VortexEnKF>VortexDynamics (line 85)
 Not enough input arguments.

 Error in VortexEnKF (line 63)
     [z_a , Gam, z_ens] = EnKF_nonlinear(num_tracers, Func_Observation,
     Observations(:,ii),vorticity,  [prev_t t], z_ens, Gam, W,R, VortexDynamics, H,
     N) ;

但VortexDynamics只要求t，z，涡度作为参数。在这两个函数中，我都以相同的方式调用ode45，但只有通过EnKFnonlinear调用它时才会得到错误。我不明白它怎么会没有正确的参数数调用。有什么想法吗？

Answer 1

愚蠢的我，我忘了在ode45调用中包括一个[]作为精确选项的占位符。

ode45(dynFunc, tspan, z0_est(:,ii), [], vorticity);

将涡度参数传递给dynFunc，运行良好。