【matlab基础知识代码】（六）极限计算&区间极限与多变量极限

>>  syms x ;f=sin(x)/x;limit(f,x,0)
 
ans =
 
1

>> syms x b a;f=sin(b/x)*x*(1+a/x)^x;limit(f,x,inf)
 
ans =
 
b*exp(a)

>> syms x;limit((exp(x^3)-1)/(1-cos(sqrt(x-sin(x)))),x,0,'right')
 
ans =
 
12

>> x=-0.1:0.001:0.1; y=(exp(x.^3)-1)./(1-cos(sqrt(x-sin(x)))); plot(x,y,'-',[0],[12],'o')

>>  syms t; f=tan(t); L1=limit(f,t,pi/2,'left'), 
L2=limit(f,t,pi/2,'right')
 
L1 =
 
Inf
 
 
L2 =
 
-Inf

>>  syms n;f=n^(2/3)*sin(factorial(n))/(n+1);
limit(f,n,inf)
 
ans =
 
0

>> syms x n; f=n*atan(1/(n*(x^2+1)+x))*tan(pi/4+x/2/n)^n;
 limit(f,n,inf)
 
ans =
 
exp(x)/(x^2 + 1)

 

区间极限与多变量极限

>> syms x n real; f=x^n; L1=limit(f,x,inf), L2=limit(f,n,inf)
 
L1 =
 
piecewise(n == 0, 1, 0 < n, Inf, n < 0, 0)
 
 
L2 =
 
piecewise(x == 1, 1, 1 < x, Inf, x < -1, limit(x^n, n, Inf), x in Dom::Interval(-1, 0) | x in Dom::Interval(0, 1), 0)

分段函数：条件（范围），结果。

>> syms a b positive, syms x; f=a*sin(8*x^2)+b*cos(-2*x+2); L=feval(symengine,'limit',f,'x=infinity','Intervals')
 
L =
 
[- a - b, a + b]

 >> f=piecewise('abs(x)>1.1','1.1*sign(x)','abs(x)<=1.1','x');syms x;x0=-3:0.01:3;f1=subs(f,x,x0);plot(x0,f1)

f = @(x) (abs(x) > 1.1) .* (1.1 * sign(x)) + (abs(x) <= 1.1) .* x;
x = -3:0.01:3;
f1 = f(x);
plot(x, f1);

 

>> syms x a; syms y positive; f=exp(-1/(y^2+x^2))*sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2);C=limit(limit(f,x,1/sqrt(y)),y,inf)
 
C =
 
exp(a^2)

 

>>  syms x y; f=(x*y/(x^2+y^2))^(x^2); 
L1=limit(limit(f,x,inf),y,inf), L2=limit(limit(f,y,inf),x,inf), L3=limit(limit(f,x,y^2),y,inf), L4=limit(limit(f,y,x^2),x,inf)
 
L1 =
 
0
 
 
L2 =
 
0
 
 
L3 =
 
0
 
 
L4 =
 
0
 

>> syms r x y; f=x*y/(x^2+y^2); L=limit(subs(f,y,r*x),x,0)
 
L =
 
r/(r^2 + 1)