请问以下的把高斯白噪声添加到duffing系统的程序错在哪里
function dy=Dyt(t,y)muu=0.1;
a=1;
betaa=3;
sigma=0.003;
% N=2500;
% rm=randn(N);
% t=;
% 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=;
% dt=0.01;
%系统初始条件
y0=;
options = odeset('RelTol',1e-4,'AbsTol',,'InitialStep',1e-1);
=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 编辑 ]
页:
[1]