single1.m

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

clear all;

k=4.0;
A1=-(0.25*k*k-1)^0.5;
B1=-4*(k*k+2)^0.5;
R1=2.0;

ST=0.0; 	%开始时间
ET=3.0;		%终止时间

SX=0;	%x开始
EX=40.0;	%x终止
x=SX:1:EX;	%x序列

[X,Y]=meshgrid(x,x);
w=zeros(size(X));

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

z1=-A1*X;
for t=ST:0.1:ET
	z=z1-A1*B1*t;
    
	b=tanh(z);
    b=((b-R1).^2)./((1-R1*b).^2);
    
    u=0.25*k*k-A1*A1*b;

    surf(X,Y,u);
    
    %ribbon(x,u);
    colorbar;
    pause(1);
	%F(j)=getframe;
	%j=j+1;
end
%pause
%movie(F)