• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

matlab的符号计算

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

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

 

 

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

>> 

 

 

 

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


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
delphiTServerSocket阻塞线程单元实例发布时间:2022-07-18
下一篇:
DELPHIlistbox去除重复项发布时间:2022-07-18
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap