码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 泰勒公式理解



    文章目录

        • 1:一元泰勒展开公式
        • 2:二元泰勒展开公式
        • 3:二元函数的黑塞矩阵
        • 4:多元函数的黑塞矩阵
        • 5:多元函数的雅可比矩阵(Jacobian矩阵)
          • 参考文献


    个人笔记:

    1:一元泰勒展开公式

    在这里插入图片描述举例:f(x) = 3x² + 2x + 5 在x=0或x=1处的泰勒展开
    在这里插入图片描述
    在这里插入图片描述当x=0时:
    在这里插入图片描述
    当x=1时:

    在这里插入图片描述

    不论Xk等于多少,最后展开得公式相加都是等于f(x) = 3x² + 2x + 5

    2:二元泰勒展开公式

    x 和 y在k处的泰勒展开:
    在这里插入图片描述
    简化:
    在这里插入图片描述

    简化:
    在这里插入图片描述①
    f x x ′ ′ f''_{xx} fxx′′​是对 x 求两次导。

    ②
    f x y ′ ′ f''_{xy} fxy′′​是先对x求一次导,然后再对y求一次导。

    ③
    f y x ′ ′ f''_{yx} fyx′′​是先对y求一次导,然后再对x求一次导。
    (其中③ = ②)

    ④
    f y y ′ ′ f''_{yy} fyy′′​是对 y 求两次导。

    3:二元函数的黑塞矩阵

    二元函数点 f ( x 1 , x 2 ) f(x_1,x_2) f(x1​,x2​) 在 X ( k ) ( x 1 ( k ) , x 2 ( k ) ) X^{(k)}(x_1^{(k)},x_2^{(k)}) X(k)(x1(k)​,x2(k)​)处的泰勒展开式为:
    在这里插入图片描述

    其中 Δ x 1 Δ x_1 Δx1​ = x 1 x_1 x1​ − x 1 ( k ) x_1^{(k)} x1(k)​ , Δ x 2 Δ x_2 Δx2​ = x 2 x_2 x2​ − x 2 ( k ) x_2^{(k)} x2(k)​

    在这里插入图片描述

    即:
    在这里插入图片描述

    (1):其中
    在这里插入图片描述

    它是 f ( X ) f(X) f(X)在 X ( k ) X^{(k)} X(k)点处的梯度。

    (2): G ( X ( k ) ) G(X^{(k)}) G(X(k))是 f ( x 1 , x 2 ) f(x_1,x_2) f(x1​,x2​)在 X ( k ) X^{(k)} X(k)处的黑塞矩阵。它是由函数 f ( x 1 , x 2 ) f(x_1,x_2) f(x1​,x2​)在 X ( k ) X^{(k)} X(k)处的二阶偏导数所组成的方阵。
    在这里插入图片描述

    4:多元函数的黑塞矩阵

    1:多元函数 f ( x 1 , x 2 , . . . , x n ) f(x_1,x_2,...,x_n) f(x1​,x2​,...,xn​)在点 x ( k ) x^{(k)} x(k)处的泰勒展开式为:
    在这里插入图片描述把泰勒(Taylor)展开式写成矩阵的形式:
    在这里插入图片描述其中:
    在这里插入图片描述
    它是 f ( X ) f(X) f(X)在 X ( k ) X^{(k)} X(k)点处的梯度。
    (2): G ( X ( k ) ) G(X^{(k)}) G(X(k))是 f ( x 1 , x 2 , . . . , x n ) f(x_1,x_2,...,x_n) f(x1​,x2​,...,xn​)在 X ( k ) X^{(k)} X(k)处的黑塞矩阵。它是由函数 f ( x 1 , x 2 , . . . , x n ) f(x_1,x_2,...,x_n) f(x1​,x2​,...,xn​)在 X ( k ) X^{(k)} X(k)处的二阶偏导数所组成 n ∗ n n*n n∗n阶方阵。

    在这里插入图片描述2:
    在这里插入图片描述
    在这里插入图片描述举例:
    在这里插入图片描述

    5:多元函数的雅可比矩阵(Jacobian矩阵)

    1.概述
    设 f f f: R n R^n Rn → R m R^m Rm是一个函数,它的输入是向量 x ∈ R n x ∈ R^n x∈Rn,输出是向量 y = f ( x ) ∈ R m y = f(x)∈ R^m y=f(x)∈Rm,并且 m ≥ n m≥n m≥n是一个从欧式 n n n维空间转换到欧式 m m m维空间的函数,这个函数由 m m m个实函数组成: f 1 ( x 1 , … , x n ) f_1(x_1,…,x_n) f1​(x1​,…,xn​),…, f m ( x 1 , … , x n ) f_m(x_1,…,x_n) fm​(x1​,…,xn​)这些函数的偏导数(如果存在)可以组成一个 m ∗ n m∗n m∗n的矩阵, 这就是所谓的Jacobian矩阵:

    在这里插入图片描述
    那么雅可比矩阵是一个 m × n m×n m×n 矩阵,通常被定义为
    在这里插入图片描述

    2:举例
    来看一个实际的数据拟合过程,输入:
    自变量: x = x = x={ 1 , 2 , 4 , 5 , 8 1,2,4,5,8 1,2,4,5,8 }
    因变量: y = y = y={ 3.2939 , 4.2699 , 7.1749 , 9.3008 , 20.259 3.2939,4.2699,7.1749,9.3008,20.259 3.2939,4.2699,7.1749,9.3008,20.259 }

    目标:用函数 f = p 1 ∗ e p 2 ∗ x − y f=p_1 ∗ e^{p_2∗x}−y f=p1​∗ep2​∗x−y 进行拟合,这里自变量 x x x,因变量 y y y,参数 p 1 p_1 p1​和 p 2 p_2 p2​,对参数 p 1 p_1 p1​和 p 2 p_2 p2​进行求导:
    在这里插入图片描述
    e p ∗ x e^{p*x} ep∗x对 p p p进行求导得: x ∗ e p ∗ x x*e^{p*x} x∗ep∗x

    雅可比矩阵描述:

    Jacobian矩阵 =
    
    [   exp(p2),     p1*exp(p2)]
    [ exp(2*p2), 2*p1*exp(2*p2)]
    [ exp(4*p2), 4*p1*exp(4*p2)]
    [ exp(5*p2), 5*p1*exp(5*p2)]
    [ exp(8*p2), 8*p1*exp(8*p2)]
    
    • 1
    • 2
    • 3
    • 4
    • 5
    • 6
    • 7

    即:
    在这里插入图片描述

    参考文献

    黑森矩阵

    黑塞矩阵和雅克比矩阵

    雅克比矩阵1

    雅克比矩阵2
    雅可比矩阵直观图像理解

    在这里插入图片描述

  • 相关阅读:
    使用linux系统IO函数实现 `cp`指令的功能
    使用flink sqlserver cdc 同步数据到StarRocks
    Spark简单介绍,Windows下安装Scala+Hadoop+Spark运行环境,集成到IDEA中
    day58| 739. 每日温度、496.下一个更大元素 I
    ESP8266-Arduino编程实例-OLED显示DHT22传感器数据
    【Y 新闻】YMatrix 成立三周年,三岁的我们还真是“不简单”
    java计算机毕业设计智慧农业水果销售系统MyBatis+系统+LW文档+源码+调试部署
    haproxy使用
    案例赏析 | 土耳其开赛利:闲置屋顶坐享“阳光收益”,助力企业实现绿色低碳转型
    c++实现最大堆
  • 原文地址:https://blog.csdn.net/weixin_43763292/article/details/127930050
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号