• 大小: 147KB
    文件类型: .rar
    金币: 1
    下载: 0 次
    发布日期: 2021-05-25
  • 语言: Matlab
  • 标签:

资源简介

可以的用的matlab阵列信号处理MVDR程序,喜欢的可以试试

资源截图

代码片段和文件信息

function [MSC]=coherence_MVDR(x1x2LK);

%% This program computes the coherence function between 2 signals 
%% x1 and x2 with the MVDR method.
%% This algorithm is based on the paper by the same authors:
%% J. Benesty J. Chen and Y. Huang “A generalized MVDR spectrum“ 
%% IEEE Signal Processing letters vol. 12 pp. 827-830 Dec. 2005.

%% x1 first signal vector of length n
%% x2 second signal vector of length n
%% L is the length of MVDR filter or window length
%% K is the resolution (the higher K the better the resolution)

%initialization
xx1     = zeros(L1);
xx2     = zeros(L1);
r11     = zeros(L1);
r22     = zeros(L1);
r12     = zeros(L1);
r21     = zeros(L1);

%construction of the Fourier Matrix
F       = zeros(LK);
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);

%number of samples equal to the lenght of x1 and x2
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(r12conj(r21));
%
%for regularization
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);
%
%Magnitude squared coherence function
MSC     = real(S12.*conj(S12))./(Si11.*Si22);

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1701  2006-01-27 11:13  MVDR程序\coherence_MVDR\coherence_MVDR.m

     文件      72048  2005-09-14 08:30  MVDR程序\coherence_MVDR\icassp2006_jb.pdf

     文件       1187  2006-01-27 11:12  MVDR程序\coherence_MVDR\Illustrate.m

     文件     143512  2005-11-21 09:42  MVDR程序\coherence_MVDR\spl_dec2005.pdf

     文件          0  2006-11-28 20:46  MVDR程序\coherence_MVDR\新建 文本文档.txt

     目录          0  2007-04-10 12:35  MVDR程序\coherence_MVDR

     目录          0  2007-04-10 12:35  MVDR程序

     文件         26  2018-10-13 19:20  readme.txt

----------- ---------  ---------- -----  ----

               218474                    8


评论

共有 条评论