matlab的符号计算

#############################

01 展开:展开因式:expand(f)

--- 幂函数:则根据次数从高到低

--- 三角函数:展开角部分

--- 指数:展开指数部分
>> z=(x+y+3)*(2*x-4*y+7)+sin(x+y)+exp(x+y)+log(x*y)
 
z =
 
exp(x + y) + log(x*y) + sin(x + y) + (2*x - 4*y + 7)*(x + y + 3)
 
>> expand(z)
 
ans =
 
13*x - 5*y + log(x*y) + exp(x)*exp(y) + cos(x)*sin(y) + cos(y)*sin(x) - 2*x*y + 2*x^2 - 4*y^2 + 21
 
>> 
>> z=(x+y+3)*(2*x-4*y+7)
 
z =
 
(2*x - 4*y + 7)*(x + y + 3)
 
>> expand(z)
 
ans =
 
2*x^2 - 2*x*y + 13*x - 4*y^2 - 5*y + 21
 
>> 

02 合并:因式分解:factor(f)

>> z1=3*x^3+2*x^2+x+2
 
z1 =
 
3*x^3 + 2*x^2 + x + 2
 
>> factor(z1)
 
ans =
 
[ x + 1, 3*x^2 - x + 2]
 
>> 
>> z1=x^2+2*x*y+y^2
 
z1 =
 
x^2 + 2*x*y + y^2
 
>> factor(z1)
 
ans =
 
[ x + y, x + y]
 
>> 

03 合并同类项:collect(f)

>> z2=x^2+x*y*7-x^2+y+x-8
 
z2 =
 
x + y + 7*x*y - 8
 
>> collect(z2)
 
ans =
 
(7*y + 1)*x + y - 8
 
>> collect(z2,y)
 
ans =
 
(7*x + 1)*y + x - 8
 
>> collect(z2,x)
 
ans =
 
(7*y + 1)*x + y - 8
 
>> 

04 化简:simplify(f)

>> z3=x^2+y^2-2*x*y+sin(x)^2+cos(x)^2
 
z3 =
 
cos(x)^2 + sin(x)^2 - 2*x*y + x^2 + y^2
 
>> simplify(z3)
 
ans =
 
x^2 - 2*x*y + y^2 + 1
 
>> 

05 解方程:solve(f,x)

>> z=x^3+x-6
 
z =
 
x^3 + x - 6
 
>> solve(z,x)
 
ans =
 
 root(z^3 + z - 6, z, 1)
 root(z^3 + z - 6, z, 2)
 root(z^3 + z - 6, z, 3)
 
>> z=x^2+x-2
 
z =
 
x^2 + x - 2
 
>> solve(z,x)
 
ans =
 
 -2
  1
 
>> 

06 级数求和:symsum(f,n,a,b)

---  级数求和:symsum(f,n,a,b):f为一个级数的通项,是一个符号表达式,求自变量n为从a到b的通项和;
其中inf可表示无穷大
>> f=n
 
f =
 
n
 
>> symsum(f,n,1,n)
 
ans =
 
(n*(n + 1))/2
 
>> symsum(f,n,1,10)
 
ans =
 
55
 
>> f=1/n^2
 
f =
 
1/n^2
 
>> symsum(f,n,1,inf)
 
ans =
 
pi^2/6
 
>> 

07 求极限:limit(f,x,a)

--- 某点极限:  limit(f,x,a)
--- 某点左极限:limit(f,x,a,\'left\')
--- 某点右极限:limit(f,x,a,\'right\')

--- 无穷极限:  limit(f,x,inf)
--- 正无穷极限:limit(f,x,inf,\'right\')
--- 负无穷极限:limit(f,x,inf,\'left\')
>> f=sin(x)/x
 
f =
 
sin(x)/x
 
>> limit(f,x,0)
 
ans =
 
1
 
>> limit(f,x,0,\'right\')
 
ans =
 
1
 
>> limit(f,x,0,\'left\')
 
ans =
 
1
 
>> limit(f,x,inf)
 
ans =
 
0
 
>> limit(f,x,inf,\'right\')
 
ans =
 
0
 
>> limit(f,x,inf,\'left\')
 
ans =
 
0
 
>> 
>> help limit
--- sym/limit 的帮助 ---

 limit    Limit of an expression.
    limit(F,x,a) takes the limit of the symbolic expression F as x -> a.
    limit(F,a) uses symvar(F) as the independent variable.
    limit(F) uses a = 0 as the limit point.
    limit(F,x,a,\'right\') or limit(F,x,a,\'left\') specify the direction
    of a one-sided limit.
 
    Examples:
      syms x a t h;
 
      limit(sin(x)/x)                 returns   1
      limit((x-2)/(x^2-4),2)          returns   1/4
      limit((1+2*t/x)^(3*x),x,inf)    returns   exp(6*t)
      limit(1/x,x,0,\'right\')          returns   inf
      limit(1/x,x,0,\'left\')           returns   -inf
      limit((sin(x+h)-sin(x))/h,h,0)  returns   cos(x)
      v = [(1 + a/x)^x, exp(-x)];
      limit(v,x,inf,\'left\')           returns   [exp(a),  0]

08:求导数:diff(f,x,n)

--- diff(f,x,n) 表示函数f对自变量x求n阶导数
>> f=x^2+exp(x)+log(x)+sin(x)+cos(x)
 
f =
 
cos(x) + exp(x) + log(x) + sin(x) + x^2
 
>> diff(f,x,2)
 
ans =
 
exp(x) - cos(x) - sin(x) - 1/x^2 + 2
 
>> 

09 泰勒展开:taylor(f,x,a,\'Order\',n)

--- taylor(f,x,a,\'Order\',n) 表示函数f在自变量x=a处的泰勒展开式,n为展开的阶数
>> f=exp(x)
 
f =
 
exp(x)
 
>> taylor(f,x,0,\'Order\',3)
 
ans =
 
x^2/2 + x + 1
 
>> taylor(f,x,1,\'Order\',3)
 
ans =
 
exp(1) + exp(1)*(x - 1) + (exp(1)*(x - 1)^2)/2
 
>> 

10 求积分:int(f,x,a,b)

--- int(f,x)     表示函数f对自变量x的不定积分
--- int(f,x,a,b) 表示函数f对自变量x从a到b的定积分,a和b可以为正负inf
>> f=sin(x)
 
f =
 
sin(x)
 
>> int(f,x,0,pi)
 
ans =
 
2
 
>> int(f,x)
 
ans =
 
-cos(x)
 
>> 

美化符号表达式:pretty(f)

>> f=sin(pi/3)

f =

    0.8660

>> pretty(f)
未定义与 \'double\' 类型的输入参数相对应的函数 \'pretty\'。
 
>> f=sin(sym(pi/3))
 
f =
 
3^(1/2)/2
 
>> pretty(f)
sqrt(3)
-------
   2

>> 

##########################