1.表达式中的置换操作
subexpr
首先这是一个很坑的函数。搞了半天没搞出来什么意思,发现实际操作情况和网上书上不一样
r=subexpr(s);从表达式中提取公因式,并附给r,但是实际结果是这样 的:
>> h=(t+2)^2+sqrt(t+2)+(t+2)*4;
>> l=subexpr(h)
l =
(t+2)^2+(t+2)^(1/2)+4*t+8 ///有任何鸡儿变化吗???
但是这样就可以了,不知道为什么
subs
subs(es,new), subs(es,old,new)
pretty
>> f(1)
ans =
1/2/b*(-c+(c^2-4*b*a*t^3-4*b*d)^(1/2))
>> subexpr(f(1))
ans =
1/2/b*(-c+(c^2-4*b*a*t^3-4*b*d)^(1/2))
>> pretty(f(1))
2 3 1/2
-c + (c - 4 b a t - 4 b d)
1/2 -------------------------------
b
2.符号微积分
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]
/**************************/
>> y
y =
x^3+8*x-6
>> limit(y,0)
ans =
-6
>> limit(y,2)
ans =
18
>> limit(power(2,x),x,-inf)
ans =
0
>> limit(power(1/2,x),x,-inf)
ans =
Inf
/************************************/
微分diff
diff(f,x,n) f-表达式 x-对谁微分 n-开几重微分
y =
x^3+8*x-6
>> diff(y,x,1)
ans =
3*x^2+8
/**************/
jacobian(f,n) 雅克比矩阵
taylor(g) 泰勒展开式
就不研究了,时间紧迫,看到再说
3.见例题2.3-1