微积分入门

摘要:
本文将会更加深入地研究加法,尤其是连续加法(integral)。简单探讨无限的意义。

格式约定:

这是正文。

这是引用,导言,注释或拓展内容

斜体只是调侃,活跃气氛,没有实际意义


一种朴素的极限思想?

一个概率小问题

一些连续的求和问题,与面积无关……

(可以跳过这一小节,对全文完整性没有影响。)

例题2.1 给定一条长为 1 的线段,在上面等均匀随机地(uniform)取 2 个点,这两个点之间距离的数学期望是多少?

数学期望

数学期望(mathematic expectation)是试验中每次可能结果的概率乘以其结果 x 的总和,记作 E(x)。通俗地讲,就是将随机事件进行足够多次并取每次结果的平均值。

例如,在 3,5,8 三个数中随机选一个数 x,则 x 的数学期望为

E(x)=33+53+83=163

如果抽出 3,5 的概率都是 25%,抽出 8 的概率为 50%,则 x 的数学期望为

E(x)=3 25%+5 25%+8 50%=6

再例如,掷一枚公平的六面骰子,其每次“点数”的期望值是 72,计算如下:

E(x)=116+216+316+416+516+616=72

值得注意的是,72 虽然是“点数”的期望值,但却不属于可能结果中的任一个,没有掷出此点数的可能。

在不产生歧义的前提下,数学期望也可以简称为期望(mean)。

:不妨称所选两点分别为点 A 和点 B。

在一条线段任意的取点似乎无从下手,因此我们不妨先考虑有限的弱化情况:假如我们在线段上均匀的取两个点(其实就是线段的两个端点,称之为点 1 和点 2),并要求所选两点必须在点 1 或点 2 上

/ B 在点 1 B 在点 2
A 在点 1 0 1
A 在点 2 1 0

则距离(称为 d)有 2 种不同情况,其中距离为 012 的概率,距离为 112 的概率。此时,

E2(d)=012+112=12

(我们使用角标 2 表示这是钦定 2 个点的弱化情况)

假如我们在线段上均匀的取三个点(线段的两个端点和线段的中点,称之为点 1、点 2 和点 3),并要求所选两点必须在点 1、点 2 或点 3 上

/ B 在点 1 B 在点 2 B 在点 3
A 在点 1 0 12 1
A 在点 2 12 0 12
A 在点 3 1 12 0

则距离)有 3 种不同情况,其中距离为 013 的概率;距离为 1249 的概率;距离为 129 的概率。此时,

E3(d)=013+1249+129=49

通过类似的方法,还可以计算出

E4(d)=512E5(d)=615

我们注意到,随着钦定的点数 n 增大,答案 En(x) 逐渐趋近于 13。我们知道,钦定的点越多,算出的值越接近在连续的线段上选点的真实答案,因此你可能会猜测

limnEn(d)=13

lim 记号

lim 是英语单词 limit 的缩写,它指的是变量在一定的变化过程中,从总的来说逐渐稳定的这样一种变化趋势以及所趋向的值。例如

limx1x=0

表示当 x 足够大时,1x 趋近于 0

x 1 10 100 1000 10000
1x 1 0.1 0.01 0.001 0.0001

这是因为,对于不存在一个确定的数 n>0,使得无论当 x 多大时都有 1x>n。(换句话说就是,对于任意确定的数 n>0,总能找到一个数 x,使得 1x<n,并且所能选的所有 x 的最小值随着 n 减小而减小)

但是这如何证明呢?我们不妨假定在线段上均匀地取 n 个点(当 n 越来越大时答案会越来越接近真实值),称之为 D1,D2,D3,,Dn,并要求所选两点必须在 D1,D2,D3,,Dn 中的某个点上。

显然此时原线段被 Di 分割为 n1 段,每段长度都为 1n1。那么距离有 n 种不同情况,即对于 k{0,1,2,,n1},距离可能是 kn1。我们不妨记距离为 kn1 的概率为 Pd(k)。那么所求的期望为

En(d)=0Pd(0)+1n1Pd(1)+2n1Pd(2)++1Pd(n1)

我们知道,概率是发生该事件的情况数除以总情况数。总情况数显然有 n2 种。观察上表,我们注意到每种距离的情况都是在 2 条斜线上的(除 0 外,只有中间 1 条斜线):

