- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
用MVDR方法估计相干函数
用MVDR方法估计相干函数
一、MVDR谱
用MVDR方法做信号的谱估计是基于滤波器分解。
设输入X(n)是一组均值为0的随机数,作为长度为N,的K阶滤波器的输入。
滤波器的输出为,输出信号的能量为:
其中,。
考虑这样一个L*K的矩阵
并且有,k=0,1,2…K-1,K=L。矩阵叫做傅里叶矩阵并且是酉矩阵,则。在MVDR谱中为了使滤波器的输出方差最小,所以滤波器系数必须满足下面的约束条件:
(1)
在这个约束条件下,输入X(n)在频率Wk处无失真,在其他频率处被衰减。这也相当于是使下面的公式最小:
(2)
为拉格朗日乘数。使(2)式最小的解为:
(3)
定义X(n)在Wk的谱为
(4)
把3式代入4式可得
(5)
由3式和5式可得出:
考虑所有的向量fk,总的表达式变为:
,
,为对角阵。
二、MVDR交叉谱
这里假设有两个0均值的稳态随机信号和,他们各自的谱分别为和。由前面的公式可得出两个滤波器为
(6)
和在Wk出的谱为:
(7)
其中:
(8)
设和分别为和的输出。定义和的交叉谱:
(9)
由公式9可以的出
(10)
其中
是和的互相关矩阵。把6式代入10式中可得:
(11)
三、MVDR相关函数
定义X1(n)和X2(n)幅值均方相关函数(MSC)为:
(12)
通过公式11可以得出:
(13)
把公式13和11代入公式12中,MSC变为:
(14)
并且有
四、代码和仿真结果:
function [MSC]=coherence_MVDR(x1,x2,L,K);
xx1 = zeros(L,1);
xx2 = zeros(L,1);
r11 = zeros(L,1);
r22 = zeros(L,1);
r12 = zeros(L,1);
r21 = zeros(L,1);
F= zeros(L,K);
l= [0:L-1];
f= exp(2*pi*l*j/K);
for k = 0:K-1
F(:,k+1) = f.^k;
end
F = F/sqrt(L);
n = length(x1);
for i = 1:n
xx1 = [x1(i);xx1(1:L-1)];
xx2 = [x2(i);xx2(1:L-1)];
r11 = r11 + xx1*conj(xx1(1));
r22 = r22 + xx2*conj(xx2(1));
r12 = r12 + xx1*conj(xx2(1));
r21 = r21 + xx2*conj(xx1(1));
end
r11 = r11/n;
r22 = r22/n;
r12 = r12/n;
r21 = r21/n;
R11 = toeplitz(r11);
R22 = toeplitz(r22);
R12 = toeplitz(r12,conj(r21));
Dt1 = 0.01*r11(1)*diag(diag(ones(L)));
Dt2 = 0.01*r22(1)*diag(diag(ones(L)));
Ri11 = inv(R11 + Dt1);
Ri22 = inv(R22 + Dt2);
Rn12 = Ri11*R12*Ri22;
Si11 = real(diag(F*Ri11*F));
Si22 = real(diag(F*Ri22*F));
S12 = diag(F*Rn12*F);
MSC = real(S12.*conj(S12))./(Si11.*Si22);
clc
clear all
close all
n = 1024;
nT= [0:n-1];
Nf= 5;
f = zeros(Nf,1);
f(1) = 0.05; f(2)
您可能关注的文档
最近下载
- 中央空调系统清洗维修保养合同协议书(详细版).doc
- (ppt)苦菜花.ppt
- 免疫接种必须与选修答案-2024年全国疾控系统“大学习”活动.docx VIP
- (统编2024版)历史七年级上册全册知识点.docx
- 2023-2024在线网课《信息检索与科技写作( 理大)》单元测试考核答案.pdf
- 中医院重点专科儿科三年发展规划.docx
- 天津市河西区2023-2024学年七年级上学期期中语文试题(含答案).docx VIP
- 第九讲种质资源的保存ppt课件-教材.ppt VIP
- YDT 3762-2020 大数据 数据挖掘平台技术要求与测试方法.docx
- 2.1 网络改变世界 作业 初中道德与法治 人教部编版 八年级上册 (2021年).docx
文档评论(0)