doubleu.m

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

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=1.0;
ST=-60.0; 	%开始时间
ET=60.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);
r1=exp(k2*X+w2*Y+z20+p2*ST);
deltapt1=exp(p1*deltat);
deltapt2=exp(p2*deltat);

one=ones(size(X));
qA=A*one;

mov = avifile('example1.avi','compression','None','fps',10)

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;
    mov = addframe(mov,getframe);
    j=j+1;  
    q1=q1*deltapt1;
    r1=r1*deltapt2;
end
%pause;
%movie(F,2);
mov = close(mov);
return;