Python怎么计算NumPy数组的傅里叶变换_使用fft函数频域分析

冬瑶吖_1495

冬瑶吖_1495

2026-03-27

478人浏览

原创

NumPy 的 fft 函数默认对数组的最后一维(axis=-1)做变换;fftfreq 需传入采样间隔 d 才得真实频率;fft 输出等长复数数组,含完整频谱;fftshift 仅重排索引以居中零频。

python怎么计算numpy数组的傅里叶变换_使用fft函数频域分析

fft 函数默认对哪一维做变换?

NumPy 的 fft 系列函数(比如 fft.fft)默认只对数组的**最后一维**(即 axis=-1)做一维傅里叶变换。这不是“自动识别信号维度”,而是硬编码规则——哪怕你传入的是二维图像数据,它也不会默认沿时间轴或空间轴处理,除非你显式指定 axis。

常见错误现象:np.fft.fft(img_2d) 返回形状和输入一致的复数数组,但变换其实是逐行进行的(因为行是最后一维),不是你想做的“整张图二维频谱”。

  • 若处理时间序列(shape 为 (N,)),不用管 axis,默认就对它变
  • 若处理多通道信号(shape 为 (N, C)),想对每条通道独立变,得写 np.fft.fft(x, axis=0)(让变换沿第 0 维发生,即每列一个 FFT)
  • 二维图像想做真正的二维 FFT,必须用 np.fft.fft2,而不是套 fft 两次

fft 和 fftfreq 配合时,采样率怎么传才不翻车?

np.fft.fftfreq 生成的频率数组,默认单位是“周期/采样点”,也就是归一化频率。要得到真实物理频率(Hz),必须传入实际采样间隔 d(单位:秒),即 fs = 1/d。漏传或传错 d,画出来的横轴就全偏了。

使用场景:你有一段 1 秒长、采样率 1000 Hz 的音频,共 1000 个点。此时 d = 0.001,调用应为 np.fft.fftfreq(1000, d=0.001),结果才是 0, 1, 2, ..., 499, -500, ..., -1(Hz)。

  • 如果误写成 fftfreq(1000, d=1),得到的是 0, 0.001, 0.002, ... —— 单位变成 kHz,数值小了 1000 倍
  • 如果数据是等时间间隔但起始时间不为 0(比如从 t=2.5s 开始),d 仍只取决于间隔,和起始无关
  • 对非均匀采样数据,fft 不适用;得先插值或换用 Lombscargle 等方法

为什么 fft 结果是复数,且长度和输入一样?

np.fft.fft 输出复数数组,实部是余弦分量幅值,虚部是正弦分量幅值。长度和输入相同,是因为它计算的是离散傅里叶变换(DFT)的标准定义,不是“只返回正频率”——后半段负频率是共轭对称的,但 NumPy 保留全部,方便逆变换精确还原。

macOS UI 自动化技能包 —— 用 Python 控制你的 Mac
macOS UI 自动化技能包 —— 用 Python 控制你的 Mac

让 macOS 自动化像说话一样简单 —— 用代码控制你的 Mac,无需学习 AppleScript 系统要求 macOS 14+ (Sonoma) Python 3.9+ 需要开启:辅助功能 + 屏幕录制权限

下载

性能影响:保留全部点意味着内存和计算量都是 O(N),没做任何截断或压缩。如果你只关心幅度谱,记得用 np.abs(fft_result);如果只想要单边谱(如画图),要手动取前半(含直流)+ 合并负频共轭(或直接用 np.fft.rfft)。

  • rfft 对实数输入更省:输出长度是 N//2 + 1,只含非负频率,且是实数数组(最后一个是 Nyquist 频率,可能为实)
  • 用 ifft 还原时,必须用和 fft 完全同构的复数数组,不能只喂一半过去
  • 注意浮点误差:即使原始数组是实数,ifft(fft(x)) 的虚部通常不严格为 0,要用 .real 或 np.real_if_close

fftshift 怎么用才不把频谱中心搞错?

np.fft.fftshift 的作用是把零频分量移到频谱中心——但它**不改变数值,只重排顺序**。很多人以为它“做了某种变换”,其实只是索引平移。错用位置(比如在 fft 前调用,或对错误维度调用)会导致频谱完全错位。

典型场景:画二维图像的频谱图,希望低频在中间、高频在四周。必须按顺序来:f2d = fft2(img); f2d_shifted = fftshift(f2d),再用 abs 或 log 显示。

  • 对一维结果,fftshift 把 [0,1,...,N-1] 变成 [N//2, ..., N-1, 0, ..., N//2-1](偶数 N)
  • 对多维数组,默认对所有轴 shift;如只想 shift 第一维,写 fftshift(arr, axes=0)
  • 如果之后还要做 ifft2,别忘了先 ifftshift 回去——fftshift 和 ifftshift 对偶数长度互为逆操作,但对奇数长度也安全

真正容易被忽略的是:FFT 结果的相位信息极其敏感,哪怕输入数组末尾多一个零、少一个零,整个相位谱都可能跳变。做对比实验时,务必检查数组长度、填充方式、是否去均值——这些细节比选哪个函数影响更大。

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

相关专题

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

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

2023.07.20

1831

4

python能做什么
python能做什么

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

2023.07.25

4884

7

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

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

2023.07.31

1789

3

python教程
python教程

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

2023.08.03

29197

23

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

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

2023.08.04

3487

5

python eval
python eval

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

2023.08.04

3507

5

scratch和python区别
scratch和python区别

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

2023.08.11

1223

5

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

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

2023.08.10

616

4

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

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

2023.08.11

2723

5

热门下载

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

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习