找回密码
 注册
关于网站域名变更的通知
查看: 765|回复: 1
打印 上一主题 下一主题

matlab求解极限问题(limit函数的用法)

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-1-22 18:38 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

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

该用户从未签到

2#
发表于 2021-1-22 18:59 | 只看该作者
matlab求解极限问题(limit函数的用法)
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-11-24 01:48 , Processed in 0.156250 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表