[求助]Van der pol方程求解

Candyy · 发表于 2006-5-1 11:02

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

x''+a2(x2-1)x'+x=ax(t)
其中x(t)为均值为0的高斯白噪声
求其概率密度函数

Candyy · 发表于 2006-5-5 13:43

没有人做过Vanderlpol振子的求解吗？急求

zhaolei · 发表于 2006-5-6 16:21

x0=[-0.2;-0.7]; t_final=20; mu=1;[t1,y1]=ode45('vdp_eq',[0,t_final],x0,[],mu); mu=2;[t2,y2]=ode45('vdp_eq',[0,t_final],x0,[],mu); plot(t1,y1,t2,y2,':') figure; plot(y1(:,1),y1(:,2),y2(:,1),y2(:,2))
function y=vdp_eq(t,x,flag,mu) y=[x(2);-mu*(x(1)^2-1)*x(2)-x(1)]; 是不是你要的那？

账号		自动登录	找回密码
密码			我要加入