clear?all;
clc;
n?=10000;?????????????????????%采樣點(diǎn)數(shù)
pfa?=?0.01;???????????????????%虛警概率
snr1=-25:1:0;???????????????%信噪比
k??=?length（snr1）;
for?i?=?1:k?????????????????%不同信噪比下仿真
????s?=?0;??????????????????%1個(gè)數(shù)初值設(shè)為0，每次循環(huán)都清零
????snr（i）?=?power（10snr1（i）/10）;%db轉(zhuǎn)化十進(jìn)制
????for?kk?=?1:1?????????%每種信噪比下仿真次數(shù)?
????????t?=?1/n:1/n:1;??????????????%離散點(diǎn)1000
????????x?=?sin（pi*t）;
????????x?=?（?x?>?0）;????????????????%大于0判為1，否則0，方波產(chǎn)生方法
????????x??=2*x-1;???????????????????%信源數(shù)據(jù)產(chǎn)生
????????noise?=?randn（1length（t））;??%均值0方差1的噪聲
????????pn?=?（std（noise））.^2;????????%噪聲功率
????????a?=?sqrt（2*snr（i））;???????%信號(hào)幅值
????????xx?=?a.*sin（2*pi*4*t）;???????%調(diào)制信號(hào)????
????????y??=?x.*xx;??????????????????%bpsk調(diào)制信號(hào)
????????ps?=?mean（abs（y）.^2）;????????%信號(hào)平均功率
????????yy?=?y+noise;????????????????%輸出信號(hào)
???????
????????r?=?（n+sqrt（2*n）*sqrt（2）*erfcinv（2*pfa））;
????????%r?=?gaminv（1-pafn/21）*2;
????????sum1?=?sum（abs（yy）.^2）;??????%輸出信號(hào)的累積功率
????????if?sum1>r
????????????s?=?s+1;?????????????????%輸出大于閾值及判為1，否則0，統(tǒng)計(jì)1的個(gè)數(shù)
????????end
????end
?p（i）?=?s/100;?????????????%檢測(cè)概率：1的個(gè)數(shù)/n
pd（i）?=?0.5*erfc（（erfcinv（2*pfa）-snr（i）*sqrt（n/2））/sqrt（2+4*snr（i）））;
%awgn?channels
%pd（i）?=?marcumq（sqrt（snr（i）*2）sqrt（r）n/2）;
%pd（i）=pfa+exp（1/（（2*n）*（snr（i）^2））-sqrt（2）*erfcinv（2*pfa）/（sqrt（n）*snr（i）））
%rayleigh?fading（pfa要很小《0.01）
end
plot（snr1pd‘o-‘）;
hold?on;

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????1517??2007-10-23?16:21??energydetection.m

-----------?---------??----------?-----??----

?????????????????1517????????????????????1