Linux最常用的20条命令

知识经验 16758浏览 20评论

2.3 .1符号微积分

 求极限 

 limit(f,x,a) 求f(x)中x趋近于a的极限值

例如:

>> clear >> syms k x >> lim_t=limit((1-1/x)^(k*x),x,inf)   lim_t =   1/exp(k)

求导数

diff(f,x,n) 求f(x)的n阶导数

>> clear >> syms x >> dfdx=diff(x^3,x,1)   dfdx =   3*x^2

注意上式中x 缺省时会自动调用findsym确认,n缺省时默认n=1

还有几个函数。。

jacobian(f,v)  求多元向量函数 f(v) 的雅克比矩阵(高级)

taylor(f,n,x,a) 把f(x) 在x=a 处展开为幂级数 (泰勒级数)

 

2.3.2 序列/级数的符号求和

s=symsum(f,x,a,b) 求f(x) 从a到b的和

x缺省时findsym自动确认,a,b可同时缺省 默认为 [0,x-1]

f是矩阵时,求和将对元素逐个进行

>> clear >> syms t k >> f=[t,k^3]; >> s=simple(symsum(f))   s =   [ t^2/2 - t/2, k^3*t]

2.3.3 符号积分

int(f,x) 求f(x)dx 的不定积分

int(f,x,a,b) 从a到b的定积分

同理,x缺省findsym自动确认 a,b可以任意值或符号表达式

>> clear >> syms x >> f=sqrt((x+1)/x)/x; >> s=int(f,x),s=simple(s)   s =   - 2*(1/x + 1)^(1/2) - 2*atan((1/x + 1)^(1/2)*i)*i     s =   - 2*(1/x + 1)^(1/2) - 2*atan((1/x + 1)^(1/2)*i)*i

 。。结果怎么整成复数了。。

矩阵求积分

>> clear >> syms a b x >> f=[a*x,b*x^2;1/x,sin(x)]; >> disp('the ans is'); the ans is >> int(f)   ans =   [ (a*x^2)/2, (b*x^3)/3] [    log(x),   -cos(x)]   >> pretty(int(f))     +-                 -+    |      2        3   |    |   a x      b x    |    |   ----,    ----   |    |    2        3     |    |                   |    |  log(x), -cos(x)  |    +-                 -+

感觉pretty()一下更丑了。。

还能求多重积分。。

>> clear >> syms x y >> int(int(x*y,y,2,3),x,1,2)   ans =   15/4

发表我的评论
取消评论

表情

Hi,您需要填写昵称和邮箱!

  • 昵称 (必填)
  • 邮箱 (必填)
  • 网址

网友最新评论 (20)

  1. 我喜欢评论。
    乐心湖 4年前 (2014-08-16) 回复 编辑
  2. 这个我喜欢,不错不错。
    鸿念网 4年前 (2014-08-17) 回复 编辑
  3. 非常感谢站长分享。
    悦色 4年前 (2014-08-19) 回复 编辑
  4. 惭愧。。。从没没有玩过linux。。。
    空城旧梦 4年前 (2014-11-22) 回复 编辑
  5. 收藏起来备用
    鲜活优惠码 4年前 (2015-03-16) 回复 编辑
  6. 很详细
    东北黑木耳价格 4年前 (2015-06-27) 回复 编辑
  7. 支持
    杭州网站设计 4年前 (2015-06-27) 回复 编辑
  8. 不错的文章 顶一个
    COD试剂 4年前 (2015-07-10) 回复 编辑
  9. 玩过linux
    广州seo 4年前 (2015-07-12) 回复 编辑
  10. 不错的
    北京男科医院 4年前 (2015-07-20) 回复 编辑
  11. 不错的
    北京男科医院 4年前 (2015-07-20) 回复 编辑
  12. 感觉挺厉害的
    松露视频 3年前 (2015-08-15) 回复 编辑
  13. 长知识了
    小旺财 3年前 (2015-08-19) 回复 编辑
  14. 我是门外汉 看一下 多了解些
    卡丝面膜 3年前 (2015-08-22) 回复 编辑
  15. 我是同性恋 3年前 (2015-10-04) 回复 编辑
  16. 我是同性恋 3年前 (2015-10-04) 回复 编辑
  17. 我是同性恋 3年前 (2015-10-04) 回复 编辑
  18. 深圳罗湖的糗友,有木有黄贝岭。1月12号下午5点多,有木有糗友在黄贝岭靖轩小学门口捡到钱包啊!身份证银行卡深圳通押金条,还有现金,现金可以不要。证件给我就行,都说糗友力量大。我相信还是好人多的!!!!!!
    乐儿 3年前 (2016-01-29) 回复 编辑
  19. [电影] 新手报到,认真学习
    家具电商 3年前 (2016-03-06) 回复 编辑