% Comparison between Element Music and MPM Using Angle% the sensor start with 0, not central symmetricclear all;M=20;  % number of elementsN=M;d=1;   % number of signalstheta=72;f=10e6;w=2*pi*f;nsample=1;    % snapshotsB=exp(j*w*[0:nsample-1]);m=1:M;  % the array is central symmetricu=exp(j*pi*sin(theta/180*pi)*m);%%%E=1;           % engergy of signalrmse_music=[];rmse_mpm=[];for snr_db=0:20    snr=10^(snr_db/10);    No=E/snr;    noise_rms=sqrt(No/2);    theta_music=[];    theta_mpm=[];    for k=1:100               % numbers of  samples of the estimate of theta        x=u.'*B+noise_rms*(randn(size(u.'*B))+j*randn(size(u.'*B)));        theta_music=[theta_music;music(x,d,N)];        theta_mpm=[theta_mpm;mpm(x,d,N)];    end;    sum_music=( theta_music - theta ).^2;    s_music=sqrt(1/100*sum(sum_music));    sum_mpm=( theta_mpm - theta ).^2;    s_mpm=sqrt(1/100*sum(sum_mpm));        rmse_music=[rmse_music s_music];    rmse_mpm=[rmse_mpm s_mpm];end;figure(1)snr_db=0:20;semilogy(snr_db,rmse_music,'b',snr_db,rmse_mpm,'r--',snr_db,rmse_music,'bo',snr_db,rmse_mpm,'rx')grid onxlabel('SNR(dB)');ylabel('RMS Error(degree)');%title('Comparison between Element space Music and Beamspace Music Using Angle');legend('Element space Music','MPM','Location','Northeast')%text(-2,1.5,'The number of antenna elements is 21')