Skip to content
Go back

导数的应用

Edit page

导数的应用

极值

绝对极值(最值)(最大、最小)

相对极值(局部)(f(c)=0\rightarrow f'(c)=0

临界点(f(c)=0f(c)=0\Leftrightarrow f'(c)=0或f(c)=0

极值相关定理

微分中值定理

我们可以看到,罗尔定理其实是微分中值定理的特殊情况。

图形的形状

增减函数检验

  1. 定义检验
  2. 一阶导数正负检验

局部极值的检验
一阶导:

  1. 局部极小值f(x)f'(x)从负变正
  2. 局部极大值f(x)f'(x)从正变负
  3. 没有局部极值,f(x)f'(x)两边同号

二阶导:

  1. f(c)=0f'(c)=0f(c)>0f''(c) \gt 0,那么ffx=cx=c处取得极小值
  2. f(c)=0f'(c)=0f(c)<0f''(c) \lt 0,那么ffx=cx=c处取得极大值

凹性

  1. 凹向上y>0y'' \gt 0
  2. 凹向下y<0y'' \lt 0

拐点

凹性改变的点,叫做拐点。 二阶导函数等于零,且一阶导函数有极大值或极小值的点是拐点。
但二阶导数等于零也不一定是拐点,二阶导数不存在也可能是拐点。

自治微分方程的图形解

导数只是关于 y 的函数,我们称该方程为自治微分方程。

使自治微分方程dydx=g(y)\dfrac{dy}{dx}=g(y)等于零的 y 值,称之为平衡点或静止点。

是一条用来表示函数的高阶导数的正负性的数轴。

借助相直线,我们可以得到自洽微分方程的图形解(没有表达式,但知道解的函数大致图像)。

三个案例

冷却

假设环境摄氏温度为PP,汤的温度为HH,经过的时间为tt,我们还假设环境足够大,使得汤的温度对环境的影响忽略不计。根据牛顿的冷却定律有:

dHdt=k(HP)\frac{dH}{dt} = -k(H - P)

H>PH \gt P时,dHdt<0\dfrac{dH}{dt} \lt 0

经过求导之后,画出相直线。
picture 41

微分方程的图形解如图:
picture 42

图中的 15 是一个稳定的平衡点,上下的两条直线都稳定朝 15 靠近。

分析有阻力时的落体

按照动量定理,我们有:

F=ddt(mv)=mdvdt+vdmdt=mdvdtF=\frac{d}{dt} (mv) = m \frac{dv}{dt} + v \frac{dm}{dt} = m \frac{dv}{dt}

当物体在介质中下落时,往往会受到阻力,阻力往往与速度相关,这时得到的公式为:

F=FgFv=m(gkv),a=gkvm.F = F_g - F_v = m (g - kv), \newline a = g - \frac{kv}{m}.

当物体受到的合力F=0F=0时,速度便稳定下来,故稳定速度为v=mgkv = \frac{mg}{k} 我们考察速度的变化过程:

d2vdt2=ddt(gkmv)=kmdvdt\frac{d^2v}{dt^2} = \frac{d}{dt} (g - \frac{k}{m} v) = - \frac{k}{m} \frac{dv}{dt}

画出相直线:
picture 43

画出图形解:
picture 44

阻力图形解

logistic 增长

假设P=P(t)P=P(t)为某个种群的数目,在时间Δt\Delta t内,种群的增长率为kk,则种群数目为: ΔPΔt=kP(t)\frac{\Delta P}{\Delta t}=kP(t) k>0k \gt 0是种群的增长率,对增长率进行假设:在自然界中,由于有限的资源,我们通常会假设一个极限种群承载容量MM,资源变得短缺了,种族的增长率kk便会降低,

k=r(MP)k = r(M-P)

rr代表了增长率变化的速度,对kk进行代换得到:

dPdt=r(MP)P=rMPrP2\frac{dP}{dt} = r(M-P)P=rMP-rP^2

该模型称之为logistic 增长模型

考察该模型的函数性态,我们对上述模型两边求导:

d2Pdt2=rMdPdt2rPdPdt=rdPdt(M2P)\frac{d^2P}{dt^2} = rM \frac{dP}{dt} - 2rP \frac{dP}{dt}=r \frac{dP}{dt} (M - 2P)

画出相直线得:

picture 45

logistic 增长相直线

画出图形解:
picture 46

logistic 增长图形解

我们可以看到,拐点在M2\frac{M}{2}位置处。

建模和最优化

书中使用了 Fermat 原理和 Snell 定律来进行建模举例。

Fermat 原理:光在介质中传播永远以时间最短的路径传播。
Snell 定律:在介质之间相互传播的时候,入射角和折射角满足一定的关系:

sinθ1c1=sinθ2c2\frac{\sin \theta_1}{c_1} = \frac{\sin \theta_2}{c_2}

线性化和微分

如果ffx=ax = a可微,那么近似函数

L(x)=f(a)+f(a)(xa)L(x) = f(a) + f'(a) (x - a)

就是ffaa线性化 近似f(x)L(x)f(x) \approx L(x)ffaa标准线性近似,点aa就是线性中心。

y=f(x)y = f(x)是一个可微函数,微分dxdx是一个自变量。微分dydy

dy=f(x)dxdy = f'(x)dx

变量dydy永远是一个因变量,它既依赖于xx,又依赖dxdx.

变化的敏感度

df=f(x)dxdf = f'(x)dx,我们知道f(x)f'(x)可以衡量ffdxdx的敏感程度。(这一点可能在多变量的情形下起作用)

微分变化估计实际变化的误差

Δfdf=f(a+Δx)f(a)f(a)Δx=Δx(f(a+Δx)f(a)Δxf(a))\begin{aligned} \Delta f - df & = f(a + \Delta x) - f(a) - f'(a) \Delta x \\ & = \Delta x (\dfrac{f(a + \Delta x) - f(a)}{\Delta x} - f'(a)) \end{aligned}

f(a+Δx)f(a)Δxf(a)=ϵ,\dfrac{f(a + \Delta x) - f(a)}{\Delta x} - f'(a) = \epsilon, 则上式 =Δxϵ= \Delta x \epsilonΔx0\Delta x \rightarrow 0 时,ϵ0\epsilon \rightarrow 0.

Newton 法

本节介绍一个求解根的方法:牛顿(Newton)法

通过切线的零点,确定下一条切线,再根据下一条切线的零点,确定再下一条切线…如此循环往复,最后得到要找的根。如下图:
picture 47
Newton 法 由线性化可知,牛顿法代数上计算就是根据给定的x0x_0求得f(x0)f(x_0)f(x0)f'(x_0)再代入:

x=xnf(xn)f(xn)x = x_n - \frac{f(x_n)}{f'(x_n)}

然后循环以上步骤,直至xnx_n收敛于某个固定的值(也可能不收敛,此时f(x)f(x)无零点。

Newton 法失效的情况:

求得的根不是想要的根的情况:

picture 48

Newton 法不是想要的根

我们使用 Newton 法解 z61=0z^6 - 1 = 0,我们根据不同的起始点,经过多次迭代之后得到如图:
分形盆
这一个个盆形区域被称为分形盆,这六个分形盆他们分别导向这个复数方程的六个解1,1,±3/2,±1/21,-1, \pm \sqrt{3} / 2, \pm 1/2


Edit page
Share this post on:

Previous Post
导数
Next Post
Linux命令杂记