n = 3 的情况

容易看出,距离为 0 的情况有 n 种,1n12(n1) 种,2n12(n2) 种,……,即距离为 kn1 的有 2(nk) 种情况。由此可得

Pd(k)={1nifk=02(nk)n2if1kn1

进一步,代入到 En(d)

En(d)=0Pd(0)+1n1Pd(1)++n1n1Pd(n1)=01n+1n12(n1)n2++12(n(n1))n2=2n2(n1)(1(n1)+2(n2)++(n1)(n(n1)))=2n2(n1)(16(n3n))=n3n3n2(n1)=13+13n

那么当 n 足够大,我们有

limnEn(d)=13

由此,我们确信,所求之期望就是 13w5

路程 速度 加速度

汽车速度表上的速度代指是什么?

  1. 假如我们有一个匀速直线运动的物体,0 时刻时位置 s=0,其速度 v=5t 个时刻后物体的速度、位置分别是多少?

    t 0 1 2 3 4 5
    v 5 5 5 5 5 5
    s 0 5 10 15 20 25

    学过函数的同学,很容易就能看出,速度和位置满足

    v(t)=5s(t)=5t

    函数图像

  2. 假如我们有一个匀加速直线运动的物体,0 时刻时位置 s=0,速度 v=0,其加速度 a=2t 个时刻后物体的速度、位置分别是多少?

    t 0 1 2 3 4 5
    a 2 2 2 2 2 2
    v 0 2 4 6 8 10
    s 0 1 4 9 16 25

    学过物理的同学,很容易就能看出,加速度、速度和位置满足

    a(t)=2v(t)=2ts(t)=t2

    函数图像

导数

以上例子中,函数 a(t),v(t),s(t) 之间的关系是什么?

不难发现,加速度 a(t) 本质上是速度 v(t) 的速度,而速度 v(t) 是路程 s(t) 的速度。因此 a(t)v(t) 是这样一种函数,它指示另外一个函数的增长速度

指示原函数 f(x) 在自变量 x 极小的变化时因变量 f(x) 的变化速度的函数,称为导函数(derivative function),简称导数(derivative),记作 f(x)Df(x)。形式化的

f(x)=limh0f(x+h)f(x)h

觉得这个定义不专业,不严谨?你可以看一看数学分析教材,那里的定义要多专业有多专业,除了看不懂之外没有任何缺点。

瞬时速度

导数是函数的瞬间变化率(instant change rate)。这听起来很矛盾,为什么既有“瞬间”,又有“变化”呢?我的意思是变化是在时间中体现的,当在一个时刻里,物体就没有变化的空间了。

我们知道时刻 a 至时刻 b 的速度的定义是

v=sbsaba

但是我们刚才定义的速度函数 v(t) 就很奇怪了,它只需要一个参数时刻参数 t 就可以得到一个速度。但是我给一辆行使中的小汽车拍一张照片,你不可能说出他的速度吧?

汽车的速度仪表盘本质上是计算了时刻 t 到时刻 t+dt 的速度,其中 dt 是一个非常小的值,例如 0.01 秒。

例如当 s(t)=t2,t=3 时,若取 dt=0.01v(3)=s(3+0.01)s(3)0.01=6.01,当 dt 越来越小时,v(t) 就会越来越接近 s(t) 的瞬间变化率:

dt 1 0.1 0.01 0.001 0.0001
v(3) 7 6.1 6.01 6.001 6.0001

我们注意到,随着 dt 越来越小,v(3) 的值越来越接近 6,也就是说

v(3)=s(3+dt)s(3)dt=6

由此,我们对“瞬间变化率”有了初步的认识,更进一步的,

v(t)=s(t+dt)s(t)dt=(t+dt)2t2dt=t2+2t dt+(dt)2t2dt=2t+dt

由于 dt 是一个极小的,趋近于 0 的数,因此可以忽略不计,即

v(t)=2t

上面那个飞矢不动的悖论吗?人家古希腊几千年前就研究过啦

记号 dx

dx 在微积分中很常用,表示 x 的一个极小变化量。也就是说 dx 默认趋近于 0。因此,下面两种写法是等价的:

f(x)=limh0f(x+h)f(x)h

f(x)=f(x+dx)f(x)dx

不定积分

之前的例子中,函数 a(t),v(t),s(t) 之间的关系是什么?

不难发现,速度 v(t) 本质上是加速度 a(t) 的某种连续的和,而路程 s(t) 是速度 v(t) 的某种连续的和。因此 v(t)s(t) 是这样一种函数,它指示另外一个函数的连续和

几何上看,这种连续的和是函数下方区域的面积,下图为 v(t)=2ts(3) 所对应的面积。

面积示意图

关于函数 f(x) 下方的面积的函数,称为函数 f(x)不定积分(indefinite integral),记作 f(x)dx

求导

函数千千万万,怎么计算每一个的导数?

单项式的导数

我们刚才已经知道,对于 f(x)=x2,其导数 f(x)=2x,那么对于任意一个单项式 f(x)=xn 呢?

f(x)=xn(n0),则其导数

f(x)=f(x+dx)f(x)dx=(x+dx)nxndx

(x+dx)n 展开后的 xn 项与后面的 xn 抵消,nxn1dx 项保留下来了,除去分母之后变为 nxn1,而剩下的众多项都至少包含 (dx)2,除去分母之后任然包含 dx,因此忽略不计。所以我们得到

f(x)=nxn1

定理2.2 上面这条规律(不规范地)记作

D(xn)=nxn1(n0)

sin 和 cos 的导数

让那些被“解析”、“代数”夺去的直观意义,在这里补回来吧!

f(x)=sin(x)f(x)=cos(x),那么函数的导数是什么呢?

例题2.3 求函数 f(x)=sin(x) 的导数。

我们考虑使用几何方法,在单位圆(unit circle)上讨论这个问题。

单位圆

以原点为圆心,半径长度为 1 的圆称为单位圆。单位圆在三角学、复数乃至数论等众多方向有重要的应用。

单位圆的解析式为

x2+y2=1

我们学习三角函数时就已经知道

sin2(θ)+cos2(θ)=1

由此我们就知道了单位圆的一个重要性质:单位圆上的点都可以写成 (cos(θ),sin(θ)) 的形式,同样的,形如 (cos(θ),sin(θ)) 的点都是在单位圆上。

解: 现在我们想知道,

f(x)=limh0sin(x+h)sin(x)h

图下所示,点 C 在单位圆 A 上,使得 FAC=xAx 轴交于点 B(1,0),过点 CCFAB 于点 F,做 A 的切线 CD,使得 CAD=h,过 DDECF

几何图示

由于 FAC=x,容易知道 C(cos(x),sin(x))。当 h 足够小时 D 几乎在 A 上,所以我们是要求出 CE 的长度,即

sin(x+h)sin(x)=CE

通过圆的基础知识,我们知道连接切点和圆心的半径与切线相互垂直,即 ACCD,所以

AFCCED

根据相似三角形的性质可知,

CE=AFCDAC=cos(x)sin(h)1=cos(x)sin(h)

代入上式可得

f(x)=limh0sin(x+h)sin(x)h=limh0cos(x)sin(h)h=cos(x)limh0sin(h)h

所以现在只要求出 limh0sin(h)h 的值,即可知道 f(x)=sin(x) 的导数。

绘制一个中心角为 h 弧度、半径为 1扇形(sector) S。这将用于建立可用于缩放定理的不等式:

扇形图

通过扇形的面积公式,我们知道扇形的面积是 S=h2

如下图,画三角形 TS 的角相同,其底面为 1,高度为 sin(h)。这意味着 T 的面积是 T=sin(h)2。由于 T 包含在 S 中,其面积必然更小:

扇形图2,包含三角形 T

如下图,画半径为 cos(h) 的扇形 S,其的面积为hcos2(h)2。因为 S 包含在三角形 T 中,所以它的面积必然更小:

扇形图3,包含扇形 S'

由此我们得到不等式

12h(cos2(h))sin(h)2h2

每一项同时乘以 2h 得到,

cos2(h)sin(h)h1

h=0 时,有

cos2(h)=(cos(0))2=1

1sin(h)h1

由夹逼定理可得

limh0sin(h)h=1

代入之前的式子,得到

f(x)=cos(x)

至此,我们就求出了 sin(x) 的导数(即 cos(x))。w5

试一试2.4 尝试仿造上面的方法(或自己设计新方法)求出函数 f(x)=cos(x) 的导数。

答案: f(x)=sin(x)

这两条规律记作(分被称为定理2.3定理2.4

D(sin(x))=cos(x)D(cos(x))=sin(x)

自然常数

传说中有一个人,叫做 FISHER_,拥有无穷无尽的神力和高深莫测的魔法。他掌管这世界上一切吉祥的和不祥的事情的发生和停止……

FISHER_ 在 3 月 5 日(神渔标准时间)的开始时刻重量为 1 个单位,从此之后他的体重以每天翻一倍的速度增长,那么他的体重关于自 3 月 5 日以来的天数的函数 f(x) 是什么?

注意到,这里的增长速度就是函数 f(x) 的导数。不假思索地,你可能会回答

f(x)=2x

我们来验证一下:由于体重每天翻一倍,函数 f(x) 的导数应该和自己一样,即

f(x)=f(x)

对于函数 f(x)=2x,它的导数我们目前还不知道怎么求,但是从数值上拟合 x=5 处的导数,可以发现似乎 f(x)=2x 的导数并不是自己(看起来这个数趋近于 22.1807097779

x 5
2x 32
dx 0.00001
D(2x) 22.1808
dt 0.1 0.01 0.001 0.0001 0.00001
D(2x) 22.9675 22.2578 22.1884 22.1815 22.1808

我们尝试对 f(x)=2x 求导

f(x)=f(x+dx)f(x)dx=2x+dx2xdx=2x2dx2xdx=2x(2dx1dx)

注意上面式子里括号中的一部分,2dx1dxx 的取值无关,数值计算得到它是一个大约为 0.69 的常数。因此 f(x)=2x 的导数并不是他它自身,而是自身的约 0.7 倍。

这么说来,f(x)=2x 是一个不错的尝试,我们已经很接近正确答案了。

那么那个函数的导数是它自己呢?或者说以什么数为底数能使得后面那部分的常数恰好等于 1 呢?不妨设未知数 e,使得函数 f(x)=ex 满足 f(x)=ex,即

D(ex)=ex

在解出 e 的值之前,我们最好先搞清楚为什么 f(x)=2x 的导数不是它自己,难道第二天的重量翻倍不是变为 2 个单位吗?

如果“天”为时间的最小颗粒单位,那么 1 天之后,重量确实为 2 个单位;

但是如果时间的最小颗粒单位为 12 天,那么每 12 天,重量变为原来的 (1+12) 倍,所以 1 天后 FISHER_ 的重量为 (1+12)2=49,比 2 稍大;

相似地,那么如果时间的最小单位为 1n 天,那么每 1n 天,重量变为原来的 (1+1n) 倍,所以 1 天后 FISHER_ 的重量为 (1+1n)n,在本问题中,时间是连续的,因此 1 天之后,FISHER 真正的重量为

f(1)=limn(1+1n)n

由于 f(x)=ex,所以 f(1)=e,解得

e=limn(1+1n)n

所以 e 等于多少呢,这很难说,简单的数值拟合可以发现 e(2,3),而且离 3 比较接近。事实上,e 是一个无理数,我们只知道它的值大约是 2.71828limn(1+1n)n 被称为自然常数(natural constant),记作 e

这里的常数 e 和上面未知数的 e 还用了不同的字体?搞这么多花里胡哨的玩意儿!

还有一个问题没有搞清楚:上面函数 f(x)=2x 的导数到底是多少呢?

我们要搞清楚下面式子,

limh02h1h

就可以就出 f(x)=2x 的导数。

例题2.5 若函数 f(x)=2x,就函数 f(x) 的导数 f(x)

解: 由我们之前的推导已经知道

f(x)=2xlimh02h1h

k=ah1,则 h=log2(k1),换元得

f(x)=2xlimh0klog2(k1)

注意到当 h 趋近于 0k 也趋近于 0,因此可以转化为

f(x)=2xlimk011klog2(k1)

应用定理1.10,

f(x)=2xlimk01log2(k1)1k

转化 lim 的位置

f(x)=2x1log2limk0(k1)1k

现在只需要知道 limk0(k1)1k 即可,我们根据 e 的定义

limk(1+1k)k=e

1k 替换 k

lim1k(1+k)1k=e

由于当 1k 时,k0,最终得到

limk0(1+k)1k=e

代回上面式子

f(x)=2x1log2e

代入 1=log22,应用换底公式(定理1.12),

f(x)=2xln2

自然对数

在应用中,以 e 为底对数十分常见,称为自然对数(natural logarithm),记作 ln。即

alna=loge2

至此,我们就求出了 2x 的导数(即 2xln2)。w5

将上述过程中的 2 替换成任意其他常数,可以得到

定理2.6 对于任意常数 c,函数 f(x)=cx 的导数 f(x)=cxlnc

那么假如指数不是 x,而是 2x 或者其他的什么玩意儿呢?

例题2.7 若函数 f(x)=e2x,求函数 f(x) 的导数 f(x)

这个问题相当简单,用现实意义来解释,就是 FISHER_ 达到某个重量只用原来一半的时间,此时 FISHER12 天就能使重量翻一倍。这个问题用已经研究过的知识可以轻松解决,请读者自己动手试一试。

:根据乘方的基本性质,有

f(x)=(e2)x

c=e2 代入定理2.6,

f(x)=(e2)xln(e2)

化简即得答案 f(x)=2e2xw5

函数和

x2+x 的导数是多少?

我们知道 D(x2)=2x,D(x)=1,那么 D(x2+x) 是否就是 D(x2)+D(x)=2x+1

从这一节开始,我们将讨论通用的求导法则。我们再看一下导数的定义

f(x)=f(x+dx)f(x)dx

导数就是因变量的变化量除以自变量的变化量,我们把自变量的变化量记作 dx,那么我们不妨把因变量的变化量记作 df。所以导数就是

f(x)=dfdx

假设我们有函数 f(x)g(x),想要求函数 (f+g)(x) 的导数 (f+g)(x)

函数加减运算

为了方便,我们记

(f+g)(x)=f(x)+g(x)(fg)(x)=f(x)g(x)

(f+g)(x)=(f+g)(x+dx)(f+g)(x)dx=f(x+dx)+g(x+dx)f(x)g(x)dx=df+dgdx=f(x)+g(x)

由此我们得出,两个函数的和的导数等于这两个函数的导数的和,即

定理2.8 对于两个函数 f(x),g(x),有 (f+g)(x)=f(x)+g(x)

函数积

x2ex 的导数是多少?

我们知道 D(x2)=2x,D(ex)=ex,那么 D(x2ex) 是否就是 D(x2)D(ex)=2xex

现在我们有函数 f(x)g(x),想要求函数 h(x)=f(x)g(x) 的导数 h(x)

无论什么时候,把代数中的乘法与几何中图形的面积联系起来,往往都是很好的。绘制一个长方形,其长宽分别为 f(x)g(x)

长方形图

现在我们 x 增加一点点 dx,则 f(x) 增加了 dfg(x) 增加了 df。总的来说,面积增加了 f(x)dg+g(x)df+df dg。其中 df dg 是极小的,可以忽略不计。

长方形增加长宽图

由此我们知道了

定理2.9 (乘积求导法则)若我们有函数 f(x)g(x),则函数 h(x)=f(x)g(x) 的导数是

h(x)=f(x)g(x)+f(x)g(x)

函数复合

函数的复合运算

函数复合就是将一个函数的输出作为下一个函数的输入的过程,即

(fg)(x)=f(g(x))(fgh)(x)=f(g(h(x)))

假如现在我们有两个函数 f(x)g(x),求 (fg)(x)

由于有

f(x)=f(x+dx)f(x)dx

x=g(x),并代入,得到

(fg)(x)=f(g(x)+dg)f(g(x))dx

x=g(x) 后代入 f(x)=f(x+dx)f(x)dx

(fg)(x)=f(g(x))dgdx

代入 g(x)=dgdx

(fg)(x)=f(g(x))g(x)dxdx

分子分母约去 dx 就得到了答案,即

(fg)(x)=f(g(x))g(x)

定理2.10 (链式求导法则)对于两个函数 f(x)g(x),其复合函数 (fg)(x) 的导数 (fg)(x)=f(g(x))g(x)

小结

对于如何求导的研究至此告一段落了,通过定理2.2 到定理2.10,我们知道了

常见函数的导数:

原函数 f(x) 导数 f(x) (限制条件)
xn nxn1 /
sin(x) cos(x) /
cos(x) sin(x) /
cx cxlnc c>0

求导法则:

组合方式 导数
f(x)+g(x) f(x)+g(x)
f(x)g(x) f(x)g(x)+f(x)g(x)
f(g(x)) f(g(x))g(x)

熟练运用上面的规律,已经可以求出大量函数的导数了,请读者自己动手算一算下面例题

例题2.11 已知函数 f(x)=2xx2+x,求该函数的导数 f(x)

解: 对函数 g1(x)=x2 应用定理2.2,

g1(x)=2x

对函数 g2(x)=x 应用定理2.2,

g2(x)=1

对函数 g3(x)=g1(x)+g2(x) 应用定理2.8,

g3(x)=2x+1

对函数 g4(x)=x 应用定理2.2(x=x12),

g4(x)=x122

整理得

g4(x)=12x

对函数 g5(x)=(g3g4)(x) 应用定理2.10并整理,

g5(x)=2x+12x2+x

对函数 g6(x)=g2(x)g5(x) 应用定理2.9并整理,

g6(x)=x(4x+3)2x(x+1)

对函数 g7(x)=2x 应用定理2.6,

g7(x)=2xlnx

对函数 f(x)=(g6g7)(x) 应用定理2.10并整理,

f(x)=2xx(x+1)1x(4x+3)ln2x(x+1)

即为函数 f(x) 的导数。w5

此题计算量确实较大(我用了将近 30 分钟),但是思路其实是相当简单的:反复运用这些法则,把复杂的函数转化为基本的函数。

积分

加法和减法,乘法和除法……

回顾先前的例子,我们知道函数 v(t)s(t) 的导数,反过来函数 s(t)v(t) 的(不定)积分。这使得我们意识到,求导和积分互为逆变换。也就是说,

定义2.12 对于给定的函数 f(x),如果存在一类函数 F(x) 使得 F(x)=f(x),则称所有 F(x) 中常数项为 0 的那个函数为 f(x)不定积分

例如,对于函数 f(x)=x2,其不定积分为 F(x)=13x3

定积分

考虑这样一个问题:

例题2.13 求函数 f(x)=x2 图像下方从 x=1x=3 区域(换句话说就是 y=x2,x=1,y=0,x=3 四条(曲)线围成的图形的面积)的面积。

x^2下方的面积图像

解: 我们使用与最开始那个概率问题相似的方法,不妨先用几个矩形拟合一下 y=x2 的曲线:

2个矩形的拟合图像

4个矩形的拟合图像

16个矩形的拟合图像

50个矩形的拟合图像

100个矩形的拟合图像

容易发现,所使用的的矩形越多,拟合结果越精准。不妨假设现在用 n 个矩形进行拟合,这 n 个矩形的面积记作 Sn。每个矩形的宽为 31n=2n,第 k 个矩形的高为

f(1+2kn)=(1+2kn)2=1+4kn+4k2n2=1+4nk+4n2k2

所以我们要求的就是

Sn=2n(f(1+2n)+f(1+4n)+f(1+6n)++f(1+2nn))=2n(1+1+1++1)+8n2(1+2+3++n)+8n3(12+22+32++n2)=2nn+8n2n(n+1)2+8n3n(n+1)(2n+1)6=263+10n+43n2

那么当 n 足够大时,

limnSn=263

故该面积为 263w5

这种有明确的上下界的积分(上面例题中上下界为 13),叫做定积分(definite integral)。记作 abf(x)dx。例如例题2.13就是要求

13x2dx

例题解答中所使用的定积分求法较为复杂,一个简单的办法是使用牛顿-莱布尼茨公式(Newton-Leibniz formula,又名微积分基本定理)。

定理2.14 (牛顿-莱布尼茨公式)一个连续函数在区间 [a,b] 上的定积分等于它的任意一个原函数(不定积分)在区间 [a,b] 上的增量。即如果有 F(x)=f(x),则

abf(x)dx=F(b)F(a)

使用这个方法求定积分变得十分简单,还是以 13x2dx 为例:注意到

D(13x3)=x2

固有

13x2dx=13321312=263

极限

y=sin(x)x 的图像是什么样子的?在 x=0 时应该如何画?

类似于 f(x)=sin(x)x 这样的函数在一些特殊点无法直接取值,但是其值任然是有意义的。

sin(x)/x的图像

我们注意图中红圈标记处,直接计算会得到 sin00 这样没有意义的结果,但是我们从图像中可以很清晰的得出 f(0)=1 的结论,这是因为

limx0sin(x)x=1

这就使用了极限(limitation)的概念。

极限的标准定义

f:DRR 是一个定义在实数上的函数。 并在某个开区间 x>Ax<A 上有定义。 L 是一个给定的实数。 c 是一个实数,并且函数 fc 的某个去心邻域上有定义。如果对任意的正实数 ϵ ,都存在一个正实数 δ ,使得对任意的实数 x ,只要 f 在点 x 处有定义,并且 xc 的某个 δ -(去心)邻域中(即 |xc|δ ),就有 |f(x)L|ϵ ,那么就称 L 是函数 fx 趋于 c 时的极限,或简称 Lfc 的极限,记为 limxcf(x)=L 。反之则称 L 不是 fx 趋于 c 时的极限。

来自 维基教科书(啥都看不懂系列

注意,极限要求从函数两边趋近能达到同样的结果。举个例子,函数 f(x)=1xx=0 时的值是未定义的,(如下图红线所示)因为从小于 0 出发回趋向于 而从大于 0 出发回趋向于 +

y=1/x的极限示意图

洛必达法则

微积分需要连续连续需要无穷小,可是谁知道无穷小长什么样子?
———— 伯特兰·罗素

之前我们尝试求三角函数的导数时,曾经用相当巧妙的几何方法求出过

limh0sin(h)h=1

那个方法虽然巧妙,但是相当复杂,而且不具备通用性,有较高的思维难度。下面介绍洛必达法则(L'Hospital's rule),此方法能够较为简单地求出大多数函数的极限。

定理2.15 (洛必达法则 描述简化版)对于一个形如

limxaf(x)F(x)

的极限,如果它满足

  1. x 趋向于常数 a 时,函数 f(x)F(x) 都趋向于 0
  2. a 附近(可以理解为 x=adxx=a+dx),f(x)F(x) 都有定义,且 F(x)0
  3. limx0f(x)F(x) 有定义。

则有

limxaf(x)F(x)=limxaf(x)F(x)

证明:(通俗不严谨版,某种程度上说,仅供感性理解)

由于在 a 附近,f(x)F(x) 都有定义,我们知道,所求的极限也是有定义的。

由于当 x=0 时,f(x)=0,F(x)=0,所以我们知道

limxaf(x+dx)=f(x)limxaF(x+dx)=F(x)

x 的变化量足够小时(即变化量为极小的 dx),应该有

f(x)=limh0f(x+h)F(x)=limh0F(x+h)

由此得出

f(x)=f(x+dx)F(x)=F(x+dx)

代入之前的式子,

limxaf(x)=f(x)limxaF(x)=F(x+dx)

由此相除得到

limxaf(x)F(x)=limxaf(x)F(x)

即为所求证的等式。

严谨的证明请见 参考文献 - 洛必达法则的证明。

泰勒展开

没有其他外力,没有初速度的条件下,质点只能静止在原地。毫无自由可言。

给质点一个初速度,我们可以使质点单向匀速运动;

若再给定一个加速度,我们可以使速度均匀变化,从而产生拐弯运动;

若再给定加速度的变化率,我们使加速度均匀变化,速度拐弯变化,产生可转向拐弯运动;

……如果一开始就设定好质点的初速度,加速度,加加速度,加加加…加速度的话,正如用一只无形的手调控着它的命运,那么无论想让它何时拐,往何处拐,如何拐……就全都在初始条件的设计之中了。

  1. 考虑函数 f(x)=sin(x),它的导数为

    f(x)=cos(x)

    x=0 时,显然有 f(0)=0,f(0)=1。但是满足上述条件的函数不止一个,至少还有

    f(x)=x

  2. 考虑函数 f(x)=cos(x),它的导数为

    f(x)=sin(x)

    x=0 时,显然有 f(0)=0,f(0)=1,f(x)=0。但是满足上述条件的函数不止一个,至少还有

    f(x)=x

  3. 考虑函数 f(x)=sin(x),它的导数为

    f(x)=cos(x)

    x=0 时,显然有 f(0)=0,f(0)=1,f(x)=0,f(x)=1。但是满足上述条件的函数不止一个,至少还有

    f(x)=x16x3

  4. 考虑函数 f(x)=cos(x),它的导数为

    f(4)(x)=sin(x)

    高阶导数记号 f(n)(x)

    对原函数求 1 次导数,我们记作 f(x);以此类推求 2 次导数,就记作 f(x);但是当求很多次导数时,这种记号阅读和书写都不方便,所以我们用 f(n)(x) 表示对原函数 f(x)n 次导数后的函数。

    x=0 时,显然有 f(0)=0,f(0)=1,f(x)=0,f(x)=1,f(4)(x)=0。但是满足上述条件的函数不止一个,至少还有

    f(x)=x16x3

以此类推,我们不停的对 f(x)=sin(x) 求导,会得到一个循环,即

f(n)(x)={sin(x)ifnmod4=0cos(x)ifnmod4=1sin(x)ifnmod4=2cos(x)ifnmod4=3

而且我们总能找到一个多项式函数,使得当 x=0 时这个多项式的前 n 阶导数与 sin(x) 的导数相同,例如多项式

x5120x36+x

sin(x) 的前 6 项导数相同,而多项式

x1139916800+x9362880x75040+x5120x36+x

sin(x) 的前 12 项导数相同,这使得它们的函数图像在 x=0 附近极为接近(如下图)。

泰勒展开值与真实值的对比图

我们之所以总能够构造出这样的多项式,是因为多项式的每求导一次,最高次项的次数减 1,并且其导数是完全可控的(由多项式的系数决定)。这些系数的分母是为了抵消求导导致的系数增加(下式中加粗的 n

D(xn)=nxn1

每次求导时,指数 n 减一,导致系数一层层一共乘大了 n(n1)(n2)1,即 n!,所以分子要除去 n!

综合上面描述,我们把这种多项式称为原函数的泰勒级数(Taylor's series),把这种展开过程称为泰勒展开(Taylor's expand)。

定理2.16 (泰勒公式)函数 f(x)x=x0 处的泰勒级数为

f(x)=f(x0)0!+f(x0)1!(xx0)+f(x0)2!(xx0)2+f(x0)3!(xx0)3+

特别地,当取 x0=0 的级数尤为常用,可以称为麦克劳林级数0 处的泰勒级数。此时有

f(x)=f(0)0!+f(0)1!x+f(0)2!x2+f(0)3!x3+

如果我们求 f(x)=ex 的麦克劳林级数,并代入 x=1,可以得到 e 的一个渐进的表达形式,极为有趣

e=10!+11!+12!+13!+

习题

  1. 已知二次函数 y=x22x8 与直线 l 交于 A(4,0),B(0,8),在二次函数上有一动点 P,且点 P 在直线 l 上方运动,求 PAB 面积的最大值。
  2. [0,1] 区间内等概率地随机选取一个数 x,求 x2 的数学期望 E(x2)
  3. f(x)=tan(x) 的导数 f(x)
  4. f(x)=x2+2x+x2+2+x 的导数 f(x)
  5. f(x)=x2+2x+x2+2+x4x4 时的取值范围。
  6. 求出 f(x)=cosx 的麦克劳林级数的前 5 项。
  7. limx1sin(πx)x21 的值。
  8. f(x)=arcsin(x) 的导数。
  9. 求半径为 r 的球的表面积公式(用含 r 的式子表示)。
  10. 平面上的解析式 y2sin(x)=x 上有一点 A(m,n),求 A 点处该图像的切线的斜率。
  11. 求函数 f(x)=lnx 的导数 f(x)
  12. 在边长为 1 的正方形内等概率随机选取两点,它们之间的距离的数学期望是多少?

参考文献

posted @   方而静  阅读(531)  评论(1编辑  收藏  举报
相关博文:
阅读排行:
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· 记一次.NET内存居高不下排查解决与启示
点击右上角即可分享
微信分享提示