[信号与系统]傅里叶变换、卷积定理、和为什么时域的卷积等于频域相乘。

本文主要是介绍[信号与系统]傅里叶变换、卷积定理、和为什么时域的卷积等于频域相乘。,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

前言

最近学习以下IIR滤波器和FIR滤波器

前置

1. 时域和频域

时域和频域代表着频率和时间与振幅的一一对应关系
在这里插入图片描述

在这里插入图片描述

2. 卷积运算

关于卷积的定义,详情请看 这篇文章能让你明白卷积

卷积运算是一种数学运算,广泛应用于信号处理、图像处理、控制系统和概率论等领域。卷积运算可以看作是两个函数之间的一种积分操作,用于描述一个函数在另一个函数上的“滑动”效果。

连续卷积:

对于连续函数 f ( t ) f(t) f(t) g ( t ) g(t) g(t)它们的卷积定义为:

( f ∗ g ) ( t ) = ∫ − ∞ ∞ f ( τ ) g ( t − τ ) d τ (f * g)(t) = \int_{-\infty}^{\infty} f(\tau) g(t - \tau) \, d\tau (fg)(t)=f(τ)g(tτ)dτ

离散卷积

对于离散函数 f [ n ] f[n] f[n] g [ n ] g[n] g[n],他们的卷积定义为:

( f ∗ g ) [ n ] = ∑ k = − ∞ ∞ f [ k ] g [ n − k ] (f * g)[n] = \sum_{k=-\infty}^{\infty} f[k] g[n - k] (fg)[n]=k=f[k]g[nk]

卷积运算的性质

1. 交换律
f ∗ g = g ∗ f f * g = g * f fg=gf
这意味着两个函数的卷积不受顺序影响。

2. 结合律(结合性)
( f ∗ g ) ∗ h = f ∗ ( g ∗ h ) (f * g) * h = f * (g * h) (fg)h=f(gh)
这意味着多重卷积的计算顺序可以随意改变。

3.分配律
f ∗ ( g + h ) = ( f ∗ g ) + ( f ∗ h ) f * (g + h) = (f * g) + (f * h) f(g+h)=(fg)+(fh)

这意味着卷积运算对加法是分配的。

4. 与冲激函数的卷积(单位冲激函数)

对于单位冲激函数 δ ( t ) \delta(t) δ(t)有:

f ∗ δ = f f * \delta = f fδ=f
这意味着任何函数与单位冲激函数的卷积等于该函数本身。

5.平移性

f ( t − t 0 ) ∗ g ( t ) = ( f ∗ g ) ( t − t 0 ) f(t - t_0) * g(t) = (f * g)(t - t_0) f(tt0)g(t)=(fg)(tt0)
这意味着函数的平移在卷积后仍然保留。

卷积定理

在傅里叶变换域中,卷积运算可以转化为点乘运算。具体来说,如果 F ( ω ) F(\omega) F(ω) G ( ω ) G(\omega) G(ω)分别是 f ( t ) f(t) f(t) g ( t ) g(t) g(t) 的傅里叶变换,那么:

F { f ∗ g } = F { f } ⋅ F { g } \mathcal{F}\{f * g\} = \mathcal{F}\{f\} \cdot \mathcal{F}\{g\} F{fg}=F{f}F{g}

反之亦然,即傅里叶变换的点乘可以通过逆傅里叶变换转化为时域的卷积运算。

3. 傅里叶变换

傅里叶变换将一个时域信号转换到频域,使得可以分析信号的频率成分。对于非周期信号,傅里叶变换定义为:

F ( w ) = ∫ − ∞ ∞ f ( t ) e − j π t d t F(w) = \int_{-\infty}^{\infty} f(t) e^{-j\pi t} \, dt F(w)=f(t)etdt

其中 F ( w ) F(w) F(w)是频域表示,称为频谱。 f ( t ) f(t) f(t)是时域信号,w是角频率

逆傅里叶变换可以将频域信号转换回时域:

f ( t ) = 1 2 π ∫ − ∞ ∞ F ( w ) e j π t d w f(t) = \frac{1}{2\pi}\int_{-\infty}^{\infty} F(w) e^{j\pi t} \, dw f(t)=2π1F(w)etdw

