鹤啸九天 自律更自由,平凡不平庸 Less is More

数学基础:微积分系列

2020-06-08
鹤啸九天
阅读量

Notes(温馨提示):

  1. ★ 首次阅读建议浏览:导航指南
  2. 右上角工具条搜索文章,右下角二维码关注微信公众号(鹤啸九天),底栏分享、赞赏、评论
  3. ★ 转载请注明文章来源,知识点积累起来不容易,水滴石穿,绳锯木断,谢谢理解
  4. ★ 如有疑问,邮件讨论,欢迎贡献优质资料


总结

【2021-11-3】平稳转动的几何形状不止圆形,比如:勒洛三角形,介绍,不能做轮子因为内有固定轴。

  • 勒洛三角形应用案例:
    • 转子发动机,也叫汪克尔引擎,汪克尔在1957年做了第一个转子发动机原型机。马自达(以前好像叫东洋工业)好像是1961年有的授权,不过马自达研究比任何车企都用心而已,奔驰以前也研究过,没结论
    • 什么形状的钻头,能钻出方孔? 答案就是勒洛三角形,这个近似三角形的几何图形,绘制起来也十分简单,以一个正三角形的三个顶点为圆心,边长为半径画圆,三条边围绕的三段弧形,即可组成勒洛三角形
  • 勒洛三角形是什么?将水杯放在三角形上推动,有趣的现象发生了

椭圆

神奇的椭圆,三种不同定义(拉伸,斜切圆锥体,画线)最终都能得到椭圆

导读


目录


微积分(Calculus)引言

微积分回忆

  • 求导公式
  • 乘积法则
  • 链式法则
  • 隐函数求导
  • 积分、微分的互逆关系
  • 泰勒级数

微积分的三个中心思想:

    1. 积分
    1. 微分
    1. 积分与微分(导数)的互逆

(几位)微积分之父

  • 发现微积分:巴罗(Barrow)、牛顿(Newton)、莱布尼茨(Leibniz)
  • 给出严格定义:柯西(Cauchy)、魏尔施特拉斯(Weierstrass)

本系列视频中,微分的概念全部包含于导数之中,所有导数都以 df/dx 的形式出现,也就是微商

推导圆的面积 - 积分的直观理解

圆的面积公式:

如何从积分的角度推导出圆的面积公式?

  • 不同的划分方法会带来不同的积分公式,下面考虑将圆划分为大量的同心圆环,这种方法保留了圆的对称性。

  • 考虑其中一个环的面积,可以将其看做一个“类矩形”

  • 虽然这不是标准的矩形,但只要dr越小,它就越接近。它的面积可表示为:

  • 于是,圆的面积可以看作是这一系列矩形面积的叠加。

  • 这部分面积的求和可以等价于求“函数y = 2πr图像在区间[0, R]下的面积”。

  • 这个推导的过程其实可以看作是对函数y = 2πr[0, R]下的积分。

积分与导数

直观来说,对函数f(x)[a, b]积分就是求函数f(x)在区间[a,b]下的图像与坐标轴包围的面积。记作:

\[\int_{a}^{b}f(x)dx\]

这实际上是定积分的概念,此外还有不定积分。

如果是其他图像,比如抛物线,该怎么求这部分的面积呢?

能不能找到一个函数 A(x) 表示 0x 之间函数图像下的面积?——这个函数 A(x) 就是该函数的积分(函数)。

这里强调 0x 之间,是为了使问题具有实际意义

以抛物线 f(x)=x^2 为例。类似的,我们可以将这块区域划分成一系列细长的矩形。

x 增加 dx,增加的面积可以看做是一个长f(x)、宽dx的矩形,只要dx越小,这条窄带就越接近矩形。

把这部分面积记作 dA,表示面积的微小变化(difference in Area)

通过这个矩形,可以得到 Af(x)dx 之间的关系:

\[{\displaystyle{\begin{aligned}dA\approxf(x)dx\\\Rightarrow\frac{dA}{dx}\approxf(x)\\\Rightarrow&\frac{A(x+dx)-A(x)}{dx}\approxf(x)\end{aligned}}}\]

这里引出了微积分中另一个重要的概念——导数dA/dx 就是 “A 的导数”

更严格的说法是:”A 的导数”是“当 dx → 0 时,dA/dx 所趋向的值”。(下一节会讨论导数的定义)

