首页 > 其他分享 >傅里叶变换以及快速傅里叶变换在雷达信号处理中的应用

傅里叶变换以及快速傅里叶变换在雷达信号处理中的应用

时间:2022-10-04 18:33:48浏览次数:76  
标签:变换 信号处理 正弦波 频域 频率 傅里叶 时域

  傅里叶变换应该是在大一或者大二的时候就开始接触了,一直对其都是一知半解的状态。不是很清楚到底是干啥的,想趁着国庆假期好好学习一下(主要是算法太难了,想换换心情,算法太虐了)。本文参考了几位大佬的文章再加上一些FFT在雷达信号处理的用处以及自己的一点理解,另外本文只是简单了解傅里也变换的基础定义和内容,不涉及具体的公式推导。

1、傅里叶变换(FT)的目的

  傅里叶变换的目的是将时域(即时间域)上的信号转换到频率域上的信号。随着域的不同,对同一个事物的了解角度也就随之改变,因此有时在时域中不好处理和观察的信号,在频域上却能更好的观察和处理。具体傅里叶变换的公式和推导大家可以看一下相关书籍的推导过程。

2、为什么要进行傅里叶变换

  关于为什么要进行傅里叶变换,有一个很经典的声音的例子,如果大家感兴趣可以去B站搜《形象展示傅里叶变换》上面讲的很形象生动。

  首先发送一个每秒440拍的纯A音,意思是若你测量耳机或扬声器附近的气压,将其当成时间的函数,那么这个函数会在它的平衡点附近上下振荡。如下面照片所示。每秒钟会产生440次振荡。

   其次对于一个更低的音调,比如D音,与A音相比,结构是相同的,只是每秒节拍数变少。

   之后,当两个音同时播放时,最后的压强-时间图像是什么样的呢?

   在任意时刻,压强的变化就是每个音调产生的压强的总和。

  这个合成的东西看起来多少有点复杂,从图中的结果可以看出,在某些时刻A音和D音两个峰值会相互叠加产生一个很高的气压,而在另外的某些时刻,A音和D音的两个峰值相互抵消。总体来说最后得到一个波浪型的图,得到的图不是一个简单的正弦波,而是更复杂的波形。

两个音合成曲

 两峰值相互重合的情况

 两峰相互抵消

两个信号叠加动画图

  当我们加入更多的音调例如F音和C音,波形也会越来越复杂。

 多个音符合成的效果图

多个音符合成动态图

       但光从组合来看其实就是四个纯音符组合得到的结果,但从合成波形来看,却很复杂。而关于声音信号的提取装置-------麦克风来说,他在记录声音时,只能获取不同时刻的气压值,只知道最后的总和,也就是图中的合成信号。所以现在的问题就出来了,要如何通过测量的合成信号分解为其中纯音的频率呢?那就要我们的傅里叶变换登场了。傅立叶变换(的三角函数形式)的基本原理是:多个正余弦波叠加(蓝色)可以用来近似任何一个原始的周期函数(红色)

  可以简单地理解为,我们去菜市场买菜的时候,无论质量如何奇怪,都可以转变为“5个 1 斤的砝码,2个 1 两的砝码”来表示出来,那么上面的图我们也可以近似地想象成周期函数就是质量特别奇怪的物品,而正余弦波就是想像成成“我用了5个1号波、3个2号波”来表示这个周期函数。

  还是有点迷糊?那看看如何正弦波变成周期函数的(大佬们这个怎么画的可否评论区教我一下)。

       上面这些是不是和声音的例子结合了起来,不同声音是不同的音的组合(本质上也是波也可认为是正弦波),而周期函数是不同的正弦波组成。可能这样看还不够明显。那咱们分开看看(动态图中最后合成的波形的左上侧那个调节的n,个人认为应该是指的sin(nx)中的n)。

  这样是不是更清楚的看出来了。

