adding noise to the solution obtained from ode45/dde23

Question

Muhammad 2023-12-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2059749-adding-noise-to-the-solution-obtained-from-ode45-dde23

评论： Muhammad 2023-12-13

采纳的回答： Sam Chak

what is the best way to add noise to the solution obtained from the solver like ode45/dde23?

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

请先登录，再进行评论。

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

Answer 1

Sam Chak 2023-12-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2059749-adding-noise-to-the-solution-obtained-from-ode45-dde23#answer_1370114

在 MATLAB Online 中打开

Hi @Muhammad,

Does it work if you execute the following?

[t, y] = ode45(@odefun, tspan, y0);
noise  = mynoise(t);    % user-defined math function for noise
ynoise = y + noise;     % additive method

3 个评论
显示 1更早的评论隐藏 1更早的评论

Sam Chak 2023-12-13

在 MATLAB Online 中打开

Hi @Muhammad

In theory, the approach remains the same. You store the ode45 solution in sol as a structure array, and then you need to evaluate the structure at n points in the interval

to obtain the data y.

%% ODE
odefun = @(t, y) [y(2); - y(1) - 2*y(2)];
tspan  = [0 10];
y0     = [1; 0];
%% Direct Method 
[t, y] = ode45(odefun, tspan, y0);
whos t y 
  Name       Size            Bytes  Class     Attributes

  t         73x1               584  double              
  y         73x2              1168  double              
%% 'sol' Method 
sol    = ode45(odefun, tspan, y0);
t      = linspace(0, 10, 73)';
y      = deval(sol, t)';
whos t y sol
  Name       Size            Bytes  Class     Attributes

  sol        1x1              5178  struct              
  t         73x1               584  double              
  y         73x2              1168  double              

Muhammad 2023-12-13

Thank you for your help

请先登录，再进行评论。

Answer 2

Jalal Khan 2023-12-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2059749-adding-noise-to-the-solution-obtained-from-ode45-dde23#answer_1370109

The best way to add noise to the solution obtained from solvers like ode45 or dde23 depends on the specific requirements of your problem and the characteristics of the noise you want to introduce