一般不会刻意区分导数导数函数的区别,都统称为导数,具体含义视语境而定;同样,积分也是如此。

导数是解决积分问题的关键——积分需要还原出某个导数原本的函数——如果你能熟练的计算导数,那么你也能解决这个问题。

积分与导数之间的这种互相转化的关系,也就是“某个图像下方面积函数的导数能够还原出定义这个图像的函数”,就叫做微积分基本定理。该定理表明,“在某种意义上”,两者互为逆运算。

导数(Derivative)的意义

总结:导数、偏导数、方向导数、梯度

【2022-8-23】直观理解梯度,以及偏导数、方向导数和法向量等

总结

  • (1)导数是一元函数的变化率(斜率)。导数也是函数,是函数的变化率与位置的关系。
  • (2)如果是多元函数呢?则为偏导数
  • 偏导数是多元函数“退化”成一元函数时的导数
    • “退化”的意思是固定其他变量的值,只保留一个变量,依次保留每个变量,则𝑁元函数有𝑁个偏导数。
  • (3)如果是方向不是沿着坐标轴方向,而是任意方向呢?则为方向导数
    • 如下图所示,点𝑃位置处红色箭头方向的方向导数为黑色切线的斜率; 链接Directional Derivative (可交互体验)
  • (4)方向导数取得最大值,最大值为梯度的模;最小时,是梯度负数
    • 梯度的几何意义:
      • 当前位置的梯度方向,为函数在该位置处方向导数最大的方向,也是函数值上升最快的方向,反方向为下降最快的方向;
      • 当前位置的梯度长度(模),为最大方向导数的值

以二元函数为例,令𝑧=𝑓(𝑥,𝑦),绘制在3维坐标系如下图所示

  • 在分别固定𝑦和𝑥的取值后得到下图中的黑色曲线——“退化”为一元函数,二维坐标系中的曲线——则偏导数∂𝑧∂𝑥和∂𝑧∂𝑦分别为曲线的导数(切线斜率)
  • 一个变量对应一个坐标轴,偏导数为函数在每个位置处沿着自变量坐标轴方向上的导数(切线斜率)

目标函数的等高线图示意图,来自 Applet: Gradient and directional derivative on a mountain

  • 红点为当前位置,红色箭头为梯度,绿色箭头为其他方向,其与梯度的夹角为𝜃。将左图中𝑧=𝑓(𝑥,𝑦)曲面上的等高线投影到𝑥𝑦平面,得到右图的等高线图。梯度与等高线垂直

总结

  • 偏导数构成的向量为梯度;
  • 方向导数为梯度在该方向上的合成,系数为该方向的单位向量;
  • 梯度方向为方向导数最大的方向,梯度的模为最大的方向导数;
  • 微分的结果为梯度与微分向量的内积
  • 等高线全微分的结果为0,所以其梯度垂直于等高线,同时指向高度更高的等高线
  • 隐函数可以看成是一种等高线,其梯度为高维曲面(曲线)的法向量

“瞬时变化率”引起的歧义——导数的悖论

“瞬时变化率”的歧义:

  • 只有在不同的时间点之间,变化才能发生;
  • 而将时间限制在某个瞬间点的时候,变化也就不存在了。

考虑这个示例:一辆车从 A 点起,先加速,再减速,至 100 米外的 B 点停下,整个过程花费 10 秒。

把车速的图像加入其中,可以发现:两者存在着某种联系:其中之一改变的话,也会引起另一个发生变化。

速度的大小是如何随着距离-时间函数的变化而变化的?——“瞬时”速度的矛盾

  • 在绘制速度图像的时候,需要给每个单独的时间点关联一个速度值,但是计算速度却需要两个时间点上的距离。
  • 记时间差为 dt,距离差为 ds,那么这段时间内的速度就能用 ds/dt 表示

    当年的微积分创始人们也经历了同样的思维冲突。

  • 换言之,速度只有在一段会时间内的才有意义;“瞬时”的说法会带来矛盾
  • 实际的做法是:会选取一个很小的 dt 值,然后把 ds/dt 看做是这个瞬间的速度。

导数的定义与计算

