BPSK，QPSK调制与解调matlab代码.rar

BPSK，QPSK调制与解调matlab代码，过高斯信道的仿真结果与理论值对比，西电研究生估值大作业。

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