solvedouble.m

来自「这是我和同事门一起努力编的的孤立子程序。基本上所以孤立子的情况都考虑到了」· M 代码 · 共 43 行

clear all;
clc;

%A='-(k1-k2)^2/(k1*k2)';

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

%p1='k1^3/4';
%p2='k2^3/4';

z1='k1*x+w1*y+p1*t+z10';
z2='k2*x+w2*y+p2*t+z20';

f='1+exp(z1)+exp(z2)';
g='exp(z1)+exp(z2)+A*exp(z1+z2)';

f=subs(f,'z1','k1*x+w1*y+p1*t+z10');
f=subs(f,'z2','k2*x+w2*y+p2*t+z20');

g=subs(g,'z1','k1*x+w1*y+p1*t+z10');
g=subs(g,'z2','k2*x+w2*y+p2*t+z20');

u=log(g/f);
u=-diff(u,'x')/2;

v=log(g*f);
v=-diff(v,'x',2)/2;


u=subs(u,'k1*x+w1*y+p1*t+z10','z1');
u=subs(u,'k2*x+w2*y+p2*t+z20','z2');
u=simplify(u);

v=subs(v,'k1*x+w1*y+p1*t+z10','z1');
v=subs(v,'k2*x+w2*y+p2*t+z20','z2');
v=simplify(v);