导数的定义

  • 在纯数学领域,导数不是 dt 为某个具体值时 ds/dt 的值 ,而是当 dt 的值无限逼近 0 时这个比值的极限
    • dt 非常非常小,无限逼近0;两点连线的斜率→某一个点的斜率。
    • 左上角长方形里的内容才是导数的完全体。
    • dt 逼近0时,后面两项就能完全忽略了
  • 从图像的角度,(某一点的)导数有一个直观的含义:就是经过图像上该点切线的斜率。
  • 注意:这里的 dt 不是“无穷小”,也不是 0;它永远是一个有限小的量,接近 0 而不是 0。

    这种说法在试图规避“瞬时”带来的矛盾,使“某个时间点的变化率”有意义。

导数的计算

  • 抛开求导公式,先来看一下面对一个实际的问题,该如何求解(在某一点处的)导数。
  • 对于 s(t)=t^3t=2 处的导数,根据导数的定义,有
  • dt 趋向 0 时,后两项也会趋于 0,进而消去。代入 t=2 可以得到在该点处的导数为 12
  • 更一般的,称 s'(t) = 3*t^2s(t) = t^3导函数

对于常见的函数,有一系列总结出的求导公式可以快速计算(下一节会演示如何从几何的角度来推导这些公式)

导数的含义是“变化率的最佳近似”

  • 回顾之前关于距离-速度的示例,思考这个问题:当 t=0 时,车在不在移动?
  • 一方面,利用导函数公式可以得到 t=0 时的速度为 0——这似乎在说“车没有移动”;另一方面,如果车在 0 时刻没有移动,那么它是何时开始移动的?——关键在于这个问题本身就是没有意义的。
  • 因为导数并不是用来测量“瞬时变化”的。
  • t=0点的导数为 0 的真正含义是指“在第 0 秒附近,车速的最佳近似为 0 米/秒”——换句话说,就是当时间间隔 dt 越来越小时,表示速度的比值 ds/dt 就越趋向于 0——这并不表示车在 0 时刻就是静止的,只能说它此时的速度近似于 0.

用几何来求导

为什么导数很重要?——当需要使用微积分来解决现实中的实际问题时,需要将其抽象成各种代表性的函数来描述;而如果能掌握这些抽象函数的变化率,那你就学会了这门可以精准描述事物变化率的语言。

从几何的角度看“微小变化量”

f(x) = x^2 为例:

  • 从坐标轴上看,x^2 的图像是一条抛物线,我们已经知道,导数可以描述为切线的斜率
  • 此外,x^2 还有另一个更直接的含义:长为 x 的正方形的面积。
  • 假如给边长 x 一个微小的增量 dx,那么正方形的增量(变化量)是多少?

  • 应该时刻记住的 dx 是一个微小的量——这意味着你可以忽略所有次数高于 1 的 dx——换言之,一个微小量的平方(或更高次方)是一个可以忽略的变化量
  • 由此,可以得到:
  • 上一节给出了 f(x) = x^3 导数的代数推导过程,这里也可以作为立方体体积来用几何的方式推导。
  • 示例

负的变化量

  • 上述两个例子的变化量都是增量(正值),但实际上变化量也可能是负值,比如 f(x) = 1/x——考虑这样一个特殊的矩形,长 x,宽 1/x,它的面积恒为 1.
    • 注意:这里的变化量不再是矩形的面积了,而是矩形的高

  • 通过简单的几何知识,可知矩形的变化量为:
  • 从而得到 1/x 的导数为:

幂函数的导数

  • 以上 x^2x^3x^-1 都遵循了幂函数的求导公式:

  • 思考一下为什么这个公式也适用于 2 和 3 以外的指数——求导的一个关键点在于很大一部分项因为包含 dx 的高次幂,可以被忽略

  • 当指数小于 0 时,会更复杂一些——比如 x^-2 可以考虑一个边长为 √x 的正方形。

大多数时候,我们都有合适的求导公式来使用,但是以上利用几何求导的过程能锻炼我们借助微小变化量来考虑的导数的能力。

