x0=0;
T=1;
M=100;%number of periods
N=1000;%number of simulations
dt=T/M;%length of time step
shocks=randn(N,M);
x=zeros(N,M+1);
x(:,1)=x0;
for j=1:M
    x(:,j+1)=x(:,j)+(-x(:,j).^3+48*x(:,j)+100)*dt+sqrt(dt)*shocks(:,j);   
end
plot(mean(x))

以上是matlab codes。参数值你可以自己调。

非常感谢

很感谢你的回复，这个程序运行结果类似于指数曲线，与上下波动形曲线差距较大，能否烦请您调整一下能否出来这样的图形，万分感谢！此外，还需要补充一些内容，dx=(-x3+48x+100)dt+kdw(t),k为常数，x=(y-0.16)/1.71,x(0)=-10,或者再调整，看看能否实现以下图形。

y是什么，然后为什么x=(y-0.16)/1.71.还有就是-x3是什么意思，我昨天猜的是x的三次方。然后dt的选取影响也比较大。你如果会matlab，可以自己调的。

求报酬。。

y先不用管它，是线性关系，就是-x的三次方，关于dw你可以参考一下这个文章，报酬肯定会有。