|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
. G0 M9 x9 L- V- u& u+ s本文介绍利用MATLAB求解函数或序列的极限问题,顺便介绍limit函数的用法。内容主要包括单变量函数的极限和多变量函数的极限。# d7 [0 f0 Z" ~* x1 [, L
目录9 i( T: S2 Q- M+ P. _2 p
单变量函数的极限" O+ m; O1 H( ]( V4 S5 V! A" N: h
极限的定义% j+ C: m6 D0 r: b+ R
普通极限/ P* f7 B: r1 o5 U
左极限8 N2 Y6 _: v7 W' w0 k Q* i. U( p
右极限
2 ?* _/ Y! g _+ gmatlab实现方法
& ~8 B5 F5 J( y应用举例
' i, N- B5 g& h" @: j, A多变量函数的极限 Y- S3 n2 X& K' X/ A5 U0 D
matlab实现方法
{" l, D* t! V7 q l* t. ~+ Y应用举例
+ }/ G& G2 I6 ?单变量函数的极限; B! O4 |( g5 f1 m8 y6 k* n0 E' M
极限的定义
; _: j3 [6 W# A) D S: o; e6 f
& ?" e3 G. r: h3 Z$ P& Z6 P, v1 i8 ~
( [4 z* z5 G% G
# J- F8 J3 x2 c: D- N8 ^ r. v0 s$ f% k( z
matlab实现方法
5 J% q3 V' o' o9 i6 q3 X- L=limit(fun, x, x0) % //普通极限
- L=limit(fun, x, x0, 'left') % //左极限
- L=limit(fun, x, x0, 'right') % //右极限" s1 F) Y! m f' ] u- n
: d1 C4 g5 v+ a8 b% T! |
K/ q! ]' b5 r. M' u应用举例
% o! P1 H- b- m+ Q求解极限:/ T' D4 Y/ M: }8 i0 W9 N! X3 k
; [% R4 O7 \& H& S) ^; m, P0 T0 o
; m+ b- _4 a2 u" P3 H: w- `; S * o$ g- {; ?+ n7 U6 `; Q) W
3 T' }# W5 y" D0 p- N- syms x; f=sin(x)/x; L=limit(f, x, 0)" |& T \* D2 v
/ X& t# {5 f0 @- r w9 [7 q
+ d p$ r# R, q! n. t, ~( y求解极限:
& X0 P1 y* s q7 v$ b1 r& q5 P% E. p$ h4 X7 e' x
) {. o4 t h% A
H7 K# b- Z! ^6 u5 e- Q
- syms x a b
- f = x*(1+a/x)^x*sin(b/x)
- L = limit(f, x, inf): F0 }# g8 G' v/ V
" i; J: E S0 P. L
/ N4 j. w: g0 v% J5 q! c求解单边极限:/ M" N9 O! o2 f6 x4 u1 a6 D( B
5 V- P( L- |3 k2 N' N+ i- y$ o4 ]
5 v& S9 W) j! p0 Y
0 t% N1 L8 ?( F% C- syms x; L = limit((exp(x^3)-1)/(1-cos(sqrt(x-sin(x)))),x,0,'right')3 [& r2 w7 h$ R* E E+ L# e
) I- ?5 b6 U# O5 f; E
( {. ~2 c, C4 c) ^, ]0 E用下面的语句还可以绘制出 ( − 0.1 , 0.1 ) (-0.1,0.1) (−0.1,0.1)区间的函数曲线。% M8 H- a! j# }% U4 n
, ]1 _# V# q2 f' d
- x0=-0.1:0.001:0.1;
- y0=((exp(x0.^3)-1)./(1-cos(sqrt(x0-sin(x0)))));
- plot(x0, y0, '-', [0], [L], 'o') Z* A5 b* V% E
, g9 b9 _; u' b. v: t
6 F/ K* i& r) o2 V7 j* E函数曲线如下:# y6 S5 c1 x. Z
. l7 Q; Z: h0 v" ^: i; L( e; S" O
可见, 对这个例子来说, 即使不用单边极限也能求出函数极限值是12。: C; ]) C! o6 K8 ] I- T: W: T6 j* E
5 V4 k A/ L: I# p5 T
- L = limit((exp(x^3)-1)/(1-cos(sqrt(x-sin(x)))),x,0)4 l0 \2 W+ K. K- B! Q+ W$ Q+ a
$ K! H- S) }, }$ H& B
2 _% H/ @9 }. h1 o/ r' Z
求函数 t a n t tan t tant 在 π / 2 \pi/2 π/2 点处的左右极限。9 ^1 ` @4 m, {
- syms t; f=tan(t);
- L1=limit(f,t,pi/2,'left')
- L2=limit(f,t,pi/2,'right')" T: Y ]' l5 V. }7 p/ q0 P
- k" w, L# P& b$ g) t/ @% b
% L1 l2 A+ y( \5 W求下面序列的极限% n& i' ?( D2 `$ t; M
$ M) _2 F' ^2 w6 R. M8 T- syms n positive
- f = n^(2/3)*sin(factorial(n))/(n+1);
- F = limit(f,n,inf). R% j6 e! n8 k7 b7 {7 _2 o9 N
& U7 i( I' ?( {
$ v" U+ p* l3 q* h( \% ?/ P" K0 R r0 E* A求下面序列函数的极限
- V3 @4 h! p/ ?: b) N$ u! ~
" i; ]8 J9 H2 ]7 Y; H- syms x n
- f = n*atan(1/(n*(x^2+1)+x))*tan(pi/4+x/2/n)^n;
- F = limit(f,n,inf)
& Y5 \; x# E: _5 {0 p
) ?2 B1 U" j) b( }; |' ^$ ^+ o7 P# U2 c: s( q9 q! M
多变量函数的极限
2 r- T* @) M& m5 ]+ gmatlab实现方法* v, H% m8 K* L) p
多元函数的极限也可以同样用MATLAB中的limit()函数直接求解。/ e' W0 s9 s8 C( l7 X! i2 u& v
0 q; ~, T r, ?+ ]
假设有二元函数 f ( x , y ) f(x,y) f(x,y), 若想求出二元函数的累极限
" g# I# v4 Z" u. G3 s
; D E2 J4 o# X& M* M' b
则可以嵌套使用limit()函数。例如:, Z* _* Y. r6 K" a: V6 {
: q& l+ u5 F; M9 a/ W& U
- L1 = limit(limit(f,x, x0), y, y0)
- L2 = limit(limit(f,y, y0), x, x0); P( X: z$ e P. m: ^
4 r7 ?. k: v+ E4 [1 w+ A* s
& W0 e3 C% V4 D2 x6 Y如果 x0或y0不是确定的值, 而是另一个变量的函数, 例如 x → g ( y ) x \rightarrow g(y) x→g(y), 则上述的极限求取顺序不能交换。
5 e8 s9 T; j# r& P L& N1 O
, e9 m1 m- @& r2 {' f/ T/ E假设有二元函数 f ( x , y ) f(x,y) f(x,y), 若想求出二元函数的重极限
$ @" s% ^/ M0 D1 z0 D/ g) l0 |9 A; f5 G' K' M. G2 J3 P1 f
7 L8 B, V; \; A5 v# I# N, Z5 _! Y" ~% I+ Z/ a
理论上不易求解,只有沿所有方向得出相同的极限才可,不可能用累极限方法求解。
% `* t9 o( ^7 j; V% ~5 W1 B; q, ^ x e( u1 N
应用举例
- p3 @+ l5 T" N8 B试求出二元函数极限值+ B# [/ S9 x+ @
% L" R/ z+ t) h# p- 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);
- L = limit(limit(f, x, 1/sqrt(y)), y, inf)
' d# i4 G: m: q) m7 n! o ; g9 T/ F; u5 I! s4 x6 I
# `6 V5 v, F3 ?9 ^9 N9 R* I5 J0 {- `重极限的尝试 ,求解重极限
3 W8 V) h" _5 f# d" C2 q6 F; b
5 k* r# n, s( E: k* X4 g
$ r" B2 C$ L- }2 b; X, M- 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)
8 q, G" w. q s/ w+ K# M0 W
+ G/ m$ h$ i4 e" | }; b& [7 |6 ^
, u2 V. _2 [8 i* r( o. E判断重极限是否存在
' V7 ^* o* p0 @# X# f- e s5 h
4 M- I6 h7 s1 J D+ l: N2 Z: o
证明极限不存在比求重极限容易的多,可以沿 y = k x y=kx y=kx趋近。$ e! K0 C" J$ O' E# x$ P9 S
7 U3 I7 G2 |6 E% d- d- syms r x y
- f=x*y/(x^2+y^2);
- L=limit(subs(f,y,r*x),x,0)4 n; H% } f2 V& H, b
5 W, ~1 a8 y a' R3 M |
|