singleu.m

这是我和同事门一起努力编的的孤立子程序。基本上所以孤立子的情况都考虑到了

clear all;
k=1.0;
w=k*k/2;
p=w*k/2;

deltat=0.1; %大于0
ST=-3.0; 	%开始时间
ET=3.0;	    %终止时间

SX=-20.0;	%x开始
EX=20.0;	%x终止
x=SX:1:EX;	%x序列

SY=-20.0;	%y开始
EY=20.0;	%y终止
y=SY:1:EY;	%y序列

z0=0.0;		%z初值

j=1;
figure(1);
pause(3);

[X,Y]=meshgrid(x,y);
z=k*X+w*Y+z0+p*ST;
b=exp(z); %初值
deltapt=exp(p*deltat);%偏移
for t=ST:deltat:ET
	c=1+b;
	c=1./c;
    u=-0.5*k*c;
	surf(X,Y,u)
	F(j)=getframe;
	j=j+1;
    b=b*deltapt;%下一时刻的值
end
pause
movie(F)