doubleu.asv

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

clear all;
clc;

k1=-1.0;
k2=1.0;
A=-(k1-k2)^2/(k1*k2);

w1=k1^2/2;
w2=k2^2/2;

p1=w1*k1/2;
p2=w2*k2/2;

deltat=0.1;
ST=0.0; 	%开始时间
ET=5.0;		%终止时间

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

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

z10=0.0;	%z1初值
z20=1.0;	%z2初值

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

[X,Y]=meshgrid(x,y);

q1=exp(k1*X+w1*Y+z10+p1*ST);
q2=exp(k2*X+w2*Y+z20+p2*ST);
deltapt1=exp(p1*deltat);
deltapt2=exp(p2*deltat);
qA=A*one;

for t=ST:deltat:ET
	q2=q1.*q1;

	r2=r1.*r1;

	u=(k1*q1.*(one+qA.*(r1+r2))+k2*r1.*(one+qA.*(q1+q2)))./((q1.*(one+r1.*qA)+r1).*(one+q1+r1));

	u=-0.5*u;
	surf(X,Y,u);
    colorbar;
	xlabel('X');
	ylabel('Y');
	zlabel('Z');
    
	view([30 50]);
	axis([-40 40 -40 40 -1 1]);
	pause(0.5);
	F(j)=getframe;
	j=j+1;
    q1=q1*delta
end
pause;
movie(F,2);