一种适用于非连续可导模拟信号的采样与重构方法与流程

本发明涉及信号采样
技术领域：
，特别涉及一种适用于非连续可导模拟信号的采样与重构的时域方法。
背景技术：
：在这个数字时代，把一个模拟信号先转换成数字信号再进行处理或存储等工作几乎成了一个标准过程。在这一转换过程中信号采样是基本的一个步骤。到目前为止用来确定采样率的基本理论是经典的针对带限信号的香农采样理论(例如文献[1],[2])。该理论的核心内容可以描述如下：对任意一个带宽为w的信号x(t)，如果采样率fs至少为2w，那么x(t)就可以由其采样点和sinc函数完美地重构出来。在进行采样器设计时香农理论可以看成一个基于频域的工具，因为其分析基础是基于信号的频谱。在过去几十年中沿着这一思路的跟进工作有许多(例如文献[1],[3]–[6])。针对某些类特殊信号的采样方法也被提出(例如文献[7]–[10])。香农理论虽然在理论上简洁漂亮，在实际设计采样器的过程中却可能带来不便或困难。首先，所有需要处理的实际模拟信号在时间上都是有限的，那么这些信号就必然是非带限信号[11]。所以任何频谱的截断就必然带来采样后的信号频谱混叠现象，从而产生误差。实际上，通常以下四种误差都需要考虑：频谱混叠误差，幅度误差，截断误差，以及时间抖动误差，如文献[12],[13]。这些误差分析也被用在小波分析上，如文献[14]。傅立叶分析中涉及到的著名的gibbs现象可能产生较大的瞬时误差，如文献[15]。随着高速dsp的应用越来越广，人们经常希望在保证时域逐点(pointwise)重构误差的前提下采样尽量少的数据点。香农理论在这方面没有提供什么支持。虽然有方法被提出来用于在保证重构误差的前提下尽量减少采样点的个数(例如文献[16])，目前还是缺乏适用于工程应用的简单有效的方法。在一些应用场合，我们并不能提前得到信号的频谱信息，也就不能直接应用香农理论了。并且，重构模拟信号的电路常常使用常数或线性插值，如文献[17]，但是在香农理论中用于插值的sinc函数在实际中是不能完全得到的。虽然香农理论中用到的能量稳定性在推导中方便使用，但是在工程实践中逐点(pointwise)的稳定性是更常期望的，尤其在研究信号瞬态特征的时候，如文献[3]。虽然之前有基于时域的采样方法被提出，如文献[18]，但此类方法都要求被采样模拟信号连续可导，从而使用范围有限。技术实现要素：本发明提出一种模拟信号采样与重构的方法，是一个纯粹基于时域的方法，能够保证逐点最大重构误差的前提下进行有效的均匀采样。本方法不但对连续可导的模拟信号适用，而且对所有连续，并且除了有限个点外连续一阶可导的模拟信号都适用。本发明所采用的技术方案为：一种适用于非连续可导模拟信号的采样与重构方法，要求连续模拟信号x(t)在[0,1]区间上除了有限个点ts,0＝0＜ts,1＜...＜ts,k-1＜ts,k＝1之外连续一阶可导；待采样连续模拟信号x(t)经过一次微分电路后，对输出信号通过模拟电路求得其幅值的最大值，用以计算采样周期t，计算得到的采样周期t能够保证在每相邻两个采样点之间最多包含一个不可导的点，并以该采样周期t对延时后的待采样模拟信号x(t)以及该延时后的待采样连续模拟信号x(t)的一阶导数进行采样，在重构时，利用各个采样值组成一个分段线性函数来近似待采样模拟信号x(t)。所述采样周期设置为式中：x(t)在点ts,k不可导，k＝0,1,...,k；ε为最大允许重构误差，单位与待采样连续模拟信号x(t)的单位一致；η2为待采样模拟信号x(t)的二阶导数幅值的最大值。利用所述采样周期t进行采样，信号重构方法：对任意t∈[nt,(n+1)t]，设置y1(t)＝x(nt)+x'(nt)(t-nt)y2(t)＝x((n+1)t)+x'((n+1)t)(t-(n+1)t)得到方程的解区间[t3,t4]。如果[t3,t4]＝[nt,(n+1)t]，设置xr(t)＝y1(t)或y2(t)否则设置或t0'取为[t3,t4]区间的任意一点，以及与传统的香农采样理论相比，本发明的有益效果：1)易于分析与实现。本文提出的方法允许采样器的设计完全在时域进行。如果信号最大的变化率x”(t)已知，可以很容易地计算出有效的采样周期。如果不知，如图3所示，可以由模拟电路自动获得。但是在香农采样理论中，信号带宽无法很容易地通过电路估计得到。2)保证每一个点的重构误差精度。在经典的香农采样理论中，这个通常在实践中所需要的功能几乎是缺失的。但是在本文所提出的方法中，可以设计保证每个点的最大重构误差。如具体实施部分论述所示，本发明重构误差范围很严谨。但是在香农采样理论中时域重构误差范围很难被估计，而且可能随着信号的变化而变化较大。3)适用于非带限和非连续可导模拟信号。所有的物理信号都是时间有限的，因此它们必然是非带限信号。所以在香农采样理论的实际应用中总是涉及到近似，而且也很难去分析每个点的时域重构误差。相比之下，本发明所提出的方法能够很容易地处理非带限和非连续可导模拟信号，并且重构误差可控。而且它也能够有效地对连续分段线性信号进行完美采样重构。4)适用于未知特征的采样信号。在经典的香农采样理论中，在采样之前我们需要知道信号的带宽。但本发明提出的方法的实现结构不需要采样信号的任何信息，因为相关参数都可以在采样电路中得到。附图说明图1中a、b依次为本发明采用的两种待采样模拟信号x(t)示意图；图2中a、b依次为两种待采样连续模拟信号x(t)经过香农采样理论进行采样后的重构误差仿真结果示意图；图3为本发明的电路原理图，其中延时td与计算t的时间相当；图4中a、b依次为本发明实施例1的两种待采样模拟信号x(t)经过采样后的重构误差仿真结果示意图。具体实施方式下面结合附图对本发明及其效果进一步说明。我们先简单回顾香农采样理论。给定一个带宽为w的待采样模拟信号x(t)，先定义一个采样冲击序列信号其中采样周期t在下面确定，从而保证完美重构。s(t)的频谱可以写成其中采样率fs＝1/t。然后得到的采样序列就是xs(t)的频谱然后就可以写成其中*是卷积运算。公式(4)告诉我们为了得到完美重构，我们需要以下两个条件：1)fs≥2w(5)2)将xs(t)通过一个理想低通滤波器th(t)其中并且w＜fc＜fs-w(7)满足以上两个条件将使得重构信号xr(t)与x(t)具有相同的频谱。那么我们就得到在l2范数下xr(t)＝x(t)。滤波器的时域冲击响应h(t)是h(t)＝2fcsinc(2fct)(8)所以，利用公式(3)和(8)，xr(t)的解析表达式就是公式(9)是利用x(t)的采样点进行重构的公式。我们将以上的香农采样与重构方法归纳为算法af。在图1中我们列出两个测试信号，并且将算法af的重构效果列在图2中。图1中的两个信号全部定义在[0,1]区间上，并且依次为如下形式：x1(t)＝g(2t)+2g(2t-1)(10)其中为了有效地采样(采集尽可能少的数据点)，我们希望在保证期望的重构信号精度的前提下使用尽可能大的采样周期t。第一个困难就是估计信号的带宽w。即使有了信号的表达式，我们通常在艰苦的推导后才得到如下的频谱：其中是1[0,1](t)的频谱，并且我们使用了小-o的记号，意思是注意这两个信号中没有一个的带宽w是有限的。所以我们就把带宽w定义为频谱从此以后衰减的足够小的频率点，即对于所有f>w，本文中我们固定εf＝10-3。然后我们就可以得到信号xk(t)的带宽wk如下：为了满足公式(5)和(7)，我们选择fs＝2.1w以及fc＝fs/2。现在我们可以根据公式(9)进行信号采样与重构。但是在这之前，我们并不能对逐点的重构误差有什么保证，因为香农理论是通过l2范数建立起来的，或者如文献[3]描述的，香农理论只保证能量稳定性。如何将频域的误差上界转化为时域的误差还是一个尚未解决的问题。事实上，在4g和5g通信中广为研究的峰均功率比(papr)问题中这是一个核心难点(例如文献[19])。在本发明中我们提出一个基于时域的方法来回答如下问题：给定一个有限长度的待采样非连续可导模拟信号x(t)以及一个最大的逐点重构误差上界ε，如何有效地对x(t)进行均匀采样，使得最大的逐点重构误差不超过ε？为此，本发明采用的方法是：将待采样连续模拟信号x(t)经过两次模拟一阶微分电路后，对输出信号通过模拟电路求得其幅值的最大值，用以与ts一同计算采样周期t，并以该采样周期t对延时后的待采样模拟信号x(t)以及该延时后的待采样模拟信号x(t)的一阶导数进行采样，在重构时，利用各个采样值组成分段线性函数来近似待采样模拟信号x(t)。下面通过具体的实施例对上述方法进行验证说明。一种适用于非连续可导模拟信号的采样与重构方法，待采样连续模拟信号x(t)经过两次一阶微分电路后，对输出信号通过模拟电路求得其幅值的最大值η2，用以与ts一同计算采样周期并以该采样周期t对延时后的待采样模拟信号x(t)以及该延时后的待采样模拟信号x(t)的一阶导数进行采样，在重构时，利用各个采样值组成一个分段线性函数来近似待采样连续模拟信号x(t)。式中：x(t)在点ts,k不可导，k＝0,1,...,k；ε为最大允许重构误差，单位与待采样连续模拟信号x(t)的单位一致；η2为待采样连续模拟信号x(t)的二阶导数幅值或绝对值的最大值。利用所述采样周期t进行采样，信号重构方法：对任意t∈[nt,(n+1)t]，设置y1(t)＝x(nt)+x'(nt)(t-nt)y2(t)＝x((n+1)t)+x'((n+1)t)(t-(n+1)t)得到方程的解区间[t3,t4]。如果[t3,t4]＝[nt,(n+1)t]，设置xr(t)＝y1(t)或y2(t)否则设置或t0'取为[t3,t4]区间的任意一点，以及依据香农理论进行信号的采样与重构并不能对逐点的重构误差有什么保证，如附图2所示，两个信号的最大重构误差都在0.01数量级，本发明能够使得最大信号重构误差完全满足设计要求，如附图4所示，本发明的重构方法产生的重构误差均小于设计要求的10-3，特别地，对于信号x1(t)这样的分段线性信号，最大重构误差是0，也就是说，达到了完美重构。不但如此，如表i所示，本发明的方法所需的采样周期比香农理论所确定的采样周期高出许多，说明本发明所需的采样点数更少。表i两个测试信号由不同算法得到的采样周期talg.t1t2香农理论0.0135140.000598本发明0.1250.00151976一，证明如下：根据假设，采样周期t满足t＜ts。这将使得任一采样区间至多含有一个不可导的点。对任意采样区间[nt,(n+1)t)，没有不可导的点的情况很容易处理。所以我们处理其他情况，即在[nt，(n+1)t]区间上x(t)在t0不可导，但在[nt，(n+1)t]区间上除t0外的其他部分连续可导。在区间[nt,t0]上，将x(t)表示为在区间t∈[t0,(n+1)t]上，将x(t)表示为在区间[nt,(n+1)t]上定义函数y1(t)＝x(nt)+x'(nt)(t-nt)(48)y2(t)＝x((n+1)t)+x'((n+1)t)(t-(n+1)t)(49)如果t∈[nt,t0]，我们可以推导得到如果t∈[t0,(n+1)t]，我们可以类似得到在t＝t0处，由(50)和(51)我们可以得到事实上，对任意t满足(52)，都有|yi(t)-x(t)≤ε，i＝1,2(53)因为,例如，如果t∈[nt,t0]，则(50)式对于y1(t)成立，并且(53)对于i＝1也成立；式(50)和(52)表明(53)对于i＝2也有效。因为y1(t)和y2(t)都是线性函数，所以(52)式可以被表示为如下形式其中a和b均为常数。如果a＝0，那么(52)式对所有的t∈[nt,(n+1)t]均有效，(52)式之后的分析表明无论y1(t)或者y2(t)用来表示xr(t)均能满足|xr(t)-x(t)≤ε。如果a≠0，在区间nt和(n+1)t上解(54)式，可以得到t∈[t3,t4](55)其中其中sgn(·)是符号函数。既然t0∈[t3,t4]，所以对于任意t0'∈[t3,t4]都会满足|xr(t)-x(t)≤ε。特别地，为了进一步减少重构误差，可以取二，具体实现算法0)输入x(t)，t∈[0,1]，时域最大误差允许值ε和ts。1)估计2)设置3)获取采样点x(nt)和x'(nt)。4)对任意t∈[nt,(n+1)t)，设置y1(t)＝x(nt)+x'(nt)(t-nt)y2(t)＝x((n+1)t)+x'((n+1)t)(t-(n+1)t)得到方程的解区间[t3,t4]。如果[t3,t4]＝[nt,(n+1)t]，设置xr(t)＝y1(t)(59)否则设置和以下为本发明中所涉及的现有文献：[1]abdulj.jerri,―theshannonsamplingtheorem—itsvariousextensionsandapplications:atutorialreview,”proceedingsoftheieee,vol.65,no.11,pp.1565–1596,1977.[2]johng.proakisanddimitrisg.manolakis,digitalsignalprocessing:principles,algorithms,andapplications,4thed.,prentice-hall,2006.[3]p.p.vaidyanathan,―generalizationsofthesamplingtheorem:sevendecadesafternyquist,”ieeetransactionsoncircuitsandsystemsi:fundamentaltheoryandapplications,vol.48,no.9,pp.1094–1109,2001.[4]ahmedizayed,advancesinshannon’ssamplingtheory,crcpress,1993.[5]robertjiimarks,advancedtopicsinshannonsamplingandinterpolationtheory,springerscience&businessmedia,2012.[6]michaelunser,―sampling—50yearsaftershannon,”proceedingsoftheieee,vol.88,no.4,pp.569–587,2000.[7]rodneygvaughan,neillscott,anddrodwhite,―thetheoryofbandpasssampling,”ieeetransactionsonsignalprocessing,vol.39,no.9,pp.1973–1984,1991.[8]raymondboute,―thegeometryofbandpasssampling:asimpleandsafeapproach[lecturenotes],”ieeesignalprocessingmagazine,vol.29,no.4,pp.90–96,2012.[9]jasondmcewen,gillespuy,jean-philippethiran,pierrevandergheynst,dimitrivandeville,andyveswiaux,―sparsesignalreconstructiononthesphere:implicationsofanewsamplingtheorem,”ieeetransactionsonimageprocessing,vol.22,no.6,pp.2275–2285,2013.[10]hiromiuedaandtoshinoritsuboi,―asamplingtheoremforperiodicfunctionswithnominusfrequencycomponentanditsapplication,”incommunications(apcc),201319thasia-pacificconferenceon.ieee,2013,pp.225–230.[11]stephanemallat,awavelettourofsignalprocessing,2nded.,academicpress,1999.[12]georgecstey,―upperboundsontimejitterandsamplingrateerrors,”inieeproceedingsg-electroniccircuitsandsystems.iet,1983,vol.130(5),pp.210–212.[13]jingfanlong,peixinye,andxiuhuayuan,―truncationerrorandaliasingerrorforwhittaker-shannonsamplingexpansion,”incontrolconference(ccc),201130thchinese.ieee,2011,pp.2983–2985.[14]wenchangsunandxingweizhou,―samplingtheoremforwaveletsubspaces:errorestimateandirregularsampling,”ieeetransactionsonsignalprocessing,vol.48,no.1,pp.223–226,2000.[15]holgerbocheandullrichjmoenich,―reconstructionbehaviorofshannonsamplingserieswithoversampling-fundamentallimits,”insourceandchannelcoding(scc),20087thinternationalitgconferenceon.vde,2008,pp.1–6.[16]zhanjiesong,beiliu,yanweipang,chunpinghou,andxuelongli,―animprovednyquist–shannonirregularsamplingtheoremfromlocalaverages,”ieeetransactionsoninformationtheory,vol.58,no.9,pp.6093–6100,2012.[17]chung-hsunhuangandchao-yangchang,―anareaandpowerefficientadder-basedstepwiselinearinterpolationfordigitalsignalprocessing,”ieeetransactionsonconsumerelectronics,vol.62,no.1,pp.69–75,2016.[18]philipeluftandtimoilaakso,―adaptivecontrolofsamplingrateusingalocaltime-domainsamplingtheorem,”incircuitsandsystems,1994.iscas’94.,1994ieeeinternationalsymposiumon.ieee,1994.[19]taojiangandyiyanwu,“anoverview:peak-to-averagepowerratioreductiontechniquesforofdmsignals,”ieeetransactionsonbroadcasting,vol.54,no.2,pp.257–268,2008.当前第1页12