关键点包括:

  1. 频域与时域的对应关系:时域信号可以通过傅里叶变换转换到频域,反之亦然。这提供了分析和处理信号的新方法。

  2. 频谱:傅里叶变换的结果 F ( w ) F(w) F(w)称为信号的频谱,表示信号在不同频率成分上的分布。

  3. 正交性 : 正弦和余弦函数是正交的,这使得傅里叶级数能够分解任何周期信号,而傅里叶变换能够分解任何非周期信号。

  4. 卷积定理 : 时域中的卷积对应于频域中的乘积,这大大简化了信号处理中的卷积运算。

我们常见的其实会由这个卷积定理延申出来一个定理:

为了接下来的这个问题,我们需要用到一些傅里叶变换的性质:

1. 线性性质

傅里叶变换是线性的。对于任意两个信号 x 1 ( t ) x_1(t) x1(t) x 2 ( t ) x_2(t) x2(t),以及任意常数 a a a b b b,有:

F { a x 1 ( t ) + b x 2 ( t ) } = a X 1 ( ω ) + b X 2 ( ω ) \mathcal{F}\{a x_1(t) + b x_2(t)\} = a X_1(\omega) + b X_2(\omega) F{ax1(t)+bx2(t)}=aX1(ω)+bX2(ω)

2. 平移性质

如果一个信号 x ( t ) x(t) x(t) 在时域上平移 t 0 t_0 t0,其傅里叶变换为:

F { x ( t − t 0 ) } = X ( ω ) e − i ω t 0 \mathcal{F}\{x(t - t_0)\} = X(\omega) e^{-i\omega t_0} F{x(tt0)}=X(ω)et0

3. 调制性质

如果一个信号 x ( t ) x(t) x(t) 在时域上乘以一个复指数函数 e i ω 0 t e^{i\omega_0 t} eiω0t,其傅里叶变换为:

F { x ( t ) e i ω 0 t } = X ( ω − ω 0 ) \mathcal{F}\{x(t) e^{i\omega_0 t}\} = X(\omega - \omega_0) F{x(t)eiω0t}=X(ωω0)

4. 微分性质

如果一个信号 x ( t ) x(t) x(t) 的导数 d n x ( t ) d t n \frac{d^n x(t)}{dt^n} dtndnx(t) 存在,其傅里叶变换为:

F { d n x ( t ) d t n } = ( i ω ) n X ( ω ) \mathcal{F}\left\{\frac{d^n x(t)}{dt^n}\right\} = (i\omega)^n X(\omega) F{dtndnx(t)}=()nX(ω)

5. 卷积性质
如果两个信号 x 1 ( t ) x_1(t) x1(t) x 2 ( t ) x_2(t) x2(t) 的卷积 ( x 1 ∗ x 2 ) ( t ) (x_1 * x_2)(t) (x1x2)(t) 存在,其傅里叶变换为:

F { ( x 1 ∗ x 2 ) ( t ) } = X 1 ( ω ) ⋅ X 2 ( ω ) \mathcal{F}\{(x_1 * x_2)(t)\} = X_1(\omega) \cdot X_2(\omega) F{(x1x2)(t)}=X1(ω)X2(ω)

时域的卷积等于频域相乘

这是怎么来的呢,我们首先知道两个时域信号 x ( t ) x(t) x(t) h ( t ) h(t) h(t)的卷积:

( x ∗ h ) ( t ) = ∫ − ∞ ∞ x ( τ ) h ( t − τ ) d τ (x * h)(t) = \int_{-\infty}^{\infty} x(\tau) h(t - \tau) \, d\tau (xh)(t)=x(τ)h(tτ)dτ

我们得到X(f)和H(f)是他们分别通过傅里叶变换得到的频域表示,那么时域中的卷积 x ( t ) ∗ h ( t ) x(t) * h(t) x(t)h(t)就对应于频域中的乘积 X ( f ) ⋅ H ( f ) X(f)·H(f) X(f)H(f)
也就是说,时域中的复杂操作,我们最后用频域中的简单操作就可以指代了。

换句话说:

