singlev.asv

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

clear all;

k=-1.0;
k1=1.0;
z0=0.0;		%z初值

if k>0 
    errordlg('error!');
    return;
end
theta=log(-k/(12*(k1+k1')*(k1^3+(k1')^3)));
w=2*k1.^3;
p=-i*k1.^2;

deltat=0.2; %大于0
ST=0.0; 	%开始时间
ET=5.0;		%终止时间

SX=-10.0;	%x开始
EX=10.0;	%x终止
x=SX:0.1:EX;	%x序列

SY=-10.0;	%y开始
EY=10.0;	%y终止
y=SY:0.1:EY;	%y序列

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

[X,Y]=meshgrid(x,y);
z2=k1*X+p*Y+z0+w*ST;
z1=z2+conj(z2)+theta;
b1=exp(z1); %初值
b2=exp(z2);
deltapt1=exp(w*deltat+conj(w*deltat));  %偏移
deltapt2=exp(w*deltat);
for t=ST:deltat:ET
    c=1+b1;
    c=b2./c;
    v=real(c);
    subplot(1,2,1);
	surf(X,Y,v);
    subplot(1,2,2);
	surf(X,Y,real(v);
    b1=b1*deltapt1;
    b2=b2*deltapt2;
    pause(0.5);
end