BPSK与QPSK调制解调系统的Matlab仿真程序

发送滤波器与接收滤波器均为根升余弦滚降滤波器； 信道噪声为加性高斯白噪声； 绘制BPSK与QPSK调制下的误码率与信噪比曲线图，并与理论曲线进行对比。

function [ber]=BPSK（K）
fc=2;fs=16;fd=1;
snr=K-10*log10（fs/2）;   %SNR与Eb/N0之间的转换，具体详见实验报告
alpha=1;
N=1000000;              %取一百万点，使计算更精确
ts=1/fs;Ts=1;
t=0:ts:N*Ts-ts; 
t=t‘;
B=randint（1N）;
sym=[-1;1];
A=sym（B+1）;
H=rcosine（fdfs‘sqrt‘alpha）;        %根升余弦滤波器
A1=upsample（Afs/fd）;
A1=filter（H1A1）;
X=A1.*cos（2*pi*fc*t）;
X=awgn（Xsnr‘measured‘）;             %信道中加入高斯白噪声
R=X.*cos（2*pi*fc*t）;
R=filter（H1R）;
R2=downsample（Rfs/fd）;               %下采样，符号判决
R2=R2>0;
delay=（length（H）-1）/fs;               %计算通过滤波器的延时
for i=1:1:N-delay                        
    Sr（i）=R2（i+delay）;
end
%{
for i=N-delay+1:N     %由于延时，有一部分没被接收，此处全置0，会有一定误差
    Sr（i）=0;
end
%}
errbit=sum（Sr~=B（1:N-delay））;
ber=errbit/（N-delay）;          %误比特率

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         847  2013-12-08 12:32  BPSK.m
     文件         271  2013-12-01 14:09  mybpsk.m
     文件         245  2013-12-01 13:34  myqpsk.m
     文件        1490  2013-12-01 14:15  QPSK.m