讲师中心 微信公众号
AI工具推荐 视频效率加速

如何在Python中用NumPy实现快速傅里叶变换FFT

胖雪君_3005

胖雪君_3005

发布时间:2026-09-05 06:22:18

|

248人浏览过

|

来源于php中文网

原创

numpy.fft.fft 默认对输入一维数组沿最后一轴(axis=-1)计算离散傅里叶变换(DFT),返回等长复数数组,零频分量位于索引0处,不归一化、不补零、不预处理。

如何在python中用numpy实现快速傅里叶变换fft

fft 函数默认做的是什么变换?

numpy.fft.fft 对输入的一维数组做离散傅里叶变换(DFT),返回复数数组,长度与输入相同。它不自动归一化,也不对输入做任何预处理(比如去均值、加窗),输出的零频分量(DC)在索引 0 处,正频率从 1 到 n//2,负频率从 n//2+1 开始(对偶索引)。如果你直接画出 abs(fft(x)),会看到镜像对称——这是实信号 DFT 的共轭对称性导致的,不是 bug。

  • 实信号的 FFT 结果满足 fft[k] == conj(fft[n-k])(模长对称,相位相反)
  • 如果你只关心幅度谱,通常取前半部分(含 DC 和 Nyquist 点)就够了:np.abs(fft_result[:len(fft_result)//2 + 1])
  • 不要误以为 fft 默认返回“物理频率轴”——它只返回系数,频率坐标得自己算:freqs = np.fft.fftfreq(n, d=dt)

为什么 fftshift 后图形才看起来“正常”?

numpy.fft.fft 输出的顺序是 [0, 1, 2, ..., n/2, -n/2+1, ..., -1](对偶排列),而人眼习惯把零频放在中间。所以直接绘图时,正负频率被劈开在两端,看着别扭。

  • 调用 np.fft.fftshift 可以把负频率块移到左边、正频率移到右边,零频居中
  • 它只是重排数组顺序,不改变数值本身,也不影响逆变换
  • 注意:必须对频率轴和 FFT 结果同时用 fftshift,否则横纵坐标对不上
    freqs = np.fft.fftfreq(n, dt)
    X = np.fft.fft(x)
    freqs_shifted = np.fft.fftshift(freqs)
    X_shifted = np.fft.fftshift(X)
  • 如果你后续要用 ifft,记得先 ifftshift 回去(ifftshift 和 fftshift 对偶长数组效果一样,但对奇数长度更安全)

real vs complex 输入对 fft 结果有什么影响?

NumPy 的 fft 自动适配输入类型:实数输入返回复数结果,复数输入也返回复数结果,但内存和计算路径略有差异。

提示词大师-python版
提示词大师-python版

图片提示词生成器?不止如此。 马甲系统 —— 把脑海中的画面,翻译成AI能理解的专业表达。 用得越多,它越懂你:首次需要多问几句确认方向,用久了几乎一说就懂。 用得越多,它越快:缓存机制让后续对话越来越省。 RAG进化:成功案例持续入库,越跑越聪明。 输入「新手指南」查看完整功能介绍

下载
  • 实数输入会触发底层优化(如用 RFFT),速度略快,且结果具备共轭对称性
  • 但如果你手动传入 np.complex64 或 np.complex128,它就走通用复数路径,不启用实数优化
  • 更关键的是:若原始信号本是实数,但你无意中加了极小虚部(比如 x + 1e-15j),就会失去对称性,幅度谱不再镜像,还可能引入数值噪声
  • 检查方式很简单:np.allclose(X[1:n//2+1], np.conj(X[-1:-n//2-1:-1])) —— 若为 False,说明输入被污染或用了非实数 dtype

采样率和点数没设对,fft 频率分辨率就全错了

FFT 本身不“知道”你的采样率;fftfreq 里的 dt(采样间隔)和 n(点数)共同决定频率分辨率 df = 1/(n*dt) 和最大可观测频率 fs/2。

立即学习“Python免费学习笔记(深入)”;

  • 常见错误:用 fs = 1000 Hz 采样,但传给 fftfreq 的 dt = 1/fs 写成 1000 或 0.001 漏掉单位换算
  • 点数 n 必须是实际参与 FFT 的长度,不是原始信号长度。如果做了截断或补零,n 就得用补零后的新长度
  • 补零(zero-padding)不会提高真实分辨率(那是由采样时长决定的),但会让频谱看起来更“平滑”,便于观察峰值位置
  • 真实分辨率只取决于总采样时间 T = n * dt,即 df = 1/T。想提高分辨率,只能延长采集时间,不能靠多补零

FFT 的核心是快速算法,但真正用准它,得盯住采样参数、数据类型、频谱排列这三处。最容易漏的是 dt 单位错和忘记 fftshift 配对,一画图就发现频谱“歪了”。

热门AI工具

更多
WorkBuddy

一款AI办公效率工具,主要用于腾讯云推出的AI原生桌面智能体工作台,适合需要提升相关任务效率的用户。

UpDream
UpDream Hot

一款AI视频创作工具,主要用于哔哩哔哩推出的自研AI视频创作工具,适合需要提升相关任务效率的用户。

二狗PPT
二狗PPT Hot

一款AI演示文稿工具,主要用于专为中式职场打造的AI PPT生成工具,适合需要提升相关任务效率的用户。

DeepSeek

DeepSeek是一款面向对话、写作、编程和推理场景的AI大模型工具。

超级简历WonderCV

一款AI办公效率工具,主要用于免费求职简历模版下载制作,应届生职场人必备简历制作神器,适合需要提升相关任务效率的用户。

火山引擎

火山引擎是一款面向企业的云计算与AI服务平台。

豆包大模型

豆包大模型是一款由字节跳动推出的企业级大语言模型服务平台。

VibeKnow
VibeKnow Hot

一款AI视频创作工具,主要用于全球首个AI知识视频创作平台,文档、文章、网页,一键生成视频,适合需要提升相关任务效率的用户。

讯飞绘文

讯飞绘文是一款由科大讯飞推出的一站式 AIGC 内容运营平台。

相关专题

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

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

1611

2023.07.20

python能做什么
python能做什么

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

3924

2023.07.25

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

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

1629

2023.07.31

python教程
python教程

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

22577

2023.08.03

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

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

2767

2023.08.04

python eval
python eval

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

2807

2023.08.04

scratch和python区别
scratch和python区别

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

1123

2023.08.11

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

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

596

2023.08.10

LLVM自定义Pass怎么写
LLVM自定义Pass怎么写

本专题聚焦LLVM自定义Pass开发,整理Pass类结构、run()方法、PreservedAnalyses、CMake构建、插件注册、-load-pass-plugin加载和测试用例编写流程。

0

2026.09.30

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
关于我们 免责申明 举报中心 意见反馈 讲师合作 广告合作 最新更新
php中文网:公益在线php培训,帮助PHP学习者快速成长!
关注服务号
PHP中文网订阅号
每天精选资源文章推送

Copyright 2014-2026 https://www.php.cn/ All Rights Reserved | php.cn