singlev.m

来自「这是我和同事门一起努力编的的孤立子程序。基本上所以孤立子的情况都考虑到了」· M 代码 · 共 37 行

clear all;

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

deltat=0.1;
ST=-4.0; 	%开始时间
ET=4.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);
z=k*X+z0+p*ST;
deltapt=p*deltat;
for t=ST:deltat:ET
    v=-0.125*k*k*((sech(z/2)).^2);
    
	surf(X,Y,v);
	F(j)=getframe;
	j=j+1;
    
   	z=z+deltapt;
end
pause
movie(F)