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

各位大虾,能不能帮帮忙,我想求高斯白噪声下的杜芬振子数值解q,q',mq"+k1q+k3q³+cq'=w,w是噪声,强度以及上面的参数都知道,我求出来的数值解太小了,只有正负0.5,应该范围在正负1.5左右,是为什么,龙格库塔应该没错,噪声数值也在5左右,求求各位解答了😭😭