我正在尝试使用ode45来解决ODE问题
[eta_ode_i, fg_ode_i] = ode45(@UVWTI,eta,FI,options);
但是我的函数在参数上取决于预先计算的数组,如下所示
function [dfi]=UVWTI(t,fi)
dfi=zeros(7,1);
dfi(1)=A11(t)*fi(1)+A12(t)*fi(2)+A13(t)*fi(3)+A14(t)*fi(4)+A15(t)*fi(5)+A16(t)*fi(6)+A17(t)*fi(7);
...
dfi(7)=similar
end
A11,A12等是长度为K的数组。但是,由于我事先不知道ode45在内部将采用什么t,所以我无法知道应该使用A11的哪个元素。我可以插值吗?应该可以,但是我想知道是否有更简单的方法。