doublevk_k.asv

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

clear all;
clc;

k=1;

k1=k;
k2=-k;
w1=k1^2/2;
w2=k2^2/2;
p1=w1*k1/2;
p2=w2*k2/2;
A=-(k1-k2)^2/(k1*k2);

deltat=5;
ST=-50.0; 	%开始时间
ET=50.0;		%终止时间

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

y=x;	%y序列,和x相同

z10=-2.0;	%z1初值
z20=2.0;	%z2初值

[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);

qa=A*ones(size(X));
q2a=qa.*qa;

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

for t=ST:deltat:ET
    q2=q1.*q1;
    q3=q2.*q1;
    q4=q3.*q1;
    
    r2=r1.*r1;
    r3=r2.*r1;
    r4=r3.*r1;
    
    s=(q1+r1+q1.*r1.*qa).*(1+q1+r1);
    s=1./(s.*s);
   
    v =-1/2*k^2*s.*(11*(q2.*r3.*qa)+11*(q2.*r1)+8*(q3.*r1)+(q4.*r1.*qa)+(q2.*r1.*qa)+4*(q1.*r1)+4*(q3.*r1.*qa)+8*(q1.*r3)+4*(q3.*r3.*q2a)+(r3)+4*(q1.*r3.*qa)+(q2.*r3.*q2a)+(q1.*r4.*qa)+(q3.*r2.*q2a)+(q1.*r2.*qa)+11*(q3.*r2.*qa)+16*(q2.*r2)+(q3)+11*(q1.*r2)+8*(q2.*r2.*qa));    
    
    surf(X,Y,v);
    
	xlabel('X');
	ylabel('Y');
	zlabel('Z');
    
	view([30 50]);
	axis([-50 50 -50 50 -0 1]);
	pause(0.5);
	%F(j)=getframe;
	j=j+1;
    
    q1=q1*deltapt1;
    r1=r1*deltapt2;
end
%pause;
%movie(F,2);