MVDR方法在信号频率估计领域的应用.ppt

《MVDR方法在信号频率估计领域的应用.ppt》由会员分享，可在线阅读，更多相关《MVDR方法在信号频率估计领域的应用.ppt（22页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、12345目目 录录课题背景MVDR频率估计算法仿真仿真结果分析思考题解答MVDR算法基本原理3 课课 题题 背背 景景课题背景课题背景MVDR是什么是什么MVDR（minimum variance distortionless response）是最小方差无失真响应算法的简称，是一种有别于经典功率谱估计和参数模型估计的另一类信号频谱估计方法，最早于1969年由Capon提出。MVDR的优缺点的优缺点优优点：点：在于低阵元和低快拍数的情况下，仍然能保证较高的分辨率。缺点：缺点：当阵元数下降和分辨率要求提高时，它的解相干性能受到严重限制，甚至无法实现。（MMVDR）MVDR的应用的应用日常：日常

2、：被用于无线定位技术中，可以为移动通信网中的用户提供位置信息。军事：军事：被用于水声无线通信技术中，可以实现水面舰艇和潜艇之间、潜艇和潜艇之间的通信及水中武器的遥测遥控。5MVDR算法基本原理MVDR滤波器滤波器原理：原理：图1 M抽头的FIR滤波器 如图1所示，有M个权系数的横向滤波器（transversal filter），滤波器的输入为随机过程x（n），输出为定义输入信号向量和权向量分别为 i-nxny1-M0ii* Tx nx nx n-1x n-M1，T1-M10， *THnnnyxx 信号y(n)的平均功率可以表示为 其中矩阵 为向量输入信号向量x(n)的M维自相关矩阵 RxxHH

3、H2nnEnyEPMMCRMVDR滤波器滤波器原理：原理：MVDR算法的原理：算法的原理： （1）约束 ，这是为了使 无失真地通过滤波器。 （2）输出平均功率P= 最小，达到抑制其他频率信号和噪声的目的。MVDR算法的原理：算法的原理： 实现最小方差无失真相应的基本思路为通过调节权向量，使噪声和来自非期望信号频率的所有干扰所贡献的功率最小，同时保持期望信号频率增益恒定。故滤波器权向量 应满足：nj11eRH11Hw 更一般的情况，实际工程中常采用N个观测样本值得到相关矩阵 。对于无失真通过系统的信号频率为 的复正弦信号，其最优权向量和滤波器最小输出功率分别为11-M VDR-HRRR 11MV

4、DR-HPR 在 内改变 ,画出 曲线。当 与输入信号的频率不相等时，信号和噪声都被滤波器抑制；当 与输入信号的频率相等时，该信号可以无失真的通过，因此曲线呈现出一个峰值。, MVDRP11MVDRMVDR频率估计算法仿真频率估计算法仿真MVDR仿真实例仿真实例： 假设随机过程u(n)为其中，v(n)是零初值、方差为1 的白噪声， 和 是相互独立并在 上服从均匀分布的随机相位。请使用MVDR方法进行信号频率估计的仿真实验，画出频率估计谱线，并给出正弦信号频率的估计值。 （要求：信号样本数取1000，估计的自相关矩阵为8阶）11(n)exp(j0.5 nj)exp( j0.3 nj) v(n)+

5、u120,2 算法流程框图 我们分别取滤波器抽头系数个数为M=4及M=8，并在 范围内均匀选取2048个频率点，此时的MVDR频率估计谱线以及正弦信号频率的估计值如图2、图3所示。2wf图2 M=4时的频率估计谱线0.5,0.5f-图3 M=8时的频率估计谱线 我们分别取滤波器抽头系数个数为M=16及M=32，并在 范围内均匀选取2048个频率点，此时的MVDR频率估计谱线以及正弦信号频率的估计值如图4、图5所示。2wf图4 M=16时的频率估计谱线0.5,0.5f-图5 M=32时的频率估计谱线MVDR仿真结果分析仿真结果分析：表1 实验结果信号信号理论频率理论频率M=4时测得时测得实际频率

6、实际频率M=8时测得时测得实际频率实际频率M=16时测时测得得实际频率实际频率M=32时测时测得得实际频率实际频率1 1-0.15-0.147289-0.1502200.1502200.1502202 20.250.2484120.2503660.2498780.249878平均相平均相对误差对误差1.2%0.15%0.1%0.1% 从实验结果图中可以看出随着抽头数M的增加，MVDR算法对频谱的分辨率得到了提高。但通过对表1 的数据进行分析发现，当抽头数M小于16时，抽头数的增加使频率估计的精度得到了提高；但当抽头数M大于16时，抽头数的增加对频率估计的精度没有明显的提升效果。17思考题解答思

