singleu.m

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

clear all;
clc;
k=1.0;
p=-k*k/2;

deltat=0.1;
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+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)