music算法通过麦克风阵列估计声源方向matlab

Signal_No = 1;
Interference_No = 2;
S_No = Signal_No + Interference_No;
sensor_No = 8;
azimuth = [0*pi/180 10*pi/180 70*pi/180];
Fs = 30e6;
Fc = 10e6;
F1 = [8e6 20e6];
theta = [10 45 52 93 43]*pi/180;
RB = 2e6;
M = 2;
wavelength = 3e8/Fc;  %波长
d = wavelength/2;
K = 30;
Ps =  -30;
SNR = 10;
SIR = -20;
M_No = 40;
Data_No = M_No*Fs/RB;
t = 1/Fs:1/Fs:Data_No/Fs;
Ps_1 = 10^（Ps/10）;
Pi_1 = 10^（（Ps-SIR）/10）;
Pn_1 = 10^（（Ps-SNR）/10）;
bit = randint（1M_No）;
bitstream = [];
 for i = 1:M_No
       if bit（i）==1
             bitstream = [bitstream ones（1 Fs/RB）];
       else
             bitstream = [bitstream -ones（1 Fs/RB）];
       end
end
 Carrier_R = cos（2*pi*Fc*t）;
 S_R = Carrier_R .* bitstream;
 Carrier_I = sin（2*pi*Fc*t）;
 S_I = Carrier_I .*bitstream;
 
 b = sqrt（（S_R*S_R‘/length（S_R）））;
 
 Signal_R = sqrt（Ps_1）*S_R/sqrt（（S_R*S_R‘/length（S_R）））;
 Signal_I = sqrt（Ps_1）*S_I/sqrt（（S_I*S_I‘/length（S_I）））;
 Signal（1:） = comp