7、考题解答思考题：思考题：1.给出DFT的定义和主要性质l定义： 离散傅里叶变换（Discrete Fourier Transform，缩写为DFT），是傅里叶变换在时域和频域上都呈离散的形式。在形式上，变换两端（时域和频域上）的序列是有限长的，而实际上这两组序列都应当被认为是离散周期信号的主值序列。即使对有限长的离散信号作DFT，也应当将其看作其周期延拓的变换。在实际应用中通常采用快速傅里叶变换计算DFT。 设x(n)为M点有限长序列，则可定义x(n)的N点离散傅里叶变换为120:XknekNjnNnDFTx,(k=0,1,N-1)l性质：（1）线性（2）圆周移位性（3）圆周共轭对称性 共轭对

8、称序列： 共轭反对称序列：（4）对偶性 把离散谱序列当成时域序列进行DFT，结果是原时域序列反褶的N倍；如果原序列具有偶对称性，则DFT结果是原时域序列的N倍2(k)DFT( )DFT () R (N)(k)=ejkmNmmNNXxnx nmX*(n)( n)(n)(N),n0,1,N 1epepNNepxxRxn*(n)( n)(n)(N),n0,1,N 1-opopNNopxxRxn1212ax (n)bx (n)a( )b(k)DFTX kXl性质：（5）离散圆卷积 时域离散圆卷积 频域离散圆卷积（6）帕斯瓦尔定理 (n)y(n)y(n)(n)(k)Y(k)xxX1(n)y(n)(k)Y

9、(k)xXN2112001(n)(k)=NNnnxXN2.抽样过程为什么要先进行滤波，此滤波器应逼近什么样的指标？ 因为里面的谐波含量太高,会造成控制系统的误动作,所以,需要先经过一个低通滤波器进行过滤,将尽可能多的主要谐波过滤掉,保证整个工控系统控制的精度。 该低通滤波器的通带截止频率应小于1/2采样频率，以达到滤去高频成分，防止频率混叠的目的。实际仪器设计中,这个低通滤波器的截止频率(fc) 为：截止频率（fc）= 采样频率（fz） / 2.56进入夏天，少不了一个热字当头，电扇空调陆续登场，每逢此时，总会进入夏天，少不了一个热字当头，电扇空调陆续登场，每逢此时，总会想起那一把蒲扇。蒲扇，

10、是记忆中的农村，夏季经常用的一件物品。记想起那一把蒲扇。蒲扇，是记忆中的农村，夏季经常用的一件物品。记忆中的故乡，每逢进入夏天，集市上最常见的便是蒲扇、凉席，不论男女老忆中的故乡，每逢进入夏天，集市上最常见的便是蒲扇、凉席，不论男女老少，个个手持一把，忽闪忽闪个不停，嘴里叨叨着少，个个手持一把，忽闪忽闪个不停，嘴里叨叨着“怎么这么热怎么这么热”，于是三，于是三五成群，聚在大树下，或站着，或随即坐在石头上，手持那把扇子，边唠嗑五成群，聚在大树下，或站着，或随即坐在石头上，手持那把扇子，边唠嗑边乘凉。孩子们却在周围跑跑跳跳，热得满头大汗，不时听到边乘凉。孩子们却在周围跑跑跳跳，热得满头大汗，不时听

11、到“强子，别跑强子，别跑了，快来我给你扇扇了，快来我给你扇扇”。孩子们才不听这一套，跑个没完，直到累气喘吁吁，。孩子们才不听这一套，跑个没完，直到累气喘吁吁，这才一跑一踮地围过了，这时母亲总是，好似生气的样子，边扇边训，这才一跑一踮地围过了，这时母亲总是，好似生气的样子，边扇边训，“你你看热的，跑什么？看热的，跑什么？”此时这把蒲扇，是那么凉快，那么的温馨幸福，有母亲此时这把蒲扇，是那么凉快，那么的温馨幸福，有母亲的味道！蒲扇是中国传统工艺品，在我国已有三千年多年的历史。取材的味道！蒲扇是中国传统工艺品，在我国已有三千年多年的历史。取材于棕榈树，制作简单，方便携带，且蒲扇的表面光滑，因而，古人

12、常会在上于棕榈树，制作简单，方便携带，且蒲扇的表面光滑，因而，古人常会在上面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名，实即今日的蒲扇，江浙称之为面作画。古有棕扇、葵扇、蒲扇、蕉扇诸名，实即今日的蒲扇，江浙称之为芭蕉扇。六七十年代，人们最常用的就是这种，似圆非圆，轻巧又便宜的蒲芭蕉扇。六七十年代，人们最常用的就是这种，似圆非圆，轻巧又便宜的蒲扇。蒲扇流传至今，我的记忆中，它跨越了半个世纪，也走过了我们的扇。蒲扇流传至今，我的记忆中，它跨越了半个世纪，也走过了我们的半个人生的轨迹，携带着特有的念想，一年年，一天天，流向长长的时间隧半个人生的轨迹，携带着特有的念想，一年年，一天天，流向长长的时间隧道，袅道，袅