请问以下的把高斯白噪声添加到duffing系统的程序错在哪里

岳叶

function dy=Dyt(t,y)
muu=0.1;
a=1;
betaa=3;
sigma=0.003;
% N=2500;
% rm=randn(N);
% t=[0:0.01:100];
% dt=0.01;
m=randn(size(t));
rm=sigma*m;
% B=wgn(1,2500,sigma);
dy=zeros(2,1);
    dy(1)=y(2);
    dy(2)=-muu*y(2)+a*y(1)-betaa*y(1)^3+rm;





close all
%数据的初始化
% mu=0.1;
% a=1;
% b=1;
% sigma=0.003;
% B=randn(1,2500);
% B=wgn(1,2500,sigma);
t0=0;
tfinal=100;
%方针开始和结束时间
tspan=[t0,tfinal];
% dt=0.01;
%系统初始条件
y0=[0,0];
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4],'InitialStep',1e-1);
[tt,yy]=ode45(@Dyt,tspan,y0,options);
figure(1);
plot(tt,yy(:,1))
xlabel('t'),title('y(1)')%实践历程图
figure(2);
plot(yy(:,1),yy(:,2))%像轨迹
ylabel('y(1)'),ylabel('y(2)')
  出现的不是混沌现象，大家帮忙看一下怎么改，谢谢了

[ 本帖最后由 eight 于 2008-5-14 12:41 编辑 ]