修正的MUSIC算法的DOA估计

signal_no=1;
interference_no=1;
s_no=signal_no+interference_no;
sensor_no=8;
azimuth=[0*pi/180 5*pi/180];
fs=30e6;
fc=10e6;
f1=[8e6];
theta=[10 45 52 93 43]*pi/180;
rb=2e6;
m=2;
wavelength=3e8/fc;
d=wavelength/2;
k=300;
ps=-30;
snr=10;
sir=00;
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=[bitstreamones（1fs/rb）];
    else
        bitstream=[bitstream-ones（1fs/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;
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:）=complex（signal_rsignal_i）;
for i=2:interference_no+