untitled.m

Matlab仿真比较分析三种矩阵迹检测器的检测性能。

clear;clc;
M=3;%发射阵数；
N=3;%接收阵数；
L=20;%编码脉冲数；
sigma0=5;
% 先考虑均匀通道
% snrdb=[-20:1:20];snr=10.^(snrdb./10);
% sigmas=(snr.*sigma0^2).^(1/2);
% for i=1:1:length(snrdb)
%     for k=1:1:10000
%         Xnoise=sigma0.*(randn(N*M,L)+j.*randn(N*M,L))./sqrt(2);
%         Xsignal=sigmas(i).*(randn(N*M,L)+j.*randn(N*M,L))./sqrt(2);
%         % 有目标：
%         X1=Xsignal+Xnoise;
%         W1=X1*X1';
%         T11(i,k)=trace(W1);
%         T21(i,k)=(norm(X1,'fro'))^2;
%         % 无目标：
%         W0=Xnoise*(Xnoise)';
%         T10(i,k)=trace(W0);
%         T20(i,k)=(norm(Xnoise,'fro'))^2;
%     end
%     % 检测器
%     tmp1=abs(T10(i,:));tmp2=sort(tmp1,'descend');tmp3=tmp2(100);
%     tmp4=abs(T11(i,:));tmp5=find(tmp4>tmp3);tmp6=length(tmp5);
%     T1(i)=tmp6/length(tmp4);
%     
%     tmp1=abs(T20(i,:));tmp2=sort(tmp1,'descend');tmp3=tmp2(100);
%     tmp4=abs(T21(i,:));tmp5=find(tmp4>tmp3);tmp6=length(tmp5);
%     T2(i)=tmp6/length(tmp4);
%     
%     i,
% end
% figure(1);plot(snrdb,T1,'r',snrdb,T2,'k');hold on;legend('T1','T2');hold on;

% 再考虑非均匀通道
snrdb=[-20:1:20];snr=10.^(snrdb./10);
sigmas=(snr.*sigma0^2).^(1/2);
D0=diag([2.2310,2.2310,2.2310,2.2310,0,0,0,0,0]);%0.01,0.05,0.01,0.005,0.001]);
for i=1:1:length(snrdb)
    D=sigmas(i).*D0;
    for k=1:1:10000
        Xnoise=sigma0.*(randn(N*M,L)+j.*randn(N*M,L))./sqrt(2);
        Xsignal=D*(randn(N*M,L)+j.*randn(N*M,L))./sqrt(2);
        % 有目标：
        X1=Xsignal+Xnoise;
        W1=X1*X1';
        T31(i,k)=trace(W1);
        T41(i,k)=(trace(W1))^(N*M)/det(W1);
        W11=X1([1:1:4],:)*(X1([1:1:4],:))';
        T51(i,k)=trace(W11);
        % 无目标：
        W0=Xnoise*(Xnoise)';
        T30(i,k)=trace(W0);
        T40(i,k)=(trace(W0))^(N*M)/det(W0);
        W00=Xnoise([1:1:4],:)*(Xnoise([1:1:4],:))';
        T50(i,k)=trace(W00);
    end
    % 检测器
    tmp1=abs(T30(i,:));tmp2=sort(tmp1,'descend');tmp3=tmp2(100);
    tmp4=abs(T31(i,:));tmp5=find(tmp4>tmp3);tmp6=length(tmp5);
    T3(i)=tmp6/length(tmp4);
    
    tmp1=abs(T40(i,:));tmp2=sort(tmp1,'descend');tmp3=tmp2(100);
    tmp4=abs(T41(i,:));tmp5=find(tmp4>tmp3);tmp6=length(tmp5);
    T4(i)=tmp6/length(tmp4);
    
    tmp1=abs(T50(i,:));tmp2=sort(tmp1,'descend');tmp3=tmp2(100);
    tmp4=abs(T51(i,:));tmp5=find(tmp4>tmp3);tmp6=length(tmp5);
    T5(i)=tmp6/length(tmp4);
    
    i,
end
figure(6);plot(snrdb,T3,'r',snrdb,T4,'k',snrdb,T5,'b');hold on;legend('T3','T4','T5');hold on;