bpsk_test.m

来自「bpsk 调制解调 该仿真程序是BPSK在AWGN信道下的调制解调性能验证实例」· M 代码 · 共 37 行

%%%%% 该仿真程序是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,'-*');