MATLAB--符号与数值

%符号表达式的整理
%% 化简 
syms a b
y=(cot(a/2)-tan(a/2)*(1+tan(a)*tan(a/2)));
y=simplify(y);
%y=simple(y) 已经删除
%% 因式分解 
syms m n x
y=-24*m^2*x-16*n^2*x;
y=factor(y);
y1=m^3-n^3;
y1=factor(y1);
y2=factor(104);
%% 展开
syms a x
y=a*(x^2-a)^2+(x-2);
y=expand(y)
%% 合并 
syms x y
z=(x+y)^2*y+5*y*x-2*x^3;
z=collect(z)
z=collect(z,y)
%% 计算分子和分母 
syms x y
z=1/x*y+x/(x^2-2*y)
[z1,z2]=numden(z)%z1分子，z2分母
%有理数
[z1,z2]=numden(sym(4.57))
%% pretty
%mupad
syms a b c
M1=[cos(a) sin(b);
    -sin(c) cos(b);];
M2=[1/2*a+b c/b^2;
    (a+b)/c   c;];
M=M1*M2;
M=M(1)*M(2)*M(3)*M(4);
pretty(M)
%% 反函数 
syms x y
z=x+y
z1=finverse(z,x)
z2=finverse(z,y)
%% 函数复合 
syms x y z g h
z=x^2+y-2*x;
g=cos(h)
z=compose(z,g,'x','h')%指出z中的变量是x，g中的是h
z=x^2+2/x;
g=sin(h)
z=compose(z,g)%均为单变量
%% 嵌套 
syms x y
y=2*x^6-3*x^3-5*x^2+12;
y=horner(y)

%% 符号微分
syms x1 x2 x3
y=x1^4-5*x1^2+6;
dy=diff(y);%求一阶导数
d2y=diff(y,2);%求二阶导数

y1=x1^5*x2+x2*x3-x1^2*x3;
py1=diff(y1,x1,1);%先指明对哪个求偏导
py2=diff(y1,x1,2);
%symvar(y1)   该函数可查看未指明时求导的顺序[ x1, x2, x3]
%% 混合偏导
syms a b x y 
u=exp(a*x)*cos(b*y);
duxy=diff(u,x,y);
duyx=diff(u,y,x);

w=x*sin(x*y);
dwx=diff(w,x,3);
dwy=diff(dwx,y);
dw3xy=diff(w,x,x,x,y);%与上式等价
%% 差分
A=[4 5 6 3 2 1];
a=diff(A);
A=[4 5 6;7 4 2;5 6 2];
A1=diff(A);
A2=diff(A,2);