[整理]BPSK扩频误码率matlab仿真.pdf

------------- 实验报告 实验目的 1. 掌握 BPSK 信号调制、相干解调方法； 2. 掌握扩频、解扩的方法； 3. 掌握 BPSK 扩频信号误码率计算 实验内容 1. BPSK 信号的调制； 2. BPSK 信号扩频； 3. 不同信噪比环境下 BPSK 扩频信号误码率计算，并与理论误码率曲线对比 实验原理 BPSK 信号调制原理 1. 系统原理 双极性码脉冲成形BPFLPF抽样判决抽样时钟载波相干载波高斯白噪声图 5 BPSK 调制系统原理框图 BPSK 调制系统的原理框图如图 1 所示，其中脉冲成形的作用是抑制旁瓣，减少邻道干扰，通常选用升余弦滤波器；加性高斯白噪声模拟信道特性，这是一种简单的模拟；带通滤波器 BPF 可以滤除有效信号频带以外的噪声，提高信噪比； 在实际通信系统中相干载波需要使用锁相环从接收到的已调信号中恢复， 这一过程增加了系统的复杂度， 同时恢复的载波可能与调制时的载波存在 180 度的相位偏差，即 180°相位反转问题，这使得 BPSK 系统在实际中无法使用；低通------------- ------------- 滤波器 LPF 用于滤除高频分量，提高信噪比；抽样判决所需的同步时钟需要从接收到的信号中恢复，即码元同步，判决门限跟码元的统计特性有关，但一般情况下都为 0。

2. 参数要求 码元速率 1000 波特，载波频率 4KHz，采样频率为 16KHz BPSK 信号扩频原理 通信理论和通信技术的研究， 是围绕着通信系统的有效性和可靠性这两个基本问题展开的， 所以有效性和可靠性是设计和评价一个通信系统的主要性能指标 通信系统的有效性， 是指通信系统传输信息效率的高低 这个问题是讨论怎样以最合理、最经济的方法传输最大数量的信息在模拟通信系统中，多路复用技术可提高系统的有效性显然，信道复用程度越高，系统传输信息的有效性就越好在数字通信系统中，由于传输的是数字信号，因此传输的有效性是用传输速率来衡量的 通信系统的可靠性， 是指通信系统可靠地传输信息 由于信息在传输过程中受到干扰， 收到的信息与发出的信息并不完全相同 可靠性就是用来衡量收到信息与发出信息的符合程度 因此， 可靠性决定于系统抵抗干扰的性能， 也就是说，通信系统的可靠性决定于通信系统的抗干扰性能在模拟通信系统中，传输的可靠性是用整个系统的输出信噪比来衡量的在数字通信系统中，传输的可靠性是用信息传输的差错率来描述的 扩展频谱通信由于具有很强的抗干扰能力， 首先在军用通信系统中得到了应用近年来，扩展频谱通信技术的理论和应用发展非常迅速，在民用通信系统中也得到了广泛的应用。

