singlev.asv

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

clear all;

k=1.0;
p=-k*k/2;

ST=-2.0; 	%开始时间
ET=1.0;		%终止时间

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

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

z0=-0;		%z初值

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

[X,Y]=meshgrid(x,y);

for t=ST:0.02:ET
	z=k*X+p*t+z0;
    v=-0.125*k*k*((sech(z/2)).^2);
	surf(X,Y,v)
	F(j)=getframe;
	j=j+1;
end
pause
movie(F)