Question

该等式具有以下形式：

x'' + w.^2 x=n
w=1

和n是高斯噪音 mean = 0和standard deviation = 1。

没有高斯噪声我可以使用ODE45中的matlab来解决方程。问题是，当采用高斯噪声时，我该如何处理这个方程？考虑到了吗？

Answer 1

这实际上取决于噪音如何添加到系统中。如果您想任意向系统添加噪声，每次调用该函数时，都会将其添加到表示数据的等式中：

function dydt = solve(t,y)
dydt = [y(2); -y(1)+randn(1)];

然后致电

[t,y] = ode45(@solve, [0 10],[1 -1]);

这里的问题是，如果噪声与信号大小相比较大，则需要更多迭代，因此需要更多时间。

另一方面，如果噪声是预先确定的，您可以采样并保持，或者包含一阶保持，然后将其添加到系统