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