三角函数的导数

  • 正弦函数的定义
    • θ=0.8 表示单位圆弧长为 0.8 时对应的角度为 θ(单位圆的周长为
    • sin(θ) 表示此时该点距离 x 轴的高度(可能为负)
    • θ 增加时,sin(θ) 的值会在 -1 到 1 之间上下摆动
  • sin(θ)θ 增加时的微小变化量

    • 在近似的观点下,应该有这样的直觉——可以把那一段微小的圆弧看作是直线
  • 根据三角函数的定义,有

直观理解链式法则

求导公式概览(3)

  • 函数的和/差——加法法则
  • 函数的积/商——乘法法则
  • 复合函数的求导法则

加法法则

  • dx 是很小很小的数,为了直观画成比较宽的样子
  • 加法法则比较简单,可以在坐标轴中展示

乘法法则

“左乘右导 + 右乘左导”:

在数学中,如果你要处理两项的乘积,用面积来理解会更方便

复合函数的求导法则——链式法则

链式法则

链式法则的表达式在说:看看输出值 g 的微小变化除以 h 的微小变化是多少(h 是要带入函数 g 中的值);然后乘以 h 微小变化与 x 微小变化的比值

这些 dh 最终会被消去,结果就是输出值 g 的微小变化与输入值 x 的微小变化的比值

dh 的消去并不只是符号上的技巧,而是真实反映出在求导时,各微小变化量发生了什么(?)

了解这些求导法则 != 灵活的使用它们

指数函数的导数——自然常数 e 的定义

2^t开始

  • 根据导数的定义,有

  • dt → 0 时,h 趋向于一个常数 0.6931...
  • 也就是说,2^t 图像上各点处切线的斜率 = 该点的函数值 * 一个常数

    原视频对 2^t 赋予了一些实际意义(质量)来说明其变化量与自身成比例

  • 类似的,3^t 也存在这样一个 h=1.0986...8^t 则是 h=2.0794...
  • 注意到,8^th 大概是 2^t 的三倍,而 8^t = 2^{3t}

是否存在某个数,使 h 恰为 1?——这个数就是自然常数 e

  • 不要问为什么当 e^t 时,h=1;因为 e 就是如此定义的。

    和“为什么 π 正好等于圆的周长比直径”一样,π 就是这么定义的。

    自然常数不是完全这样发现的,但过程类似

指数函数的导数

  • 有了自然常数 e 以及链式法则,就可以求出其他指数函数的导数了
  • 根据指数函数的性质与链式法则,有

  • 事实上,在微积分中,指数函数基本都是以 e^ct 的形式出现的,很少会直接使用 c^t 的形式

为什么称 e 为自然常数?

  • 现实世界中,很多自然现象里的变化率变化量是成正比的
  • 比如:在室温环境下,热水变凉的速率与水和房间的温差成正比(或者说,温差的变化率与温差本身成正比)

  • 尽管指数函数有多种写法,但是将其表达为以 e 为底的幂函数是非常自然的。
  • 因为 e 有着非常自然的含义,它就是变化率与数量本身的比例系数。(?)

隐函数求导

什么是隐函数?

  • 如果方程 F(x,y)=0 能确定 y 是 x 的函数,那么称这种方式表示的函数是隐函数;
  • 直观来说,就是满足某种关于变量 x 和 y 的关系的、所有 (x, y) 点的集合,相应的曲线就是“隐函数曲线”。

示例 1:圆上某一点切线的斜率

  • 圆的方程就是一个隐函数,比如 x^2 + y^2 = 5^2,下面需要求圆上一点的斜率(不考虑几何方法)

  • 因为这里的曲线并不是一个函数图像,所以不能单纯对其求导——它不存在变量 x 的微小变化对函数值 y 的微小变化
  • 当然,如果我们的目标只是求 dy/dx,那么对等式两边同时求导可以解决:

  • 为什么这么做?

示例 2:相关变化率

  • 这是一个更实际的例子:开始时,梯子上端距地面 4m(y=4),下端距墙 3m(x=3)。如果梯子的上端以 1 m/s 的速度下滑,那么下端左滑的速度是多少?

  • 等式的左边可以看做是一个关于时间 t 的函数,于是对表达式左边求导的含义相当于在问“经过短暂的时间 dty 会减少一点,x 会增加一点,那么整体的变化量是多少?”,而表达式右边告诉我们这个变化量为 0
  • 带入 x(t)=3, y(t)=4, dy/dt=1,就能求出 dx/dt 了。

隐函数求导的含义

  • 从纯数学的角度看,示例 1 与示例 2 的做法没有区别。但是示例 2 求导时带有明确的物理意义——表达式随时间的变化率。但是求圆切线时似乎难以解释为什么这么做。
  • 其实隐函数的导数也是一个隐函数,它对 dxdy 的关系做了限制——在圆的例子上,它要求整体变化量 dS = 2xdx + 2ydy 为 0;也就是说为了让 S = x^2 + y^2 始终保持在 25,那么 dS 就需要为 0.

    严格来讲,这个条件实际上是保证每一步落在过该点的切线上,而不是落在圆本身。当 dxdy 足够小时,这两者才没有区别。

    在隐函数的导数中,原来的变量看做常数,比如这里 2xdx + 2ydy = 0 中的 xy

    显然,(左)的正体变化量就不满足使 S 保持不变

示例 3:sin(x)y^2 = x

  • 该函数的图像是一系列 U 型曲线

    有了这个等式,就能方便求出 dy/dx 了。

  • 对表达式的两边求导,就可以得到导数的隐函数,它对 xy 的整体微小变化量作出了限制。

示例 4:ln(x)的导数——从已有的导函数推算出其他函数的导函数

  • 根据指数函数与对数函数的关系,有

多元微积分

  • 隐函数求导是多元微积分的入门。两者的要点是一样的,需要理解这多个变量是如何联系在一起变化的。

极限

简单来说,“极限”就是逼近/趋近更“洋气”的说法。

导数的正式定义

  • 在[2.2. 导数的定义与计算]已经给出了导数的计算公式:

    导数的正式定义:

    为什么不用 dx?——dx 本身已经表达了求极限的含义——使用 dx 的表示容易将其理解成“无穷小的变化量”;更正确的解读应该是把它看做一个具体的、有限小的变化量,并时刻考虑 dx → 0 时的情况(3. 用几何来求导)。

    为什么用一个新的变量 h?——明确变化量 h 只是一个普通的数,跟“无穷小”没有任何关系。

极限的 (ϵ, δ) 定义

  • 所谓极限,指的是变量逼近 0 时的影响,而非无穷小变量的影响
  • 导数由极限定义,而极限本身由 (ϵ, δ) 定义,下面是比较书面的说法:

设函数f(x)在点x0的某一去心邻域内有定义,如果存在常数a,对于任意给定的正数ϵ,都∃δ>0,使不等式|f(x)-a| < ϵ|x-x0| ∈ (0,δ)恒成立,那么a就叫做函数f(x)x→x0时的极限,记作

  • 右边的式子就是导数的正式定义。h 和 dx 是一回事。
  • dx 可以解读为一个具体的,有限小的变化量。随时考虑 dx 逼近于 0 时的情况。
  • 讨论极限时,指当变量逼近于 0 时的影响,而非无穷小的变化量的影响

极限存在与不存在的两个例子

极限存在的前提是,你总能在距离极限点x0距离为δ(>0)的取值范围内,找到一系列取值点,使得范围内的任一取值点,它的函数值都在距离f(x0)ϵ(>0)的范围内——关键在于这种情况对于任意ϵ都成立——无论ϵ多小,你总能找到与之对应的δ

利用导数来求极限——洛必达法则

如果函数在点x0处是有定义的,那么该点的极限值 == 函数值本身。

如果函数在该点没有定义呢?比如 sin(πx)/(x^2-1),该函数在x0=±1处就没有定义。

洛必达法则

  • 洛必达法则专门用于求解 0/0 型和 ∞/∞ 的极限值。
  • 导数是由极限定义的,但反过来,也能利用导数来求极限
  • 洛必达法则利用了这样一个性质——当f(a)=0时,在a的附近,有f(a+dx)=f'(a+dx)。如图所示:

    dx 越小,这个比值就越精确

  • 洛必达法则:

    条件:1)f(x)g(x)在点a处的极限都为 0;2)f(x)g(x)在点a的某去心领域内可导,且g'(x)!=0;3)A可以为实数,也可以为±∞

    注意,只要满足以上条件,洛必达法则是可以继续进行的

  • 这么看,实际上,求解导数的过程也是在使用洛必达法则

    分子分母在 h→0 处的极限都为 0

