doublevp.m

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

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

clear all;
clc;

%k1=10;
%k2=-10;

v ='(2*A*exp(3*z1+z2)*k2^2+k1^2*exp(3*z1)+2*k1^2*exp(z1+3*z2)-4*k1*exp(3*z1+z2)*k2-4*k1*exp(z1+2*z2)*k2-4*k1*exp(2*z1+z2)*k2+2*A*exp(z1+3*z2)*k2^2+4*A*exp(2*z1+3*z2)*k1^2+2*A*exp(z1+3*z2)*k1^2+4*A*exp(2*z1+2*z2)*k1^2-4*k1*exp(z1+3*z2)*k2+A*exp(z1+4*z2)*k1^2-8*k1*exp(2*z1+2*z2)*k2+3*A*exp(3*z1+2*z2)*k1^2+A*exp(4*z1+z2)*k2^2+4*A*exp(3*z1+2*z2)*k2^2+2*k1^2*exp(3*z1+z2)+2*k2^2*exp(z1+3*z2)+4*k2^2*exp(z1+2*z2)-4*k1*exp(3*z1+2*z2)*k2*A-4*k1*exp(2*z1+3*z2)*k2*A+3*A*exp(2*z1+3*z2)*k2^2+4*k2^2*exp(2*z1+2*z2)*A+k1^2*exp(z1+z2)+k2^2*exp(3*z2)+A^2*exp(2*z1+3*z2)*k2^2+A^2*exp(3*z1+3*z2)*k1^2+A^2*exp(3*z1+3*z2)*k2^2-2*A^2*exp(3*z1+3*z2)*k1*k2-2*k1*exp(z1+z2)*k2+2*A*exp(3*z1+z2)*k1^2+A*exp(z1+2*z2)*k1^2+A*exp(2*z1+z2)*k2^2+k2^2*exp(z1+z2)+A^2*exp(3*z1+2*z2)*k1^2+3*k2^2*exp(2*z1+z2)+3*k1^2*exp(z1+2*z2)+4*k1^2*exp(2*z1+z2)+2*k2^2*exp(3*z1+z2)+4*k2^2*exp(2*z1+2*z2)+4*k1^2*exp(2*z1+2*z2))';
%v=subs(v,'k1','k');
%v=subs(v,'k2','-k');

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

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'