How to apply time varying input with dlode45?

Question

Emebet Gedlu 2022-11-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1848528-how-to-apply-time-varying-input-with-dlode45

评论： Walter Roberson 2022-11-11

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

Hi all,

I wan to solve the NODE model

function y = NodeModel(tspan,y,theta,inputftrs)
y = tanh(theta.fc1.Weights_state*y+theta.fc1.Weights_input,inputftrs(:,:,tspan.*2)+theta.fc1.Bias);
y = pagemtimes(theta.fc2.Weights,y) + theta.fc2.Bias;
end

Where inputftrs is the input matrix. I called the dlode45 function like this

function X = model(tspan,X0,inputftrs,neuralOdeParameters)
    odeModel=@(tspan,X0,neuralOdeParameters)NodeModel(tspan,X0,neuralOdeParameters,inputftrs);    
    X = dlode45(@odeModel,tspan,X0,neuralOdeParameters,DataFormat="CB");
end

However there comes an error "Unable to apply ODE function."

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Walter Roberson 2022-11-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1848528-how-to-apply-time-varying-input-with-dlode45#answer_1097043

在 MATLAB Online 中打开

X = dlode45(@odeModel,tspan,X0,neuralOdeParameters,DataFormat="CB");

remove the @ there. odeModel is already a function handle.

2 个评论
显示无隐藏无

Emebet Gedlu 2022-11-11

Thank you. I found also I have to interpolate the input matrix to the ode time step.

Walter Roberson 2022-11-11

@Emebet Gedlu

If you have to interpolate something to use it with any of the ode*() functions, then typically that means that you are violating the mathematical contraints for the ode*() functions. The mathematics of Runge-Kutta is only valid if the equations you give and their first two derivatives are continuous. In most cases, when people interpolate they use linear interpolation, but linear interpolation between lists of points has discontinuous first derivatives.

In order to use interpolation with the ode*() functions you need to use 'cubic' or 'spline' method so that the derivatives get matched up piecewise. This will, of course, only be valid for situations in which you can say that the data values are the results of processes that are continuous, that it is valid for a point between A and B to be influenced both by A and by B. You would not be able to use this kind of interpolation for impulses such as drug dosing, or for road height measurements that reflect stones and other sharp transitions.

请先登录，再进行评论。