matlab开发-最大最小特征值组合

matlab开发-最大最小特征值组合。有说服力无线电的CMME检测

clc;
clear;
close all;

tic;
%% SNR vs pd using maximum eigenvalue detection
Ns=10000;
SNR=-25:-16;
snr=10.^（SNR./10）;
L=8;
pf=0.1;
F1_inv=0.45;
num_iter=2000;
a=（（sqrt（Ns）+sqrt（L））^2）/（Ns）;
b=1+（（sqrt（Ns）+sqrt（L））^（-2/3）*F1_inv）/（（Ns*L）^（1/6））;
threshold=a*b;
%
w=waitbar（0‘calculating probability of false alarm‘）;
l=1;
for i=1:length（SNR）
count=0;
for h=1:num_iter
signal=randn（1Ns）;
noise=randn（1Ns）;
noise_power=norm（noise）^2;
signal_power=norm（signal）^2;
mult=sqrt（snr（i）*noise_power/signal_power）;
signal=mult*signal;
signal=signal+noise;
%noise=noise/std（noise）;
k=0;
Cx=zeros（1L）;
for n=0:L-1
for j=1:1:Ns-L-1
Cx（n+1）=Cx（n+1）+signal（j）*signal（j+k）;
end
k=k+1;
end
Cx=Cx/Ns;
Cx_mtx=toeplitz（Cx）;
eig_value=eig（Cx_mtx）;
max_eig=max（max（eig_value））;
min_eig=min（min（eig_value））;
ratio=max_eig/（max_eig-min_eig）
if ratio>threshold
count=count+1;
end
end
pd（i）=count/2000;
waitbar（i/length（SNR））;
end
close（w）;
plot（SNRpd‘-^‘）;

toc;
xlabel（‘SNR‘）;
ylabel（‘P_d‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1029  2016-07-23 05:36  ch2.m
     文件        1314  2016-07-23 05:36  license.txt