谱估计 MUSIC算法 matlab程序

clear all
N=256;                               %观测样本个数                  
noise=sqrt（0.01）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
                                     %噪声信号
noise2=sqrt（6）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
noise3=sqrt（6）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
noise4=sqrt（6）*（（randn（1N）+1i*randn（1N））/sqrt（2））;
f1=0.1;                              %信号归一化频率
f2=0.25;
f3=0.27;
A1=0.1;                                %信号幅度   
A2=0.1;
A3=0.1;
snr1=-15;
snr =power（10snr1/10）;%db转化十进制
pn = （std（noise））.^2;
A=sqrt（pn*snr）;
signal1=A1*exp（1i*2*pi*f1*（0:N-1））;  %有用信号
signal2=A1*exp（1i*2*pi*f2*（0:N-1））;
signal3=A1*exp（1i*2*pi*f3*（0:N-1））;
for j=1:30
    noise1=（0.5*i）*sqrt（20）*（（randn（1N）+1i*randn（1N））/sqrt（2