3、时域和频域

  看完上面的知识点,这些基本上都是在时域里进行的,而傅里叶变换则是将状态从时域变换到频域中进行解析的,那么我们有必要学习一下时域和频域。

  为什么要进行频域分析 
  1)信号不只和时间有关,还和频率有关,在不同频率下信号的响应是不一样的,所以就要知道信号随着频率是怎么变化的。

  2)还有信号的计算,在时间域内往往要解微分方程,而用傅立叶和拉普拉斯变换到复频域后就变成了代数方程,求解起来很方便。
  3)从频谱图上可以看到幅值和相位随着频率变化是如何响应的,也可以求出系统的截止频率。

  4)无需求解微分方程,图解(频率特性图)法,间接揭示系统性能并指明改进性能的方向和易于实验分析.可推广应用于某些非线性系统(如含有延迟环节的系统)以及可方便设计出能有效抑制噪声的系统。

  关于时域和频域我们先看一下另一篇文章是如何讲解的,个人觉得很形象易懂。

  标准正式版

  1、时域(时间域)——自变量是时间,即横轴是时间,纵轴是信号的变化。其动态信号x(t)是描述信号在不同时刻取值的函数。

  2、频域(频率域)——自变量是频率,即横轴是频率,纵轴是该频率信号的幅度,也就是通常说的频谱图。
下面是图文讲解:

  图1是正弦波的时域图,示出了振幅与时间的关系。在时域图中,横轴是时间,纵轴是振幅。时域图显示振幅随时间的变化,可以看出峰值振幅为5V,可以算出频率f=6 Hz。

  图2是图1中正弦波的频域图在频域图中,横轴是频率,纵轴是峰值振幅。频域图仅仅示出峰值振幅与频率,而不显示振幅随时间的变化。从频域图可以看出,正弦波的频率为6Hz,这个6Hz的正弦波的峰值振幅为5V 。频域图的优点是,从频域图中,可以一眼看出正弦波的频率和峰值振幅整个正弦波在频域图上只是一个立柱立柱的位置显示了正弦波的频率立柱的高度显示了正弦波的峰值振幅 

  非正式

  看完正式版的我们再看一个不太正式的例子解释。

    时间是永远在流动的花谢花开、潮来潮往,世界永远在不停地变化,而以时间为参照系去看待这个世界,我们就叫它时域分析。就好像心电图一样,心电图是记录心脏每一心动周期所产生的电活动变化,所以随着时间变化心电图也会变化。这就是时域。

   而频域呢,就是描述信号在频率方面特性时用到的一种坐标系,频域就是装着正弦函数的空间,自然而然的,正余弦波是频域中唯一存在的波形。

 心电图符号

  我们从时域我们可以观察到心脏随着时间变化在不停地跳动的情形,但是从频域来看,就是一个简单的心电图符号。如果时域是运动永不停止的,那么频域就是静止的(不过要是心跳加快的话心电图符号应该会变吧,因为频率改变了)。

  我们从时域我们可以观察到心脏随着时间变化在不停地跳动的情形,但是从频域来看,就是一个简单的心电图符号。如果时域是运动永不停止的,那么频域就是静止的。

  在很多领域我们都可以用到时域和频域,在时域,我们观察到钢琴的琴弦一会上一会下的摆动,就如同一支股票的走势;而在频域,只有那一个永恒的音符。

  刚刚我们讲了多个正余弦波叠加可以用来近似任何一个原始的周期函数,我们心脏不同时间、不同强度的跳动就成了我们所看到的心电图。就可以看作正余弦波叠加成的周期函数。同样的,利用对不同琴键不同力度,不同时间点的敲击,可以组合出任何一首乐曲,也可以看作余弦波叠加成的周期函数。

  而对于信号来说,信号强度随时间的变化规律就是时域特性,信号是由哪些单一频率的信号合成的就是频域特性傅里叶变换实质涉及的是频域函数和时域函数的转换。

 

  通过上面我们对时域和频域的了解,可以知道对于由多个不同频率的正弦信号组成的矩形信号来说,他的频域图如下图所示。

  一般书上到这就结束了,然后就是你猜的阶段。其实教科书只要补一张图就足够了:频域图像,也就是俗称的频谱,就是

 

   再精细一点:

  在这两幅图中,最前面的黑色的线就是所有正弦波叠加而成的总和,也就是最接近矩形波的那个图形。而后面以不同颜色排列开的分别是sin(nx),n=1,3,5,7,......对应的正弦波。而这些正弦波正是组合成矩形波的分量。这些正弦波按照频率从低到高从前向后排列开来,而每一个波的振幅都是不同的。同时每两个正弦波之间都还有一条直线,那并不是分割线,而是振幅为0的正弦波!也就是说,为了组成特殊的曲线,有些正弦波成分是不需要的。

  这里,不同频率的正弦波我们称为频率分量。

  接下来就是关键的地方

  如果我们把第一个频率最低的频率分量看作“1”,我们就有了构建频域的最基本单元

  对于我们最常见的有理数轴,数字“1”就是有理数轴的基本单元。

  时域的基本单元就是“1秒”,如果我们将一个角频率为ω0的正弦波cos(ω0t)看作基础,那么频域的基本单元就是ω0。

  有了“1”,还要有“0”才能构成世界,那么频域的“0”是什么呢?cos(0t)就是一个周期无限长的正弦波,也就是一条直线!所以在频域,0频率也被称为直流分量,在傅里叶级数的叠加中,它仅仅影响全部波形相对于数轴整体向上或是向下而不改变波的形状

  说了那么多正弦波,可能有的已经忘了它的定义,既然气氛都到这了我们不回忆一下正弦波好像有点不合适。看图。

   正弦波就是一个圆周运动在一条直线上的投影。所以频域的基本单元也可以理解为一个始终在旋转的圆(可参见本文第二节最后一个动图,更有利于理解,算了再贴一次吧,我也不喜欢往上翻)

  介绍完这些后我们回头再看看矩形波的频域变换

   在频谱中,偶数项的振幅都是0,奇数项的振幅随着频率由前到后依次减小。

  从以上的信息,通过时域到频域的变换,我们得到了我们得到了一个从侧面看的频谱,但是这个频谱并没有包含时域中全部的信息。因为频谱只代表每一个对应的正弦波的振幅是多少,而没有提到相位。基础的正弦波A.sin(wt+θ)中,振幅,频率,相位缺一不可,不同相位决定了波的位置,所以对于频域分析,仅仅有频谱(振幅谱)是不够的,我们还需要一个相位谱。