时域信号可以分解成一串不同频率正弦信号的叠加。根据卷积的分配率,两个时域信号的卷积最终可以展开成两两正弦信号的卷积的和。由于不同频率的正弦信号的卷积为0,所以最终只剩下相同频率的正弦信号的卷积。而卷积的结果就是频率不变,幅度相乘。
在频域里边就表现为直接相乘。

关于这句话,我们还可以生硬的推导一下:

推导

卷积定理的推导

我们首先知道两个时域信号 x ( t ) x(t) x(t) h ( t ) h(t) h(t) 的卷积定义为:

( x ∗ h ) ( t ) = ∫ − ∞ ∞ x ( τ ) h ( t − τ ) d τ (x * h)(t) = \int_{-\infty}^{\infty} x(\tau) h(t - \tau) \, d\tau (xh)(t)=x(τ)h(tτ)dτ

y ( t ) = ( x ∗ h ) ( t ) y(t) = (x * h)(t) y(t)=(xh)(t),则 y ( t ) y(t) y(t) x ( t ) x(t) x(t) h ( t ) h(t) h(t) 的卷积。

傅里叶变换

y ( t ) y(t) y(t) 进行傅里叶变换:

Y ( f ) = F { y ( t ) } = F { ∫ − ∞ ∞ x ( τ ) h ( t − τ ) d τ } Y(f) = \mathcal{F}\{y(t)\} = \mathcal{F}\left\{\int_{-\infty}^{\infty} x(\tau) h(t - \tau) \, d\tau\right\} Y(f)=F{y(t)}=F{x(τ)h(tτ)dτ}

根据傅里叶变换的线性性质,我们可以将积分符号放到傅里叶变换操作符的前面:

Y ( f ) = ∫ − ∞ ∞ x ( τ ) F { h ( t − τ ) } d τ Y(f) = \int_{-\infty}^{\infty} x(\tau) \mathcal{F}\{h(t - \tau)\} \, d\tau Y(f)=x(τ)F{h(tτ)}dτ

时移性质

根据傅里叶变换的时移性质,如果 h ( t − τ ) h(t - \tau) h(tτ) 的傅里叶变换为 H ( f ) H(f) H(f),则:

F { h ( t − τ ) } = H ( f ) e − i 2 π f τ \mathcal{F}\{h(t - \tau)\} = H(f) e^{-i2\pi f \tau} F{h(tτ)}=H(f)ei2πfτ

将其代入上式:

Y ( f ) = ∫ − ∞ ∞ x ( τ ) H ( f ) e − i 2 π f τ d τ Y(f) = \int_{-\infty}^{\infty} x(\tau) H(f) e^{-i2\pi f \tau} \, d\tau Y(f)=x(τ)H(f)ei2πfτdτ

分离变量

H ( f ) H(f) H(f) 从积分符号中分离出来:

Y ( f ) = H ( f ) ∫ − ∞ ∞ x ( τ ) e − i 2 π f τ d τ Y(f) = H(f) \int_{-\infty}^{\infty} x(\tau) e^{-i2\pi f \tau} \, d\tau Y(f)=H(f)x(τ)ei2πfτdτ

这里, ∫ − ∞ ∞ x ( τ ) e − i 2 π f τ d τ \int_{-\infty}^{\infty} x(\tau) e^{-i2\pi f \tau} \, d\tau x(τ)ei2πfτdτ x ( t ) x(t) x(t) 的傅里叶变换:

X ( f ) = ∫ − ∞ ∞ x ( τ ) e − i 2 π f τ d τ X(f) = \int_{-\infty}^{\infty} x(\tau) e^{-i2\pi f \tau} \, d\tau X(f)=x(τ)ei2πfτdτ

因此:

Y ( f ) = X ( f ) ⋅ H ( f ) Y(f) = X(f) \cdot H(f) Y(f)=X(f)H(f)

总结

我们得到了频域中的乘积:

Y ( f ) = X ( f ) ⋅ H ( f ) Y(f) = X(f) \cdot H(f) Y(f)=X(f)H(f)

这表明时域中的卷积对应于频域中的乘积。这就是卷积定理的内容。

这篇关于[信号与系统]傅里叶变换、卷积定理、和为什么时域的卷积等于频域相乘。的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/1080080