积分与微积分基本定理

1.2. 积分与导数就提到过,积分其实就是求导的逆运算

示例:距离是速度对时间的积分

  • 已知速度-时间的关系 v(t)=t(8-t),求这段时间驶过的距离。
  • 在[2.1. 节]介绍这个例子时,是“已知距离函数s(t),来推导速度函数v(t)”;现在“已知每个时间点的速度 v(t),来求距离s(t)
  • 从数学的角度,前者的问题等价于求“s(t)的导数”,于是现在的问题也就是求“哪个函数的导数是v(t)”——这通常被称为“求函数的‘原函数’(或‘反导数’)”

  • 类比[1.1. 节]中推导圆面积的公式,这个问题可以转化成求解曲线下方的面积

  • 这个过程可以表示为 v(t) 的“积分”:

    为什么不适用?——指的并不是具体的dt的加和,而是当dt→0时,加和逼近的值,这个逼近的值就是曲线下方的面积;当dt越趋近 0,加和就越趋近精确解。

    积分的含义就是将所有小量“累积”在一起

  • “求函数图像与横轴围成的面积”是许多不相干问题的共通之处——可以被拆成小量,然后近似为这些小量的和

积分与导数是一组互逆的运算

  • 根据以上的分析,距离实际上是速度对时间的积分,而速度是距离对时间的导数——这表明积分与导数是一组互逆的运算

