doublevp.asv

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

%根据k1和k2的值计算双孤子解v的表达式
%
%

clear all;

%k1=10;
%k2=-10;

v='-1/2*(4*exp(2*z1+3*z2+A)*k1^2+3*exp(3*z1+2*z2+A)*k1^2+2*exp(3*z1+z2+A)*k2^2+k1^2*exp(3*z1)+k2^2*exp(3*z2)+4*k2^2*exp(2*z1+2*z2)-4*k1*exp(2*z1+z2)*k2+2*k2^2*exp(z1+3*z2)-8*k1*exp(2*z1+2*z2)*k2-4*k1*exp(z1+2*z2)*k2-2*exp(3*z1+3*z2+2*A)*k1*k2+2*k1^2*exp(z1+3*z2)-4*exp(3*z1+2*z2+A)*k1*k2-4*k1*exp(z1+3*z2)*k2+4*k1^2*exp(2*z1+2*z2)+4*exp(3*z1+2*z2+A)*k2^2-4*exp(2*z1+3*z2+A)*k1*k2+3*exp(2*z1+3*z2+A)*k2^2+4*k2^2*exp(z1+2*z2)+2*exp(z1+3*z2+A)*k2^2+3*k2^2*exp(2*z1+z2)+exp(3*z1+3*z2+2*A)*k1^2+exp(4*z1+z2+A)*k2^2+3*k1^2*exp(z1+2*z2)+exp(3*z1+3*z2+2*A)*k2^2+2*exp(3*z1+z2+A)*k1^2+exp(z1+4*z2+A)*k1^2+2*k2^2*exp(3*z1+z2)+exp(3*z1+2*z2+2*A)*k1^2+4*k1^2*exp(2*z1+z2)+2*k1^2*exp(3*z1+z2)+4*exp(2*z1+2*z2+A)*k1^2+4*exp(2*z1+2*z2+A)*k2^2+exp(2*z1+z2+A)*k2^2+exp(z1+2*z2+A)*k1^2+k1^2*exp(z1+z2)+k2^2*exp(z1+z2)+2*exp(z1+3*z2+A)*k1^2-2*k1*exp(z1+z2)*k2-4*k1*exp(3*z1+z2)*k2+exp(2*z1+3*z2+2*A)*k2^2)/(exp(z1)+exp(z2)+exp(z1+z2+A))^2/(1+exp(z1)+exp(z2))^2';

v=subs(v,'k1','k');
v=subs(v,'k2','k');

v=subs(v,'2*A','');  %注意先替换'2*A'
v=subs(v,'A','');

v=subs(v,'4*z1','q4');  %注意替换的顺序
v=subs(v,'3*z1','q3');
v=subs(v,'2*z1','q2');
v=subs(v,'z1','q1');    %注意这里没有'1*'

v=subs(v,'4*z2','r4');  %注意替换的顺序
v=subs(v,'3*z2','r3');
v=subs(v,'2*z2','r2');
v=subs(v,'z2','r1');    %注意这里没有'1*'

v=simplify(v)

%后面的化简需要人工介入
%在Matlab中显示v，拷贝v的公式到记事本中，后面的操作在记事本中完成
%去掉分母，因为分母单独计算
%将指数里面的'+'换为'.*'，即点乘运算，不属于指数里面的'+'不要替换
%将‘exp'换为空，即去掉指数，因为在上一步中已经指数运算化成了点乘运算
%用s点乘v，即v=s.*v;
%拷贝公式到绘图程序中的适当位置用于计算v





%当k1=k2=k时的表达式
%v ='-1/2*k^2*(3*exp(q3+r2+qa)+3*exp(q1+r2)+exp(q3+r2+q2a)+exp(q1+r2+qa)+8*exp(q2+r2+qa)+3*exp(q2+r3+qa)+exp(q2+r3+q2a)+4*exp(q1+r3+qa)+exp(q4+r1+qa)+exp(q2+r1+qa)+4*exp(q3+r1+qa)+3*exp(q2+r1)+exp(q3)+exp(r3)+exp(q1+r4+qa))/(exp(q1)+exp(r1)+exp(q1+r1+qa))^2/(1+exp(q1)+exp(r1))^2'


%当k1=k,k2=-k时的表达式
%v
%='-1/2*k^2*(11*exp(q2+r3+qa)+8*exp(q3+r1)+exp(q4+r1+qa)+exp(q2+r1+qa)+4*exp(q1+r1)+4*exp(q3+r1+qa)+11*exp(q2+r1)+exp(q2+r3+q2a)+4*exp(q1+r3+qa)+8*exp(q1+r3)+4*exp(q3+r3+q2a)+exp(q1+r4+qa)+exp(r3)+8*exp(q2+r2+qa)+exp(q3)+exp(q3+r2+q2a)+exp(q1+r2+qa)+11*exp(q3+r2+qa)+16*exp(q2+r2)+11*exp(q1+r2))/(exp(q1)+exp(r1)+exp(q1+r1+qa))^2/(1+exp(q1)+exp(r1))^2'