蒙特卡洛模拟高斯白噪声下的杜芬振子

蒙特卡洛模拟高斯白噪声下的杜芬振子

function dydt = duffing(t,y,pt,W,xi,w0,epsilon)
Wt=interp1(pt,W,t);
dydt = [y(2);(Wt-2xiw0y(2)-w0^2y(1)-epsilonw0^2y(1)^3)];
end

白噪声激励生成之后，用interp1进行迭代运算，1W个sample要run 20秒左右，10w个Sample就卡住了，想请问一下有没有比较好的迭代白噪声激励代入方程进行循环的的方式呀？

朋友你可以看看这篇论文，里面有广义平方程求解
https://www.doc88.com/p-3357828127704.html
我们当初模拟分析的泊松白噪声用的它，