面积函数的导数等于函数本身

  • 如果固定左端点,将右端点当做一个变量 T,那么这个积分可以看做以上限为自变量的函数 s(T)——这是距离关于时间的导数,同时也是图形下方的面积函数
  • 这表明任一函数图像下方面积函数的导函数等于该函数本身

    • 从几何的角度看,矩形的面积为 ds,宽为 dT,它的高为 v(T),于是有 ds/dT = v(T)
    • 图像所表示的函数就是面积函数的导数
  • 因此,求解积分的过程就是求导的逆过程

每个函数都有无数个原函数

  • 因为常数的导数为 0,所以在原函数的基础上加上任意常数,其导数不变——这意味着每个函数的原函数有无数个
  • 所以v(t) = 8t - t^2正确的原函数应该是

  • 从图像上来看,曲线上下移动并不会影响其在每一点的斜率

  • 在不同的问题中,一般会有额外的条件帮你决定该使用哪个原函数
    • 比如,在本问题中,t=0 时,s(t)=0,于是 C=0

微积分基本定理

  • 在对任意函数求(定)积分时,你是在对 x一定范围内的所有 f(x)*dx 求和,而积分值就是 dx → 0 时,这个和趋近的值。
  • 求积分的第一步是找出原函数 F,使其导数为积分内的函数。
  • 积分值就等于原函数在上限时的值减去其在下限时的值。
  • 这个过程就是“微积分基本定理”

  • 积分异于常识的一点在于:它连续地遍历了从下限到上限中的每一个自变量的值,而我们在利用原函数求值时,只需要关注上限与下限两个自变量

负面积

  • 如果图像有部分出现在横轴的下方,那么这部分就是负面积

    有符号的面积

  • 积分计算的不是真正的面积,而是图像与横轴围城的带有正负的面积。

连续变量的平均值

结论:区间上的平均斜率等于起点和终点连线的斜率

  • 函数的平均值提供了积分和求导是互逆运算的另一个角度。
    • sin(x)(高度)的平均值就是原函数(cos(x))从 x=0 到 x=pi 所有切线斜率的平均值。

如何求连续变量的平均值?

  • 等价于求图像下方面积的平均高度
  • 如果把面积作为一个整体来看,这是一个很简单的问题——平均高度 = 积分面积 / 上下限宽度

  • 值得注意的是,(F(b)-F(a))/(b-a) 实际上也是原函数F(x)x=ax=b 两点连线的斜率
  • 这表明在某一区间上所有点处切线的平均斜率等于起点和终点连线的斜率

    根据定义,f(x)是其原函数F(x)的导函数,也就是说它给出了F(x)在每个点上切线的斜率,所以f(x)(a,b)上的平均值也就是原函数从x=ax=b上所有切线斜率的平均值。

  • 换言之,求解连续函数的平均值,可以转化为求解其原函数在各点切线的平均斜率;而两点间的平均斜率等于这两点的斜率。

适用“积分”的场景


8章和第9章给出了两个适用于积分的场景

  • 可以通过细分然后相加的方式进行估算时
  • 求解连续变量的均值时(特别在概率中)

