基于MATLAB的π

基于MATLAB的4QPSK的仿真实现（刘海亮通信12011230440103）摘要在以前的数字蜂窝系统中，往往采用FSK、ASK、PSK等调制方式。随着数字蜂窝系统的发展，对调制和数字蜂窝系统的技术要求越来越高,许多优秀的调制技术应运而生，其中π/4QPSK技术是无线通信中比较突出的一种二进制调制方法本文概述了π/4QPSK的调制解调原理及其所实现的功能。并通过MATLAB编程对系统在相同条件下，对比了加噪声和没有噪声的情况进行了比较，并画出了它的眼图，及已调制信号的时域波形和频谱。并通过眼图分析其性能。关键字：π/4QPSK；调制解调原理；MATLAB编程；眼图1引言无线通信在现代社会中起着举足轻重的作用。从日常生活到航空航天，从工商业运作到军事领域，无线通信得到了越来越广泛的应用。现代数字调制技术的发展，使得传输速率和频谱的利用率进一步得到提高，功率更加节省。在相同的码元速率下，多进制系统的信息传输速率显然比二进制系统高，但信息速率的提高是以牺牲功率为代价的。显然增大码元宽度，就会增加码元的能量，同时也减少了由于信道特性引起的码间串扰等。恒包络调制适用于限带非线性信道中，能有效地防止非线性引起的幅频效应，节省功率，提高频谱的利用率。多进制调制和恒包络调制这两种技术结合在一起能取得更好的调制效果。为了使基带信号更好的利用信道进行传输，必须使代表的原始信号经过调制，而调制技术的好坏影响频谱资源的利用和通信性能的好坏。π/4QPSK在QPSK的基础上，克服了180°相位突变，信号轨迹不经过原点，有了更好的功率效率。2π/4QPSK调制解调原理2.1π/4QPSK调制原理π/4QPSK也是在QPSK的基础上发展起来的。不同的是这里把信号的相位平面分成π/4的八种相位，八种相位相间的又分成两个相位组。如图1所示。图中带‘*’的相位为一组，带‘。’的相位为另一组。规定π/4QPSK信号相位每隔Ts=2Tb秒必须从一个组跳变到另一个组。如果当前码元的信号相位等于‘*’组4个相位的一个，那么，下一个码元的信号相位只能变成‘。’组中4个相位的一个，反之也一样。这说明，图1中，符号不同的相位分别构成一个QPSK相量图，只是两者在相位上错开π/4。图1可以看出，在相邻码元之间，信号相位共有4种，即4,43不会出现π。图1π/4QPSK信号的矢量图π/4QPSK信号的表示式可以写成ttttsckkckkckckkcksin)1sin(cos)cos(sinsincoscos)cos(1（式1）式中k是当前码元信号相位k与前一码元信号k之差。所谓差分相位编码，就是利用相位差k来携带所需传输的信息。对于π/4QPSK信号来讲，对当前码元数据取值，k的取值范围为4,43四种取值，其编码规则为：AB=00对应π/4的相位，AB=10对应3π/4的相位，AB=11对应-3π/4的相位，AB=01对应-π/4的相位。显然令信号的初始相位为0，则当前码元的相位可能有0，π，4,43和π/2这八种初始相位。如图1所示。我们令)sin(),cos(11kkkkkkYX则有kYkkXkkkkkXksin1cos1sin)1sin(cos)1cos(（式2）kkkkkkkkXYYsincossin)cos(cos)sin(1111（式3）上式说明，kX和kY完全取决于前一码元的相位及前后码元的相位差，且kX和kY的取值只有0，1，21五种。因此π/4QPSK的信号包络不是恒定的。为了获取已调π/4QPSK信号，只要获取输入的当前码元数据所对应得kX和kY取值，再将其分别与相互正交载波信号ttccsin,cos相乘，并进行加法运算即可，其组成框图如图2所示图2π/4QPSK调制的原理框图图中的成形滤波器的目的，一是为了抑制已调信号的带外功率辐射，二是去除接受端的码元串扰。2.2π/4QPSK的解调原理π/4QPSK的常用解调方式也有相干解调和非相干解调两种方式。非相干解调主要分为鉴频检测和中频差分解调两种。在这里，我们采用的是中频差分解调。图3为中频差分解调框图。图3中频差分解调框图从图中可以看出，这种电路不需要另外的振荡器产生本地正交载波；经过延时的信号)cos()(1kcktts与两个支路的信号)cos(kct和)sin(kct分别相乘，即)cos()sin()cos()cos(1)(1)(kckckkckckttVttU（式4）经虑波和取样，可得)sin(21)cos(211)(1)(kkkkkkQI（式5）根据π/4QPSK信号基带信号的编码规则，1kkk只有四种取值。编码规则为：AB=00对应4的相位，AB=10对应43的相位，AB=11对应-3π/4的相位，AB=01对应4的相位。当k为4时，;0)(,0)(kQkI当k为43时，0)(,0)(kQkI；当k为43时,;0)(,0)(kQkI当k为4时，0)(,0)(kQkI。因此，可以对)(),(kQkI的符号进行判决，进而直接判决输出发射端的原始数据。判决规则为：0)(kI时判为0，否则判为1；0)(kQ时判为0，否则判为1。3π/4QPSK信号的调制解调建模根据前面的原理分析可知，产生调制信号并不难，直接按照框图来就是，但是解调在延式那块不能用直接延时来处理，必须采用Hilbert变换滤波器来实现。因为直接延时的前提是获取准确的载波频率，显然，这种处理的方法随着载波频率的估计误差，以及移动环境下的载波频率的偏移会带来较大误差。Hilbert滤波器是一个准确的相位延时系统，Hilbert滤波器是以个全通的滤波器，会产生90度的相移，其频响为：0,0,{)(jjjeh（式6）MATLAB中提供了现成的Hilbert滤波器设计函数firm()。4MATALAB仿真π/4QPSK信号的调制解调（差分解调）过程1）仿真4QPSK信号的产生方法；2）仿真差分解调π/4QPSK信号的过程；3）符号速率Rb=1Mbps;4）基带成形滤波器滚降系数α=0.8；5）采样速率为fs=8Rb;绘制π/4QPSK信号已调信号的频谱及时域波形，绘制解调后的同相反向支路眼图4.1π/4QPSK调制过程的程序清单ps=1*10^6;%码速率为1MHza=0.8;%成形滤波器系数B=(1+a)*ps;%中频信号处理带宽Fs=8*10^6;%采样速率fc=2*10^6;%载波频率N=2000;%仿真数据的长度t=0:1/Fs:(N*Fs/ps-1)/Fs;%产生长度为N,频率为fs的时间序列s=randint(N,1,4);%产生随机四进制数据作为原始数据%将绝对码变换为相对码%xk=ones(1,N);yk=ones(1,N);fori=2:Nifs(i)==0xk(i)=xk(i-1)*cos(pi/4)-yk(i-1)*sin(pi/4);yk(i)=yk(i-1)*cos(pi/4)+xk(i-1)*sin(pi/4);elseifs(i)==1xk(i)=xk(i-1)*cos(-pi/4)-yk(i-1)*sin(-pi/4);yk(i)=yk(i-1)*cos(-pi/4)+xk(i-1)*sin(-pi/4);elseifs(i)==2xk(i)=xk(i-1)*cos(3*pi/4)-yk(i-1)*sin(3*pi/4);yk(i)=yk(i-1)*cos(3*pi/4)+xk(i-1)*sin(3*pi/4);elseifs(i)==3xk(i)=xk(i-1)*cos(-3*pi/4)-yk(i-1)*sin(-3*pi/4);yk(i)=yk(i-1)*cos(-3*pi/4)+xk(i-1)*sin(-3*pi/4);endend%对相对码数据以Fs频率采样%Ads_i=upsample(xk,Fs/ps);Ads_q=upsample(yk,Fs/ps);%加噪声%SNR=20;%Ads_i=awgn(Ads_i,SNR);%Ads_q=awgn(Ads_q,SNR);%设计平方根升余弦滤波器%n_T=[-22];rate=Fs/ps;T=1;Shape_b=rcosfir(a,n_T,rate,T,'sqrt');%对采样后的数据进行升余弦滤波%rcos_Ads_i=filter(Shape_b,1,Ads_i);rcos_Ads_q=filter(Shape_b,1,Ads_q);%产生同相正交两路载频信号f0_i=cos(2*pi*fc*t);f0_q=sin(2*pi*fc*t);%产生PI/4_QPSK已调信号piqpsk=rcos_Ads_i.*f0_i-rcos_Ads_q.*f0_q;4.2Hilbert滤波器的设计程序清单%设计Hilbert滤波器及相同阶数的普通带通滤波器fpm=[00.25133.754]*10^6*2/Fs;%firpm函数的频段向量magpm=[001100];%firpm函数的幅值向量n=30;%滤波器阶数h_bpf=firpm(n,fpm,magpm,'hilbert');%Hilbert带通滤波器bpf=firpm(n,fpm,magpm);%普通带通滤波器3)4QPSK的解调过程程序清单%完成对PI/4_QPSK信号的Hilbert滤波及普通滤波piqpsk_i=filter(bpf,1,piqpsk);piqpsk_q=filter(h_bpf,1,piqpsk);%对普通带通滤波后的数据进行一个符号周期延时处理piqpsk_di=[zeros(1,Fs/ps),piqpsk_i(1:length(piqpsk_i)-Fs/ps)];%实现差分解调demod_mult_i=piqpsk_i.*piqpsk_di;demod_mult_q=piqpsk_q.*piqpsk_di;%对乘法运算后的同相正交支路滤波demod_i=filter(Shape_b,1,demod_mult_i);demod_q=filter(Shape_b,1,demod_mult_q);5结果分析图44QPSK信号的频谱及时域波形图图5同相支路眼图图6正交支路眼图02468101214-1-0.500.51时间(us)幅度(V)00.511.522.533.54-80-60-40-200频率(MHz)幅度(dB)PI/4QPSK信号频谱PI/4QPSK时域信号波形-0.500.5-0.3-0.2-0.100.10.20.30.4TimeAmplitudeEyeDiagram-0.500.5-0.4-0.3-0.2-0.100.10.20.30.4TimeAmplitudeEyeDiagram图7加噪音之后的正交支路眼图图8加噪音之后的正交支路眼图图9Hilbert滤波器-0.500.5-0.8-0.6-0.4-0.200.20.40.6TimeAmplitudeEyeDiagram-0.500.5-0.8-0.6-0.4-0.200.20.40.6TimeAmplitudeEyeDiagram00.20.40.60.81-3000-2000-100001000NormalizedFrequency(rad/sample)Phase(degrees)00.20.40.60.81-150-100-50050NormalizedFrequency(rad/sample)Magnitude(dB)由同相支路和正交支路眼图在加噪音和没加噪音的情况下对比可知，没加噪音的眼图眼孔张得很大，眼图端正，表示码间串扰小，扫描到的波形重叠在一块，波形良好，说明具有良好的解调性能。反之，加了噪音的眼图，眼孔很小，眼线