clear?all;
clc;
n?=10000;?????????????????????%采樣點數
pfa?=?0.01;???????????????????%虛警概率
snr1=-25:1:0;???????????????%信噪比
k??=?length（snr1）;
for?i?=?1:k?????????????????%不同信噪比下仿真
????s?=?0;??????????????????%1個數初值設為0，每次循環都清零
????snr（i）?=?power（10snr1（i）/10）;%db轉化十進制
????for?kk?=?1:1?????????%每種信噪比下仿真次數?
????????t?=?1/n:1/n:1;??????????????%離散點1000
????????x?=?sin（pi*t）;
????????x?=?（?x?>?0）;????????????????%大于0判為1，否則0，方波產生方法
????????x??=2*x-1;???????????????????%信源數據產生
????????noise?=?randn（1length（t））;??%均值0方差1的噪聲
????????pn?=?（std（noise））.^2;????????%噪聲功率
????????a?=?sqrt（2*snr（i））;???????%信號幅值
????????xx?=?a.*sin（2*pi*4*t）;???????%調制信號????
????????y??=?x.*xx;??????????????????%bpsk調制信號
????????ps?=?mean（abs（y）.^2）;????????%信號平均功率
????????yy?=?y+noise;????????????????%輸出信號
???????
????????r?=?（n+sqrt（2*n）*sqrt（2）*erfcinv（2*pfa））;
????????%r?=?gaminv（1-pafn/21）*2;
????????sum1?=?sum（abs（yy）.^2）;??????%輸出信號的累積功率
????????if?sum1>r
????????????s?=?s+1;?????????????????%輸出大于閾值及判為1，否則0，統計1的個數
????????end
????end
?p（i）?=?s/100;?????????????%檢測概率：1的個數/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;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????532482??2010-12-28?16:13??資料\A+Survey+of+Spectrum+Sensing+Algorithms+for2009.pdf

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

?????文件????????883??2010-12-28?16:10??資料\energydetection\energydetection.rar

?????文件????1143296??2010-12-28?16:52??資料\CHAPTER5_Spectrum_Sensing.ppt

?????文件?????612430??2010-12-28?16:53??資料\一種高效利用資源的協作感知方法.pdf

?????文件????1035128??2010-12-28?16:55??資料\基于能量檢測的頻譜感知方法.pdf

?????目錄??????????0??2010-12-28?16:39??資料\energydetection

?????目錄??????????0??2010-12-28?14:27??資料

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

??????????????3325736????????????????????8