高中可能讲了,但是高中讲了不太可能。
莱布尼茨记法
高中一般都喜欢用拉格朗日记法来表示导数(如 f′ 这种),但是有一个弊端是不好表示自变量和因变量。而莱布尼茨记法就能更好地展现变量之间的关系。
比如 f 的关于 x 的一阶导写作:dxdf,二阶导:dx2d2f,三阶导:dx3d3f。
这里的 dxdf 其实就相当于 limΔx→0ΔxΔf 的简写,因此在运算中可以像分数那样操作。
如链式法则,将 f[g(x)] 对 x 求导:
dxdf=dgdf⋅dxdg
又比如 y=f(x),则很显然可以看出:
dydx⋅dxdy=1
这是反函数求导法则。
e
假如说你有 100 元,存到某银行,1 年之后利率 100%,那么一年之后你就有 200 元。
但是若能分 2 次结算利率,每次利率 21,则两次结算利率之后,就能得到 100×(1+21)2=225 元。
如果能分 3 次结算利率,每次利率 31,则最后能得到 100×(1+31)3≈237 元。
好像越来越多了,那么当结算利率的次数无限多,最终收益就会变成:
n→∞lim100(1+n1)n=100e
这就是 e 最早、最经典的定义就是用极限:
e=n→∞lim(1+n1)n
高中教材中曾经可能有这个定义,但是现在删掉了。
此外,我们可以用均值不等式证明数列 (1+n1)n 单调递增:
(1+n1)n=1⋅(1+n1)n<(n+11+n(1+n1))n+1=(1+n+11)n+1
高中都学过二项式定理,我们把这个 e 的表达式展开:
e=n→∞lim(1+n1)n=n→∞limi=0∑n(in)ni1=n→∞limi=0∑ni!nini1
这里 ni 是 n 的 i 次下降幂的意思。当 n→∞ 的时候,可以认为 ni=ni,于是就有:
e=i=0∑∞i!1
看到 e 不由自主想要求幂,还是类似上面的推导:
ex=n→∞lim(1+n1)nx=n→∞limi=0∑nx(inx)ni1=n→∞limi=0∑nxi!xi=i=0∑∞i!xi
证明 ex 的导数等于其自身:
(ex)′=Δx→0limΔxex+Δx−ex=exΔx→0limΔxeΔx−1
把 eΔx 展开:
exΔx→0limΔx1+1!Δx+2!Δx2+⋯−1
高端无穷小丢掉,就是:
exΔx→0limΔx1+Δx−1=ex
从而我们就证明了 (ex)′=ex。
定义 lnx 为 ex 的反函数。要给 lnx 求导,不妨令 y=lnx,也有 x=ey。
dxdy=dydx1=ey1=x1
由此就证明了 (lnx)′=x1。
e 太伟大了。现在我们可以来求 xn 的导数,但是这个其实不好搞,我们先对它取个对数再 exp。
(xn)′=(enlnx)′=(nlnx)′⋅enlnx=xn⋅xn=nxn−1
并没有用到广义二项式定理。
可以求 ax 的导数:
(ax)′=(exlna)′=(xlna)′⋅exlna=axlna
甚至还可以求 fg 的导数:
(fg)′=(eglnf)′=(g′lnf+fgf′)fg
天不生 ex,万古如长夜。
洛必达法则
由于是高中,不用讨论函数能不能导问题。
先从一个很显然的东西出发。
罗尔定理:若函数 f(x) 满足 f(a)=f(b),则 ∃c∈(a,b) s.t. f′(c)=0。
假如说 f(x) 是常函数,那么显然有导数为 0 的地方。
否则,f(x) 中间一定有最大值或者最小值,不妨假设有最大值 f(x0)。则 x0 左侧导数:
Δx→0lim−Δxf(x−Δx)−f(x)≥0
右侧导数:
Δx→0limΔxf(x+Δx)−f(x)≤0
这两个导数要相等,所以在 x0 处的导数应该为 0。
这个是 f(a)=f(b) 的情况,那么如果 f(a)=f(b) 该怎么办呢。一个简单的思路是把 f(x) 减去过 (a,f(a)),(b,f(b)) 的直线,即 y=b−af(b)−f(a)(x−a)+f(a)。这样就转化成了罗尔定理的情况。
令 g(x)=f(x)−b−af(b)−f(a)(x−a)−f(a),根据罗尔定理,这个函数在 (a,b) 之间也有导数为 0 的地方。对 g(x) 求导,得:
g′(c)=f′(c)−b−af(b)−f(a)=0
即:
f′(c)=b−af(b)−f(a)
这就是拉格朗日中值定理。
直观来看中值定理,说的就是平面上固定两点的可微曲线,一定存在一点使得这点的切线斜率和两端点之间的斜率相同。
若这个线甚至不是函数,而是类似 (f(t),g(t)) 描述的曲线,类似的结论呢?
首先我们要知道这种曲线怎么求切线斜率,模拟求导:
k=Δt→0limf(t+Δt)−f(t)g(t+Δt)−g(t)=Δt→0limΔtg(t+Δt)−g(t)⋅f(t+Δt)−f(t)Δt=f′(x)g′(x)
即我们希望证明 (a,b) 之间存在一点 c,使得 f′(c)g′(c)=f(b)−f(a)g(b)−g(a)。
模拟拉格朗日中值定理的证法,搞一条过 (f(a),g(a)) 和 (f(b),g(b)) 的直线,然后相减。即构造:
h(t)=g(t)−f(b)−f(a)g(b)−g(a)(f(t)−f(a))−g(a)
h(a) 和 h(b) 都为 0,套用罗尔定理,存在 c∈(a,b):
h′(c)=g′(c)−f(b)−f(a)g(b)−g(a)f′(c)=0
即:
f′(c)g′(c)=f(b)−f(a)g(b)−g(a)
当然,有分数不太好,因为很可能出现分母为 0 的情况,即斜率不存在。我们还是最好全部化成乘积形式:
g′(c)(f(b)−f(a))=f′(c)(g(b)−g(a))
这是柯西中值定理。
柯西中值定理的形式让人想起洛必达法则。现在我们要求 limx→ag(x)f(x),但是 f(a)=g(a)=0。
对于任意一个 b>a,根据柯西中值定理,我们知道在 (a,b) 中存在一点 c 满足:
g′(c)f′(c)=g(b)−g(a)f(b)−f(a)=g(b)f(b)
当 b→a,也有 c→a。所以就有:
x→alimg(x)f(x)=x→alimg′(x)f′(x)
泰勒级数
退役 OIer 还不会泰勒展开正常吗?但是我就是不会。
我们希望用一个多项式函数去逼近一个复杂函数。具体来说,我们设我们要逼近的函数 f(x) 可以用 (x−a) 为基底的无穷次多项式表示:
f(x)=c0+c1(x−a)+c2(x−a)2+c3(x−a)3+⋯
其中 ci 是待定的 i 次项系数。
为了逼近,首先在 x=a 处式子左右两边一定要相等:
f(a)=c0+c1(a−a)+c2(a−a)2+c3(a−a)3+⋯
这样就得到了 c0=f(a)。
继续,在 x=a 处不但要相等,我们还希望它们的一阶导数也要相等:
f′(a)=c1+2c2(a−a)+3c3(a−a)2+4c4(a−a)3+⋯
于是 c1=f′(a)。
二阶导:
f′′(a)=2c2+6c3(a−a)+12(a−a)2+20(a−a)3+⋯⇒c2=2f′′(a)
……
规律不言自明了:ci=i!f(i)(a),这里 f(i)(a) 是 f(x) 的 i 阶导在 x=a 处的值。
也可以写作:
f(x)=i=0∑∞i!f(i)(a)(x−a)i
特别地,当 a=0 的时候,是马克劳林级数:
f(x)=i=0∑∞i!f(i)(0)xi
来试试把一些常用函数给泰勒展开一下:
首先是 ex,由于其任意次导数均为 ex,而 e0=1,则可以容易得到:
ex=1+x+2x2+3!x3+⋯+n!xn+⋯
和上面我们推的吻合。
然后来试试 sinx。sinx 的导数有周期性,分别是 sinx,cosx,−sinx,−cosx,它们在 0 处的值分别为 0, 1, 0, -1。容易得到:
sinx=x−3!x3+5!x5−⋯+(−1)n(2n+1)!x2n+1+⋯
cosx 同理,有:
cosx=1−2!x2+4!x4−⋯+(−1)n(2n)!x2n+⋯
考虑将 (1+x)α 展开:
(1+x)α=1+αx+2!α(α−1)x2+3!α(α−1)(α−2)x3+⋯+n!αnxn+⋯
这个结果和广义二项式定理的一致。
多项式有个好处是可以很自然地将扩展到复数域。
eix=1+ix−2x2−3!ix3+⋯+n!(ix)n+⋯
其偶数次项正好是 (−1)k(2k)!x2k,也是 cosx 的展开项;奇数次项是 (−1)ki(2k+1)!x2k+1,是 isinx 的展开项。
于是就有了欧拉公式:
eix=cosx+isinx
微分方程
f(x)=f′(x) 的唯一解是 f(x)=Cex,其中 C 是常数。
假设有其它函数满足这个条件,假设为 g(x),那么我们考虑 h(x)=exg(x)。h′(x)=e2xg′(x)ex−g(x)ex=0,即 h(x) 是个常数,所以得到 g(x) 也属于 Cex。
由于导数之间可以线性相加,所以如果一个微分方程有若干解 f1(x),f2(x),f3(x),⋯,则它们的线性组合 C1f1(x)+C2f2(x)+C3f3(x)+⋯ 也是解。
物理中有很多物理量之间有导数的关系。比如位置 x,对时间 t 求导是速度 v,再对 t 求导是加速度 a。求解这些量之间的关系,就相当于解微分方程。
单杆问题
考虑这样的问题:
平行导轨间距为 L,连通且忽略电阻,有垂直与两导轨所成平面的匀强磁场 B。有长度为 L,电阻为 R 导体棒放置在平行导轨上,其初速度为 v0,方向沿导轨方向。问导体棒什么时候停下来。
导体产生的电动势:
E=vLB
产生的电流:
I=RE=RvLB
产生的安培力:
F=ILB=RvL2B2
牛二:
F=ma
a 是什么?a 是 v 的导数,所以我们相当于有这样一个方程:
−RvL2B2=mv˙
这里 v˙ 是牛顿的导数记号,表示对时间求导。
我们只需要求出 v 关于时间的表达式,然后解出 v=0 对应的 t 就行了。
看到这个导数形式让人不由得想起了 e,不妨假设 v=Cert,只需解出 C 和 r。 C 很简单:当 t=0 的时候 v=v0,所以 C=v0。
解 r,将 v=Cert 直接带入:
−RertL2B2=mrert
解得:
r=−mRL2B2
所以有:
v=v0e−mRL2B2t
但是指数函数永远不会变成 0,所以导体棒永远不会停下来。
有没有可能还有其它的解呢?类似前面关于 f(x)=f′(x) 的解唯一的讨论,可以证明没有其它解。
简谐振动
又来考虑简谐震动。简谐振动中有回复力:
F=−kx
还有牛二:
F=ma=mx¨
设 x=ert,带入可得:
−kert=mr2ert
解得:
r=±imk
带回去:
v1v2=cos(mkt)+isin(mkt)=cos(mkt)−isin(mkt)
为什么有两个解?速度怎么还有虚部?实际上,v1 和 v2 的所有线性组合都能满足上面的方程,但我们应该选择符合实际情况的。
假如振幅为 A,初相位为 φ。
先看怎么搞出来一个 φ,根据复数乘法规则(模长相乘,幅角相加),乘一个 cosφ+isinφ:
(cosφ+isinφ)v1=cos(mkt+φ)+isin(mkt+φ)
为了能抵消掉虚部,给 v2 乘个 cosφ−isinφ:
(cosφ−isinφ)v2=cos(mkt+φ)−isin(mkt+φ)
相加,再乘一个 2A 就能得到:
v=Acos(mkt+φ)
所以发现简谐运动的周期:
T=mk2π=2πkm
最后结果可以看出来,初相位和振幅都不影响周期。而且在前面解出 r=±imk 时,你就已经可以知道 ω=mk 了,因为后面的线性操作都没法影响到 ω。
LC 振荡电路
最后我们看到 LC 振荡电路。
LC 振荡电路由一个电容和一个线圈组成。对于电容我们有:
C=UCQ
对于线圈,我们也有:
UL=LI˙
一条回路电压和应该为 0:
UC+UL=0
同时我们还知道 I 是 Q 的导数:
I=Q˙
联立得微分方程
CQ+LQ¨=0
和简谐振动那个几乎一样,解得:
Q=e±iLC1
所以 LC 振荡电路的周期:
T=2πLC
后面不会了。