lilei2.asv

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

clear all;

lambda0=100.0;
choice=1;   %只有当lambda0<0时使用

if lambda0>=0
    k=sqrt(lambda0);
    flag=true;      %lambda0>=0
else
    k=sqrt(-lambda0);
    flag=false;     %lambda0<0
end

w=4*k^2;
p=w*k;

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

SX=-10.0;	%x开始
EX=10.0;	%x终止
x=SX:0.1:EX;	%x序列

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

figure(1);
pause(3);

[X,Y]=meshgrid(x,y);
A=2*(k*X+p*ST);
deltapt=2*p*deltat;%偏移
for t=ST:deltat:ET
    subplot(1,2,1);
    if flag     %
        u1=w*sec(A).^2;
        surf(X,Y,u1);
        subplot(1,2,2)
        u2=-w*sin(A).*(sec(A).^2);
        surf(X,Y,u2);    
    else
        
    end
    A=A+deltapt;    %下一时刻的值
    pause(0.5);
end