4、相位谱信息

  关于相位谱信息话不多说直接上图,图中是7个波叠加的图(奇数4个,偶数三个)。

  鉴于正弦波是周期的,我们需要设定一个用来标记正弦波位置的东西。在图中就是那些小红点。小红点是距离频率轴最近的波峰,而这个波峰所处的位置离频率轴有多远呢?为了看的更清楚,我们将红色的点投影到下平面,投影点我们用粉色点来表示。当然,这些粉色的点只标注了波峰距离频率轴的距离,并不是相位。(大家一定要注意一点图中的相位和频率的坐标系,里面表示的不是相位,当时看这个图一直没绕出来这一点,当时坐标轴还让我困扰了一段时间)。

  这里需要纠正一个概念:时间差并不是相位差。如果将全部周期看作2Pi或者360度的话,相位差则是时间差在一个周期中所占的比例。我们将时间差除周期再乘2Pi,就得到了相位差。

  在完整的立体图中,我们将投影得到的时间差依次除以所在频率的周期,就得到了最下面的相位谱。所以,频谱是从侧面看,相位谱是从下面看。

  注意到,相位谱中的相位除了0,就是Pi。因为cos(t+Pi)=-cos(t),所以实际上相位为Pi的波只是上下翻转了而已。对于周期方波的傅里叶级数,这样的相位谱已经是很简单的了。另外值得注意的是,由于cos(t+2Pi)=cos(t),所以相位差是周期的,pi和3pi,5pi,7pi都是相同的相位。人为定义相位谱的值域为(-pi,pi],所以图中的相位差均为Pi。

  最后来一张时域、频域和相位域的总图:(总的来说我觉得相位谱我还是有点懵逼,后面理解了再补充吧)

  到这基本上学习完了傅里叶变换的一些概念的知识。既然是搞雷达DSP的就列一些傅里叶变换在雷达信号处理中的应用的知识。方便后面的记录和相关知识的学习。

5、FFT在雷达信号处理中的应用

  FFT是DFT的快速实现算法。快速傅里叶变换(FFT)和傅里叶反变换(IFFT)一直是最常用、最普遍的雷达信号处理方法,它在雷达信号处理中的作用包括:

  1)脉冲压缩。现代雷达中的脉压比较大。由于时域脉压处理的运算量大,故大多数采用频域脉冲压缩处理。无论是采用线性或者非线性调频脉冲信号还是相位编码脉冲信号,均可以采用FFT和IFFT实现脉冲压缩处理。

  2)步进频率综合处理。针对利用步进频率实现高分辨的雷达系统,都采用IFFT实现多脉冲回波的步进频率综合处理,实现高的距离分辨率。

  3)拉伸处理。对线性调频信号的拉伸处理系统接收的信号采用FFT处理,通过频率的高分辨实现距离的高分辨。

  4)调频连续波雷达距离的测量。调频连续波雷达的目标回波信号与发射信号的副本(耦合的发射泄露信号)混频后,通过FFT将距离的测量转换为频率的测量。

  5)相干积累和脉冲数较多的MTD。对一个波位的多个脉冲在各距离单元通过FFT实现相干积累或MTD处理。

  6)数字波束形成。针对均匀线阵,可以采用FFT实现多个波位的数字波束形成。

  7)雷达成像。在雷达成像中,距离维的高分辨、合成孔径处理得到横向的高分辨大多是采用FFT实现的。

