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

如何在Python中通过NumPy实现快速傅里叶逆变换还原信号

阿静姑娘_9343

阿静姑娘_9343

发布时间:2026-11-02 11:02:03

|

623人浏览过

|

来源于php中文网

原创

np.fft.ifft返回复数是因浮点误差导致虚部微小非零,需用.real安全取实部;频谱必须满足k=0和k=N/2为实数、其余点共轭对称,否则虚部不可忽略。

如何在python中通过numpy实现快速傅里叶逆变换还原信号

为什么 np.fft.ifft 返回复数,但原始信号是实数?

直接调用 np.fft.ifft 后得到的数组类型通常是 complex128,即使输入是纯实信号的频谱(比如从 np.fft.fft 得到的),还原结果的虚部也常有微小非零值(如 1e-16j)。这不是 bug,而是浮点计算累积误差所致。实际使用中必须显式取实部,但要注意:仅当确认原始信号为实数、且频谱满足共轭对称性时,才能安全用 np.real() 或 .real。

  • 错误做法:signal = np.fft.ifft(freq_data) 后直接用于后续计算(可能触发复数参与运算)
  • 推荐做法:signal = np.real(np.fft.ifft(freq_data)) 或更稳妥地 signal = np.fft.ifft(freq_data).real
  • 若频谱不满足共轭对称(例如手动构造的非对称频谱),虚部不可忽略,强行取实会丢失信息

如何避免因长度不匹配导致的时域信号错位?

np.fft.ifft 的输出长度严格等于输入频域数组长度。若你只保留了部分频谱(比如截断高频分量做滤波),再直接 ifft,还原信号长度会变短,且相位关系被破坏——这会导致时域波形整体平移甚至失真,而非“干净去噪”。

  • 关键原则:滤波后必须保持频谱长度与原始 fft 输出一致,空缺位置补零(不是丢弃)
  • 正确流程:freq_full = np.fft.fft(original_signal) → 修改 freq_full[10:-9] = 0(保留低频)→ recon = np.real(np.fft.ifft(freq_full))
  • 错误示例:low_freq = freq_full[:10] → np.fft.ifft(low_freq)(长度突变,信号崩坏)

用 np.fft.irfft 替代 ifft 能省什么?

当你确定原始信号是实数,且手头只有「一半」频谱(即 np.fft.rfft 输出),就该用 np.fft.irfft。它专为实信号设计,自动补全共轭对称的另一半频谱,省去手动拼接,还默认返回实数数组,无需 .real。

python-flask-vue-creator
python-flask-vue-creator

创建 Flask + Vue.js 全栈项目,采用单端口架构、uv 依赖管理、结构化日志和生产级前端设计。当用户要求以下操作时触发:创建新的 Flask+Vue 项目、搭建 Web 应用、构建全栈 Python+Vue 应用、初始化 Flask Web 项目、设置单端口 Flask 前后端项目。触发词:新建项目、创建项目、搭建项目、Flask项目、Vue项目。

下载
  • 适用场景:从传感器读取的实测时间序列,经 rfft 处理后需还原
  • 参数注意:irfft 默认按 n = len(rfft_output) * 2 - 2 推断原长(偶数长度);若原始长度为奇数,必须显式传入 n= 参数
  • 性能差异:比 ifft 略快,且内存更省(不存冗余共轭部分)

为什么还原信号和原始信号在幅度上差了一个缩放因子?

np.fft.ifft 默认执行的是**未归一化**逆变换,其数学定义包含 1/N 因子,但 NumPy 把这个因子放在了 ifft 而非 fft 上——也就是说,ifft(fft(x)) 结果是 x,但中间每一步都不带归一化。可一旦你手动修改频谱(比如放大某段),再 ifft,幅度就不再匹配原始信号。

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

  • 验证方法:np.allclose(np.fft.ifft(np.fft.fft(x)), x) 应为 True(浮点误差内)
  • 若需与物理单位对齐(如电压/伏特),通常要乘以 len(x):还原后做 recon * len(x)
  • 更稳妥的做法:统一用 norm="ortho" 参数(NumPy 1.19+),即 np.fft.ifft(..., norm="ortho") 和 np.fft.fft(..., norm="ortho") 构成严格正交变换,双向缩放一致

真正容易被忽略的是频谱的对称性维护和归一化意图的显式表达——这两点不出错,ifft 就只是个可靠的工具;一旦模糊处理,还原出来的信号连基本周期都对不上。

热门AI工具

更多
创客贴AI

创客贴AI是一款AI图像与设计工具,AI辅助的智能在线设计工具。

WorkBuddy

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

Atoms
Atoms Hot

Atoms是一款AI智能体工具,第一支自动构建真实业务的 AI 团队。

豆包大模型

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

音述AI
音述AI Hot

一款AI音频处理工具,主要用于音述AI是一个以“用声音述说故事”为核心的 AI 音乐创作与声音分享社区,适合需要提升相关任务效率的用户。

火山引擎

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

Seko
Seko Hot

一款AI视频创作工具,主要用于商汤科技推出的创编一体的AI短视频创作Agent,适合需要提升相关任务效率的用户。

UP简历
UP简历 Hot

一款AI办公效率工具,主要用于基于AI技术的免费在线简历制作工具,适合需要提升相关任务效率的用户。

DeepSeek

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

相关专题

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

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

1851

2023.07.20

python能做什么
python能做什么

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

5004

2023.07.25

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

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

1789

2023.07.31

python教程
python教程

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

30097

2023.08.03

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

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

3587

2023.08.04

python eval
python eval

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

3587

2023.08.04

scratch和python区别
scratch和python区别

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

1243

2023.08.11

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

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

616

2023.08.10

Uvicorn 报错与问题解决教程
Uvicorn 报错与问题解决教程

本专题汇总了Uvicorn服务器常见问题与报错解决的教程合集,涵盖启动失败、端口占用、模块导入错误、事件循环冲突及热重载失效等典型故障的排查思路。同时整理了多进程与Worker配置、SSL/HTTPS证书报错、跨域与代理配置等实战解决方案,助你快速定位并解决Uvicorn运行中的各类问题。

0

2026.11.02

热门下载

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

精品课程

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

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