<span>matlab的符号计算</span>

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

 

 

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

>> 

 

 

 

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

全部评论

相关推荐

10-15 16:27
门头沟学院 C++
LeoMoon:建议问一下是不是你给他付钱😅😅
点赞 评论 收藏
分享
一个菜鸡罢了:哥们,感觉你的简历还是有点问题的,我提几点建议,看看能不能提供一点帮助 1. ”新余学院“别加粗,课程不清楚是否有必要写,感觉版面不如拿来写一下做过的事情,教育经历是你的弱势就尽量少写 2. “干部及社团经历”和“自我评价”删掉 3. 论文后面的“录用”和“小修”啥的都删掉,默认全录用,问了再说,反正小修毕业前肯定能发出来 4. 工作经验和研究成果没有体现你的个人贡献,着重包装一下个人贡献
点赞 评论 收藏
分享
评论
点赞
收藏
分享
牛客网
牛客企业服务