doublev_21.m

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

clear all;

k1=-2;
k2=1;

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

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

ST=-100.0; 	%开始时间
ET=100.0;		%终止时间

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

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

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

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


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

z1_1=k1*X+w1*Y+z10;
z2_1=k2*X+w2*Y+z20;

one=ones(size(z1_1));

qa=A*one;
q2a=qa.*qa;

for t=ST:5:ET
	z1=z1_1+p1*t;
	z2=z2_1+p2*t;
    
    q1=exp(z1);
    q2=q1.*q1;
    q3=q2.*q1;
    q4=q3.*q1;
    
    r1=exp(z2);
    r2=r1.*r1;
    r3=r2.*r1;
    r4=r3.*r1;
    
    s=(q1+r1+q1.*r1.*qa).*(one+q1+r1);
    s=1./(s.*s);
   
    v =-1/2*s.*(27*q2.*r3.*qa+24*q3.*r2.*qa+10*q3.*r1.*qa+4*q3+r3+36*q2.*r2+27*q2.*r1+18*q1.*r3+24*q1.*r2+9*q3.*r3.*q2a+10*q1.*r3.*qa+q4.*r1.*qa+4*q1.*r4.*qa+4*q3.*r2.*q2a+20*q2.*r2.*qa+q2.*r1.*qa+4*q1.*r2.*qa+9*q1.*r1+q2.*r3.*q2a+18*q3.*r1);    
    surf(X,Y,v);
	xlabel('X');
	ylabel('Y');
	zlabel('Z');
    
	view([20 30]);
	%axis([-10 10 -10 40 -1 1]);
	pause(0.5);
	%F(j)=getframe;
	j=j+1;
end
%pause;
%movie(F,2);