clear?all
close?all
clc
N?=?10^2;??%信號(hào)采樣個(gè)數(shù)
Pf?=?0.001;??%設(shè)定虛警概率
K?=?3;??%感知節(jié)點(diǎn)個(gè)數(shù)
SNR?=?-15:1:5;??%設(shè)置不同的信噪比
SNR1?=?power（10SNR/10）;??%將db轉(zhuǎn)化為十進(jìn)制
lambda?=?norminv（1-PfK*N*1sqrt（2*K*N）*1）;??%門(mén)限值???
%-------------------理論結(jié)果
Pd_ana?=?1-normcdf（lambdaK*N*（1+SNR1）sqrt（2*K*N）*（1+SNR1））;?%檢測(cè)概率
%-------------------仿真結(jié)果
Np=10^3;?%統(tǒng)計(jì)循環(huán)次數(shù)
Pd_sim?=?zeros（1length（SNR））;
for?i?=?1:length（SNR）
????d?=?0;
????for?n?=?1:Np
????????s?=?sqrt（SNR1（i））*randn（KN）;
????????noise?=?randn（KN）;
????????r?=?s?+?noise;
????????r_power?=?sum（sum（abs（r）.^2））;
????????if?（r_power?>?lambda）
????????????d?=?d+1;
????????end
????end
????Pd_sim（i）?=?d/Np;
end


plot（SNRPd_ana‘r-+‘SNRPd_sim‘r-s‘）;
hold?on;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????1554??2013-10-16?16:16??single_node_multi_node.m