扩频通信是扩展频谱通信的简称我们知道，频谱是电信号的频域描述承载各种信息（如语音、图象、数据等）的信号一般都是以时域来表示的，即信息信号可表示为一个时间的函数)(tf 信号的时域表示式)(tf可以用傅立叶变换得到其频域表示式)( fF频域和时域的关系由式(1-1)确定： tetffFftjd)()(π2 ------------- ------------- fefFtfftjd)()(π2 (1-1) 函数)(tf的傅立叶变换存在的充分条件是)(tf满足狄里赫莱(Dirichlet)条件，或在区间(-∞，+∞)内绝对可积，即ttfd)(必须为有限值 扩展频谱通信系统是指待传输信息信号的频谱用某个特定的扩频函数 （与待传输的信息信号)(tf无关）扩展后成为宽频带信号，然后送入信道中传输；在接收端再利用相应的技术或手段将其扩展了的频谱压缩， 恢复为原来待传输信息信号的带宽， 从而到达传输信息目的的通信系统 也就是说在传输同样信息信号时所需要的射频带宽， 远远超过被传输信息信号所必需的最小的带宽 扩展频谱后射频信号的带宽至少是信息信号带宽的几百倍、 几千倍甚至几万倍。

信息已不再是决定射频信号带宽的一个重要因素，射频信号的带宽主要由扩频函数来决定 由此可见，扩频通信系统有以下两个特点： (1) 传输信号的带宽远远大于被传输的原始信息信号的带宽； (2) 传输信号的带宽主要由扩频函数决定，此扩频函数通常是伪随机（伪噪声）编码信号 以上两个特点有时也称为判断扩频通信系统的准则 扩频通信系统最大的特点是其具有很强的抗人为干扰、抗窄带干扰、抗多径干扰的能力 这里我们先定性地说明一下扩频通信系统具有抗干扰能力的理论依据 扩频通信的基本理论根据是信息理论中香农(C·E·Shannon)的信道容量公式 NSBC1log2 (1-2) 式中： C——信道容量，b/s； B——信道带宽，Hz； S——信号功率，W； N——噪声功率，W 香农公式表明了一个信道无差错地传输信息的能力同存在于信道中的信噪比以及用于传输信息的信道带宽之间的关系 ------------- ------------- 令 C 是希望具有的信道容量，即要求的信息速率，对(1-2)式进行变换 NSBC1ln44. 1 (1-3) 对于干扰环境中的典型情况，当1NS时，用幂级数展开(1-3)式，并略去高次项得 NSBC44. 1 (1-4) 或 SNCB7 . 0 (1-5) 由式(1-4)和(1-5)可看出，对于任意给定的噪声信号功率比SN /，只要增加用于传输信息的带宽B，就可以增加在信道中无差错地传输信息的速率 C。

或者说在信道中当传输系统的信号噪声功率比NS /下降时，可以用增加系统传输带宽 B 的办法来保持信道容量 C 不变或者说对于任意给定的信号噪声功率比NS /，可以用增大系统的传输带宽来获得较低的信息差错率 若100/SN(20dB)，kb/s3C，则当kHz21031007 . 0B时，就可以正常的传送信息，进行可靠的通信了 这就说明了增加信道带宽 B，可以在低的信噪比的情况下，信道仍可在相同的容量下传送信息甚至在信号被噪声淹没的情况下，只要相应的增加信号带宽也能保持可靠的通信如系统工作在干扰噪声比信号大 100 倍的信道上，信息速率 R＝C＝3kb/s，则信息必须在kHz210B带宽下传输，才能保证可靠的通信 扩频通信系统正是利用这一原理， 用高速率的扩频码来扩展待传输信息信号带宽的手段， 来达到提高系统抗干扰能力的目的 扩频通信系统的带宽比常规通信系统的带宽大几百倍乃至几万倍， 所以在相同信息传输速率和相同信号功率的条件下，具有较强的抗干扰的能力 香农在其文章中指出， 在高斯噪声的干扰情况下， 在受限平均功率的信道上，实现有效和可靠通信的最佳信号是具有白噪声统计特性的信号。

这是因为高斯白------------- ------------- 噪声信号具有理想的自相关特性，其功率谱密度函数为 2)(0NfS -∞< f <∞ (1-6) 对应的自相关函数为 )(2d)()(0π2τδNfefSτRfτj (1-7) 其中：为时延，)(定义为 000)(τττδ (1-8) 白噪声的自相关函数具有)(函数的特点，说明它具有尖锐的自相关特性但是对于白噪声信号的产生、 加工和复制， 迄今为止仍存在着许多技术问题和困难 然而人们已经找到了一些易于产生又便于加工和控制的伪噪声码序列， 它们的统计特性近似于或逼近于高斯白噪声的统计特性 伪噪声序列的理论在本书以后的章节中要专门讲述， 这里仅简略引用其统计特性，借以说明扩频通信系统的实质 通常伪噪声序列是一周期序列假设某种伪噪声序列的周期(长度)为 N，且码元ic都是二元域1, 1上的元素一个周期（或称长度）为 N，码元为ic的伪噪声二元序列 ic的归一化自相关函是一周期为 N 的周期函数，可以表示为 kckNRR)()()( (1-9) 其中)(cR为伪噪声二元序列 ic一个周期内的表示式 01011)(1τNτccNτRNiiic (1-10) 式中0，1，2，3，…N。

当伪噪声序列周期(长度)N 取足够长或 N→∞时，式(1-10)可简化为 00101)(τNττRc (1-11) 比较式(1-7)和式(1-11)，看出它们比较接近，当序列周期(长度)足够长时，式(1-11)就逼近式(1-7)式(1-10)是自相关函数归一化的形式，乘周期 N 后就是一般表达式，在一般表达式中NR)0(）所以伪噪声序列具有和白噪声相类似的统计特性， 也就是说它很接近于高斯信道要求的最佳信号形式 因此用伪噪声码扩展待传输信息信号频谱的扩频通信系统，优于常规通信系统 ------------- ------------- 哈尔凯维奇(А·А·Харкевич)早在上世纪 50 年代，就已从理论上证明：要克服多径衰落干扰的影响， 信道中传输的最佳信号形式应该是具有白噪声统计特性的信号形式 采用伪噪声码的扩频函数很接近白噪声的统计特性， 因而扩频通信系统又具有抗多径干扰的能力 下面我们以直接序列扩频通信系统为例，来研究扩频通信系统的基本原理图 1-1 给出了直接序列扩频通信系统的简化原理方框图 由信源产生的信息流 na通过编码器变换为二进制数字信号)(td。

二进制数字信号中所包含的两个符号的先验概率相同，均为2/1，且两个符号相互独立，其波形图如图1-2(a)所示，二进制数字信号)(td与一个高速率的二进制伪噪声码)(tc的波形（如图 1-2(b)所示，伪噪声码作为系统的扩频码序列）相乘，得到如图 1-2(c)所示的复合信号)()(tctd，这就扩展了传输信号的带宽一般伪噪声码的速率ccTR/1是 Mb/s 的量级，有的甚至达到几百 Mb/s而待传输的信息流 na经编码器编码后的二进制数字信号的码速率bbTR/1较低，如数字话音信号一般为16 kb/s ~32kb/s ，这就扩展了传输信号的带宽 图 1 扩展频谱通信系统模型 (a) 发射系统；(b) 接收系统 频谱扩展后的复合信号)()(tctd对载波)π2cos(0tf（0f为载波频率） 进行调制（直接序列扩频一般采用PSK 调制） ， 然后通过发射机和天线送入信道中传输发射机输出的扩频信号用)(ts表示，其示意图如图 1-2(d)所示扩频信号)(ts的带宽取决于伪噪声码)(tc的码速率cR在 PSK 调制的情况下，射频信号的带宽等于伪噪声码速率的 2 倍， 即cRR2RF， 而几乎与数字信号)(td的码速率无关。

s(t) 数据源 编码器 m 序列 发生器 发射机 射频 振荡器 {an} d(t) c(t) cos(2f0t+) 中频 滤波器 VCO m 序列 发生器 射频 滤波器 R(t) r(t) v(t) 至数据 检测器 fˆ ˆdTˆ )ˆ(drTtc2cos[2(f0+fIF+dfˆ)t+] (b) (a) ------------- ------------- 以上对待传输信号)(td的处理过程就是对信号)(td的频谱进行扩展的过程 经过上述过程的处理，达到了对)(td扩展频谱的目的 图 2 理想扩展频谱系统波形示意图 在接收端用一个和发射端同步的参考伪噪声码)ˆ(drTtc所调制的本地参考振荡信号] ˆ)ˆ(π2cos[2IF0tfffd（IFf为中频频率），与接收到的)(ts进行相关处理相关处理是将两个信号相乘，然后求其数学期望（均值），或求两个信号瞬时值相乘的积分 当两个信号完全相同时（或相关性很好），得到最大的相关峰值，经数据检测器恢复出发射端的信号)(td 若信道中存在着干扰，这些干扰包括窄带干扰、人为瞄准式干扰、 单频干扰、 多径干扰和码分多址干扰等等， 它们和有用信号)(1ts同时进入接收机，如图 1-3(a)所示。

图 1-3 中，cR为伪噪声码速率，0f为载波频率，IFf为中频频率 图 3 扩频接收机中各点信号的频谱示意图 (a) 接收机输入；(b) 混频器输出；(c) 中频滤波器输出 (a) (b) (c) 窄带干扰 有用信号 多径干扰 白噪声 f f0 BRF=2Rc 窄带干扰 有用信号 多径干扰 白噪声 Bb=2Rb f fIF 2Rc 4Rc 窄带干扰 有用信号 多径干扰 白噪声 f fIF Bb=2Rb Tb (a) d(t) +1 -1 Tb (e) d’ (t) +1 -1 s(t)=Ad(t)c(t)cos[2f0t+(t)](d) s(t) (b) c(t) +1 -1 Tc (c) d(t)c(t) +1 -1 A -A ------------- ------------- 由于窄带噪声和多径干扰与本地参考扩频信号不相关， 所以在进行相关处理时被削弱，实际上干扰信号和本地参考扩频信号相关处理后，其频带被扩展，也就是干扰信号的能量被扩展到整个传输频带之内，降低了干扰信号的电平( 单位频率内的能量或功率) ，如图 1-3(b)所示。

由于有用信号和本地参考扩频信号有良好的相关性，在通过相关处理后被压缩到带宽为bbRB2的频带内，因为相关器后的中频滤波器通频带很窄，通常为bbRB2，所以中频滤波器只输出被基带信号)(td调制的中频信号和落在滤波器通频带内的那部分干扰信号和噪声，而绝大部分的干扰信号和噪声的能量 （功率）被中频滤波器滤除， 这样就大大地改善了系统的输出信噪比，如图 1-3(c)所示关于这一特性，将在扩频通信系统的性能分析一章中作进一步分析为了对扩频通信系统的这一特性有一初步了解，我们以解扩前后信号功率谱密度示意图来说明这一问题 假设有用信号的功率为01PP ，码分多址干扰信号的功率02PP ，多径干扰信号的功率03PP ，其他进入接收机的干扰和噪声信号功率0PN 再假设所有信号的功率谱是均匀分布在cRB2RF的带宽之内解扩前的信号功率谱见图 1-4 中的(a)，图中各部分的面积均为0P解扩后的信号功率谱见图 1-4 中的(b)，各部分的面积保持不变通过相关解扩后，有用信号的频带被压缩在很窄的带宽内，能无失真的通过中频滤波器（滤波器的带宽为bbRB2）其他信号和本地参考扩频码无关， 频带没有被压缩反而被展宽了， 进入中频滤波器的能量很少，大部分能量落在中频滤波器的通频带之外， 被中频滤波器滤除了。

我们可以定性的看出，解扩前后的信噪比发生了显著的改变 图解扩前后信号功率谱密度示意图(a) 解扩前；(b) 解扩后 BPSK 信号解调原理 BPSK 信号的解调方法是相干解调法 由于 PSK 信号本身就是利用相位传递2Rc 2Rc Bb (a) (b) 有用信号 多址信号 多径信号 噪声 有 用 信 号 多址信号 多径信号 噪声 4Rc f0 fIF ------------- ------------- 信息的，所以在接收端必须利用信号的相位信息来解调信号图 5 中给出了一种2PSK 信号相干接收设备的原理框图图中经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决判决器是按极性来判决的即正抽样值判为 1，负抽样值判为 0 带通滤波器相乘器低通滤波器抽样判决器( )BPSKetacde定时脉冲bcos()cw t输出 图 5 BPSK 解调系统原理框图 BPSK 信号误码率 在 AWGN 信道下 BPSK 信号相干解调的理论误码率为：1()2ePerfcr，其中r为信噪比222nar 在大信噪比(1)r条件下，上式可近似为：12rePer 实验结果与分析 BPSK 基带信号产生 产生二进制随机码流作为传输信息，假定每个码元中有 16 个码片，产生BPSK信号基带信号。

Walsh 码产生 利用 Hadamard 矩阵进行迭代运算产生 16 阶 Walsh 矩阵， 得到所需的 Walsh码本次实验中使用的是[ -1,1,1,-1,-1,1,1,-1,-1,1,1,-1, -1,1,1,-1]Walsh 码 ------------- ------------- 扩频 将源码与扩频因子模二加，即可得到扩频之后的码流 解扩 将接受到的码字与扩频因子相乘，即可恢复出源码再进行抽样判决可得到最终信息 误码率计算 在扩频码流中添加不同的信噪比的噪声，进行调制接受接收后进行解调、解扩并且计算误码率与理论结果对比，如图 6 所示 可以看出在进行 16 阶扩频后，仿真误码率与理论误码率曲线相差很小，结果较为理想 图 6 仿真误码率与理论误码率 ------------- ------------- 实验代码 walsh 矩阵产生代码如下： % 产生 Walsh 函数通用函数 % 参数 N 表示 Walsh 函数阶数,当 N 不是 2 的幂时,通过向无穷大取整使得所得Walsh 阶数为 2 的幂 function [walsh]=walsh(N) M=ceil(log2(N)); hn=0; for i=1:M h2n=[hn,hn;hn,~hn]; hn=h2n; end walsh=hn; 扩频、误码率计算程序如下： %% walsh 码扩频、解扩、误码率仿真等 clear all;clc;close all; %383 %% 产生随机比特流作为传输信息 bit_len=100000; bits=randint(1,bit_len); %% 初始各码元对应的基带波形 s(t) RB=1000; %符号速率 Fs=16000; %采用率 Fd=Fs/RB; %每个码元的采样点数 %-------------16 阶 walsh 码产生--------------% B=walsh(16); %定义 walsh 函数产生 16 阶 walsh 矩阵 B wal=B(4,:); %walsh 码调制 wal=sign(wal-0.5); %-------------walsh 码与原码相加-------------% sig_len=Fd*bit_len; s=zeros(1,sig_len); %基带信号 for i=1:bit_len if bits(i)==1 s((i-1)*Fd+1:i*Fd)=ones(1,Fd).*wal; else s((i-1)*Fd+1:i*Fd)=-1*ones(1,Fd).*wal; end end %% BPSK 调制 ------------- ------------- wumalv=zeros(1,length(-10:0.5:13)); lilun=zeros(1,length(-10:0.5:13)); for j=-10:0.5:13 sig=awgn(s,j); %解调信号 z=sign(sig); %抽样判决 解扩 for i=1:bit_len z((i-1)*Fd+1:i*Fd)= z((i-1)*Fd+1:i*Fd).*wal; end z(z>=0)=1; %判决 z(z<0)=0; z=z(Fd/2:Fd:sig_len); %抽样 wuma=z-bits; c=find(wuma~=0); wumalv((j+10)/0.5+1)=length(c)/bit_len; k=sqrt(10^(j/10)); lilun((j+10)/0.5+1)=0.5*erfc(k); end %%由于噪声不是加载在调制信号上，而是加载基带信号上， %因此抽样判决的信噪比比带通滤波器后的信噪比大三个 dB figure semilogy(-13:0.5:10,wumalv,'r*');hold on; semilogy(-10:0.5:10,lilun(1:length(-10:0.5:10)));hold off; grid on;title('误码率');legend('仿真误码率','理论误码率') 。