相关文章

使用Python构建一个高效的日志处理系统

《使用Python构建一个高效的日志处理系统》这篇文章主要为大家详细讲解了如何使用Python开发一个专业的日志分析工具,能够自动化处理、分析和可视化各类日志文件,大幅提升运维效率,需要的可以了解下... 目录环境准备工具功能概述完整代码实现代码深度解析1. 类设计与初始化2. 日志解析核心逻辑3. 文件处

golang程序打包成脚本部署到Linux系统方式

《golang程序打包成脚本部署到Linux系统方式》Golang程序通过本地编译(设置GOOS为linux生成无后缀二进制文件),上传至Linux服务器后赋权执行,使用nohup命令实现后台运行,完... 目录本地编译golang程序上传Golang二进制文件到linux服务器总结本地编译Golang程序

Linux系统性能检测命令详解

《Linux系统性能检测命令详解》本文介绍了Linux系统常用的监控命令(如top、vmstat、iostat、htop等)及其参数功能,涵盖进程状态、内存使用、磁盘I/O、系统负载等多维度资源监控,... 目录toppsuptimevmstatIOStatiotopslabtophtopdstatnmon

linux重启命令有哪些? 7个实用的Linux系统重启命令汇总

《linux重启命令有哪些?7个实用的Linux系统重启命令汇总》Linux系统提供了多种重启命令,常用的包括shutdown-r、reboot、init6等,不同命令适用于不同场景,本文将详细... 在管理和维护 linux 服务器时,完成系统更新、故障排查或日常维护后,重启系统往往是必不可少的步骤。本文

Mac系统下卸载JAVA和JDK的步骤

《Mac系统下卸载JAVA和JDK的步骤》JDK是Java语言的软件开发工具包,它提供了开发和运行Java应用程序所需的工具、库和资源,:本文主要介绍Mac系统下卸载JAVA和JDK的相关资料,需... 目录1. 卸载系统自带的 Java 版本检查当前 Java 版本通过命令卸载系统 Java2. 卸载自定

基于Python实现一个简单的题库与在线考试系统

《基于Python实现一个简单的题库与在线考试系统》在当今信息化教育时代,在线学习与考试系统已成为教育技术领域的重要组成部分,本文就来介绍一下如何使用Python和PyQt5框架开发一个名为白泽题库系... 目录概述功能特点界面展示系统架构设计类结构图Excel题库填写格式模板题库题目填写格式表核心数据结构

全屋WiFi 7无死角! 华硕 RP-BE58无线信号放大器体验测评

《全屋WiFi7无死角!华硕RP-BE58无线信号放大器体验测评》家里网络总是有很多死角没有网,我决定入手一台支持Mesh组网的WiFi7路由系统以彻底解决网络覆盖问题,最终选择了一款功能非常... 自2023年WiFi 7技术标准(IEEE 802.11be)正式落地以来,这项第七代无线网络技术就以超高速

Linux系统中的firewall-offline-cmd详解(收藏版)

《Linux系统中的firewall-offline-cmd详解(收藏版)》firewall-offline-cmd是firewalld的一个命令行工具,专门设计用于在没有运行firewalld服务的... 目录主要用途基本语法选项1. 状态管理2. 区域管理3. 服务管理4. 端口管理5. ICMP 阻断

Windows 系统下 Nginx 的配置步骤详解

《Windows系统下Nginx的配置步骤详解》Nginx是一款功能强大的软件,在互联网领域有广泛应用,简单来说,它就像一个聪明的交通指挥员,能让网站运行得更高效、更稳定,:本文主要介绍W... 目录一、为什么要用 Nginx二、Windows 系统下 Nginx 的配置步骤1. 下载 Nginx2. 解压

如何确定哪些软件是Mac系统自带的? Mac系统内置应用查看技巧

《如何确定哪些软件是Mac系统自带的?Mac系统内置应用查看技巧》如何确定哪些软件是Mac系统自带的?mac系统中有很多自带的应用,想要看看哪些是系统自带,该怎么查看呢?下面我们就来看看Mac系统内... 在MAC电脑上,可以使用以下方法来确定哪些软件是系统自带的:1.应用程序文件夹打开应用程序文件夹