matlab有趣的短代码,MATLAB小代码

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

function schur_36=schur36(fname)

syms x y z sumcyc_

g_61=sumcyc_*(x^4*(x-y)*(x-z));

g_62=sumcyc_*(x^3*(y+z)*(x-y)*(x-z));

g_63=(x-y)^2*(y-z)^2*(z-x)^2;

g_64=sumcyc_*((y*z)^2*(x-y)*(x-z));

g_65=x*y*z*sumcyc_*(x*(x-y)*(x-z));

g_66=x*y*z*sumcyc_*((y+z)*(x-y)*(x-z));

g_67=(x*y*z)^2;

a=subs(fname,{x,y,z},{1,0,0});

n=subs(fname,{x,y,z},{1,1,1});

d=subs(fname,{x,y,z},{1,1,0});

b=(-subs(fname,{x,y,z},{0,-1,1})-2*1i*subs(fname,{x,y,z},{0,1,1i})+8*a+d)/8;

c=(subs(fname,{x,y,z},{0,-1,1})-2*1i*subs(fname,{x,y,z},{0,1,1i})+3*d)/8;

e=(6*a-8*b+16*c-2*d-subs(fname,{x,y,z},{-1,1,1})-subs(fname,{x,y,z},{-1,1,1i}))/2;

m=(16*a-24*b+32*c+n-3*subs(fname,{x,y,z},{-1,1,1})-2*subs(fname,{x,y,z},{-1,1,1i}))/8;

schur_36=a*g_61+b*g_62+c*g_63+d*g_64+e*g_65+m*g_66+n*g_67;

end

%三元齐六次不等式的schur分拆

%例如命令窗口输入：schur36('4*(x*y+y*z+z*x)*((y+z)^2*(z+x)^2+(y+x)^2*(z+x)^2+(y+z)^2*(y+x)^2)-9*(x+y)^2*(y+z)^2*(z+x)^2')

%ans =

%3*(x - y)^2*(x - z)^2*(y - z)^2 + 4*sumcyc_*x^3*(y + z)*(x - y)*(x - z) +

%16*sumcyc_*x^2*y*z*(x - y)*(x - z)+ 4*sumcyc_*x*y*z*(y + z)*(x - y)*(x - z)