write a code to do 3d beamforming in the y direction in matlab

main.m
numSources = 3;      % number of sources
numMics = 4;         % number of microphones
micSpacing = 0.1;    % distance between microphones (in meters)
150 chars
4 lines

main.m
yGrid = -1:0.1:1;   % grid in the y-direction (from -1 to 1 with a step size of 0.1)
85 chars
2 lines

main.m
angleGrid = linspace(-pi/2, pi/2, length(yGrid));   % angles from -pi/2 to pi/2
steeringMatrix = exp(1j*2*pi*micSpacing*sin(angleGrid(:))*(0:numMics-1));
154 chars
3 lines

main.m
sourcePositions = zeros(numSources, 3);   % source positions in 3D space
sourcePositions(:, 2) = [0.5, -0.4, 0.8]; % positions along the y-direction

receivedSignal = zeros(length(yGrid), numMics);   % received signal at each microphone

for s = 1:numSources
    distance = sqrt((sourcePositions(s, 1) - micSpacing*(0:numMics-1)).^2 ...
        + (sourcePositions(s, 2)).^2 + (sourcePositions(s, 3)).^2);
    receivedSignal = receivedSignal + exp(1j*2*pi*distance) * randn(1);
end
481 chars
11 lines

main.m
estimatedDOA = abs(steeringMatrix' * receivedSignal);
54 chars
2 lines