total.m

Matrix Pencil Method（矩阵束） 和MUSIC 算法的比较

% Comparison between Element Music and MPM Using Angle
% the sensor start with 0, not central symmetric
clear all;
M=20;  % number of elements
N=M;
d=1;   % number of signals
theta=72;
f=10e6;
w=2*pi*f;
nsample=1;    % snapshots
B=exp(j*w*[0:nsample-1]);
m=1:M;  % the array is central symmetric
u=exp(j*pi*sin(theta/180*pi)*m);
%%%
E=1;           % engergy of signal
rmse_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 on
xlabel('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')