Python中如何使用TensorFlow计算梯度和雅可比矩阵

轻静吖_3686

轻静吖_3686

2026-09-30

371人浏览

原创

tf.gradienttape用于计算标量函数对张量的梯度,要求输入为tf.variable或显式watch的张量,目标必须是标量,否则报错;支持嵌套求高阶导、vectorized_map批量算雅可比,需注意形状匹配与persistent内存管理。

python中如何使用tensorflow计算梯度和雅可比矩阵

用 tf.GradientTape 计算标量函数对张量的梯度

TensorFlow 2.x 默认启用 Eager Execution,必须用 tf.GradientTape 显式记录运算才能求导。不包在 tape 里,tf.gradients 或 tape.gradient 会返回 None。

常见错误是忘了把输入变量设为可训练(trainable=True)或没用 tf.Variable 包裹——纯 tf.constant 或 NumPy 数组不会被 tape 追踪。

  • 输入必须是 tf.Variable 或设置了 watch 的张量:tape.watch(x)
  • 被求导目标必须是标量(shape=()),否则 tape.gradient(y, x) 报错 ValueError: Cannot take gradient of non-scalar tensor
  • 示例:计算 y = x**2 + 2*x 在 x=3 处的导数
import tensorflow as tf
x = tf.Variable(3.0)
with tf.GradientTape() as tape:
    y = x**2 + 2*x
dy_dx = tape.gradient(y, x)  # 返回 <tf.tensor: shape="()," dtype="float32," numpy="8.0"></tf.tensor:>

用 tf.vectorized_map + tf.GradientTape 批量算雅可比矩阵

雅可比是向量函数对向量输入的一阶偏导矩阵,TensorFlow 没有开箱即用的 tf.jacobian(v2.11+ 实验性支持但不稳定),主流做法是循环或向量化调用 gradient。

直接 for 循环效率低;用 tf.vectorized_map 可并行计算每行雅可比,前提是输出函数能写成 batch-aware 形式。

  • 若 f: R^n → R^m,先确保 f 接受 shape (batch, n) 输入,返回 shape (batch, m)
  • 对每个输出分量 f_i 单独 tape 求梯度,再 stack 成矩阵
  • 更稳的做法:用 tf.stack([tape.gradient(f[i], x) for i in range(m)], axis=0),但需保证 f[i] 是标量
x = tf.Variable([[1.0, 2.0]])  # shape (1, 2)
with tf.GradientTape(persistent=True) as tape:
    f0 = x[0, 0]**2 + x[0, 1]
    f1 = x[0, 0] * x[0, 1]
# 分别求两个输出对标量输入的梯度,拼成 2x2 雅可比
jacobian_row0 = tape.gradient(f0, x)  # [2., 1.]
jacobian_row1 = tape.gradient(f1, x)  # [2., 1.]
J = tf.stack([jacobian_row0, jacobian_row1], axis=0)  # shape (2, 2)

处理高维输入/输出时的形状陷阱

雅可比维度容易搞反:对 f: R^N → R^M,正确雅可比是 M × N 矩阵,但 tape.gradient(f_i, x) 返回的梯度 shape 和 x 一致,不是自动展平的。

Python Use Agent
Python Use Agent

智能执行Python任务,自动生成、执行代码并反馈结果,无需额外配置,兼容旧命令。

下载

比如 x 是 (3, 4),tape.gradient(f0, x) 也返回 (3, 4),想得到 1 × 12 行向量就得手动 tf.reshape;否则 stack 后维度混乱。

  • 始终检查 x.shape 和 grad.shape 是否一致,不一致说明没正确 watch 或输入类型不对
  • 批量输入时,避免用 tf.expand_dims 错误增加维度导致 tape 无法匹配
  • 如果输出是 rank > 1(如图像 loss),先 tf.reduce_sum 或 tf.reduce_mean 降为标量再求梯度,否则不是雅可比而是高阶导数

替代方案:用 tf.function 加速但注意 tape 生命周期

tf.function 能加速梯度计算,但 tape 必须在 function 内部创建——把 tape 提到外面再传入会失效,因为 graph mode 下 tape 不记录。

另外,persistent=True 的 tape 在 tf.function 中要小心内存,不用完记得 del tape,否则可能 OOM。

  • 不要写 @tf.function 包裹外层函数再调用已定义的 tape —— tape 必须在被装饰函数内部实例化
  • 调试时先关掉 @tf.function,确认梯度逻辑正确后再加,否则报错信息难读
  • 涉及控制流(if/while)时,tape 仍有效,但所有分支都必须产生可微路径,否则某分支返回常量会导致梯度截断

实际用雅可比的地方不多,多数场景只需梯度;真需要完整雅可比时,得反复确认输入展平方式、输出索引顺序和最终 reshape 是否匹配数学定义——这里最容易漏掉转置或维度错位。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1611

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

3924

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1629

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

22597

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2767

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

2827

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1123

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

596

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2183

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程