singleu.m

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

clear all;
clc;

deltat=1.0; %时间分辨率
ST=-10.0; 	%开始时间
ET=10.0;	%终止时间

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

y=x;	%y序列,和x相同
[X,Y]=meshgrid(x,y);

lambda=1;%参数
r=-11;
beta=sqrt(1+lambda*lambda);

xi=beta*(X+lambda*ST);
deltapt=beta*lambda*deltat;

for t=ST:deltat:ET
    sigma=beta*(tanh(xi)-r)./(1-r*tanh(xi))+lambda;
    u=(1+2*lambda*sigma-sigma.*sigma)./(6*sigma.*(sigma-2*lambda));
    
    surf(X,Y,u);
    title(t);
    pause(0.5);
    
    xi=xi+deltapt;
end