专利名称：：频率估计的多段采样信号融合处理方法
技术领域：
：本发明涉及信号频率估计，具体的说，涉及多段采样信号的处理方法。技术背景多段釆样信号指的是，对同一被测对象，在不同时间段或不同采样参数下，通过数次釆样过程，得到了数段采样信号。各次被采样对象在该次釆样过程中的频率不变，各段采样所得信号(以下简称为"釆样信号")的频率之间的相对关系已知，相邻段采样之间的时间间隔可以未知，在相邻次采样之间的时间间隔内被采样信号可以发生未知变化，各次采样过程的釆样频率可以不同。多段釆样信号尤其是多段同频信号处理问题十分常见，例如，在一些测量设备的校准过程中，需要对已知对象进行多次测量，然后综合评估设备的工作状况；还例如，在LFMCW雷达测距过程中，多个回波信号经过延迟-混频、低通滤波后，表现为同一单频信号的多次釆样；在电子对抗过程中，对相雷达进行数字侦收时，也将获得多段(多个脉冲)釆样信号。众所周知，相比于单段信号，多段同频信号应含有更多的有用信息。但一直以来对多段同频信号进行深度信息融合式的综合处理技术不受重视。其原因在于一方面，该融合处理过程会产生较大的运算量，当时运算处理器的性能相对不足；另一方面，以往的研究中，多将待处理信号近似为某一段时间内的平稳信号且背景噪声的幅度相对较低，从而得到采样持续时间较长的平稳信号，即单段采样数据中包含有足量的、易于提取的信息，已能满足当时的测量要求。在此情况下，虽然多段采样信号十分常见，但其融合处理技术没有得到应有的发展。近年来，非平稳信号、低信噪比信号的研究日益受到重视，由于被测频率保持相对稳定的时间非常短暂，相当于用一个持续时间很短的窗函数去乘该信号，且采样结果中包含大量噪声，这给信号处理带来了很大的困难。例如，为达到测量精度要求，需要一段信噪比为5必、持续时间超过0.1毫秒的平稳信号，而实际采样过程只能得到4段信噪比为-5d5、持续时间仅为0.05毫秒的平稳信号。为此，除了进一步改进现有的单段采样信号处理方法和提高硬件性能外，另一种重要的途径就是尽可能利用多段釆样信号中包含的信息。近几年来，多段釆样信号处理技术逐渐得到关注，已有研究人员开展了初步研究，目前已有的处理方法可分为三种。(1)直接累加法对各段釆样信号分别进行如DTFT等频域分析，然后将分析结果累加。该方法计算量小，能够应用于各种情况，但是精度较低。(2)相位关联法一种对多个相位相参的分段信号进行频谱估计的方法，利用各段信号之间的相位关系，通过相位联系对由单段信号得到的频谱估计做进一步的修正。该方法对多段采样信号提出了较多的限制条件，且推导过程中采用了近似处理较粗略，难以实际应用。(3)相位积累法包括直接相位积累法和旋转相位积累法，将分段采样信号进行相位关联，以实现相参积累和提高频率估计精度。直接相位积累法运算量大，并要求各段采样之间的空闲间隔时间已知且在该空闲间隔内被测信号保持稳定不变，应用范围较小。旋转相位积累法是目前较先进的一种方法，但是受噪声影响大，且要求各段采样信号的初相已知，严重制约了该方法的应用与推广。综上所述，对多段采样信号处理方法的研究具有重要意义，本专利要解决的技术问题是提出一种高精度方法，具有普遍适用性、能够处理各种类型的多段采样信号，以较小的实时计算量获得较高的精度。发明的内容由于本专利方法是线性的，具有可加性，因此分析的结果同样适用于多段多频信号；由于本专利方法能够应用于多段不同频信号，因此分析的结果同样适用于多段不同采样频率情况;为便于叙述，本文均以被测信号为单频信号、釆样频率相同为例。1.三类多段釆样信号本专利方法将现有的多段采样信号阿题划分为三类(1)在噪声和干扰背景下，对被测信号进行重复釆样，被测信号的频率和各段釆样信号的初相未知。(2)在噪声和干扰背景下，对被测信号进行多次采样得到多段采样信号，各次釆样所得信号的频率、初相未知且频率并不相等，但各次釆样所得信号频率之间的频率差已知。(3)在噪声和干扰背景下，对被测信号进行多次釆样得到多段采样信号，各次采样所得信号的频率、初相未知且频率并不相等，但各次采样所得信号频率之间的倍数(可为任意小数)已知。将以上三类情况依次简称为多段同频信号、多段降频信号和多段分频信号，其中后两类可统称为多段变频信号处理。可以看出，本专利方法将多段釆样信号处理研究从单纯的重复测量即多段同频信号处理拓展到非同频信号处理。这是因为在多段变频信号中，虽然各段信号的频率各不相同，但是由于它们的频率之间的相对关系已知，所以本质上仍然可以将它们视作同频信号的重复测量。这种新型的重复测量关系是本专利方法的一个特色，将本专利方法的应用领域从平稳信号扩展到非平稳信号，扩展了本专利方法的应用范围。2.通用处理流程由于本专利方法对以上三类多段采样信号进行频率估计的构思基本相同，为便于说明本专利方法的核心内容，特归纳出本专利方法的通用处理流程如附图所示，并给出具体说明如下(为便于叙述，以下数学表达式，若未经特别说明，则均参照MATLAB7函数规定解释)设有^/段采样信号，首段信号频率为力，各段信号频率间的相对关系已知，釆样频率为y;.,各段初相为代w),各段长度为"(m)，前(w-1)段的长度和为5(w)。(1)根据首段采样信号频率的估计值/2和最大相对估计误差C确定首段待比较频率序列5F0(k)，;te[l，"/2](2)根据各段釆样信号之间的相对关系，得出各段釆样信号的待比较频率序列5F(附,k)。设Z)/(m)为各段信号降频幅度，p(附)为各段信号分频比，则5FO，A:)=-￡)/(/)]。(3)将各段信号进行丑尸(m，k)频率处的频域分析，得到Z。(w，A:)。(4)根据z。(附，w利用直接累加法得到首段釆样信号频率的直接累加法估计值y;。。(5)设2'sin^[,—^^)]//s},其中sin{/7i[胸(/)-綱糊/X/￡>M}sin(/!'[5F0(/)—S尸0(/t)]/X}调出预存的的极大值表s("，y;(o，w，该表为函数i^c/;(o，;a)的巧个极大值序列，"e[l,nj。设频率估计区间Gi^[g/j，gf2]--512[",_/;(/),/〗|。取&("，W中较大值对应的力序列&(/)，取S2("，y;(f)，/i)中Ae/io(/)部分，得到S4("，y;(OA(/》。因为一般情况下，、-i即可满足精度要求，为便于说明，取《!=1,即&(",yi(a^o(/))-S4(y;(r)，w/》。.(7)设&的长度为",，确定搜索频率序列SF0(/)-5F0(/^)),/e[l，"y;](8)部分情况下可以直接根据经验给出代m)的估计值《(/w)。如果发现该经验值的误差过大，例如《(加)-e(/w)〉/7/，其中《O0,代7w)e，则可以通过近似计算得到对应于不同信号段的初相估计值《(/n),方法如下>找到^>a)中最接近,。的频率点5F(7/a。(/))。>得到《(w)=a"g/e[Z。(m，&。(/))],其中a"g/e(f)表示复数f的角度。在接下来的说明中以较复杂的《(w)为例，若可以采用《(w)计算，则只霜作相应简化即可，不再赘述。(9)根据《(m)、SFO(it)和SFO(/)确定相参补偿因子矩阵-(/，Ara)。其中，针对多段同频信号和多段降频信号，有附，A:)=《(/)-《(1)-2[S尸0(/)—SFO(A:)]/乂；针对多段分频信号，有J"w,A:)=《(附)一6>2(1)—2/w';5(w)[S尸0(/)-SFO(A;)]/乂/5(/w)=s鹏(/D(1:w-1))(10)计算如下<formula>formulaseeoriginaldocumentpage6</formula>(13)若/3位于频率区间(尸的左端处即/3-^;^厶(A根据实际情况取值)，则左移频率区间GF和相应修改&，然后重复上述流程中的第(12)步。若/;位于频率区间GF的右端处^/2-/3^厶(A根据实际情况取值)，则右移频率区间GF和相应修改&，然后重复上述流程中的第(12)步。(14)重复上述流程中的第(13)步，直至/3不再位于频率区间GF的左右端，该频率/3即为本专利方法的频率估计结果。3.理论证明为更好说明本专利方法的优越性，特给出理论证明过程如下(1)针对多段同頻信号的理论证明凡是对同一平稳信号进行多次测量的场合都会产生多段同频信号，其一般形式如下设有W段同频信号，信号频率为y;，釆样频率为乂，各段长度为Z)(m),各段初相为e(/t7),即<formula>formulaseeoriginaldocumentpage6</formula><formula>formulaseeoriginaldocumentpage7</formula>证毕。(2)针对多段降類信号的理论证明产生多段降频信号的情况也比较多，例如LFMCW雷达测距过程中，对同一目标使用不同延时长度的延迟线组合则可以得到多段降频信号，其一般形式如下设有#段降频信号，首段信号频率为/i，各段信号降频幅度为Z)/(w),采样频率为,，各段长度为Z)(w)，各段初相为代m)，即=0，D/(1)=O，fle[l，D(w)],w€[l，M]<=D(附)，e[2，M+1]x0[乖)+a]=cos刚+2W[_/;-D/(附)]("-1)//s}根据上述通用处理流程，可知<formula>formulaseeoriginaldocumentpage7</formula>证毕。(3)针对多段分類信号的理论证明产生多段分频信号的情况也比较多，例如LFMCW雷达测距过程中，对同一目标使用不同电压调频斜率的压控振荡器则可以得到多段分频信号，其一般形式如下设有^段分频信号，各段信号的频率为/(加)^/K附)，y;，;7(l)y，p(m+l)^p(m),釆样频率为乂.，各段长度为Z)(附)，各段初相为代/w)，即5(1)=0,ae[1,_D(/m)]，we[1，M]50)=UD(附),se[2,+1]+fl]=cos[e(w)+2-y;("—1)/乂]根据上述通用处理流程，可知g2[/-WOW]//,}m='expL/讽/n)--(/，附，"+;;'L/;-WO(O][Z)(w)+1]//J]荟:".^、"V、i/人ex力(代1)+"[SF0(/)-5F0(A:)]{2"in{"[/i-5F0(A:)]//J、一(/)+;7乖)+1]//"/.[/;-柳(/)].[;^(/)+1]/,}〗_sin{/n'-gFO(A:)]/^}sin{//[S尸0(/)-丑尸O(A:)]/,M}eXp[>^(i)+b5(M+i)+i]，p〃[SF0(/)—5F0W]/y;+(^+i),/7/.L/;—wo(/)]//s}]证毕。4.实验证明本专利方法给出了理论证明，而且通过大量实验对本专利方法的精度和适应性进行了验证，这一点在后文"具体实施方式"中的两个应用实例的实验结果分析中得到了体现。图本专利方法的具体处理流程具体实施方式为进一步说明本专利方法的处理流程，特给出如下应用实例。1.多段同频信号处理的一个实例及实验结果分析设有四段同频信号，已知各段采样信号的频率相同，首段釆样信号的频率的初始估计值/2二3.15M/fe，其最大相对估计误差为0.1,即首段釆样信号频率y;e[2.863636MHz，3.5MHz]。已知釆样信号的频率为8MHz，各段采样信号的长度依次为[50444844]。.为精确测量首段采样信号频率，根据本专利方法，可釆用如下处理过程(1)确定首段待比较频率序列丑FO=linspace(2.863636MHz,3.5MHz，207)。设定该序列长度为207以减小后续的Chirp—Z变换的运算量。(2)确定各段采样信号的待比较频率序列BF(/n,k)=APO(k)。(3)将各段同频信号进行5F0(k)频率处的Chirp-Z变换，得到Z。(w,W。(4)根据Zp(附，矢)利用求和取平均的方法得到首段釆样信号频率的直接累加法估计值乂o。计算过程如下<formula>formulaseeoriginaldocumentpage9</formula>兩,=8一飼,.万}(5)设wl刈—、ys,，其中d二47，sin{/i[腦(/)-講糊/y>￡)M}"sin{/7"[SF0(/)-5F0(it)〗，5//s},(0e[2.863636MHz,3.5MHz]，；ae[1，207]。调出预存的的极大值表S(y;(0，/0,该表为函数Fsc/;(a/a)的最大值序列。取该函数表中y;we(o.99、y;。,i.oo"y;。)部分，得到&(/(0,/0。(6)计算&(/0=1;|52[/;(/+1)，/1]一52[1/;(/),/2]1。取&(w中最大的io个值所对应的力序列&(/)，取&c/;(o，w中^/^(/)部分，得到&(y;(0A。)。(7)确定搜索频率序列<formula>formulaseeoriginaldocumentpage9</formula>(8)计算<formula>formulaseeoriginaldocumentpage9</formula>(9)确定<formula>formulaseeoriginaldocumentpage9</formula>(10)计算如下<formula>formulaseeoriginaldocumentpage9</formula>(11)取Z2(/，;t)中的谱峰频率位置Z3(/)-5F0(M4X[Z2(/,"])。'Z4"(,),/HZ3(/)—s4"(OA(/))l(12)根据Z3和&计算得到频率/3。j兵，取Zs的最小值对应y;(/。)，得到/3=/;(。)；(13)若/3位于频率区间G尸的左端处即/3-_g/;^0.01M//z,则左移频率区间GF和相应修改&，然后重复上述流程中的第(12)步。若,位于频率区间c^的右端处gf2_/3S0.01Mfe,则右移频率区间GF和相应修改&,然后重复上述流程中的第(12)步。(14)重复上述流程中的第(13)步，直至《不再位于频率区间GF的左右端，该频率/3即为本专利方法的频率估计结果。接下来，结合该例，对本专利方法和目前适应性最好的直接累加法、处理精度最高的旋转相位积累法(为满足旋转相位积累法的特殊应用条件，专为该方法提供其所需特别信息，其他两种方法不使用这些特别信息)进行了全面对比，给出了定量结果。仿真实验中使用的参数设置如下信噪比^、^=-5^/万；噪声类型为高斯白噪声；各段采样信号的初相均满足幅度为2，p/的高斯分布。(1)抗噪性仿真该仿真过程包括500次实验，每次实验的参数除随机函数和信噪比外完全相同，每进行50次实验后信噪比改变一次，得到实验结果如表1所示。表1中的信噪比SNR单位为d5,标准差单位为幻&。可以看出，在多段同频信号、相同信噪比情况下，本专利方法的精度大约提高到其他方法的4~5倍，表1多段同频信号情况下三种方法的抗噪性,比标\\SNR\^\-2.5-3.1-3.6-4.2-4.7-5.3-S.8-6.4-6.9-7.5标准差平均值直接象加法l加91114l幼卯143146148l卯124旋转梱位私累法78抑卯105108102114卯110122咖1227181721182737274125(2)初始估计值/2的相对误差对处理结果的影响仿真该仿真过程包括500次实验，每次实验的参数除随机函数和初始估计值外完全相同，每进行50次实验后初始估计值改变一次，得到实验结果如表2所示。表2中的初始估计值/2的相对误差。=(/2-,)/,，其单位为千分之一，标准差单位为iOfe。可以看出，在多段同频信号、相同初始估计值慵况下，本专利方法的精度大约提髙到其他方法的3~4倍。表2多段同频信号情况下初始估计误差对三种方法的精度影响对比<table>tableseeoriginaldocumentpage11</column></row><table>(3)对不同首段采样信号频率y;的适应性仿真该仿真过程包括500次实验，每次实验的参数除随机函数和首段采样信号频率外完全相同，每进行50次实验后首段采样信号频率改变一次，得到实验结果如表3所示。表3中的首段釆样信号频率_/;的单位为^//&,标准差单位为尤wz。可以看出，在多段同频信号、相同首段釆样信号频率情况下，本专利方法的精度大约提髙到其他方法的4~5倍表3多段同頻信号情况下首段采样信号频率对三种方法的精度影响对比<table>tableseeoriginaldocumentpage11</column></row><table>(4)对不同采样频率y;的适应性该仿真过程包括500次实验，每次实验的参数除随机函数和采样频率外完全相同，每进行50次实验后釆样频率改变一次，得到实验结果如表4所示。表4中的采样频率,的单位为Mife，标准差单位为《/fe。可以看出，在多段同频信号、相同采样频率情况下，本专利方法的精度大约提高到其他方法的4~5倍。表4多段同频信号情况下采样频率对三种方法的精度影响对比<table>tableseeoriginaldocumentpage11</column></row><table>(5)对不同单段信号长度平均值万的适应性该仿真过程包括500次实验，每次实验的参数除随机函数和单段信号长度平均值外完全相同，每进行50次实验后单段信号长度平均值改变一次，得到实验结果如表5所示。表5中的单段信号长度平均值5的单位为点，标准差单位为夂/fe。可以看出，在多段同频信号、相同单段信号长度平均值情况下，本专利方法的精度大约提高到其他方法的4~5倍。表5多段同频信号情况下单段信号长度平均值对三种方法的精度影响对比<table>tableseeoriginaldocumentpage12</column></row><table>(6)对不同信号总段数的适应性该仿真过程包括500次实验，每次实验的参数除随机函数和信号总段数外完全相同，每进行50次实验后信号总段数改变一次，得到实验结果如表5所示。表5中的信号总段数M的单位为段，标准差单位为《/fe。可以看出，在多段同频信号、相同信号总段数情况下，本专利方法的精度大约提高到其他方法的45倍。表6多段同频信号情况下信号总段数对三种方法的精度影响对比<table>tableseeoriginaldocumentpage12</column></row><table>(7)实时计算量对比直接累加法实时计算量"A/.("log2(丄)+丄+l)，其中丄》"/2+max(/)(m))。旋转相位积累法实时计算量"^，万*"/2本专利方法实时计算量4*M*(i:.log2(￡)+Z+l)+4*(M+w/>M"/2)+"*M"_/;*0.02/0.0005,其中，",一般取30即可；"表示重复上述流程中第(13)步骤的次数，一般取2即可；0.02表示上述流程中第(13)步骤中的频率搜索区域的长度；0.0005表示频率搜索中的的频率步进幅度。在本例多段同频信号的一般处理情况下，三种方法的实时计算量依次为36880、16S600和184256，本专利方法的实时计算量增幅很小。2.多段分频信号处理的一个实例及实验结果分析设有四段分频信号，已知首段釆样信号的频率的初始估计值/2=3.15^/&，其最大估计误差为0.1,即首段采样信号频率yje[2.863636MHz,3.5MHz]。已知各段采样信号频率与首本发明提供的用于大型集装箱辐射检查系统的主被动结合扫描方法包括以下步骤在扫描通道中设置手动控制按钮，将手动控制按钮和检测车内的控制系统相关联；集装箱卡车司机将卡车开到指定待检位置；司机自主选择按下通知按钮通知搡作人员启动检测的被动扫5描方式或者选择离开检测区域等待检测的主动扫描方式；辐射检查系统扫描检测待检集装箱后，司机驾车离开。其中，主动扫描时，扫描车能够沿扫描通道方向作前后平移运动，以便对待检车辆进行扫描；被动扫描时，集装箱卡车低速通过扫描通道，进行扫描。io其中，手动控制按钮被设置于一支架上，支架设置在竖探测器臂前方，辐射危险区域之外。当待检车辆驾驶室与所述手动控制按钮支架处于垂直于扫描通道的同一直线位置时，司机处在安全位置而待检集装箱处于扫描区本发明的另一目的是提供一种用于大型集装箱辐射检查系统的15主被动结合扫描装置，包括辐射检查系统和手动控制按钮支架两部分，辐射检查系统和手动控制按钮支架通过数据信号线电连接。进一步，在手动控制按钮支架上装有用于输入启动扫描的信号的手动控制按钮。进一步，辐射检查系统包括扫描车，用于根据控制信号对待检20车辆进行扫描；扫描通道，用于使待检车辆通过扫描区域；其中所述扫描车能够沿扫描通道方向作前后平移运动，以便对待检车辆进行扫描。优选的，手动控制按钮支架安装在所述竖探测器臂远离扫描区域的一侧的前方，辐射危险区域之外，当待检车辆驾驶室与所述手动控25制按钮支架处于垂直于扫描通道的同一直线位置时，司机处在安全位置而待检集装箱处于扫描区。扫描车内具有控制室，所述手动控制按钮产生的启动信号传输至所书控制室，从而提示控制室内搡作人员发(9)确定<formula>formulaseeoriginaldocumentpage14</formula>(10)计算如下<formula>formulaseeoriginaldocumentpage14</formula>(11)取22(/，6)中的谱峰频率位置<formula>formulaseeoriginaldocumentpage14</formula>(12)根据23和&计算得到频率/3<formula>formulaseeoriginaldocumentpage14</formula>取^的最小值对应,(/。)，得到/3-,(。；(13)若,位于频率区间GF的左端处即/3-g/^0.01M/fe,则左移频率区间GF和相应修改&,然后重复上述流程中的第(12)步。若厶位于频率区间GF的右端处即#2-/3《0.01M/fe，则右移频率区间GF和相应修改&，然后重复上述流程中的第(12)步(14)重复上述流程中的第(13)步，直至/3不再位于频率区间(尸的左右端。接下来，结合该例，对本专利方法和目前适应性最好的直接累加法进行了全面的精度对比，给出了定量结果。仿真实验中使用的参数设置如下噪声类型为髙斯白噪声，信噪比SA^:-5^5;各段采样信号的初相均满足幅度为2，;7/的髙斯分布。(1)抗噪性仿真该仿真过程包括500次实验，每次实验的参数除随机函数和信噪比外完全相同，每进行50次实验后信噪比改变一次，得到实验结果如表7所示。表7中的信噪比SNR单位为标准差单位为K/fe。可以看出，在多段分频信号、相同信噪比情况下，本专利方法大约能将精度提高到现有方法的6倍。表7多段分频信号情况下两种方法的抗噪性能对比<table>tableseeoriginaldocumentpage14</column></row><table>(2)初始估计值/2的相对误差对处理结果的影响仿真该仿真过程包括500次实验，每次实验的参数除随机函数和初始估计值外完全相同，每进行50次实验后初始估计值改变一次，得到实验结果如表8所示。表8中的初始估计值/2的相对误差c-(y;-其单位为千分之一，标准差单位为K/fe。可以看出，在多段分频信号、相同初始估计值情况下，本专利方法大约能将精度提高到现有方法的6倍。表8多段分频信号情况下初始估计误差对两种方法的精度影响对比<table>tableseeoriginaldocumentpage15</column></row><table>(3)对不同首段采样信号频率乂的适应性仿真该仿真过程包括500次实验，每次实验的参数除随机函数和首段采样信号频率外完全相同，每进行50次实验后首段采样信号频率改变一次，得到实验结果如表9所示。表9中的首段釆样信号频率乂的单位为M/fe，标准差单位为/C/fe。可以看出，在多段分频信号、相同首段采样信号频率情况下，本专利方法大约能将精度提髙到现有方法的6倍。表9多段分親信号情况下首段采样信号频率对两种方法的精度影响对比<table>tableseeoriginaldocumentpage15</column></row><table>(4)对不同采样频率义的适应性该仿真过程包括500次实验，每次实验的参数除随机函数和采样频率外完全相同，每进行50次实验后采样频率改变一次，得到实验结果如表10所示。表IO中的釆样频率X的单位为M/fe,标准差单位为K/fe。可以看出，在多段分频信号、相同釆样频率情况下，本专利方法大约能将精度提髙到现有方法的6倍。表10多段同频信号情况下采样频率对两种方法的精度影响对比<table>tableseeoriginaldocumentpage15</column></row><table>(5)对不同单段信号长度平均值D的适应性该仿真过程包括500次实验，每次实验的参数除随机函数和单段信号长度平均值外完全相同，每进行50次实验后单段信号长度平均值改变一次，得到实验结果如表11所示。表ll中的单段信号长度平均值万的单位为点，标准差单位为iOfe。可以看出，在多段分频信号、相同单段信号长度平均值情况下，本专利方法大约能将精度提高到现有方法的6~7倍。表11多段分频信号情况下单段信号长度平均值对两种方法的精度影响对比<table>tableseeoriginaldocumentpage16</column></row><table>(6)对不同信号总段数的适应性该仿真过程包括500次实验，每次实验的参数除随机函数和信号总段数外完全相同，每进行50次实验后信号总段数改变一次，得到实验结果如表12所示。表12中的信号总段数M的单位为段，标准差单位为尺/fe。可以看出，在多段分频信号、相同信号总段数情况下，本专利方法大约能将精度提髙到现有方法的6倍。表12多段分频信号情况下信号总段数对两种方法的精度影响对比<table>tableseeoriginaldocumentpage16</column></row><table>(7)实时计算量对比直接累加法实时计算量4争M"丄，log2(丄)+丄+l)，其中丄2w/2+max(D(w))。旋转相位积累法实时计算量4*M*D*w/2本专利方法实时计算量4*A/*(￡.log2(Z)+Z+l)+4，(M+",*il/"/2)+"A/"y>0.02/0.0005,其中，",一般取30即可；"表示重复上述流程中第(13)步骤的次数，一般取2即可；0.02表示上述流程中第(13)步骤中的频率搜索区域的长度；0.0005表示频率搜索中的频率步进幅度。在本例多段分频信号的一般处理情况下，三种方法的实时计算量依次为36880、165600和184256,本专利方法的实时计算量增幅很小。权利要求1.一种信号频率估计方法，其特征在于，给出频域分析参数矩阵、相参补偿因子序列和搜索频率序列的计算公式，实现对多段同频信号、多段降频信号和多段分频信号的频率估计。全文摘要多段采样信号融合处理是一种适用于信号频率估计的频域分析方法，通过为各段采样信号设定不同的频域分析参数和相参补偿因子，本方法能补偿由于各段被测信号之间的频率不等、相位不连续对信号处理结果的影响，实现相关信息融合和有用信息提取，以获得更高的频率估计精度。本方法的应用范围广泛，适用于各段采样信号的初相未知、相邻采样段之间的时间间隔未知、在相邻采样段之间的时间间隔内被采样信号发生未知变化、各次采样过程的采样频率不同等情况，能够实现对多段不同频信号的高精度处理，满足多段同频信号、多段降频信号、多段分频信号等情况下的频率估计需要。本方法给出了理论推导过程和验证实验结果与对比，进行了实时计算量分析。文档编号G01R23/02GK101231315SQ200710078138公开日2008年7月30日申请日期2007年1月24日优先权日2007年1月24日发明者刘良兵,涂亚庆申请人:涂亚庆;刘良兵