6、参考文献:

 1、 什么叫时域和频域?

    2、 傅里叶分析之掐死教程(完整版)

       作 者:韩 昊    

  知 乎:Heinrich   微 博:

  @花生油工人    

  知乎专栏:与时间无关的故事    

  谨以此文献给大连海事大学的吴楠老师,柳晓鸣老师,王新年老师以及张晶泊老师。

 3、傅里叶变换(一)——认识傅里叶变换

 4、傅里叶变换看不懂,5分钟教你快速理解!

 5、(四十九)通俗易懂理解——傅里叶变换

 6、现代雷达系统分析与设计

    7、 官方双语】形象展示傅里叶变换

 

标签:变换,信号处理,正弦波,频域,频率,傅里叶,时域
From: https://www.cnblogs.com/sbb-first-blog/p/16750937.html

相关文章

  • 三--图像处理中的正交变换
    3.1傅里叶变换3.2离散余弦变换3.3沃尔什变换3.4哈尔函数及哈尔变换3.5斜矩阵与斜变换3.6小波变换 数字图像处理的方法分为:空间域处理法(空域法)频域法(变换......
  • 6. Z 字形变换
    将一个给定字符串s根据给定的行数numRows,以从上往下、从左到右进行 Z字形排列。比如输入字符串为"PAYPALISHIRING" 行数为3时,排列如下:PAHNAPLSIIGY......
  • 用于语义线检测的深度霍夫变换
    文章目录​​摘要​​​​论文的贡献​​​​主要的改进​​​​参考​​摘要我们专注于检测有意义的线结构的基本任务,也就是自然场景中的语义线。以往的许多方法都将此问题......
  • 雷达信号处理的基本方法------------雷达
    1、什么是雷达信号处理雷达信号处理是指对观测信号进行分析、变换、综合等处理,抑制杂波、干扰等非期望信号,增强有用信号,并估计有用信号的特征参数,或是将信号变成某种......
  • 【Numpy总结】第二节:Numpy 的属性与形状变换
    一、最基本的属性在NumPy中,每一个线性的数组称为是一个轴(axis),也就是维度(dimensions),维度的数量称为秩(rank);比如说,二维数组相当于是两个一维数组,即二维数组有两个轴,秩为2。......
  • 跟我学Python图像处理丨傅里叶变换之高通滤波和低通滤波
    摘要:本文讲解基于傅里叶变换的高通滤波和低通滤波。本文分享自华为云社区《[Python图像处理]二十三.傅里叶变换之高通滤波和低通滤波》,作者:eastmount。一.高通滤波傅......
  • 跟我学Python图像处理丨傅里叶变换之高通滤波和低通滤波
    摘要:本文讲解基于傅里叶变换的高通滤波和低通滤波。本文分享自华为云社区《​​[Python图像处理]二十三.傅里叶变换之高通滤波和低通滤波​​》,作者:eastmount。一.高通滤......
  • 13、OpenCV4.4 仿射变换
    Parameters:src–Coordinatesoftriangleverticesinthesourceimage.dst–Coordinatesofthecorrespondingtriangleverticesinthedestinationimage.Thefunc......
  • 信号处理(1) --常用信号平滑去噪的方法
    前言:最近研究汽车碰撞的加速度信号,在信号的采集过程中难免遇到噪音,导致信号偏差,为了更好的反映系统情况,故常需要信号去噪,本文分享一些常用信号平滑去噪的方法。关键字:信号;去......
  • 坐标系变换——“旋转矩阵/欧拉角/四元数”
    向量的旋转一共有三种表示方法:旋转矩阵、欧拉角和四元数,接下来我们介绍一下每种旋转方法的原理以及相互转换方式。旋转矩阵坐标变换的作用在一个机器人系统中,每个测量元......