课设基于MATLAB的BPSK调制.完本

.wd摘 要 本次课程设计为基于MATLAB的BPSK原理电路仿真。本次课设着重介绍了算法的实现，并采用MATLAB程序仿真测试了BPSK过程中单极性不归零编码、脉冲成形、PSK调制、信号通过AWGN信道、载波恢复、解调、解码等过程。关键词：BPSK；2PSK；MATLAB；数字频带通信；目 录绪论11 BPSK数字调制原理21.1数字带通传输分类21.2 BPSK调制原理分析21.2.1调制原理分析21.2.2解调原理分析42 MATLAB软件62.1 MATLAB软件介绍63基于的MATLAB的BPSK调制分析和仿真73.1基于MATLAB的BPSK调制系统总述73.2编码过程的MATLAB实现83.3 BPSK调制的MATLAB的实现113.4 AWGN信道MATLAB的实现133.5载波恢复的MATLAB实现163.5.1接收端带通滤波器163.5.2通过FFT实现载波的直接频率估计203.5.3自适应迭代算法验证恢复频率223.6 BPSK解调254总结30附录31致谢32参考文献33绪论数字信号传输方式分为数字带通传输和数字基带传输。对于本次课程设计二进制相移键控BPSKBinary Phase Shift Key是利用载波的相位变化来传递数字信息，而振幅和频率保持不变的一种数字带通调制方式。在实际应用中，PSK具有恒包络特性，频带利用率比FSK高，在一样信噪比的条件下误码率也较低，同时PSK调制实现相对简单，故卫星通信，遥测遥控中用得最多的是BPSK方式调制。1 BPSK数字调制原理1.1数字带通传输分类数字带通传输中一般利用数字信号的离散取值特点通过开关键控载波，从而实现数字调制，比方对载波的振幅、频率和相位进展键控可获得振幅键控ASK、频移键控FSK和相移键控PSK。1.2 BPSK调制原理分析1.2.1调制原理分析相移键控是利用载波的相位变化来传递数字信息，而振幅和频率保持不变。在BPSK中，通常用初始相位0或分别表示二进制的0和1。因此，BPSK信号的时域表达式为： 1-1其中，表示第n个符号的绝对相位： 1-2 因此，式1-1可以改写为：1-3由于表示信号的两种码元的波形一样，极性相反，故BPSK信号一般可以表述为一个双极性全占空100% duty ratio矩形脉冲序列于一个正弦载波的相乘，即 1-4其中， 1-5这里,g(t)脉冲宽度为的单个矩形脉冲，而得统计特性为 1-6即发送二进制符号位0时取+1，取0相位；发送二进制符号为1时取-1，取相位。这种以载波的不同相位直接去表示响应的二进制数字信号的调制方式，称为二进制绝对相移调制BPSK。调制原理框图如图1.1所示：图1.1 BPSK调制原理图典型波形如图1.2所示:图1.2 发送码元为1 0 0 1 1的BPSK波形1.2.2解调原理分析BPSK信号的解调方法是相干解调法。由于PSK信号本身就是用相位传递信息的，所以在承受端必须利用信号的相位信息，故采用相干解调法来解调信号。BPSK解调原理框图如1.3所示：抽样判决低通滤波器相乘器 带通滤波器 图1.3 BPSK解调原理框图给出了一种BPSK信号想干承受设备原理框图。图中经过带通滤波的信号在相乘器与本地载波相乘，在相干解调中，如何得到与接收的BPSK信号同频同相的相干载波是关键，然后用低通滤波器去除高频分量，再进展抽样判决，判决器是按极性进展判决，得到最终的二进制信息。图中，假设相干载波的基准相位于BPSK信号的调制载波的基准相位一致。但是，由于在BPSK信号的载波恢复过程中存在的相位迷糊phase ambiguity，即恢复的本地载波与所需的相干载波可能同相，也可能反相，这种相位关系的不确定性将会造成解调出数字基带信号与发送的数字基带信号正好相反，即1变为0，0变为1，判决器输出数字信号全部出错。这种现象称为BPSK方式的倒现象。具体波形如图1.4所示图1.4 解调信号示意图2 MATLAB软件2.1 MATLAB软件介绍MATLAB 简介 MATLAB 的名称源自Matrix Laboratory,它是一种科学计算软件，专门 以矩阵的形式处理数据。MATLAB将高性能的数值计算和可视化集成在一起， 并提供了大量的内置函数，从而被广泛地应用于科学计算、控制系统、信息处理 等领域的分析、仿真和设计工作，而且利用MATLAB产品的开放式构造，可以 非常容易地对MATLAB的功能进展扩大，从而在不断深化对问题认识的同时，不断完善MATLAB产品以提高产品自身的竞争能力。在通信仿真中，欲将寄载消息的信号保质保量的传送给接收者，就必须对它们进展加工处理，对相应的部件进展分析和研究，而这种研讨过程中会用到大量的数学计算。它是以矩阵，数组为 基本处理对象，矩阵的维数不需要预先指定，且随时可变，具有丰富的矩阵运算功能。仿真功能我们采用MATLAB语言编程，以Windows和MATLAB为软件运行环境。将通信原理中的知识点用MATLAB来实现。通过灵活的菜单和界面设计使 某些参数可调。同时运用Simulink软件包，使波形在仿真的同时显示系统。更加有利于在课堂上使用。MATLAB界面：图2.1 MATLAB界面3基于的MATLAB的BPSK调制分析和仿真3.1基于MATLAB的BPSK调制系统总述基于MATLAB的BPSK调制解调系统流程框图如3.1所示：道信调制成形冲脉编码 文本 信息决判样抽载波恢复 解码解调图3.1 BPSK调制解调系统流程框图通信过程分析：整个发送与接收过程仿真了实际中的通信过程，输入一段信息如ph经过编码过程，变换为一串二进制字符，再加载矩形窗进展脉冲成形，生成原始的数字基带信号，为了发送信息，通过BPSK调制成数字带通信号发射，发射信号进入模拟的AWGN信道，被接收机所承受，BPSK解调采用的是相干解调，故需要进展载波恢复，恢复载波的频率，用已恢复的载波对承受信号进展相乘，再对其进展抽样判决，恢复最原始的二进制字符，再进展解码恢复原始的信息内容。仿真过程的要求：本次课设使用MATLAB进展BPSK通信过程仿真，要求不仅能将原始信息转换成二进制字符信息，解调后也能很好的恢复成原来的信息，过程中对信号观察不仅需要从时域上进展观察分析，也需要再频域上进展观察分析，以分析通过AWGN信道对信号的影响。3.2编码过程的MATLAB实现对于输入的文本信息，如果要进展调制发送，必然而且也可以将其转换成由0与1二进制字符，这个过程是信源编码过程。算法原理：对于输入的文本信息，在MATLAB中有对于的ASCII编码，MATLAB会自动的任意一行文档表示为一列ASCII码。再将十进制的ASCII码转换为二进制的的字符。对于数字形式的信息，必须转化成模拟形式，也就是脉冲成形过程，它将数字形式的每个字符转换成适宜的模拟脉冲，经过传输后，承受端可以从承受到的信号中恢复出原来的二进制字符。本次课设采用矩形窗加载，即为矩形脉冲。算法原理框图如3.2所示： 基带信号 脉冲成形二进制信息文本信息 图3.2 算法框图MATLAB代码：%编码和脉冲成形%采样率100decimalvalue=real(str);%使用自带real函数读取字符串的ASCII码binaryvalue=dec2base(decimalvalue,2,8);% dec2base函数转成8位二进制形式字符%因为使用dec2base函数转换出的二进制字符是使用char型的数组%进展将char型转换成int型二进制字符matrix1=size(binaryvalue);a=matrix1(1,1);b=matrix1(1,2);bv=;for i=1:1:a; for j=1:1:b; bv1=binaryvalue(i,j); bv1=bv1*1; bv1=bv1-48; bv=bv bv1; endend%为了能使编码后的二进制可以进展脉冲成形需要对其进展过采样M=100;N=length(bv);bsignal=zeros(1,N*M);bsignal(1:M:end)=bv;%使用矩形窗进展脉冲成形p=ones(1,M);s=conv(p,bsignal);%因为使用卷积，加载后的信号需要剔除后续多余局部k1=length(s);k1=k1-99;signal=zeros(1,k1);for i=1:1:k1; signal(i)=s(i);end%画图语句k=1;%第一幅图x=signal;Ts=1/100;plotspec(x,Ts,k);%调用自编子函数运行程序输入str=ph，输出二进制字符的基带信号signal，时域图如3.3所示： 图3.3 基带信号时域图基带信号频域图如3.4所示： 图3.4 基带信号频域图由MATLAB的工作空间可知ph转换成的二进制字符为：011000001101000.3.3 BPSK调制的MATLAB的实现数字基带信号一般不直接进展发射，需经过调制，本次课设要求使用BPSK调制方式进展调制。算法原理：在MATLAB中，可以采取脉冲成形的方式对原始信号进展调制，先将信号进展过采样插值，将1s内取100k个有一样幅值的点，以便使用载波进展脉冲成形，过采样后直接使用载波对信号进展相乘，以到达调制的目的。算法原理框图如3.5所示： 图3.5 算法原理框图MATLAB代码：%BPSK调制%采样率100kHZt2=0:(2*pi)/100000:2*pi-0.00001;%载波一周期内采样点数fc=300;%载波频率m1=;c1=;b1=;%进展过采样插值for n=1:1:length(bv); if bv(n)=0; m=-ones(1,100000); b=zeros(1,100000); else bv(n)=1; m=ones(1,100000); b=ones(1,100000); end c=cos(fc*t2);%生成相应载波%过采样后插值的信号m1=m1,m；c1=c1,c;b1=b1,b;endpsk=c1.*m1;%生成BPSK信号Ts=1/100000;%采样周期x=psk; %画图k=2;%第二幅图plotspec(x,Ts,k);%调用自编函数文件运行程序得到BPSK已调信号psk，时域图如3.6所示：图3.6 BPSK信号时域图局部BPSK的频域图如3.7所示：图3.7 BPSK信号频域图从图观察，在HZ处附近有值，说明已用频率300HZ的载波调制数字基带信号。3.4 AWGN信道MATLAB的实现加性高斯白噪声 AWGN(Additive White Gaussian Noise) 是最 基本的噪声与干扰模型。它的幅度分布服从高斯分布，而功率谱密度是均匀分布的，它意味着除了加性高斯白噪声外，r(t)与s(t)没有任何失真。即H(f)失真的加性噪声：叠加在信号上的一种噪声，通常记为n(t)，而且无论有无信号，噪声n(t)都是始终存在的。因此通常称它为加性噪声或者加性干扰。白噪声：噪声的功率谱密度在所有的频率上均为一常数，则称这样的噪声为白噪声。如果白噪声取值的概率分布服从高斯分布，则称这样的噪声为高斯白噪声。在MATLAB中有模拟AWGN信道可直接使用awgn函数，y = awgn(x,SNR)在信号x中参加高斯白噪声。信噪比SNR以dB为单位。x的强度假定为0dBW。如果x是复数，就参加复噪声。y = awgn(x,SNR,SIGPOWER)如果SIGPOWER是数值，则其代表以dBW为单位的信号强度；如果SIGPOWER为measured，则函数将在参加噪声之前测定信号强度。y = awgn(x,SNR,SIGPOWER,STATE)重置RANDN的状态。y = awgn(,POWERTYPE)指定SNR和SIGPOWER的单位。POWERTYPE可以是dB或linear。如果POWERTYPE是dB，那么SNR以dB为单位，而SIGPOWER以dBW为单位。如果POWERTYPE是linear，那么SNR作为比值来度量，而SIGPOWER以瓦特为单位。MATLAB代码：%模拟AWGN信道e2psk=awgn(psk,100);%调用awgn%画图x=e2psk;Ts=1/100000;%采样率k=3;plotspec(x,Ts,k);运行程序，经过AWGN信道后的信号e2psk，其时域图如3.8所示：图3.8经过AWGN信道后的信号e2psk信号局部放大图如3.9所示：图3.9 e2psk信号时域图局部可见信号经过AWGN信道后，受到噪声污染，时域上有失真。e2psk信号频谱图如3.10所示：图 3.10 e2psk信号频谱图300HZ左右处放大图如3.11所示： 图3.11 e2psk信号300Hz左右处频谱分析图可知，由于信号经过了AWGN信号，收到噪声干扰，频谱图在整个频域上有加性白噪声污染，都有不同程度的失真，不利于信号的解调。3.5载波恢复的MATLAB实现从接收端接收到信号后，为了进展相干解调，需要从本地恢复载波。对本次课设BPSK调制方式是一种没有抑制载波的调制方式，但是在通常情况下，载波很容易的被隐藏在接收信号中，而且需要很多额外的处理过程来恢复载波。对于刚接收到的信号可先通过FFT实现载波的频率估测，得出一个频率之后，再使用自适应算法的方式，准确定位和验证载波频率值。3.5.1接收端带通滤波器在接收器承受到调制信号后，由于调制信号被AWGN信道“污染，在整个频带上增加了加性高斯白噪声，故为了过滤信号需要在承受端加上一个带通滤波器，以滤除带外信号。本次课设调用remez函数以使用雷米兹算法，雷米兹算法提供了求解切比雪夫逼近问题的有效方法。切比雪夫滤波器是在通带或阻带上频率响应幅度等波纹波动的滤波器。MATLAB代码：%带通滤波器生成及其频率响应f1=500;%滤波器阶数fs=1000；%滤波器显示频率最大值的两倍ff=0 .51 .52 .68 .69 1;%归一化的截止频率fa=0 0 1 1 0 0;h=remez(f1,ff,fa);%生成带通滤波器figure4，freqz(h,1,512,fs);%绘出滤波器的频率响应r=filter(h,1,e2psk);%信号通过带通滤波器k=5;plotspec(r,Ts,k)带通滤波器的幅频响应如图3.12所示：图3.12 带通滤波器幅频响应带通滤波器的相频响应如图3.13所示：图3.13带通滤波器的相频响应分析可知在通带或阻带上频率响应幅度等波纹波动的切比雪夫滤波器，带通滤波器的频率响应可知带通滤波器的中心频率为300Hz，通带为260Hz至340Hz，群延时为179.85度。信号通过带通滤波器后时域图如3.14所示：图3.14 信号通过带通滤波器通过带通滤波后信号的局部图如3.15所示：图3.15 接收信号通过带通滤波局部可知经过带通滤波后，时域图中信号有一定成度的改善，波形突刺减少，一定程度的恢复成原来的波形。信号通过带通滤波后的幅频响应如3.16所示：图3.16 承受信号通过带通滤波后的幅频响应分析频域图，在频带中多余的加性白噪声得到一定的抑制，而很好的保存了载频附近的频率信息。3.5.2通过FFT实现载波的直接频率估计对于没有抑制载波的调制信号来说，使用FFT变换后通常载波频率幅度是最大值，可以通过检测已调信号的整个频率范围最大值来确定最大值。但是，一般情况下，由于已调信号经过信道后，频谱受到损伤，出现一定程度的失真，降低了用最大值来确定载波频率的方法的可靠性，此时需要使用另外的方法来恢复载波频率。载波恢复流程框图如3.17所示： 图3.17 载波恢复流程图如图3.17所示，一般可以在带通滤波器后使用非线性模块平方器，是承受信号，由基带信号和载波相乘组成，平方器的输出项为 3.1 又根据等式 ，改写3.1为 3.2将用其平均值和偏移量的和来表示 3.3因此 3.4以为中心的窄带带通滤波器只能让中的纯余弦项通过，而且还会抑制当中的直流成分，低频局部和向上边频了的也不能通过。带通滤波器的输出约为 3.5式中， 是相位偏移，是带通滤波器引入的。对于承受机来说，值是的，故可用来估算出载波的频率。 MATLAB代码：rsc=r; q=rsc.2; % 通过平方器fl=500; ff=0 .0100 .0104 .0136 .0140 1; fa=0 0 1 1 0 0; h=remez(fl,ff,fa)； % 生成带通滤波器rp=filter(h,1,q); % 再通过带通滤波fftrBPF=fft(rp); % 对rp信号进展FFT变换m,imax=max(abs(fftrBPF(1:end/2); % 找到最大值，返回最大值和最大值所在数组点ssf=(0:length(rp)/(Ts*length(rp); %建设储存频率的数组freqS=ssf(imax); %找出存储频率数组对应最大幅度的频率运行程序结果：freqS = 600由结果可知估测频率fc=300Hz。3.5.3自适应迭代算法验证恢复频率对于前面已经估测处了载波的频率，此时可以由freqS为假设值，利用自适应算法对恢复的载波频率进展再估测和验证。在通信中的许多问题能够被设计成最优化问题要解决这种问题需要三个 基本的步骤：1设定一个目标-选择一个目标函数2选择一个实现目标的方法-最大化或最大化这个目标函数3检验以确定这个方法能到达预期效果“设定一个目标函数通常包括寻找一个可以最大化或可以最小化的函数，定位最大和最小值为所要解决问题提供了有用的信息。这个函数的选择必须使它和它的导数的值能基于一些信息计算出来。有很多的方法可以用来实现最大化或最小化过程。有的时候是可以使用直接的方法，如，如果为了找到多项式到达最小值的点，只要直接找到它的导数并置零就可以得到解答。但是很多时候，不可能有如此简单的解决方法，递归方法自适应方法通常在信号运算中有很好的应用。自适应算法通常以最陡下降为目的，也可以称之为最陡下降法，最陡下降自适应算法从最初猜测的位置开场，算出从估计值开场哪个是最陡下降的方向，并沿着这条下降方向做出最新的估算，为了运用最陡下降法来最小化多项式，假设在时间处的的当前估算是可以得到的，表示为在时间的估算可以用下式：3.6式中，是很小正数，称为步长，而的梯度是当前点计算出来的，当重复计算，依次迭代运算，最终接近最小值对于承受信号为，设freqS为，合理的算法为尝试最小化下式：3.7利用梯度算法重写3.6式，得：3.8将平方项展开整理为：3.9此式即为其误差外表函数，可以用来分析性能函数。MATLAB代码：%自适应算法估测频率%预先恢复频率2fc=freqS=600Hz%本次算法初始值为300Hz，采用步长mu=0.003Ts=1/100000; time=5; t=0:Ts:time-Ts; % 建设时间向量fl=300;ff=0 .005 .01 1;fa=1 1 0 0;h=remez(fl,ff,fa); % 低通滤波生成mu=.003; % 步长fc=300; % freqS估测出频率fest=zeros(1,length(t); fest(1)=fc+3; % 起始运算的频率303Hzz=zeros(1,fl+1); % 进展循环运算 迭代运算有3.8式得出 for k=1:length(t)-1 z=(r(k)-cos(2*pi*fest(k)*t(k)*sin(2*pi*fest(k)*t(k), z(1:fl); update=fliplr(h)*z; fest(k+1)=fest(k)-mu*update; end%画图figure(6),plot(t,fest);title(Frequency Tracking via the Direct method)xlabel(time);ylabel(phase offset)； 运行程序后，自适应算法结果如3.18所示：图3.18 自适应算法结果图由图可知，在一段区域内由假定的303Hz趋近于300Hz，与之前直接频率估计的结果相符，故可知载波的频率为300Hz。3.6 BPSK解调对于BPSK解调，本次课设采用匹配滤波接收，匹配滤波器可使输出信噪比达最大，有利于判决抽样。设承受滤波器传输函数为，冲击函数为，滤波器输入码元的持续时间为匹配滤波器的特性可以用其冲激响应函数来描述：3.10由式可知匹配滤波器的冲激响应就是的镜像，但在时间轴上平移了。一个实际的匹配滤波器应该是物理可实现的，其冲激响应必须符合因果关系，在输入冲激脉冲参加前不应有冲激响应出现，即必须有：3.11既要求满足条件：3.12或满足条件：3.13式3.13的条件说明，承受滤波器输入端的信号码元在抽样时刻之后必须为零。一般不希望在码元完毕之后很久才抽样，故通常选择在码元末尾抽样，即选。故匹配滤波器冲激响应可以写为： 3.14这时，假设匹配滤波器的输入电压为，则输出信号码元的波形，可以求出为：3.15式3.15说明，匹配滤波器输出信号码元波形是输入信号码元波形的自相关函数。为任意一常数，通常取。MATLAB 代码：%匹配滤波tm=0:1/100000:1-(1/100000);ps=cos(2*pi*300*tm);dm=conv(ps,r);%画图figure(7),plot(dm);grid on运行程序图如3.19所示：图3.19 匹配滤波后信号信号局部放大图如3.20所示：图3.20 信号通过匹配滤波后局部方法图对信号匹配滤波后，则可以进展抽样判决以恢复原始二进制信息。对于匹配滤波后的信号采样点为nTs。MATLAB 代码：%抽样赋值1和-1dmbv=;for i=1:1:length(r)/100000; dmbv1=dm(i*100000); dmbv=dmbv dmbv1;end%将双极性不归零码变为单极性不归零码，赋值0和1for i=1:1:length(dmbv); if dmbv(i)0;dmbv(i)=0; else dmbv(i)=1; endend%画图tm1=0:1:length(dmbv)-1;figure(8),stem(tm1,dmbv);grid on;axis(-0.5 15.5 -1 2);xlabel(时间);ylabel(幅度)title(解调后的单极性不归零二进制信号)运行程序后得结果如3.21所示：图3.21解调后的单极性不归零二进制信号由图可读出恢复数据为0111000001101000MATLAB代码：%解码dmbv1=dmbv(1:end/2);%8位dmbv2=dmbv(end/2+1:end);dmstr1=num2str(dmbv1);%数值型变成字符型dmstr2=num2str(dmbv2);dms1=bin2dec(dmstr1);%由二进制转成十进制dms2=bin2dec(dmstr2);运行程序后，由工作空间可知都dms1为112，dms2为104。再有代码：char112；char104；得结果为：ph故解调出了原始的文本信息。4总结本次课设基于MATLAB仿真了BPSK调制方式的通信传输，发送端将文本信息转换成二进制基带信号，进展BPSK调制发送，信道为AWGN信道产生加性高斯白噪声，接收端接收信号，进展载波恢复，将恢复的载波对承受信号进展匹配滤波，并解调，再进展抽样判决，恢复出原始的二进制字符，在解码成文本信息。本次课设进展了相对实际的BPSK调制通信仿真，初步验证了BPSK调制的可行性，并对通信系统中遇到的问题进展了一定的探索研究。通过本次课设，使我对已掌握专业知识有了更深一层次的把握，也学到了很多新的知识并很好的掌握了这些知识，同时整个过程充分理解了理论指导实践，实践验证理论这句话。我不仅书本上得理论知识有了进一步的理解，也同时熟练的掌握了MATLAB的仿真，使自己脑海中的理论不再只是理论而得到了实践的验证。这次课设使我能很好的综合运用自己所学的知识解决一些问题，在面对自己不懂的问题时能逐步分析并最后解决这些问题，对我以后的学习起到了抛砖引玉的作用，促使以后能不断的进步。附录%为调用的plotspec函数function plotspec(x,Ts,k)%函数功能是画出信号的时域图和频域图%x为输入信号，Ts为抽样率，k为图的序号 N=length(x);t=Ts*(1:N);ssf=(-N/2:N/2-1)/(Ts*N);fx=fft(x(1:N);fxs=fftshift(fx);figure(k)subplot(2,1,1),plot(t,x);grid on,xlabel(秒);ylabel(幅度);subplot(2,1,2),plot(ssf,abs(fxs);grid on,xlabel(频率);ylabel(幅度)end致 谢感谢这次课程设计的指导教师李宏达教师、刘庆泉教师和伍彩云教师，也感谢指导过我的教研室教师和同学们。参考文献1 C.Richard Johnson Jr.等著.潘甦译.软件无线电 机械工业出版社 2 樊昌信等编著.通信原理.国防工业出版社3 郭文彬等编.通信原理-基于MATLAB的计算机仿真.北京邮电大学出版社4 赵鸿图等编.通信原理MATLAB仿真教程.人民邮电出版社