压缩感知与奈奎斯特采样频率比较报告
从奈奎斯特采样到压缩感知拓展教学方法
从奈奎斯特采样到压缩感知拓展教学方法
盛志超;方勇;徐强荣;余鸿文;黄知雨
【期刊名称】《电气电子教学学报》
【年(卷),期】2024(46)1
【摘要】从“信号与系统”到“数字信号处理”,采样定理都是重要的教学内容。
但是在工程应用中,产生大量数据造成存储空间的极大浪费,而压缩感知突破奈奎斯特采样定理的限制,能够实现远低于奈奎斯特频率的采样。
为适应新工科背景下的教学改革,让学生接触前沿研究成果,压缩感知被引入作为传统奈奎斯特采样定理教学的补充和拓展,取得良好的教学效果。
【总页数】6页(P164-169)
【作者】盛志超;方勇;徐强荣;余鸿文;黄知雨
【作者单位】上海大学通信与信息工程学院
【正文语种】中文
【中图分类】G642
【相关文献】
1.基于采样值随机压缩矩阵核空间的亚奈奎斯特采样重构算法
2.亚奈奎斯特采样雷达的运动目标回波信号的快速重构
3.高速中高精度奈奎斯特采样ADC结构综述
4.基于非正交波形的超奈奎斯特采样
因版权原因,仅展示原文概要,查看原文内容请购买。
浅谈压缩感知
浅谈压缩感知根据传统的奈奎斯特采样定律,采样速率必须大于原始信号最高频率的两倍才能保证完全重建原始信号,但是最近十几年信号的带宽和最高频率都有了比较大的变化,这样一来就要求采样速率和处理速度要更高,如此一般,对于高分辨率的信号数据的采样、传输、存储就是一个比较大的问题。
这个问题先放着,我们看另一个问题,我们都知道信息论可以指导我们对数据进行压缩,压缩的前提是数据的信息之中存在着冗余,所谓信息的冗余,在信息论当中指的就是可以确定,或者可以根据其他信息推测出的数据,如果能将这种数据全部去除,只保留无法根据其他信息确定的信息,那么就实现了数据的压缩。
于是有人就在考虑,高速采样之后进行数据压缩,太浪费系统资源了,不如我们先处理一下这个原始的高速宽带信号,在保证信息熵无损或可以接受的范围的情况下,建立一个新的信号,之后对新的信号进行低速采样,同时还能重建原始信号。
后来这种想法经过发展,就成为了目前的压缩感知,或者更通俗的说法,就是压缩采样。
那现在我们来看一下,实现压缩感知需要的步骤和要求是什么。
前提,信号要有稀疏性。
首先,需要将原始信号进行一定的变换,得到新的信号,暂且称之为预变换。
新的信号速率不能太高,通俗的说,这是一个稀疏信号,并且这个稀疏信号携带的信息量,不能比原始信号低多少。
之后是对稀疏信号的采样,并将稀疏信号还原为原始信号,暂且称之为后处理。
可以看出,压缩感知虽然降低了采样速率,但实际上因为预变换和后处理,增加了实现的计算复杂度,这体现了一个世界的基本道理:凡事都是有代价的,有多大的优势,就要付出多大的努力。
我们继续回到信息论,如果以信息熵和符号的角度去衡量数据压缩的过程,实际上就是信息熵在符号上的再分配,并且这种分配方式的方向是朝着符号平均信息量变大,并且接近某一平均值的过程。
这过程在压缩感知上的表现,就是对信号在损失信息熵可接受的程度上进行某个变换域处理,并且变换之后的信号是稀疏的。
那么压缩感知第一步需要做的,就是找到这样一个稀疏域,而找到稀疏域过程中最为关键的一点是找到或者构建适合某类信号的正交基底来表示原始信号,对于多种不同类型的原始信号来说,就是找出一本能够根据信号类型选择合适正交基底的字典。
奈奎斯特采样频率求解
奈奎斯特采样频率求解摘要本文将介绍奈奎斯特采样频率的概念以及如何进行求解。
我们首先会解释为什么需要奈奎斯特采样频率,在此基础上提供了一种简单的计算方法。
同时,我们还会探讨一些与奈奎斯特采样频率相关的重要概念和实际应用。
希望通过这篇文档,您能够更好地理解奈奎斯特采样频率的原理和计算方法。
1.引言在信号处理和通信系统中,采样是一个非常重要的过程。
奈奎斯特采样频率是指在数字信号处理中,为了能够完美地重构原始模拟信号,需要对模拟信号进行采样的最小频率。
本文将详细介绍奈奎斯特采样频率的定义和计算方法。
2.奈奎斯特采样频率的背景在进行模拟信号的数字化处理时,我们需要将连续的模拟信号转化为离散的数字信号进行处理。
采样是这个过程中的第一步,它将连续的信号在时间上进行离散化。
然而,如果采样频率过低,将会导致采样结果中丢失了一些信号的信息。
为了在数字信号中完美地重构原始模拟信号,我们需要满足一定的采样频率。
3.奈奎斯特采样频率的定义奈奎斯特采样频率就是在理论上最低有效采样频率。
根据奈奎斯特定理,为了保证完美重构,采样频率必须是信号带宽的两倍以上。
因此,奈奎斯特采样频率的定义可以表达为:奈奎斯特采样频率=2×信号带宽4.奈奎斯特采样频率的计算方法为了计算奈奎斯特采样频率,我们需要知道信号的带宽。
信号的带宽是指信号的最高频率成分与最低频率成分之间的差异。
根据信号的具体情况,我们可以通过以下几种方法计算信号的带宽:-如果信号是理想低通滤波器的输出,那么信号的带宽就是滤波器的截止频率。
-如果信号是多个频率成分的叠加,那么信号的带宽就是最高频率成分与最低频率成分之间的差异。
在得到信号的带宽后,我们可以根据奈奎斯特采样频率的定义计算得到奈奎斯特采样频率。
以下是奈奎斯特采样频率的计算公式:奈奎斯特采样频率=2×信号带宽5.奈奎斯特采样频率的应用奈奎斯特采样频率在信号处理和通信系统中有着广泛的应用。
以下是一些常见的应用场景:-音频信号处理:在数字音频系统中,为了能够完美地重构原始音频信号,需要使用至少符合奈奎斯特采样频率的采样频率。
感知压缩理论
• 图像还原。
TITLE 01
• 医学成像
TITLE 03
• 线性编码
ADD YOUR TEXT HERE
感知压缩应用
TITLE 02
• 天文学
TITLE 04
• 图像还原
感知压缩应用
Perceptual compression applications
最典型的例子就是医学图像成像,例如断层扫描 (CT)技术和核磁共振(MRI)技术。这两种成像 技术中,仪器所采集到的都不是直接的图像像素,而 是图像经历过全局傅立叶变换后的数据。也就是说, 每一个单独的数据都在某种程度上包含了全图像的信 息。
TITLE 03
Lorem ipsum dolor sit amet, consectetuer adipiscing elit, sed diam nonummy nibh euismod tincidunt ut laoreet dolore magna aliquam erat volutpat.
T EMP L A T E
Perceptual compression theory and its application
感知压缩
1 引言 2 感知压缩理论 3 感知压缩应用
ACADEMIC DESIGN
前言 PREFACE
引言
传统的信号采样以奈奎斯特采样定理为基础,在获取信号时,为了 不丢失信号的信息,采样频率必须大于信号中最高频率的两倍,才能精确 重构信号。事实上,大部分采样数据是不重要的。在信号或图像处理 的过程中,只保留了一些重要的数据,丢弃了大量的数据。残差数据、 重构信号或图像不会引起视觉差异。
具有直接信息采样特性。由于从理论上讲任
何信号都具有可压缩性,只要能找到其相应的 稀疏T表ITL示E 空01 间,就可以有效地进行压缩采样 ,这 一理论必将给信号采样方法带来一次新的革 Lorem ipsum dolor sit amet,
基于压缩感知的图像处理
基于压缩感知的图像处理基于压缩感知的图像处理一、压缩感知在过去的几十年里,人们获取数据的能力不断提高,需要处理的数据量也越来越大,因此信号的带宽也越来越大,所以对信号处理的速度和采样速率的要求也随之提高。
众所周知,奈奎斯特采样定理要求采样率不得低于信号带宽的两倍,这对目前的信号处理能力提出了巨大的挑战。
所以人们试图找到一种新的信号处理技术。
近年来提出了一种新的信号处理理论——压缩感知理论。
压缩感知理论表明:如果信号是稀疏的或者是可压缩的,就可以通过一个测量矩阵将其投影到一个低维的空间上,得到的低维信号成为测量信号,然后将这个测量信号进行传输,在接收端通过接收到的信号和已知的测量矩阵来重构出原始的信号。
理论上指出任何信号经过一定处理后都可以转化为稀疏信号,这也为压缩感知理论在各个领域的广泛使用提供了保障。
1、压缩感知理论传统的信号处理过程包括信号的采样、压缩、传输和重构四个部分,根据奈奎斯特采样定理,信号的采样速率不能低于信号最大带宽的两倍,只有以满足这一要求的采样速率进行采样,才能保证信息不丢失,但是在很多情况下,奈奎斯特采样速率显得很高,实现起来比较困难。
压缩感知是一种新的信号获取的方法,它突破了奈奎斯特采样定理的瓶颈,它将对信号的压缩和采样合并进行,使得测量数据量远远小于传统的采样方法所得的数据量。
压缩感知主要包括三个方面的内容:信号的稀疏表示、信号的压缩采样和信号的重构。
2、信号的稀疏表示前面提到,压缩感知理论只能直接应用于稀疏信号。
如果需要处理的信号是稀疏的,那就不需要稀疏表示这一部分,直接进行压缩采样就行了,但是就目前来看,我们所要处理的大多数信号都不是稀疏信号,这就需要将其转换为稀疏信号。
假设ψ=[ψ1, ψ2, ψ3, , ψN ]为R 空间上的一组基,Ψi (i=1,2,3…N)是N一个N*1的列向量,考虑x =[x 1, x 2, x 3, , x N ]T ,它是一个实值有限长的ψ线性表示:N x ∈R 一维离散信号,。
奈奎斯特采样和压缩感知
奈奎斯特采样和压缩感知奈奎斯特采样和压缩感知:从理论到应用的探究引言在信息处理领域,信号的采样和压缩是两个关键的概念。
奈奎斯特采样理论和压缩感知是两种常用的方法,它们在传感器网络、通信系统、图像处理等领域都得到了广泛的应用。
本文将深入探讨奈奎斯特采样和压缩感知的原理、应用以及个人观点。
1. 奈奎斯特采样的原理和应用奈奎斯特采样是用于从连续时间信号中获取离散时间采样的方法,它基于奈奎斯特——香农采样定理。
根据这个定理,为了完全恢复原始信号,采样频率必须大于信号的最高频率的两倍。
奈奎斯特采样的原理可以简化为“至少两倍采样频率”。
采样频率低于此阈值会导致信号失真,无法完全还原。
奈奎斯特采样在实际应用中有着广泛的用途。
在通信系统中,奈奎斯特采样保证了信号的信息不会丢失。
在图像处理中,奈奎斯特采样确保图像的每个像素都得到准确的采样。
这种采样方法在模拟信号转换为数字信号时起着至关重要的作用。
2. 压缩感知的原理和应用压缩感知是一种通过从稀疏信号中获取少量线性投影来重构信号的技术。
相比于传统的采样方法,压缩感知可以实现更高效的信号采样和信号重构,从而极大地减少数据传输和存储的需求。
压缩感知的原理基于两个重要的概念:稀疏表示和随机投影。
稀疏表示指的是信号可以用较少的非零系数表示。
随机投影是指通过在信号上进行线性投影来得到一组稀疏的测量结果。
通过这种方式,压缩感知能够仅使用较少的测量结果来还原信号,从而实现高效的信号处理。
压缩感知在许多领域都有重要的应用。
在无线传感器网络中,压缩感知可以减少传感器数据的传输量,延长网络寿命。
在医学影像处理中,压缩感知能够减少医学影像数据的存储需求,提高图像传输速度。
3. 个人观点和理解奈奎斯特采样和压缩感知作为信号处理领域的两个重要概念,具有各自的优势和应用场景。
奈奎斯特采样保证了信号的完整性和准确性,适用于连续时间信号的离散化处理。
而压缩感知则通过提取信号的稀疏表示,实现高效的信号采样和处理,适用于稀疏信号的重构和压缩。
奈奎斯特采样率与压缩感知学习报告
数字信号处理第一次大作业奈奎斯特采样定理与信号稀疏采样学习报告专业:信息对抗技术学生姓名:石星宇02123010指导教师:吕雁目录奈奎斯特采样定理与信号稀疏采样学习报告 (1)一、奈奎斯特采样定理 (1)1、奈奎斯特采样定理说明 (1)2、信号的采样与恢复 (1)3、相关代码 (3)4、关于奈奎斯特采样定理的一些问题 (5)二、信号稀疏采样 (5)1、为什么要提出信号的稀疏采样 (5)2、压缩感知概述 (6)3、压缩感知基本概念 (6)4、压缩感知仿真 (7)5、压缩感知仿真程序 (8)三、总结 (9)四、参考资料 (10)奈奎斯特采样定理与信号稀疏采样学习报告一、奈奎斯特采样定理1、奈奎斯特采样定理说明采样过程所应遵循的规律,称为取样(采样)定理、抽样定理。
采样定理说明采样频率与信号频率之间的关系,是连续信号离散化的基本依据。
在进行模拟/数字信号的转换过程中,当采样频率s f 大于等于信号中最高频率c f 的2倍时,采样之后的数字信号完整地保留了原始信号中的信息,可由采样得到的数字信号恢复原来的模拟信号。
一般实际应用中保证采样频率为信号最高频率的5~10倍。
采样定理又称奈奎斯特采样定理。
将c s f f 2=称为奈奎斯特频率。
2、信号的采样与恢复结合实例,说明奈奎斯特采样定理与内插恢复的应用。
假设有模拟信号()t f t f t x a 212cos 2cos ππ+=,其中Hz f Hz f 50,2021==。
该信号波形及频谱如下图所示:对信号()t f t f t x a 212cos 2cos ππ+=以采样频率为Hz f f s 10022==进行采样,得到如下所示的离散时间信号,即序列()s s nT f nT f nT x 212cos 2cos ππ+=,其中s s f T /1=。
该序列的频谱如下:由此可见,采样过程对原始信号的频谱有一定的影响。
但是随着采样频率的逐渐增加,会使得采样信号的频谱与原始信号的频谱逐渐接近。
浅谈压缩感知(十二):压缩感知与奈奎斯特采样定理
浅谈压缩感知(十二):压缩感知与奈奎斯特采样定理奈奎斯特采样定理:定理:为了不失真地恢复模拟信号,离散信号系统的采样频率不小于模拟信号频谱中最高频率的2倍。
在时域上,频带为F的连续信号f(t)可用一系列离散的采样值f(t1),f(t1+Δt),f(t1+2Δt)…来表示,只要这些采样点的时间间隔Δt <= 1/2F,便可根据各采样值完全恢复原始信号。
在频域上,当时间信号函数f(t)的最高频率分量为fmax时,f(t)的值可由一系列采样间隔小于或等于1/2fmax的采样值来确定,即采样点的重复频率为fs >= 2fmax。
采样定理指出,只要离散系统的奈奎斯特频率高于采样信号的最高频率或带宽,就可以避免混叠现象。
从理论上说,即使奈奎斯特频率恰好大于信号带宽,也足以通过信号的采样重建原信号。
但是,重建信号的过程需要以一个低通滤波器或者带通滤波器将在奈奎斯特频率之上的高频分量全部滤除,同时还要保证原信号中频率在奈奎斯特频率以下的分量不发生畸变,而这是不可能实现的。
在实际应用中,为了保证抗混叠滤波器的性能,接近奈奎斯特频率的分量在采样和信号重建的过程中可能会发生畸变。
因此信号带宽通常会略小于奈奎斯特频率,具体的情况要看所使用的滤波器的性能。
需要注意的是,奈奎斯特频率必须严格大于信号包含的最高频率。
如果信号中包含的最高频率恰好为奈奎斯特频率,那么在这个频率分量上的采样会因为相位模糊而有无穷多种该频率的正弦波对应于离散采样,因此不足以重建为原来的连续时间信号。
压缩感知:压缩感知:作为一个新的采样理论,通过利用信号的稀疏特性,在远小于Nyquist采样率的条件下,用随机采样获取信号的离散样本,然后通过非线性重建算法完美重建信号。
提出背景:众所周知,在奈奎斯特采样定理为基础的传统数字信号处理框架下,若要从采样得到的离散信号中无失真地恢复模拟信号,采样速率必须至少是信号带宽的两倍。
然而,随着当前信息需求量的日益增加,信号带宽越来越宽,在信息获取中对采样速率和处理速度等提出越来越高的要求。
压缩感知在数字电视信号频谱检测中的应用
o e o sa p i ai n , EE 0 .2 h so e e eVHF UHF T b n o n t er d o u e s I r e t iet eT b n r p r n fi p lc t s I E8 2 2 a p n d t t o h / V a d t c g i v a i s r . n o d r ou i z V a d p o e- o i t l h
i g n p l ain , 0 2 4 ( ) 1 6 1 8 n dA pi t s 2 1 , 8 3 : 2 —2 . a c o
Ab t a t Ast e s e tu r s u c si r a h r g u rn l, o n t e r d o h se r e o e c e t s e i l p cr m. sr c : p cr m e o r e i g e t o t ec re t c g i v a i a me g d t f i n l u e t e s e t h n s a y i i y h d u As
ld ot tT esmpe a frdt n l y ussmpigihg , hc s s p c n me Sn e V s a i sas e a yi r n. h l rt o a io a N q i a l ih w i wat aeadt . ic g ls p r i t mp a a e t i t n s h es i T i n e nh
C m u r ni ei d p lai s o p t gn rga A pi tn 计算 机工程 与应用 eE e n n e o
压缩感知在数字 电视信 号频谱检测 中的应用
信号处理发展方向
个人认为信号处理方向有以下几个可发展方向:1.压缩感知(Compressed sensing, CS)CS是近五年来信号处理领域最为火热的研究点,其核心思想是使用低于奈奎斯特采样率的频率采样,然后尽可能完全恢复原始信号。
已经在数学上证明,一般信号都能在某一个映射域内具有稀疏性,因此可以在该映射域采样少量的样本点发射出去,然后在接收端通过反变换恢复原始信号。
目前该领域发展迅速,最早应用在图像处理(如单像素相机),现已应用至通信信号、雷达信号、语音信号、天文信息处理等各方面。
应用在空天:可应用在飞行器的海量数据处理上。
飞行器意味着快速与海量,需求信号采集量小、传输带宽小、传输速度极快、抗干扰能力强,CS 正是解决这些问题的好办法。
2.干涉合成孔径雷达现在较火热的雷达有干涉合成孔径雷达(Interferometric Synthetic Aperture Radar,InSAR),已有十几年的发展。
该雷达可机载(机载叫逆合成孔径,原理相似),也可星载。
星载InSAR可成像三维地表,不仅有地面二维图像,还具有高度维信息(原理是仿照人的双眼,一只眼睛只能成二维像,两只有一定距离的眼睛便可成三维像)。
InSAR在遥感遥测领域是一个研究热点。
应用在空天:可应用在飞行器上做干涉的三维成像,获取地表的多种信息,如地貌、植被、建筑等。
由于具备高度维信息,比传统的地面成像更具有价值。
3.认知无线电(Cognitive radio, CR)CR是下一代通信5G的频谱管理技术。
目前美国由于频谱资源紧张得已经迫在眉睫,而已分配的频谱利用率又低下,于是开发CR技术允许非授权用户接入授权用户频段进行通信。
FCC在2008年已经开放空白电视频段用于CR 技术的研究与实验。
当前美国各领域都在呼吁全面开放频段,迎接5G通信。
因此现在CR技术由纯理论发展到可商用的阶段,当前还有很多实际问题并未解决,存在很多可研究的话题。
应用在空天:可应用在空天信号传输时对所占频段采用CR的非授权接入方式,如绕隙通信。
奈奎斯特采样率和稀疏采样学习报告
奈奎斯特采样率和稀疏采样学习报告1.采样定理数字信号处理系统的基本组成(1)前置滤波器将输入信号xa(t)中高于某一频率(称折叠频率,等于抽样频率的一半)的分量加以滤除。
(2)A/D变换器在A/D变换器中每隔T秒(抽样周期)取出一次x a(t)的幅度,采样后的信号称为离散信号。
在进行A/D信号的转换过程中,当采样频率fs.max大于信号中最高频率fmax的2倍时(fs.max>2fmax),采样之后的数字信号完整地保留了原始信号中的信息,一般实际应用中保证采样频率为信号最高频率的5~10倍;采样定理又称奈奎斯特定理。
1.1 在时域频带为F的连续信号 f(t)可用一系列离散的采样值f(t1),f(t1±Δt),f(t1±2Δt),...来表示,只要这些采样点的时间间隔Δt ≤1/2F,便可根据各采样值完全恢复原始信号。
1.2 在频域当时间信号函数f(t)的最高频率分量为fmax时,f(t)的值可由一系列采样间隔小于或等于1/2fo的采样值来确定,即采样点的重复频率fs ≥2fmax。
2.奈奎斯特采样频率2.1 概述奈奎斯特采样定理:要使连续信号采样后能够不失真还原,采样频率必须大于信号最高频率的两倍(即奈奎斯特频率)。
奈奎斯特频率(Nyquist frequency)是离散信号系统采样频率的一半,因哈里·奈奎斯特(Harry Nyquist)或奈奎斯特-香农采样定理得名。
采样定理指出,只要离散系统的奈奎斯特频率高于被采样信号的最高频率或带宽,就可以真实的还原被测信号。
反之,会因为频谱混叠而不能真实还原被测信号。
采样定理指出,只要离散系统的奈奎斯特频率高于采样信号的最高频率或带宽,就可以避免混叠现象。
从理论上说,即使奈奎斯特频率恰好大于信号带宽,也足以通过信号的采样重建原信号。
但是,重建信号的过程需要以一个低通滤波器或者带通滤波器将在奈奎斯特频率之上的高频分量全部滤除,同时还要保证原信号中频率在奈奎斯特频率以下的分量不发生畸变,而这是不可能实现的。
压缩感知——精选推荐
压缩感知⼀、什么是压缩感知(CS)?compressed sensing⼜称compressed sampling,CS是⼀个针对信号采样的技术,它通过⼀些⼿段,实现了“压缩的采样”,准确说是在采样过程中完成了数据压缩的过程。
因此我们⾸先要从信号采样讲起:1. 我们知道,将模拟信号转换为计算机能够处理的数字信号,必然要经过采样的过程。
问题在于,应该⽤多⼤的采样频率,即采样点应该多密多疏,才能完整保留原始信号中的信息呢?------------------------------------------------------------------------------------------------------------------------------------------------------------2. 奈奎斯特给出了答案——信号最⾼频率的两倍。
⼀直以来,奈奎斯特采样定律被视为数字信号处理领域的⾦科⽟律。
-----------------------------------------------------------------------------------------------------------------------------------------------------------3. ⾄于为什么是两倍,学过信号处理的同学应该都知道,时域以τ为间隔进⾏采样,频域会以1/τ为周期发⽣周期延拓。
那么如果采样频率低于两倍的信号最⾼频率,信号在频域频谱搬移后就会发⽣混叠。
-----------------------------------------------------------------------------------------------------------------------------------------------------------4. 然⽽这看似不容置疑的定律却受到了⼏位⼤神的挑战。
压缩感知理论及其OMP算法FPGA实现研究
韩 林 王正彦 孟 南南 ( 青岛大学自 动化工程学谠, 山东 青岛 2 6 6 0 7 1 )
摘 要
奈 奎 斯 特 定 理要 求采 样 频 率 不 得低 于信 号 最 高频 率 的 2倍 , 这使 得 高频 信 号 的 硬件 采样 实现 变得 较 为 困难 。 压缩 感 知
( C o mp r e s s i v e S e n s i n g, CS ) 理 论 从 研 究信 号 的 稀 疏 性 出发 , 指 出在 一 定条 件 下 可 以用 低 于奈 奎 斯 特 定 理 的 频 率 对 信 号
进行采样。 介 绍 了压 缩 感 知 理 论及 其 OMP 重构 算 法 , 设 计 了 OMP重 构 算 法 的 F P G A 实现 的 总 体 框 图 和各 模 块 框 图 , 编 写 了V e r i l o g HD L程 序 代 码 , 并 给 出 了在 Qu a r t u s I I 中的仿 真 结果 , 和 Ma t l a b仿 真 结 果对 比 , 压 缩 重构 效 果 比 较 理 想 。
关键词 : 压 缩感 知 , F P G A, O MP
Ab s t r a c t Ny qu i s t t h e or e m r e qu i r e s t h at t h e s ampl i n g f r equ e n cy i s n ot l es s t h a n 2 t i me s t h e hi gh es t f r equ e nc y s i gn a l , whi ch mak e s t he hi gh- f r e qu en c y s i gn al s ampl i n g h ar dwa r e i mpl e me n t a t i o n a n d r ap i d i n f or ma t i o n pr oc e s si n g f a ci n g en or mo u s c h al l e nge s .
奈奎斯特频率和采样频率的关系
奈奎斯特频率和采样频率的关系一、引言奈奎斯特频率和采样频率是数字信号处理中非常重要的概念,它们直接影响着数字信号的采样、重构和滤波等过程。
本文将从基础概念、数学推导和实际应用等方面,全面介绍奈奎斯特频率和采样频率的关系。
二、基础概念1. 奈奎斯特定理奈奎斯特定理是数字信号处理中最基本的定理之一,它指出:如果一个连续时间信号的最高频率为fmax,那么在进行采样时,采样频率fs 必须大于2*fmax才能完全还原原始信号。
2. 奈奎斯特频率奈奎斯特频率也称为折叠频率或Nyquist折叠频率,是指当采样频率fs固定时,能够被完全还原的最高模拟信号的频率。
其计算公式为:f_nyquist = fs / 23. 采样频率采样频率是指对连续时间信号进行离散化时所使用的每秒采样次数。
在数字信号处理中,通常使用赫兹(Hz)作为单位。
其计算公式为:fs = 1 / T其中T为采样间隔时间。
三、数学推导1. 采样定理根据奈奎斯特定理,为了完全还原原始信号,采样频率必须大于等于2倍的最高模拟信号频率。
即:fs >= 2*fmax这个条件称为采样定理。
2. 折叠现象当采样频率小于2倍的最高模拟信号频率时,就会出现折叠现象。
折叠现象是指在重构过程中,高于奈奎斯特频率的信号被错误地重构成低于奈奎斯特频率的信号,从而导致信息丢失和失真。
3. 数学推导设原始模拟信号为x(t),其傅里叶变换为X(f)。
将x(t)进行离散化得到序列x[n],其傅里叶变换为X(e^jw)。
其中w=2*pi*f/fs。
根据采样定理可得:fs >= 2*fmax即:w <= pi因此,在重构过程中,只需要保留-w到w范围内的分量即可还原原始信号。
但是,当w>pi时,由于周期性扩展的存在,会出现折叠现象。
四、实际应用1. 音频采样在音频采样中,通常使用的采样频率为44.1kHz或48kHz,这是因为人耳听觉范围的最高频率为20kHz左右。
基于压缩感知的线性调频信号参数估计
基于压缩感知的线性调频信号参数估计闫浩;董春曦;赵国庆【摘要】提出一种基于压缩感知(Compressed Sensing,CS)理论的线性调频(Linear Frequency Modulation,LFM)信号参数估计算法.考虑到LFM信号在最佳分数阶Fourier变换(Fractional Fourier Transform,FrFT)域中是稀疏信号,对变换阶次进行粗搜索与精搜索,利用CS恢复信号在各个阶次FrFT矩阵中的系数向量,通过二维搜索.得到最佳变换阶次,进而得到信号的调频斜率和起始频率.在窄带干扰条件下,将形态学成分分析应用于算法中,提高了算法的抗干扰性能.实验结果表明,在随机采样点数远低于奈奎斯特采样点数的情况下,该算法能够准确估计信号参数,并且对高斯白噪声和强窄带干扰不敏感.【期刊名称】《电波科学学报》【年(卷),期】2015(030)003【总页数】8页(P449-456)【关键词】压缩感知;分数阶Fourier变换;线性调频信号;参数估计【作者】闫浩;董春曦;赵国庆【作者单位】西安电子科技大学电子信息攻防对抗与仿真技术教育部重点实验室,陕西西安710071;西安电子科技大学电子信息攻防对抗与仿真技术教育部重点实验室,陕西西安710071;西安电子科技大学电子信息攻防对抗与仿真技术教育部重点实验室,陕西西安710071【正文语种】中文【中图分类】TN971+.1引言线性调频(Linear Frequency Modulation, LFM)信号广泛应用于雷达、声纳和通信等系统中,LFM信号参数估计是电子战信号处理领域中一个重要的问题.在传统奈奎斯特采样框架下,国内外学者提出了很多LFM信号参数估计的方法,有基于最大似然[1]、短时Fourier变换和小波变换[2]、Wigner-Ville分布[3]、Radon-Wigner变换[4]、Randon-Ambiguity变换[5]和分数阶Fourier变换[6](Fractional Fourier Transform, FrFT)等方法.以上的方法都存在一个问题,就是在奈奎斯特采样框架下,随着LFM信号带宽不断增大,对信号的采样频率也越来越高,这给战场中用于信号采集、传输和处理的硬件系统造成极大压力.虽然国内学者提出的欠采样结合解线性调频的方法[7],对采样频率的要求降低了,但存在频率模糊问题,需要在后续进行解模糊的处理[8].因此,如何寻找一种新的信号参数估计算法来降低LFM信号带宽过宽对采样系统造成的压力是目前亟待解决的问题.压缩感知(Compressed Sensing, CS)[9-11]理论的提出给LFM信号处理问题带来了新的思路,以往在应用CS理论进行信号处理时,大多是以精确重构信号为目的,然而在应用CS理论进行信号检测与参数估计时,可以在不完全重构信号的情况下,达到信号检测与参数估计的目的[12-15].文献[16]构造波形匹配字典,达到LFM信号检测的目的,但是在没有先验知识的情况下,构造字典的参数不易选取,而且字典原子数目过于庞大,在利用优化算法恢复信号时,每次都要在所有原子中遍历,算法计算量过大,并且在强窄带干扰情况下,信号检测成功率不高.文献[17]提出的方法可以在强窄带干扰下估计LFM信号调频斜率,但同样是根据信号波形构造冗余字典,而且没有考虑调频斜率与起始频率的联合估计问题.文献[18]利用LFM信号在时频域的稀疏性,通过重构信号的短时Fourier变换,对信号调频斜率进行估计,但同样不能估计信号的起始频率.文中利用LFM信号在FrFT域的稀疏性,提出了一种基于CS的FrFT域LFM信号参数估计算法.首先构造Fourier变换矩阵,作为信号稀疏分解的字典,对窄带干扰信号进行抑制;然后利用LFM信号在最佳FrFT域中具有能量聚集的特性,对变换阶次p进行搜索,构造各个阶次的FrFT矩阵,通过恢复算法得到信号在各个阶次FrFT矩阵中的系数向量;最后通过二维搜索,得到最佳变换阶次,进而得到LFM信号的调频斜率和起始频率的估计.该算法用于参数估计的压缩采样点数远低于奈奎斯特框架下的采样点数,能够减轻战场中信号采样系统的负担,降低功耗和成本,缓解数据存储和传输的压力,提高战场信息传输的实时性.由于在低信噪比(Signal-to-Noise Ratio, SNR)情况下,LFM信号在FrFT域中仍然表现为稀疏信号,相比利用信号波形构造字典的算法,该算法能在更低的SNR情况下,获得更高的信号参数估计成功概率.1 压缩感知理论假设复信号x∈CN是N维的列向量,矩阵Ψ≡[ψ1,ψ2,…,ψZ]的列向量ψi(i=1,…,Z)构成复向量空间CN的一个字典,该字典可以是一个标准正交基,也可以是一个冗余框架,其中Z为字典中原子的个数.信号x可以由字典Ψ展开为(1)式中,Θ=[θ1,…,θZ]T为信号x在Ψ域中的Z×1维复系数向量,T表示向量转置.若系数向量Θ中非零元素的个数K≪N,或者当将Θ中的元素按从大到小的顺序排列,且其按照一定量级幂次速度衰减,若大系数[11]个数K≪N,则认为信号x在Ψ域中是稀疏信号或可压缩信号,且稀疏度为K.CS理论的思想是如果N维信号x在某个变换基Ψ下是稀疏的或可压缩的,且稀疏度为K,则可以利用一个与变换基Ψ不相关的M×N(K<M≪N)维观测矩阵Φ,得到原信号x的压缩采样信号:y=Φx=ΦΨΘ.(2)利用信号x的少量观测值y、变换基Ψ和观测矩阵Φ,通过求解一个优化问题,就能够以高概率精确恢复原始信号x[9-10].具体的压缩采样过程可以通过模拟信息转换器(Analog to Information Converter, AIC)[19]来完成,典型的AIC如图1所示.对信号的处理过程主要由信号调制、滤波和均匀采样三个步骤完成.首先用伪随机序列pc(t)对输入信号x(t)进行调制,它的取值必须满足奈奎斯特采样定理,调制的目的是为信号重构提供必要的随机性;然后用模拟低通滤波器h(t)对调制信号滤波,保证信号不失真;最后通过低速模数转换器(Analog to Digital Converter, ADC)采样得到压缩采样数据y[n].图1 AIC框架最直观的恢复信号的方法是在Ψ中找到信号x的稀疏表示,即找到系数向量Θ非零元素最少的解,也即求解复系数向量Θ的最小l0范数解:min‖Θ‖0s.t. y=ΦΨΘ,(3)‖·‖0表示向量中非零元素的个数.由于求解式(3)是一个NP-hard问题,常用的解决方法是用l1范数代替l0范数对复系数向量Θ进行约束,即min‖Θ‖1 s.t. y=ΦΨΘ.(4)在得到复系数向量Θ后,通过式(1)便可获得原始信号x.另外一种常用的算法是贪婪算法,典型的算法包括匹配追踪(Matching Pursuit, MP)算法和改进型的正交匹配追踪[13](Orthogonal Matching Pursuit, OMP)算法.2 LFM信号参数估计2.1 FrFT的定义从线性积分变换的角度给出FrFT的基本定义:函数f(t)的p阶FrFT是一个线性积分运算[6],且有fp(u)≡Kp(u,t)f(t)dt.(5)式中,Kp(u,t)=Aαexp[jπ(u2cot α-2utcsc α+t2cot α)],称为FrFT的核函数,是整数.当p=4n(α=2nπ)时,Kp(u,t)≡δ(u-t);当p=4n±2(α=(2n±1)π)时,Kp(u,t)≡δ(u+t).进一步地,为计算方便,经变量代换,fp(u)可以表示为fp(u) ={Fp[f(t)]}(u)=f(t)Kα(u,t)dt,0<|p|<2(0<|α|<π)(6)由式(6)看出,当p=1时,α=π/2,Aα=1,此时有f1(u)=e-j2πutf(t)dt.(7)可见f1(u)就是f(t)的普通Fourier变换.同样,可以看出f-1(u)是f(t)的普通Fourier 逆变换.2.2 基于FrFt的LFM信号参数估计单分量LFM信号x(t)可表示为x(t)=Aex p(j2πf0t+jπk0t2),(8)式中: A为信号幅度; k0为信号调频斜率; f0为信号起始频率; Δt为信号持续时间.将式(8)代入式(5),得Xα(u)= AAαexp{jπ[u2cot α-2(u csc α-f0)t+(cot α+k0)t2]}dt.(9)当cot α=-k0,即α=arccot(-k0)时,式(9)变为Xα(u)= AAαexp(jπu2cot α)·δ[2π(u csc α-f0)].(10)即一个LFM信号在旋转角度α=arccot(-k0)的FrFT域内表现为一个冲击函数,具有能量聚集的特性,该FrFT域被称为最佳FrFT域,相应的阶次p=α/(π/2)称为最佳FrFT阶次.对上述信号的调频斜率k0和起始频率f0的估计过程可以描述为[20]:(11)(12)实际工程中,需要计算离散FrFT,在经过量纲归一化[21]后信号的调频斜率k0和起始频率f0的估计式变为(13)式中: fs为采样频率; N为采样点数.3 基于CS的LFM信号参数估计高斯白噪声中LFM信号参数估计的信号模型为x=s+n.(14)式中: s为LFM信号; n为加性高斯白噪声.利用CS对LFM信号进行参数估计的目标是通过直接处理信号的压缩采样值,提取信号参数信息.由上节讨论可知,对于一个给定的LFM信号,存在一个最佳FrFT域对该LFM信号具有最好的能量聚集特性,即在最佳FrFT域中LFM信号表现为一个冲击函数,也就是说在最佳FrFT域中LFM信号是一个稀疏信号,且稀疏度K=1.这恰好满足CS理论对于信号在某个变换域是稀疏信号的要求,使利用CS理论在FrFT域内对LFM信号进行参数估计成为了可能.受上节基于FrFT的LFM信号参数估计算法启发,根据参数估计精度要求对FrFT的变换阶次p在[0,2)区间内进行搜索,分别构造各个阶次p对应的FrFT矩阵ΨFrFT_p作为信号的稀疏变换基,此时有s=ΨFrFT_pΘs_p.由于高斯白噪声在变换域内不稀疏[13],利用观测矩阵Φ对信号x压缩采样,得y=Φs+Φn=ΦΨFrFT_pΘs_p+Φn.(15)利用恢复算法求解信号在各个FrFT矩阵ΨFrFT_p中的系数向量Θs_p,即Θs(p,u).对Θs(p,u)进行二维搜索,此时式(11)变为(16)根据α=pπ/2得到旋转角度的估计值进而由式(13)得到信号调频斜率k0与起始频率f0的估计.实际中,LFM信号参数估计是一个复杂的问题,接收到的信号中常常伴随着噪声与窄带干扰的存在.此时信号模型为x=s+J+n.(17)式中: s为LFM信号; J为复正弦信号; n为加性高斯白噪声.对阶次p在[0,2)区间内进行搜索得到ΨFrFT_p,再利用观测矩阵Φ对信号x压缩采样,此时将得到y =Φ(s+J)+Φn=ΦΨFrFT_pΘ(s+J)-p+Φn.(18)由式(18)看出,通过恢复算法得到的信号系数向量Θ(s+J)_p中的元素将会受到干扰信号J的影响,因此通过对Θ(s+J)(p,u)进行二维搜索估计信号参数时,将会造成错误估计.在这里我们考虑强干扰的情况,即窄带干扰信号功率与信号功率相近或远大于信号的功率.正弦信号在Fourier变换域是稀疏信号,因此将文献[16]中的形态学成分分析的思想引入算法中.首先构造Fourier变换矩阵ΨFFT作为信号的稀疏变换基;然后利用OMP算法[13],控制算法迭代次数,迭代次数由干扰信号数目确定,进而得到信号x在ΨFFT中占主要成份的系数向量ΘJ,由于是在强干扰情况下,因此由OMP 算法得到的系数即为干扰信号J对应的成份,利用式(19)可以将干扰信号成份从y 中剔除:y=y-ΨFFTΘJ.(19)在剔除干扰成份后,对余下的信号利用恢复算法求解信号在FrFT矩阵ΨFrFT_p中的系数向量Θs-p,然后根据式(16)求得和进而根据α=pπ/2得到旋转角度由式(13)得到信号调频斜率k0与起始频率f0的估计.为了减少算法的计算量,可以对p采取两步搜索的策略,即粗搜索和精搜索.粗搜索是将p在[0,2)范围内进行间隔较大的搜索确定峰值的大致位置pcoarse;精搜索是在粗搜索得到的pcoarse附近进行间隔较小的搜索,得到峰值的精确位置pfine,进而得到角度综上所述,在高斯白噪声和强窄带干扰背景下,基于CS的LFM参数估计算法的步骤为:1) 获得信号压缩采样值y=Φx;2) 构造Fourier变换矩阵ΨFFT;3) 利用OMP算法进行有限次数迭代,求解信号在ΨFFT中的系数向量ΘFFT,利用式(19)去除窄带干扰信号;4) 以精度0.02构造阶次p∈[0,2)的FrFT矩阵5) 利用BP算法求解信号在中的系数向量即Θcoarse(p,u),进而得到6) 根据能量重心原理,在[pcoarse-0.01, pcoarse+0.01]范围内,以精度0.001构造FrFT矩阵7) 利用BP算法求解信号在中的系数向量即Θfine(p,u),然后通过式(16)搜索得到与8) 利用等式和式(13)计算得到信号调频斜率k0与起始频率f0的估计值.4 仿真实验及分析实验中,设信号长度N=512,幅度A=1,用fs=80 MHz的采样频率对信号进行模拟.利用快速Fourier变换(Fast Fourier Transform,FFT)算法直接构造N×N的Fourier变换矩阵ΨFFT作为窄带干扰信号的稀疏基.对信号的压缩采样过程是通过一个M×N维随机矩阵Φ完成的,其中Φ的元素满足高斯分布.规定检测成功率高于95%检测有意义.实验1 对单分量LFM信号参数估计进行数值实验.在信号中加入高斯白噪声和强窄带干扰信号,假设有一个复正弦信号作为干扰信号,信干比(Signal-to-Interference Ratio, SIR)为0 dB,频率f=20 MHz,此时‖ΘJ‖0=1,因此步骤3中OMP算法执行一次迭代.在SNR为5 dB的情况下,信号的调频斜率k0与起始频率f0分别取不同值时,考察参数估计的相对误差ER:(20)(21)信号参数估计结果如表1所示.表1 LFM信号参数估计结果k0/(Hz/s)^k0/(Hz/s)f0/Hz^f0/HzERk0ERf02e1220039e1215e614983e619e-312e-35e1250044e1210e610025e688e-425e-310e12100034e125e649847e634e-431e-3三组实验中,窄带干扰信号的频率均在LFM信号带宽内.由表1可以看出,在高斯白噪声和强窄带干扰条件下,当信号调频斜率与起始频率分别取不同值时,该算法的二次搜索过程可以准确估计信号参数,参数估计的相对误差能达到10-3量级,且算法不受干扰信号频带范围影响.实验2 考察文中算法在不同SNR条件下,参数估计成功概率随压缩采样点的变化情况.设输入信号起始频率f0=10 MHz,调频斜率k0=5×1012 Hz/s.在信号中加入窄带干扰信号,信号参数同实验1,SIR为0 dB.压缩采样点数M从40到140以步进20变化,在每个采样点下,进行200次实验,参数估计成功率为参数估计成功次数在200次实验中所占的比例,参数估计成功标准为估计值与实际值的相对误差小于0.01.在SNR为0、3、5和10 dB的情况下,分别考察信号参数估计的成功率.实验结果如图2所示.图2 不同SNR下参数估计成功概率随采样点数变化曲线由图2可以看出,信号参数估计成功概率随SNR和采样点数的增大而提高.当SNR 为3 dB时,采样点数M达到80,信号检测成功率即可以达到95%以上;而当SNR逐渐提高时,可以在更少的采样点数时,达到更高的检测成功率,获得更高的压缩比.由实验2可以看到,当SNR为3 dB,压缩采样点数低于奈奎斯特采样点数1/4时,信号参数估计成功率仍可以达到95%以上.实验3 考察文中算法在不同SIR条件下,参数估计成功概率随压缩采样点的变化情况.信号参数设置和参数估计成功概率准则同实验2.假设有两个复正弦信号作为干扰信号,频率分别取f1=15 MHz和f2=35 MHz,干扰信号和LFM信号频带重叠.此时满足‖ΘJ‖0=2,即干扰信号J在ΨFFT中是稀疏度为K=2的稀疏信号,因此步骤3中OMP算法执行两次迭代.在SNR为3 dB的情况下,SIR分别取0、-10和-20 dB,分别考察信号参数估计成功率随压缩采样点数变换情况.实验结果如图3所示. 由图3可以看出,在强干扰条件下,即当SIR远小于0 dB时,信号参数估计成功率基本不变.当压缩采样点数大于等于60时,信号参数估计成功率能够达到95%以上,当压缩采样点数小于60时,信号参数估计成功率迅速下降,同样是由于压缩采样点数过少时,采样值中包含原信号的信息太少,对恢复算法影响较大的原因.实验3说明文中算法对强窄带干扰不敏感,可以在强窄带干扰中准确估计LFM信号参数.图3 不同SIR下参数估计成功概率随采样点数变化曲线实验4 通过蒙特卡罗实验对文中算法与文献[17]算法进行对比实验,考察两种算法参数估计成功率随SNR和压缩采样点数M的变化情况.设信号的起始频率f0=10 MHz,带宽B=30 MHz,则此时对应的信号持续时间Δt=6.39 μs,调频斜率k0=4.696 7×1012 Hz/s.按照文献[17]算法,根据调频斜率k0和起始频率f0构造波形匹配字典,当压缩采样点数M=150,SNR从-10 dB到10 dB步进为2变化时,蒙特卡罗实验200次,两种算法信号参数估计成功率比较如图4所示.图4 参数估计成功概率随SNR变化曲线从图4可以看出:在受到噪声影响条件下,当检测成功率高于95%时,与文献[17]算法相比,由于LFM信号在FrFT域具有良好的稀疏性,相同参数估计成功概率下,文中算法允许的SNR低,抗噪声能力较强;在相同SNR条件下,文中算法能获得更高的参数估计成功概率.在SNR为-3 dB情况下,压缩采样点数从100到200步进为20变化,蒙特卡罗实验200次,考察两种算法参数估计成功率,结果如图5所示.图5 参数估计成功概率随采样点数变化曲线由图5可以看出:与文献[17]算法相比,在低SNR条件下,参数估计成功概率相同时,文中算法需要的压缩采样点数更少,具有更高的压缩率;在相同采样点数条件下,文中算法可以获得更高的参数估计成功概率.5 结论文中利用LFM信号在FrFT域具有稀疏性的特点,提出一种基于CS的LFM信号参数估计新方法.在没有LFM信号先验知识的情况下,该算法对高斯白噪声和强窄带干扰不敏感,可以对待估计LFM信号的调频斜率和起始频率准确估计.该方法将分数阶Fourier矩阵作为LFM信号的稀疏变换基,在达到相同参数估计成功率条件下,较波形匹配字典的方法需要的SNR更低、压缩采样点数更少.文中算法需要对变换阶次p搜索,计算量较大,如何寻找一种方法来降低搜索带来的计算量,还有待深入研究.参考文献[1] DJURIC P M, KAY S M. Parameter estimation of chirp signals[J]. IEEE Transactions on Acoustics, Speech and Signal Processing, 1990, 38(12): 2118-2126.[2] HAIMOVICH A M, PECKHAM C D, TETI Jr J G. SAR imagery of moving targets: application of time-frequency distributions for estimating motion parameters[C]//SPIE’s International Symposium on Opti cal Engineeringand Photonics in Aerospace Sensing. 1994: 238-247.[3] RAO P, TAYLOR F J. Estimation of instantaneous frequency using the discrete Wigner distribution[J]. Electronics letters, 1990, 26(4): 246-248. [4] BARBAROSSA S. Analysis of multicomponent LFM signals by a combined Wigner-Hough transform[J]. IEEE Transactions on Signal Processing, 1995, 43(6): 1511-1515.[5] WANG M, CHAN A K, CHUI C K. Linear frequency-modulated signal detection using Radon-ambiguity transform[J]. IEEE Transactions on Signal Processing, 1998, 46(3): 571-586.[6] SHEN Liran, YIN Qinngbo, LU Mingyu, et al. Linear FM signal parameter estimation using STFT and FRFT[J]. Chinese Journal of Electronics, 2013,22(2): 301-307.[7] 黄佑勇, 王激扬, 陈天麒. 基于欠采样的宽频段信号频率估计技术[J]. 电波科学学报, 2001, 16(2): 275-279.HUANG Youyong, WANG Jiyang, CHEN Tianqi. Wide-band frequency estimation with sub-Nyquist sampling[J]. Chinese Journal of Radio Science, 2001, 16(2): 275-279. (in Chinese)[8] 沈显祥, 叶瑞青, 唐斌, 等. 基于欠采样的宽带线性调频信号参数估计[J]. 电波科学学报, 2007, 22(1): 43-46.SHEN Xianxiang, YE Ruiqing, TANG Bin, et al. An algorithm for estimation of wideband LFM signal parameters based on subsampling[J]. Chinese Journal of Radio Science, 2007, 22(1): 43-46. (in Chinese)[9] DONOHO D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289-1306.[10] CANDS E J, ROMBERG J, TAO T. Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information[J]. IEEE Transactions on Information Theory, 2006, 52(2): 489-509.[11] 焦李成, 杨淑媛, 刘芳, 等. 压缩感知回顾与展望[J]. 电子学报, 2011, 39(7): 1651-1662.JIAO Licheng, YANG Shyuan, LIU Fang, et al. Development and prospect of compressive sensing[J]. Acta Electronica Sinica, 2011, 39(7): 1651-1662. [12] DUARTE M F, DAVENPORT M A, WAKIN M B, et al. Sparse signal detection from incoherent projections[C]//IEEE International Conference on Acoustics, Speech and Signal Processing. Toulouse,IEEE, 2006, 3: 305-308.[13] 刘冰, 付平, 孟升卫. 基于正交匹配追踪的压缩感知信号检测算法[J]. 仪器仪表学报, 2010,31(9): 1959-1964.LIU Bing, FU Ping, MENG Shengwei. Compressive sensing signal detection algorithm based on orthogonal matching pursuit[J]. Chinese Journal of Scientific Instrument, 2010, 31(9): 1959-1964. (in Chinese)[14] XU J, PI Y, CAO Z. UWB LFM echo signal detection and time-delay estimation based on compressive sensing[C]//2010 IEEE 10th International Conference on Signal Processing (ICSP). IEEE, 2010: 2435-2438.[15] 黄传禄, 晁坤, 毛云志. 基于压缩感知的空间谱估计[J]. 电波科学学报, 2014, 29(1): 150-157.HUANG Chuanlu, CHAO Kun, MAO Yunzhi. The spatial spectrum estimation based on compressive sensing[J]. Chinese Journal of Radio Science, 2014, 29(1): 150-157. (in Chinese)[16] SHI G, LIN J, CHEN X, et al. UWB echo signal detection with ultra-low rate sampling based on compressed sensing[J]. IEEE Transactions on Circuits and Systems II: Express Briefs, 2008, 55(4): 379-383.[17] LIU B, FU P, XU C, et al. Parameter estimation of LFM signal with compressive measurements[J]. Journal of Convergence Information Technology, 2011, 6(3): 303-310.[18] ZHA Song, LIU Peiguo, HUANG Jijun. Parameter estimation of LFM signal via compressive sensing[C]//IET International Radar Conference. IET, 2013: 1-5.[19] TROPP J A, LASKA J N, DUARTE M F, et al. Beyond Nyquist: efficient sampling of sparse bandlimited signals[J]. IEEE Transactions on Information Theory, 2010, 56(1): 520-544.[20] QI L, TAO R, ZHOU S, et al. Detection and parameter estimation of multicomponent LFM signal based on the fractional Fourier transform[J]. Science in China Series F: Information Sciences, 2004, 47(2): 184-198. [21] 赵兴浩, 邓兵, 陶然. 分数阶傅里叶变换数值计算中的量纲归一化[J]. 北京理工大学学报, 2005, 25(4): 360-364.ZHAO Xinghao, DENG Bing, TAO Ran. Dimensional normalization in the digital computation of the fractional Fourier transform[J]. Transactions of Beijing Institute Technology, 2005, 25(4): 360-364.。
12-第十二讲-频率响应和奈奎斯特图
tan
1
1
2 n
2
2 n
0
M 1
0
Im Re
1/2
1
T
M
1
2
90
M 0 180
g. GHs
1
2 n
s
1
2 n
s2
Im
M
s 1 1 e j
s 1 s 1
j
1
e j tan1
2 2 1
根据传递函数绘制奈奎斯特曲线
GH s
1
ss 1
G1sG2 s
GH j G1 j G2 j M1e j1 M 2e j2 M1 M 2 e j2
??s?????90111??????nj?msghnnn????s????????????????????????1???????????1????????2n22?22n222n224?12?11nmssgh????f12?reim????????????2n211tann????1?00???m?????????2??1901mt??0180?????m?????s??????????????????????????????????1????????????????????22n212n22?22n222n112tan4?12?1nnmssgh????????g1imre??1?00???m?????2??901???mt???180???m???12222100011tan100901111tan1090119009090kkkktsttsttt???????????????????传递函数幅值m相位?渐近值00mm????表
基于压缩感知的数据压缩理论及其重构算法对比研究
基于压缩感知的数据压缩理论及其重构算法对比研究李逸川;于峻川;徐红燕【摘要】The compressive sensing is a new theory ,which attracted wide attentions from the world once it was proposed ,In this study ,we evaluated the advantages and disadvantages among various reconstruction algorithm through theoretical summary and simulation experiment ,aiming at providing theoretical support for the research and application in the future ,First of all ,this thesis systematically summarized the theoretical framework and the main components of Compressive Sensing ,and then carried out one-dimensional and two-dimensional simulation experiment by using OMP reconstruction algorithm , The result shows that the Compressive Sensing algorithm can reconstruct the original signal in a high probability even under a low sampling rate ,In the case of the sampling rate is 0 .5 ,the compressing rate was achieved to 53%-60% ,Finally ,in order to estimate the time spent during reconstructing and reconstruction precision of various reconstruction algorithm ,we carried out another simulation experiment through standard test image based on the brief summary of the characteristics of these algorithms ,The study shows that IRLS algorithm can provide a higher reconstructionaccuracy ,while the GPSR algorithm costs the minimum time to reconstruct the image .%压缩感知作为一种全新的信号采样理论 ,一经提出便引起广泛关注.本文拟在前人研究的基础上 ,通过理论研究及仿真实验对常见重构算法进行评价为后续理论研究及应用提供科学依据.首先对压缩感知的理论基础和主要构成进行阐述 ,以贪婪算法中的OM P算法分别对一维信号及不同类型的二维图像进行仿真实验 ,实验结果表明压缩感知算法可以在较低采样率下实现对一维或二维信号的高效重构 ,在采样率在0 .5的情况下 ,其数据的压缩率达53% ~60% .在系统总结几种常见重构算法特点的基础上 ,以标准测试影像为对象构建仿真实验 ,分别从重构算法的运算效率和重构质量两个方面对实验结果进行评价 ,结果显示IRLS算法重构精度较高 ,而GPSR算法的运算耗时较短.【期刊名称】《中国矿业》【年(卷),期】2015(024)012【总页数】6页(P159-164)【关键词】压缩感知;图像压缩;重构算法;稀疏性【作者】李逸川;于峻川;徐红燕【作者单位】中国国土资源航空物探遥感中心 ,北京100083;地质过程与矿产资源国家重点实验室,中国地质大学 (北京) 地球科学与资源学院,北京100083;中国国土资源航空物探遥感中心 ,北京100083;中国地质调查局地学文献中心 ,北京100083【正文语种】中文【中图分类】TN911.72信号采样是我们获取数字信息的第一手段,而传统的采样理论要求采样率必须为信号带宽的两倍。
压缩采样介绍
压缩采样介绍一个与传统数据采集不同的传感、采样范例Emmanuel J. Candès and Michael B. WakinIEEE信号处理杂志2008年3月对信号和图像采样的传统方法遵循著名的香农定理:采样速率必须至少是当前信号最大频率的两倍(也称奈奎斯特速率)。
事实上,这个定理几乎是所有信号采集方式的基础,被大规模应用在消费性音频视频电子设备、医学成像设备、收音机等等。
(对于某些信号,例如并非是原始的天然影像,采样率就不是由香农定理所规定,而是取决于需要的空间或时间分辨率。
然而,这些系统经常在采样前使用反锯齿低通滤波器保持信号,使得香农定理扮演一个幕后的角色。
)例如在数据转换领域,标准模数转换器(ADC)技术使用量化香农定理:信号一律以大于或等于奈奎斯特速率来采样。
这篇论文研究了压缩采样的理论,该理论亦被称为压缩传感或CS,是一项新颖的与传统数据采集不同的传感、采样技术。
CS理论声称可以比传统方法使用更少的采样信号和测量来恢复特定信号和图形。
为了实现这一点,CS依靠两个准则:稀疏性,不连贯性。
◆稀疏性,表示连续时间信号的信息率可能比由带宽决定的还要小,或者离散时间信号自由度远小于其有限的长度。
更准确的说,CS方法发现许多自然信号是十分稀疏并可压缩的,当用适当的基础Ψ表述时,他们就能有简明的表达。
◆传感形式不连贯延续了时间和频率的二元性,并表明目标在Ψ有一个稀少的特征并一定会在他们已得的范围内扩展,正如时域里面的狄拉克或冲击信号在频域也可展开表示。
另外,不连贯表明它与信号特征不同,采样、传感的波形在Ψ有一个非常密集的表示。
最重要的是我们可以设计有效的传感、采样规则,来获取有用的信息并将其嵌入到稀疏的信号中,加入到一小段数据中。
这些规则仅十分简单的要求把信号同一些固定的波形相关联,这些波形也都没什么根据。
关于这些采样规则,最不同寻常的是它们可以使用一个传感器,从稀疏的信号中高效的获取信息,并且也不需要分析这些信号。
奈奎斯特频率总结
奈奎斯特频率总结奈奎斯特频率(Nyquist frequency)是信号处理领域中一个重要的概念,它是在采样过程中最高能够准确表示的频率。
奈奎斯特频率的理论基础是奈奎斯特采样定理(Nyquist sampling theorem),它提供了在通信和数字信号处理中正确采样信号的最低要求。
奈奎斯特采样定理奈奎斯特采样定理是由美国电气工程师哈里·S·奈奎斯特(Harry Nyquist)在20世纪20年代提出的。
该定理指出,若要准确还原一个信号,需要对其进行至少两倍于信号最高频率的采样频率。
具体来说,若信号的最高频率为ff,则采样频率必须大于等于2ff才能避免出现混叠(aliasing)现象。
混叠指的是高于奈奎斯特频率的频率成分被错误地还原到低于奈奎斯特频率的频率范围内,在数字信号中表现为频谱的重叠。
奈奎斯特采样定理为我们提供了正确采样信号的准则,保证了信号还原的准确性和完整性。
当采样频率满足奈奎斯特定理的要求时,我们可以通过恢复滤波器(reconstruction filter)对采样信号进行重建,还原出原始信号。
奈奎斯特频率的计算在实际应用中,我们需要根据信号的特性来计算奈奎斯特频率。
一般而言,信号的奈奎斯特频率是信号最高频率的两倍。
例如,对于一个音频信号,它的频率范围一般在20 Hz至20 kHz之间。
因此,它的奈奎斯特频率为40 kHz。
为了能够准确采样并还原这个音频信号,我们需要以至少40 kHz的采样频率进行采样。
同样地,对于图像信号,奈奎斯特频率是图像中最高频率分量的两倍。
如果一幅图像中包含最高频率为100线对/英寸的细节,那么它的奈奎斯特频率将是200线对/英寸。
在实际应用中,我们通常会选择稍高于奈奎斯特频率的采样频率,以避免混叠和重建误差。
因此,在设计采样系统时,我们需要仔细确认信号的奈奎斯特频率,并选择合适的采样频率。
奈奎斯特频率和信号重建奈奎斯特频率不仅在信号采样中起到重要作用,也在信号重建中扮演关键角色。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
奈奎斯特采样频率与压缩感知比较报告学生张**年级2010级班级0210** 班学号021012**专业电子信息工程学院电子工程学院西安电子科技大学2013年5月压缩感知与奈奎斯特采样频率比较报告张**摘要:经典的采样定理认为,不失真的恢复模拟信号,采样频率应该不小于奈奎斯特采样频率(模拟信号最高频率的两倍)。
但是这种方法在使采集到的数据有很大的冗杂性。
Dohono等人提出的压缩感知理论运用了大部分信号可以在预知的一组基上面稀疏表示的原理,利用随机投影实现了在低于奈奎斯特采样频率下实现了信号的采集。
本文介绍了压缩感知的一些基本理论以及,并将其与香农采样定理进行了比较。
最后讨论了压缩感知的一些信息获取算法以及压缩感知理论的应用前景。
关键词:香农采样定理奈奎斯特采样频率压缩感知引言当前大部分数据采集系统都是基于传统的香农采样定理来设计,按照这种方式采集的数据能够充分表示原始信号,但是它们存在较大的冗余。
因此,这些方法往往导致采集数据的泛滥和传感器的浪费。
研究如何根据信号的一些特征来实现低于奈奎斯特采样频率的采集,以减少所需采集的数据量具有重要的意义。
起源于对具有有限新息率信号(即单位时间内具有有限自由度的信号)进行采集的研究,利用固定的结构性基函数以两倍于新息率而不是两倍于奈奎斯特采样频率对连续信号进行采集,Donoho等人提出的压缩感知方法则提供一种可以广泛应用于可压缩信号的采集方法。
该方法所需要的传感器数目大大减少,采集到的数据也具有更小的冗余度。
因此,该理论提出后立即吸引了众多科学家的关注,目前我国关于压缩感知方法的研究也已经开始起步,相信不久将有更多的人加入到关于压缩感知的研究行列。
压缩感知采集方法并不是对数据直接进行采集,而是通过一组特定波形去感知信号,即将信号投影到给定波形上面(衡量与给定波形的相关度),感知到一组压缩数据。
最后利用最优化的方法实现对压缩数据解密,估计出原始信号的重要信息。
压缩感知关键的问题是如何给定用来感知信号的波形才能有效地恢复出原始信号的重要信息。
涉及的关键因素在于给定的波形要与可以用来压缩原始信号的波形组均不相干,并且不相干程度越高,感知数据包含的信息量越大,为准确获取重建原始信号所需的感知数据量就越少。
第一章奈奎斯特采样原理奈奎斯特频率是离散信号系统采样频率的一半,因奈奎斯特-香农采样定理得名。
采样定理指出,只要离散系统的奈奎斯特频率高于被采样信号的最高频率或带宽,就可以避免混叠现象。
从理论上说,即使奈奎斯特频率恰好大于信号带宽,也足以通过信号的采样重建原信号。
但是,重建信号的过程需要以一个低通滤波器或者带通滤波器将在奈奎斯特频率之上的高频分量全部滤除,同时还要保证原信号中频率在奈奎斯特频率以下的分量不发生畸变,而这是不可能实现的。
在实际应用中,为了保证抗混叠滤波器的性能,接近奈奎斯特频率的分量在采样和信号重建的过程中可能会发生畸变。
因此信号带宽通常会略小于奈奎斯特频率,具体的情况要看所使用的滤波器的性能。
需要注意的是,奈奎斯特频率必须严格大于信号包含的最高频率。
如果信号中包含的最高频率恰好为奈奎斯特频率,那么在这个频率分量上的采样会因为相位模糊而有无穷多种该频率的正弦波对应于离散采样,因此不足以重建为原来的连续时间信号。
第二章压缩感知理论本节对压缩感知做一个简单介绍。
待采集信号f只在k个时刻非零(k即为稀疏程度)。
为采集f中的信息,将其投影到了给定的一组感知波形φ上(也可以说,用一组给定的波形f进行感知),得到了一组远小于信号原始长度的测度数y=y fφ压缩感知方法的目的是通过远小于采集信号数据量的测度y 恢复出原始信号f 的全部信息。
从式中求解f 是一个欠定的问题,但从另一方面看,信号只有k 个未知位置上的未知变量,即信号只有1k +个自由度,因此,在一定条件下当测度数只要超过信号自由度的时候就是可以通过一些非线性的方法进行恢复。
显然,当任意选取感知矩阵φ的1k +列均线性无关的时候(即感知矩阵能够区分任意两个不同的均为k 稀疏的信号),在满足y f φ=的所有情况中找到的具有最稀疏特性的信号即为所求,即求解如下的最优化问题:min (0)..f P s t y f φ⎧⎨=⎩其中f 代表f 的零范数,即其中非零元素的个数。
此外给定一些其它约束条件以后,信号也可以通过求解如下的最优化问题来实现:min..f s t y f φ⎧⎨=⎩ 其中01p <≤,()1p p p f f =∑然而对(0)P 问题的求解只能通过对所有可能的稀疏情况进行求解后才能找到最稀疏的形式,这是一个NP 难问题.而对01p <<时式所示最优化问题的求解也存在一定困难。
幸运的是,用1p =时上式的解来估计可压缩信号,即接近最优,并且对它的求解能够等价为一个线性规划问题,从而便于利用现有方法进行求解。
当f 具有可压缩性时,上述的方法同样可以估计出f 中的具有较大幅值的项。
在实际应用中,待采集的信号f (如一副图像)本身通常不是稀疏的,但是在某个基ψ上的变换系数是稀疏的,或者是可压缩的,即由f ψθ=确定的变换系数θ是稀疏的,或者极少的系数包含了几乎全部的能量。
结合上式可以得到感知数据与变换系数之间的关系为y f φψϕθ==若令a φψ= (测度矩阵),则y A θ=由于上面两式有相同的形式和假设,因而同样可以按照上述的最优化方法估计出式中的变换系数,进而经过变换估计出待采集信号。
实际上f 本身稀疏可以看成是,在单位阵上稀疏,此时的测度矩阵即为感知矩阵。
当感知矩阵与变换矩阵不相干时,矩阵A 就能够很好地满足能够通过极小化变换系数的1l 范数来估计原始信号的约束条件。
并且随机矩阵即与任何一个固定的变换矩阵不相干的概率非常大。
也即通过收集信号随机投影系数来采集未知信号是最优的,这样的一个数据采集方式在实际应用中显然是比较容易实现的。
另外,信号的稀疏程度越高,利用上式恢复出原始信号所需的感知数据量也将越少。
第三章 奈奎斯特采样频率与压缩感知对比下面将通过一个信号采集的实例来说明压缩感知方法的优越性。
下图给出了一个周期为0.01T s =的待采集信号()f t 。
在实际工作中,需要将连续信号()f t 进行抽样以便于在计算机上实现处理,该抽样过程可以看作是原始信号和脉冲串函数的乘积:()()()f nT f t p t =其中,()()f T T nT σ=-∑,()t σ为脉冲函数。
根据香农一奈奎斯特采样定理,如果该信号为连续信号且抽样频率两倍于信号的最高截止频率(奈奎斯特采样频率),那么,可由采集信号()f nT 恢复出原始信号。
根据信号的连续性和最高截止频率有界,可以得到信号的重建方法为: sin(()/)()()()/t nT T f t f nT t nT Tππ-=-∑即利用sin c 函数插值出未采集时刻的信息。
当该信号为非带限信号或者采样频率低于奈奎斯特采样频率时重建出的信号较原始信号将存在一定的误差。
如上图所示当我们对以采样速率为5KHz 从上图(a )中获取的采样数据进行两倍的上采样时,按照重建方法获取到的估计信号与原始信号存在一定的误差。
结果如上图(b )所示与传统重建方法所依据的连续性和有限带宽相比,压缩感知重建方法利用的是信号在某个域上具有稀疏性和可压缩性,通过给定信号更好的表征域(换言之,利用一些其它的先验知识),通过同样的数据可以得到更加好的重建效果。
上图(c )结合该信号的变化是稀疏的(即只在部分时刻存在变化)这一先验知识,利用最优化方法重构的结果,该方法完全重构出了原始信号。
另一方面在压缩感知中对数据的获取方面也采用了更加灵活的方式以利于感知信号中包含的重要信息,传统的利用脉冲串函数进行采集的方法可以看作是压缩感知中感知波形为脉冲函数时的一种特殊情况。
由于采用压缩感知重建方法所需的数据量依赖于感知波形和表征波形之间的不相干性和信号的稀疏程度,而采用随机信号感知的方法能够确保与任给的一个表征域具有良好的不相干性,因此在对包含一些高频分量的信号采集时,若已知信号在某个表征域稀疏或可压缩,采用压缩感知的方法将具有更大的优越性。
实际上大部分信号均具有可压缩性,因此,压缩感知方法具有更广泛的优越性和普适性,具有广阔的应用前景。
第四章 压缩感知信息获取方法这一章主要介绍几个从感知到的数据中估计原始信号的几种常见实用方法:基追踪算法、贪婪算法● 1 基追踪算法首先需要指出的是基追踪算法并不是一个最优化原则。
其原理是上述讨论的给定一些限制条件后,通过极小化1l 范数可以获得最稀疏的解。
等价的标准线性规划问题为min....,0v s t U b ααα⎧⎨=≥⎩其中,α是重新定义的一个新变量,2m n ⇔;(,)U A A ⇔-;b y ⇔上述线性规划问题可以通过单纯形法、内点法或对数障碍法来进行求解。
● 2 贪婪算法贪婪算法的基本步骤为1。
给定初始的估计00θ=;2。
在每次迭代中根据00()A A A θθθθ-=-确定出0θθ-的估计值∆ 3。
仅保留∆中较大值的项,将其它项置为零。
更新00θθ=+∆。
关于贪婪算法的具体实现方法有很多种,主要有正交贪婪算法、规整化正交贪婪算法、分段正交贪婪算法和梯度贪婪算法等。
● 3 其他方法极小化1l 范数的方法能够有效解决压缩感知中的恢复问题,但是当结合其它的一些先验知识后,该问题可以被更加有效地解决。
在此,我们仅简单介绍贝叶斯压缩感知方法(,和基于模型的压缩感知方法。
Ji 等人提出的BCS 借助传统的贝叶斯方法和机器学习中的主动学习方法,通过将关于稀疏性的先验信息用垂直先验分布来建模,提出了自适应的感知方法以及相应的恢复方法。
而Baraniuk 等人提出的针对基于模型可压缩信号的压缩感知方法中利用小波树模型和块稀疏模型,仅需要与稀疏程度相当的测度数即可实现信号的鲁棒性恢复。
、第五章压缩感知技术的应用从压缩感知的采样机制可以看出:同一个投影方法(硬件结构)可以完成大部分可压缩信号的采集;每个采集到的数据(测度数)包含的信息量是相同的,即CS方法对测度数的丢失鲁棒性比较强;感知数据的获取方法仅需通过简单的投影来完成,需要大量运算的估计方法放在拥有更加强大计算能力的解码端进行;此外,由于信号是通过随机投影的方式获取的,感知到的数据具有一定的保密性。
由于压缩感知方法具有上述优点,可被广泛应用于以下领域:(1)数据采集莱斯大学的研究组利用压缩感知原理成功设计出了单像素相机。
该相机利用透镜和数字可控镜片阵列随机地将从物体发射来的部分光线汇聚到一点,从而利用一个像素即可实现感知数据的采集。
尽管该相机还存在一些问题,但是它为未来相机的发展提供了一种新的思路和发展方向。
此外压缩感知方法被应用到了多跳网络和无线传感器网络中具有稀疏性或可压缩性网络数据的重构上面。