FFT based adaptive MVDR beamforming

10 次查看（过去 30 天）
zozo
zozo 2012-9-22
Hello,
I have a small question on FFT based adaptive beamforming based on Spectral Matrix Inversion technique.
I have estimated the Spectral Matrix (Sxx) from my data by dividing it into blocks and averaging them for each bin, to get sufficient rank on the matrix.
Now my question is :
To produce the MVDR weighting vector (8x1 vector in my case since I have 8 microphones) for each bin, I multiply the inverse of Sxx for each bin (say bin#1) with the steering vector 'a' given by exp(i*2*pi*D*omega), where D is the set of time delays(8x1 vector) pointing towards look direction(say +10°). Here what should be 'omega'? Isit equal to the first frequency bin of data from first channel and so on for all the bins? please clarify.
  1 个评论
Gu
Gu 2012-11-22
I want to know how do you generate your signal? narrowband or wide.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（3 个）

Honglei Chen
Honglei Chen 2012-9-24
Looks like you are doing subband MVDR. It should be the frequency for each corresponding band.
  11 个评论
显示 9更早的评论
Honglei Chen
Honglei Chen 2012-9-26
For an arbitrary shaped array (sounds like your situation), I don't think there are a lot of theoretical results. It's a design problem you need to solve.
zozo
zozo 2012-9-27
编辑：zozo 2012-9-28
okay. And once more, I have tried estimating the spectral matrix from the data and adaptively place notches, but somehow it does not work. :-( I know that with increasing the data length and therefore the length of averaging, the estimated spectral matrix should converge to the true spectral matrix and so should the corresponding beampatterns. However, with true spectral matrix, i get perfect suppression and undistorted signal of interest. Just to sum up the main points of matrix estimation:
1) consider a data matrix 8x102400 [8 channels and 102400 samples(signal of duration 10seconds) with fs=10240].
2) divide the signal into 100 non-overlapping blocks each of size 8x1024.
3) consider the 1st block, take FFT of each channel to get 8x1024 matrix of complex numbers resulting due to FFT, in which 2:512 are positive freqs and 514:1024 are negative frequencies(complex conjugates of positive freqs) 1DC component and 1 Nyquist freq component at fs/2
4) Repeat the FFT for all 100 blocks(channel-wise).
5) Take the first column vector(8x1) from first block, first column(8x1) from second block and so on.. first column from 100th block to form an 8x100 matrix (say Sxx(1)).
6) Compute 1/100*(Sxx(1)*Sxx(1)') to get the spectral matrix for first bin.
7) multiply the resulting 8x8 spectral matrix with the steering vector 'a' (8x1) computed by exp(i*2*pi*D*omega(1)), to get the 8x1 weighting vector for first bin.
8) repeat the steps 5-7 for all 1024 bins to get 8x1024 steering vectors.
9) Finally, take FFT of data block again and apply the above computed steering vectors to the block of 1024 data length, take IFFT to get the beamformed signals in time domain.
10) Repeat the above steps for all look directions
Am I missing a trick?

请先登录，再进行评论。

zohar
zohar 2012-10-9
Hi zozo,
You allready posted a question
http://www.mathworks.com/matlabcentral/answers/39458-covariance-matrix-stimation-in-mvdr-beamforming
Combining the two questions ,evrything looks OK!
1) What is the shape of the array ? In the previous question you mentiond that it's spherical shape.
2) It's seems that in step 9 you want to go back to time domain, I do not understand what you doing ! can you supply the relevant matlab code ?
3) Insted of steps 7-9 calc P - the narrowband power spectrum, where P is:
P = zeros(Nbins,Nbeam)
for k = 1:Nbins
% for each bin
% calc A - sterring matrix.
% calc IRf - Inverse spatial covariance matrix.
P(k,:) = real(Nbeam./sum(A'*IRf*A));
end
Now calc
Plog = 20*log10(P);
mesh(Plog);
And let me know what the result is...
  3 个评论
显示 1更早的评论
zohar
zohar 2012-10-24
Nbeam is number of beams acorrding to sterring matrix. Let me help you do step (3) and supply matlab code I think yor are missing something befor ifft .
Have fun
zozo
zozo 2012-10-25
编辑：zozo 2012-10-25
thank you for offering help. I have mailed you the function I have written to estimate the covariance(spectral) matrix. Please have a look and correct me if I am wrong somewhere.

请先登录，再进行评论。

Gu
Gu 2012-11-22
How to generate wideband signal for wideband time domain beamforming.
  1 个评论
zozo
zozo 2012-11-22
use randn() to produce a gaussian distributed random sequence. and band pass filter it to have your desired band of frequencies. u can check it using magnitude plot of fft()

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Beamforming 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by