BPSK的matlab仿真原理加源代码

BPSK 误码率的 Matlab 仿真实验作者：黄准南京航空航天大学，江苏省 南京市 430062摘要：对于刚刚从事通信领域讨论性工作的讨论生来说,很有必要通过 Matlab 强大的仿真功能来强化理解通信信号处理的过程。本文正是通过对 BPSK 的 Matlab 仿真来提供一个比较好的实例.关键词：BPSK,高斯白噪声，带通滤波器，载波，抽样判决1 BPSK 的理论分析1。1 基本原理图 1.1 BPSK 相干解调原理框图输入信号 A（t）可表示为: (1。1-1）式中为矩形波，而表达式如下：对 A(t)进行傅里叶变换，得到 A(t）的功率谱密度（式中 Ts 为码元时间长度）：得到 A（t）的频谱图如 1。2 所示(图中 Tb 为码元时间长度)加法器s(t)相乘器BPF相乘器LPF抽样判决A(t)WGN定时脉冲A(t)图 1.2 矩形基带信号频谱A(t）与载波相乘：其频谱图如图 1.3 所示:图 1.3图中 fc 为载波频率，fs=1/Tb=Rb（码速率）,故带通滤波器的中心频率为 fc，带宽至少为 2fs.经过带通滤波器的高斯白噪声变成了窄带高斯信号,可以表示为：经过带通滤波器后的信号为 y（t）y(t）=s(t）+n（t）与载波相乘后再经过低通滤波器滤除高频重量得到 x(t)由 x(t)的功率谱密度图（图 1.4 所示)可知，最佳判决门限为 0；误码率 Pe=0。5＊erfc(sqrt（r))；式中信噪比 r=A^2/2＊σ^2,σ^2 为窄带高斯信号的方差,由于其均值为 0，故噪声功率就是其方差。图 1.4附录:clear allclcnum = 5000； ％码元数Rb = 1000； ％码速率 2*Rb〈=fstnum = 500; ％没一个码元用 tnum 个点代替fs = Rb*tnum; %采样频率 fs〉2＊(fp+Rb)fp = 100000; ％载波频率t = 0：1/fs：num/Rb-1/fs; ％时间 t 的长度与间隙。panjue1 = zeros(1,num）; ％基带传输的判决输出panjue2 = zeros（1，num）； %频带传输的判决输出pe1 = zeros（1,15）； ％基带传输的误码率矩阵pe2 = zeros(1，15）; ％频带传输的误码率矩阵a = 2*（fp—2*Rb)/fs; %带通滤波器的第一截止频率b = 2*(fp+2＊Rb)/fs; %带通滤波器的第二截止频率s = 1; %循环累加数％［a1,b1] = butter(2，[a,b］)； %采纳 butter 滤波器，得到 butter 滤波器的 H（z）的分子分母参数noise = randn（1，num＊tnum); %产生正态分布参数nwgnnoise = filter（a1,b1,noise); ％对 WGN 进行带通滤波stdnwgn = std（nwgnnoise); ％为了减少运算量。[a2，b2］ = butter（2，b)； %低通滤...