clc
clear?all
%產生信源s
ber=zeros（1011）;%對誤碼率進行初始化
for?t=1:10%蒙特卡洛循環次數為10
????num=10000;%信號點數為1000000
%產生信源s
s=randint（1num）;%用randint函數產生01比特流
%將信源進行BPSK調制信息存儲在x中
x=[];
for?i=1:length（s）
if（s（i）==0）
????x（i）=-1;
else
?????x（i）=1;
end
end
%將信號通過高斯信道，信噪比為0：10dB
?snr=0:10;
?for?i=1:length（snr）
N（i）=power（10-snr（i）/10）/2;%計算噪聲功率
n=sqrt（N（i））*randn（1length（s））;
r=x+n;
%?r=awgn（xsnr（i）‘measured‘）;
%采用最大似然檢測算法對接收信號r進行判決
for?j=1:length（s）
if（r（j）<0）
????y（j）=0;
else
????y（j）=1;
end
end
%計算仿真誤碼率
for?j=1:length（s）
???if（s（j）~=y（j））?;
???????ber（ti）=ber（ti）+1;
???end
end
ber（ti）=ber（ti）/num;
?end
end
bern=sum（ber）;
bern=bern/10;
snrn=power（10snr/10）;
q=qfunc（sqrt（2*snrn））;%計算理論誤碼率
%畫圖
semilogy（snrbern‘-‘snrq‘+‘）
grid?on
xlabel（‘Es/No（dB）‘）;?
ylabel（‘誤碼率‘）;?
title（‘BPSK調制信號在AWGN下的蒙特卡羅仿真‘）;?
legend（‘仿真誤碼率‘‘理論誤碼率‘）;?

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

?????文件??????34842??2017-07-11?17:17??大作業\1.fig

?????文件??????28361??2017-06-20?19:04??大作業\4psk.fig

?????文件??????31723??2017-07-10?12:14??大作業\bpsk.fig

?????文件???????1017??2017-11-16?10:26??大作業\bpsk.m

?????文件???????1257??2017-06-22?15:29??大作業\fpsk.m

?????文件??????28472??2017-06-22?14:47??大作業\qpsk.fig

?????文件???????1711??2017-11-16?10:29??大作業\qpsk.m

?????文件??????30843??2017-06-26?22:07??大作業\three.fig

?????文件????????437??2017-07-10?11:59??大作業\three.m

?????文件????????250??2017-10-18?21:01??大作業\LDPC.m

?????目錄??????????0??2017-07-11?17:17??大作業

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

???????????????158913????????????????????11