delay and sum beamforming(in frequency domain)
조회 수: 20 (최근 30일)
이전 댓글 표시
I have 8 microphones (excluding the origin[0;0]) arranged in circular geometry. And I have a sound source placed at a certain location (say 2mts) away from the array. The source should produce signal(say sinusoid) and it should impinge on microphones with certain delays due to its geometry.
Now i need to simulate this idea in MATLAB and plot the beam pattern by performing beamforming for the microphone array.
Here is my code so far:
clc
clear all
close all
Antenna_X_Pos=[-0.0025; 0.1406; -0.0685;-0.1134;0.1160;0.0640;-0.1471;-0.0013];
Antenna_Y_Pos=[0.0546;-0.1323;-0.2277;-0.0053;-0.0069;-0.2296;-0.1289;-0.0949];
Antenna_Z_Pos=[0.3446;0.2975;0.2605;0.2179;0.1770;0.1361; 0.0928 ; 0.2177];
scatter3(Antenna_X_Pos,Antenna_Y_Pos,Antenna_Z_Pos,'filled'),view(-60,60)
for j=1:360
theta(j)=degtorad(j);
[SP_X(j) SP_Y(j)] = pol2cart(theta(j),10);
for i=1:numel(Antenna_X_Pos)
Distance(i,j) = sqrt((Antenna_X_Pos(i) - SP_X(j))^2 +(Antenna_Y_Pos(i) - SP_Y(j))^2);
end
end
%plot(Distance)
Dist_m=Distance(1,:);
mMinusM1 = bsxfun(@minus,Distance,Dist_m);
speed_sound = 340; %speed of sound in meters
Time_delay = (mMinusM1/speed_sound)*10e3; %time delay in milliseconds
max((max(Time_delay)));
plot(Time_delay); % time delay plot for 8 microphones(including the origin of array)
Now how can I place the source(say a sinusoid) at the particular location and perform beamforming using the microphone array from that particular direction where the source is placed?
My final goal is to simulate conventional beamformer (Delay and Sum) algorithm.
am I missing some idea/trick? please help!!
댓글 수: 0
채택된 답변
Honglei Chen
2011년 12월 15일
Hi Zoro,
What you have so far is calculating the delays between the elements. The next steps will be
(1) simulate the received signal at each element. Here you need to use the delay you just calculated to delay the signal at each element accordingly.
(2) You can then use delay and sum beamformer to align the received together and see the combined signal.
(1) and (2) are essentially inverse operations. Because you are dealing with microphone signals, the signal by nature is wideband. So there are two ways you can do it, either you do the delay-and-sum in the time domain or you have to break the signal into narrowbands, do the delay-and-sum in frequency domain at each band and then combine result at each narrowband together to get the final result.
추가 답변 (1개)
Wayne King
2011년 12월 14일
편집: John Kelly
2014년 5월 27일
Hi Zozo, Are you able to use the Phased Array System Toolbox?
The Phased Array System Toolbox has a delay-sum narrowband beamformer (as well as a number of other beamforming algorithms).
See Narrowband (Phase Shift) Beamformer with a ULA for an example.
You can simulate arrays and directional of arrival easily.
댓글 수: 3
Wayne King
2011년 12월 14일
The Phased Array Toolbox is all written in the MATLAB language. If you own that toolbox, you are able to see the MATLAB code.
참고 항목
카테고리
Help Center 및 File Exchange에서 Array Geometries and Analysis에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!