微分方程matlab编程

大家好！我的大论文中需要实现以下微分方程的仿真图形，dx=(-x3+48x+100)dt+dw(t)，其中，dw(t)是通常意义上的高斯白噪声。
请matlab 高手给予指导支持，最好能编写程序，也可以私下交流QQ:541456653，非常感谢大家能够积极参与，期待您的解答。谢谢！

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。参数值你可以自己调。

非常感谢

soar1120 发表于 2014-6-24 21:33
x0=0;
T=1;
M=100;%number of periods

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

tcyrdjjlt 发表于 2014-6-25 09:44
很感谢你的回复，这个程序运行结果类似于指数曲线，与上下波动形曲线差距较大，能否烦请您调整一下能否出 ...

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

求报酬。。

soar1120 发表于 2014-6-25 15:24
y是什么，然后为什么x=(y-0.16)/1.71.还有就是-x3是什么意思，我昨天猜的是x的三次方。然后dt的选取影响也 ...

An Algorithmic Introduction to.pdf (542.21 KB)

2014-6-26 08:27:27 上传

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

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

在我的程序里，sqrt(dt)*shocks(:,j)就是dw。我用向量形式写的，你的那篇文章是用矢量写的而已。
另外你没有给k的值。
调的话，你把x0=-10就行了。我试过了，结果不好，不知道为什么。可能还是和k值有关。