%%%%% 該仿真程序是BPSK在AWGN信道下的調制解調性能驗證實例 %%%%%%%%%%%%
%%%%% 該BPSK調制解調是忽略了載波調制與發送、接收基帶成形的 %%%%%%%%%
%%%%%% 根據程序框架補充完程序并仿真 %%%%%%%%%%%%%%%%%%%%%%%
clc
clear all;
EbN0dB=1:10;%設定信噪比
errLimit=50;%設定每個信噪比下循環終止統計錯誤數
N=100;      %每次循環仿真矢量長度
for ii=1:length(EbN0dB)
errCount=0;
totalN=0;
sigma2(ii)=1/(2*(10^(EbN0dB(ii)/10)));%AWGN信道噪聲方差與信噪比的關系
while errCount<errLimit
    source=rand(1,N)>0;                        %產生比特信源
    trans= 1-2*source;              %BPSK調制映射
    receiver= trans+sqrt(sigma2(ii))*randn(1,N);           %加信道噪聲,如果是衰落信道，信號模型為y=k*x+n
    demodata=  receiver<0;          %BPSK解調，根據調制映射得到判決條件
    temp=sum(abs(source-demodata));
    errCount=errCount+temp;
    totalN=totalN+N;
end
ber(ii)=errCount/totalN;
end
ber
%仿真性能曲線
figure(1)
semilogy(EbN0dB,ber,'-o');
hold on;
grid on;
%理論性能曲線
figure(2)
berthe=berawgn(EbN0dB,'psk',2,'nodiff');
semilogy(EbN0dB,berthe,'-*');