级数

高阶导数

这里介绍高阶导数的目的是帮助得到函数的近似——泰勒级数

二阶导数描述的是曲线的弯曲程度

  • 根据导数的定义,二阶导数也就是导数的微小变化率,即斜率的变化率——直观来看,也就是曲线的弯曲程度

    • 考虑函数的一个取值 x,然后向右连续的增加两个小量 dx
    • 第一个增量是函数产生了第一个变化量 df1,第二个同理,记 df2
    • 这两个变化量的差也就是函数值变化量的变化量,记 d(df). 它和 (dx)^2 成正比
    • 所谓二阶导数,就是 d(df)(dx)^2dx → 0 时比值的极限
  • 用符号表示为
  • 其实中间的写法才是最正确的,但为了书写方便,通常写成最右侧的形式

  • 一个理解二阶导数的现实示例就是加速度
  • 二阶导数为正,说明车在加速;反之,为负,说明在减速

级数类别

【2023-6-23】学数列极限找规律,都是为了让你学级数的

级数和微积分作为数学两大主要方向

无穷级数两种重要级数

  • 幂级数:常规函数,趋近任意精度,用于工程计算
  • 傅里叶级数:周期函数,拟合任意形状

幂级数和三角级数各有所长,不能相互替代。

幂级数

工程计算,比如:

  • 不用计算器计算任意精度的 $\sqrt{2}$

幂级数是多项式的推广。

  • 无穷次多项式,收敛域非常特殊,以某一点为中心的区间,在收敛区间内,求和函数是无限可微的。

使用泰勒公式将对数、指数、三角、反三角等拆解成无穷项

  • 每项都是常系数与N次幂构成
  • 后面的项绝对值趋近于0

幂级数的收敛性

  • 1)阿贝尔第一定理:如果幂级数收敛于x1,则 当 x 属于 |x1| 时,幂级数绝对收敛; 如果级数在x2处发散,则对于所有x>|x2|,所有幂级数发散
  • 2) 收敛半径下定理:在幂级数收敛半径R下,(-R, R)收敛 绝对,并且都在 [-R, R] 之外发散。 正或负 R 可能收敛或发散
  • 3) 给定一个幂级数 Σan(x-x0)^n,如果 lim |(an 1)/an| = p, 或 lim 3√|an | = p,则收敛半径R = 1/p

幂级数和函数的解析性质:

  • 1)阿贝尔第二定理:关于收敛和一致收敛的定理
  • 2) 连续性
  • 3) 逐项推导和逐项积分

幂级数的运算:

  • 1) 两个幂级数的概念相等,同次幂项的系数也相等
  • 2) 函数为奇函数,不会有偶次幂项,反之亦然
  • 3)两个幂级数的和与乘积

函数幂级数的展开:

  • 1)泰勒展开
  • 2)McLaughlin展开
  • 3) Lagrange型余数
  • 4) Cauchy型余数
  • 5) 几种常用初等的幂级数展开 功能

傅里叶级数

周期信号函数的展开

三角级数是三角函数无穷和的一种特殊形式。 三角级数所表达的函数可以是不可微的,甚至是不连续的,因此表达的函数范围很广。

三角级数:法国傅里叶提出,又称傅里叶级数

  • 随着三角函数累加次数增加,趋近于脉冲信号、锯齿波形

【2023-6-23】幂级数、傅里叶级数(notes)

泰勒多项式与泰勒级数

泰勒级数的作用——函数近似工具。比如,在 x=0 附近,有:

泰勒级数利用多项式函数去近似复杂抽象的函数,从而化简问题

  • 多项式函数的优势:易计算、易求导、易积分
  • 深度学习中的梯度下降(一阶)、牛顿法(二阶)也是利用的泰勒级数

