《一种基于多观测向量稀疏表示的复正弦波频率估计方法.pdf》由会员分享，可在线阅读，更多相关《一种基于多观测向量稀疏表示的复正弦波频率估计方法.pdf（12页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104166795A43申请公布日20141126CN104166795A21申请号201410395022222申请日20140812G06F19/0020110171申请人电子科技大学地址611731四川省成都市高新区（西区）西源大道2006号72发明人费晓超张浩罗晓宇甘露廖红舒74专利代理机构成都宏顺专利代理事务所普通合伙51227代理人李玉兴54发明名称一种基于多观测向量稀疏表示的复正弦波频率估计方法57摘要本发明属于信号处理领域，尤其涉及复正弦波频率估计。本发明提供一种高精度复正弦波频率估计方法。由于复正弦波信号在频域具有稀疏性，基于观测数据稀疏表示的模型，首先针对。
2、网格失配的情况，通过一阶泰勒展开来逼近真实的频率向量，得到多观测向量下的稀疏表示逼近模型。然后在该模型的基础上用正交匹配追踪算法计算得到网格上离真实频率值最近的点，并且利用最小二乘方法计算频率的修正值。最后得到修正后复正弦波频率的估计值。利用本发明的估计方法，修正后的估计性能有明显的提高，即使在粗糙的网格，也有很高的估计精度。本发明的估计方法与其他估计方法相比，估计性能有明显的提高。51INTCL权利要求书2页说明书6页附图3页19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书6页附图3页10申请公布号CN104166795ACN104166795A1/2页21一种基于多观测。
3、向量稀疏表示的复正弦波频率估计方法，其特征在于，包括如下步骤S1、由长度为P的复正弦信号XI加白噪声得到接收数据其中，I1，P，K为复正弦波信号的个数，为第K个信号的复幅度，FK0，1为第K个信号的归一化频率，K1，2，K，初始相位是在0，2均匀分布的随机变量，VI是零均值、方差为的高斯白噪声，当IK时，和相互独立，S2、S1所述接收数据YI的向量形式为YAFS0V，其中，频率向量幅度向量S3、将S2所述接收数据YI的向量形式Y写成稀疏表示的模型S4、在S3所述稀疏表示的模型Y的基础上，用一阶泰勒展开来逼近真实的频率向量其中，存在第K个真实频率且是网格上距离FK最近的点，是的一阶导数，S5、根。
4、据S4所述AFK构造过完备字典记DIAG，当时，则当则N0，因此，观测模型可以写成S6、利用S5所述观测模型Y通过正交匹配追踪算法和最小二乘方法求解出真实频率值的估计值，具体为S61、用正交匹配追踪算法求解出网格上离真实频率值最近的点，具体如下S611、初始化记迭代次数T1，残差矩阵RY，支撑集其中0NL表示维数为NL的全零矩阵；S612、计算CAHR，ARGMAX|CI，|2，其中H表示矩阵的共轭转置，|CI，|2表示矩阵C第I行的L2范数，ARGMAX表示最大值的索引值；S62、用最小二乘方法计算S和Q在支撑集上元素的估计值和并更新残差矩阵更新迭代次数TT1；权利要求书CN10416679。
5、5A2/2页3S63、当迭代次数TK时，返回步骤S612，依次循环迭代；一旦TK，停止迭代，输出和把和的对应元素相除，得到维数为KL的矩阵求矩阵每一行的平均值，得到向量在支撑集上元素的估计值S64、设S63所述在网格上的索引值为N1，NK，则修正后的频率估计值为2根据权利要求1所述的一种基于多观测向量稀疏表示的复正弦波频率估计方法，其特征在于S3所述将接收数据YI的向量形式Y写成稀疏表示的模型具体如下S31、由S2所述接收数据YI，得到信号向量从而得到数据矩阵其中，LM1P，M为信号向量的长度，L为信号向量的个数；S32、把S1所述归一化频率FK在0，1的范围上过完备化为一个离散的网格N代表网。
6、格数且满足NK；S33、将S31所述信号向量YI写成多观测向量的稀疏表示形式其中，过完备矩阵频率向量是一个K行稀疏矩阵，即S中只有K行元素为非零，其他行的元素都为零，记所述K行元素的索引值为支撑集，权利要求书CN104166795A1/6页4一种基于多观测向量稀疏表示的复正弦波频率估计方法技术领域0001本发明属于信号处理领域，尤其涉及复正弦波频率估计。背景技术0002正弦波频率估计算法广泛应用于雷达、通信以及电子对抗等信号处理领域。其中经典的算法有多重信号分类MULTIPLESIGNALCLASSICATION，MUSIC算法、旋转不变子空间ESTIMATIONOFSIGNALPARAMET。
7、ERSVIAROTATIONALINVARIANCETECHNIQUE，ESPRIT算法等子空间类算法和最大似然估计类算法MAXIMUMLIKELIHOOD，ML等。然而，基于子空间理论的频率估计方法虽然实现了超分辨侧向，但是一旦接收数据不足或信噪比比较低时，这类方法不能有效地区分信号子空间和噪声子空间，其性能会急剧下降。而最大似然估计类算法由于要进行复杂的多维搜索而不具有实用性。0003近年来，随着压缩感知和稀疏表示技术的发展，许多基于稀疏表示的频率估计方法被提出。最具代表性的为L1SVD算法，它利用L1范数来重构稀疏信号，并且在多观测向量的条件下通过奇异值分解SINGULARVALUEDE。
8、COMPOSITION，SVD来减小数据矩阵的规模以及降低噪声的影响。然而，这些基于稀疏表示的方法通常假设频率位于离散化网格上，一旦真实的频率值不在预先设定的网格上，即在网格失配的情况下，其估计性能会急剧恶化。另一方面，虽然更密集的网格理论上可以减小重构误差，但是太过密集的离散化网格会使得过完备化字典原子间高度相关。在其他领域，通用一阶泰勒展开来逼近真实字典来解决字典失配的问题，例如，在DOA估计中，利用稀疏贝叶斯推论OFFGRIDSPARSEBAYESIANINFERENCE，OGSBI来重构稀疏信号，可以估计不在网格上的来波方向角。类似地，该方法也可以用于复正弦波频率估计，但是这种方法的缺。
9、点是估计结果容易受初值的影响，且计算量较大。发明内容0004本发明的目的在于，提出一种基于多观测向量稀疏表示的复正弦波频率估计方法，在频率网格失配的情况下，能够对频率值进行修正，提高估计精度。0005本发明基于观测数据稀疏表示的模型，首先针对网格失配的情况，通过一阶泰勒展开来逼近真实的频率向量，得到多观测向量下的稀疏表示逼近模型，然后用正交匹配追踪算法和最小二乘法计算复正弦波频率的估计值。0006本发明的目的通过如下步骤实现0007S1、由长度为P的复正弦信号XI加白噪声得到接收数据其中，I1，P，K为复正弦波信号的个数，为第K个信号的复幅度，FK0，1为第K个信号的归一化频率，K1，2，K，。
10、初始相位是在0，2均匀分布的随机变量，VI是零均值、方差为的高斯白噪声，当IK时，和相互独立，说明书CN104166795A2/6页50008S2、S1所述接收数据YI的向量形式为YAFS0V，其中，YY1Y2YPT，频率向量幅度向量S012KT，VV1V2VPT；0009S3、将S2所述接收数据YI的向量形式Y写成稀疏表示的模型Y，具体如下0010S31、由S2所述接收数据YI，得到信号向量YIYIYI1YIM1T，从而得到数据矩阵其中，LM1P，M为信号向量的长度，L为信号向量的个数；0011S32、把S1所述归一化频率FK在0，1的范围上过完备化为一个离散的网格N代表网格数且满足NK；0。
11、012S33、将S31所述信号向量YI写成多观测向量的稀疏表示形式其中，过完备矩阵频率向量是一个K行稀疏矩阵，即S中只有K行元素为非零，其他行的元素都为零，记所述K行元素的索引值为支撑集，VIVIVI1VIM1T；0013S4、在S3所述稀疏表示的模型Y的基础上，用一阶泰勒展开来逼近真实的频率向量其中，存在第K个真实频率且是网格上距离FK最近的点，是的一阶导数，NK1，N；0014S5、根据S4所述AFK构造过完备字典记DIAG，12NT，0015当时，则0016当则N0，0017因此，观测模型可以写成0018S6、利用S5所述观测模型Y通过正交匹配追踪ORTHOGONALMATCHINGPU。
12、RSUIT，OMP算法和最小二乘方法求解出真实频率值的估计值，具体为0019S61、用正交匹配追踪算法求解出网格上离真实频率值最近的点，具体如下说明书CN104166795A3/6页60020S611、初始化记迭代次数T1，残差矩阵RY，支撑集其中0NL表示维数为NL的全零矩阵；0021S612、计算CAHR，ARGMAX|CI，|2，其中H表示矩阵的共轭转置，|CI，|2表示矩阵C第I行的L2范数，ARGMAX表示最大值的索引值；0022S62、用最小二乘方法计算S和Q在支撑集上元素的估计值和并更新残差矩阵更新迭代次数TT1；0023S63、当迭代次数TK时，返回步骤S612，依次循环迭代；。
13、一旦TK，停止迭代，输出和把和的对应元素相除，得到维数为KL的矩阵求矩阵每一行的平均值，得到向量在支撑集上元素的估计值0024S64、设S63所述在网格上的索引值为N1，NK，则修正后的频率估计值为0025本发明的有益效果是0026本发明可在粗糙的网格上进行复正弦波频率值的精确估计，避免了密集的网格带来的高计算量，提高了估计精度。利用本发明的估计方法，修正后的估计性能有明显的提高，即使在粗糙的网格，也有很高的估计精度。本发明的估计方法与其他估计方法相比，估计性能有明显的提高。附图说明0027图1是本发明方法的流程图。0028图2是本发明方法与其他频率估计方法的均方根误差随信噪比变化曲线图。00。
14、29图3是本发明方法与其他频率估计方法的均方根误差随信号向量个数L变化曲线图。具体实施方式0030下面结合实施例和附图，详细说明本发明的技术方案。0031如图1所示，一种基于多观测向量稀疏表示的高精度频率估计方法包括如下步骤0032S1、由长度为P的复正弦信号XI加白噪声得到接收数据其中，I1，P，K为复正弦波信号的个数，为第K个信号的复幅度，FK0，1为第K个信号的归一化频率，K1，2，K，初始相位是在0，2均匀分布的随机变量，VI是零均值、方差为的高斯白噪说明书CN104166795A4/6页7声，所述VI与信号相互独立，当IK时，和相互独立，0033S2、S1所述接收数据YI的向量形式为。
15、YAFS0V，其中，YY1Y2YPT，频率向量幅度向量S012KT，VV1V2VPT，进行频率估计的目的是根据接收的P1维数据，估计K个信号的归一化频率FK，K1，2，K；0034S3、将S2所述接收数据YI的向量形式Y写成稀疏表示的模型Y，具体如下0035S31、由S2所述接收数据YI，得到信号向量YIYIYI1YIM1T，从而得到数据矩阵其中，LM1P，M为信号向量的长度，L为信号向量的个数；0036S32、把S1所述归一化频率FK在0，1的范围上过完备化为一个离散的网格N代表网格数且满足NK；0037S33、将S31所述信号向量YI写成多观测向量的稀疏表示形式其中，过完备矩阵频率向量是一。
16、个K行稀疏矩阵，即S中只有K行元素为非零，其他行的元素都为零，记所述K行元素的索引值为支撑集，VIVIVI1VIM1T；0038S4、在S3所述稀疏表示的模型Y的基础上，用一阶泰勒展开来逼近真实的频率向量其中，是网格上距离FK最近的点，K1，K且是的一阶导数；0039S5、根据S4所述AFK构造过完备字典记DIAG，12NT，0040当时，则0041当则N0，0042因此，观测模型可以写成说明书CN104166795A5/6页80043S6、利用S5所述观测模型Y通过正交匹配追踪ORTHOGONALMATCHINGPURSUIT，OMP算法和最小二乘方法求解出真实频率值的估计值，具体为0044。
17、S61、用正交匹配追踪算法求解出网格上离真实频率值最近的点，具体如下0045S611、初始化记迭代次数T1，残差矩阵RY，支撑集其中0NL表示维数为NL的全零矩阵；0046S612、计算CAHR，ARGMAX|CI，|2，其中H表示矩阵的共轭转置，|CI，|2表示矩阵C第I行的L2范数，ARGMAX表示最大值的索引值；0047S62、用最小二乘方法计算S和Q在支撑集上元素的估计值和并更新残差矩阵更新迭代次数TT1；0048S63、当迭代次数TK时，返回步骤S612，依次循环迭代；一旦TK，停止迭代，输出和把和的对应元素相除，得到维数为KL的矩阵求矩阵每一行的平均值，得到向量在支撑集上元素的估计。
18、值0049S64、设S63所述在网格上的索引值为N1，NK，则修正后的频率估计值为0050下面结合具体实施例来进行说明。0051实施例1、估计值的均方根误差随信噪比变化0052采用接收数据的长度为P207，信号向量长度为M8。两个相同功率的复正弦波的归一化频率为0253，0418，即K2。为了使真实的复正弦波频率不落在网格上，取离散化网格为0，001，099，间隔001。参考信噪比SNR从5DB到20DB变化，间隔为5DB，每个信噪比进行1000次蒙特卡洛实验。0053根据不同信噪比下的接收数据YI得到稀疏表示模型YASV；0054利用泰勒一阶展开对真实的频率向量进行逼近，进而得到多观测向量下。
19、的稀疏表示逼近模型0055通过正交匹配追踪算法计算得到网格上离真实频率值最近的点0056通过最小二乘方法求解上一步骤波达方向角的修正值0057得到修正后的频率估计值0058按照本发明的方法估计得到的频率值的均方根误差随信噪比变化的曲线如图2所示。图2可以看到，利用本发明的估计方法，修正后的估计性能有明显的提高，即使在粗糙的网格，且信噪比为0DB时都能达到0005以内的估计精度。本发明的估计方法与其他估计方法相比，估计性能有明显的提高，说明了本发明的估计方法是有效的。而且本发明说明书CN104166795A6/6页9方法所需的计算时间大约为0005秒左右，而其他估计方法的计算时间均在1秒左右，说。
20、明了本发明估计方法的高效性。用均方根误差RMSE来评估各算法的性能，其定义为其中，MON为蒙特卡洛实验次数，和FK分别代表第M次蒙特卡洛实验估计得到的第K个频率值和第K个真实频率。0059实施例2、估计值均方根误差随信号向量个数L变化0060信号向量长度M8，信噪比固定为10DB。两个相同功率的复正弦波的归一化频率为0253，0418，即K2。为了使真实的复正弦波频率不落在网格上，取离散化网格为0，001，099，间隔001。信号向量个数L从20到200变化，间隔为20，向量个数L进行1000次蒙特卡洛实验。用均方根误差RMSE来评估各算法的性能，其定义为其中，MON为蒙特卡洛实验次数，和FK。
21、分别代表第M次蒙特卡洛实验估计得到的第K个频率值和第K个真实频率。0061根据不同信号向量个数L下的接收数据YI得到稀疏表示模型YASV；0062利用泰勒一阶展开对真实的频率向量进行逼近，进而得到多观测向量下的稀疏表示逼近模型0063通过正交匹配追踪算法计算得到网格上离真实频率值最近的点0064通过最小二乘方法求解上一步骤波达方向角的修正值0065得到修正后的频率估计值0066按照本发明的方法估计得到的频率值的均方根误差随信号向量个数L变化的曲线如图3所示。图3可以看到，利用本发明的估计方法，修正后的估计性能有明显的提高，即使在粗糙的网格，且信号向量个数L40时都能达到0001以内的估计精度。本发明的估计方法与其他估计方法相比，估计性能有明显的提高，说明了本发明的估计方法是有效的。而且本发明方法所需的计算时间大约为0005秒左右，而其他估计方法的计算时间均在1秒左右，说明了本发明估计方法的高效性。说明书CN104166795A1/3页10图1说明书附图CN104166795A102/3页11图2说明书附图CN104166795A113/3页12图3说明书附图CN104166795A12。