各位大神，我想用蒙特卡洛法求高斯激励下的杜芬振子模型，matlab能求吗？运动方程式：mq''+0.5cq'＋0.5k1×q＋0.25k3×q³=w，w是高斯白噪声，强度d=0.5，m,c,k1,k3,是已知的，想求q,q'。我自己求出来的结果不对，算出来的（q,q')概率密度像个螺旋状。就是先求高斯白噪声，再用四阶龙格库塔，但我就是做不出各位能不能帮帮忙，2周了都没弄出来，心态崩了呀😢🤧😭