示例:cos(x)x=0 附近如何用二次多项式近似?

  • 问题等价于在所有可能的 c0 + c1*x + c2*x^2 中确定系数使其在 x=0 附近最近似 cos(x)
  • 不考虑任何先验知识,看看只凭直觉应该怎么确定这三个系数
    • 首先,cos(x)x=0 处等于 1,那么至少多项式这一点要满足,即 c0=1
    • 如果多项式在 x=0 处切线的斜率也与 cos(x) 相同,那么近似程度应该会更高,于是 c1=0
    • 二阶导数描述的是曲线的弯曲程度,那么让两个图像在 x=0 处的弯曲程度相同应该会更近似,cos(x) 的二阶导是 -cos(x),于是 2c2=-cos(0)=-1 -> c2=-1/2
    • 还可以用更高次的多项式去近似,比如

    • 在非 0 点,有类似的结果,如 x=π

      带入 x=π 就能消去所有无关项

  • 泰勒多项式”小结
    • 在对高阶多项式求导时,自然而然的出现了阶乘的形式
    • 当然,真正的系数需要除以这个阶乘

    • 在向近似多项式中添加更高次的项时,不会影响低次项
    • 因此,多项式任意 n 阶的导数在 x=0 时的值都由唯一的系数控制

    • 这样的多项式就称为“泰勒多项式”

泰勒多项式的本质

  • 泰勒多项式实际上是利用函数在某点处的高阶导数,来近似该点附近的函数值。其原因是多项式比别的函数友好。好计算,好求导,好积分。用导数的方法去找近似:
  • 宏观来讲:泰勒级数把某一点处高阶导数的信息转化成了在那一点附近的函数值信息

  • x=0
    • 它的常数项能让它与 f(0) 的值相等
    • 它的一次项让两者的斜率相等
    • 二次项让两者斜率的变化率相同
  • 以此类推,项数越多,近似就越精确

几个要点:

  • 阶乘的形式是自然而然出现的。
  • 往近似多项式中添加更高次项时,低阶的项并不会因此而改变。多项式任意 n 阶的导数在 x=0 时的值都由唯一的一个系数控制。

  • 累加有限多项的式子叫做 “泰勒多项式”,累加无限多项就是 “泰勒级数”。
  • 收敛:一个级数加的越多,和越接近某个确定的数值,这个级数 “收敛” 到那个值。也可以说,累加无限多项的和,即这个级数就 “等于” 它收敛到的值。
  • 发散:能让多项式和收敛的最大取值范围叫做这个泰勒级数的收敛半径。

从几何角度看二阶泰勒多项式

  • 考虑如何近似曲线下的面积函数 f_area(x)

    • dx → 0 时,我们可以将增长的部分近似为矩形
    • 但如果想将面积的变化近似得更准确,就需要考虑那块近似三角形的部分
  • 假设已知面积函数 fa 点的导数信息,想要近似在 x 时的面积

    三角形的高 = 斜率 * 底,斜率 = 曲线在 a 点的导数 = f''(a)

任意函数的泰勒多项式

  • 函数f(x)x=0处的泰勒多项式

  • 函数f(x)x=a处的泰勒多项式

级数的概念

级数(series)的定义——无限项的和

  • 当泰勒多项式无限累加下去,就成了泰勒级数

级数的一些重要概念

  • 级数的收敛(Converges)与发散(diverges)、收敛半径

    无限累加下去的和就“等于”级数收敛到的值;这里 x=1

    • 一些级数无论你带入什么值,它的级数都会收敛,比如 e^xx=0 处的泰勒级数

      无论 x 取何值,e^x 都等于该泰勒级数

    • 但有些级数就只会在一定取值范围内才会收敛,比如 ln(x)x=1 处的泰勒级数

      • 只有当 x∈(0,2]时,该级数才会收敛
      • 换言之,在 x=1 处取得的导数信息无法拓展到更广的取值范围
      • 这个泰勒级数的收敛半径1,即x∈(0, 2]

常见的泰勒级数

泰勒级数_百度百科 https://baike.baidu.com/item/泰勒级数/7289427?fr=aladdin#6

泰勒级数的更多内容

  • 拉格朗日余项
  • 判断级数收敛——审敛法
  • 求收敛半径

本系列的视频目的是帮你建立关于微积分的基本直觉,而微积分的道路才刚刚开始…

资料

结束


支付宝打赏 微信打赏

~ 海内存知已,天涯若比邻 ~

Share

Related Posts

标题:联邦学习-Federated Learning

摘要:如何保护用户隐私的同时也能充分利用各类数据?

标题:论文阅读-Paper-Club

摘要:最新技术咨询梳理,sota(the state of the art), 论文资料汇总

Comments

--disqus--

    Content
    My Moment ( 微信公众号 )
    欢迎关注鹤啸九天