BPSK和QPSK调制解调原理及MATLAB程序

1、PSK调制方式PSK®理介22(以2-PSK为例)移相键控(PSK)又称为数字相位调制,二进制移相键控记作2PSK绝对相移是利用载波的相位(指初相)直接表示数字信号的相移方式。二进制相移键控中，通常用相位0和冗来分别表示“0”或“1"。2PSK已调信号的时域表达式为s2psk(t)=s(t)cosct,2PSK移相键控中的基带信号与频移键控和幅度键控是有区别的，频移键控和幅度键控为单极性非归零矩形脉冲序列，移相键控为为双极性数字基带信号，就模拟调制法而言，与产生2ASK信号的方法比较，只是对s(t)要求不同，因此2PSK信号可以看作是双极性基带信号作用下的DSB调幅信号。在

2、二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控(2PSK)信号。通常用已调信号载波的0°和1800分别表示二进制数字基带信号的1和0。二进制移相键控信号的时域表达式为e2PSK(t)=ang(t-nTs)coswctn其中，an与2ASKf口2FSK寸的不同，在2PSKS制中，an应选择双极性。l1,发送概率为Pan=YM,发送概率为1-P若g(t)是脉宽为Ts,高度为1的矩形脉冲时，则有Lcosct,发送概率为Pe2PSK(t)=V-coLct,发送概率为1-P由上式可看出，当发送二进制符号1时，已调信号e2PSK(t)取0°相位，

3、发送二进制符号0时，e2PSK(t)取180°相位。若用小n表示第n个符号的绝对相位，则有.00,发送1符号小n=Y18。10,发送0符号由于在2PSK信号的载波恢复过程中存在着180°的相位模糊，所以2PSK信号的相干解调存在随机的“倒冗”现象，从而使得2PSK方式在实际中很少采用。为了解决2PSK信号解调过程的反向工作问题，提出了二进制差分相位键控(2DPSK),这里不再详述。2-PSK调制解调二进制移相键控信号的调制原理：如图9所示。其中图(a)是采用模拟调制的方法产生2PSK信号,图(b)是采用数字键控的方法产生2PSK信号。解调器原理：如图10所示。2PSK信号的

4、解调通常都是采用相干解调，在相干解调过程中需要用到与接收的2PSK信号同频同相的相干载波。双极性不归零(a)模拟调制产生2PSK信号开关电路(b)数字键控的方法产生2PSK信号图92PSK信号的调制原理图图102PSK信号的解调原理图11100100_v_v_图112PSK信号相干解调各点时间波形就键控法来说，用数字基带信号S(t)控制开关电路，选择不同相位的载波输出,这时s(t)为单极性NRZIE双极性NRZ脉冲序列信号均可。当基带信号为0时候，连通开关0,产生无差别的载波，当所发出的信号为1时，既连通开关口改变载波的相位。在移相键控中还有一种差分移相键，他和普通的移相键控区别在与，差分移相

5、键只有在当前传输的码元和上次传输的码元产生差别时才会产生相位的变化。移相键控相对与幅度键控和移频键控有着更好的抗干扰性，也更适合于在信道中传输。QPSKS制QPSK号可以看作两个载波正交2PSK言号的合成。用调相法产生QPS喇制器框图如图12所示，QPSK勺调制器可以看作是由两个BPSKM制器构成，输入的串行二进制信息序列经过串并变换，变成两路速率减半的序列，电平发生器分别产生双极性的二电平信号I(t)和Q(t),然后对cosAt和sinAt进行调制，相加后即可得到QPSK言号二进制信息Q(t)载波友生命移相90度Asin(wt)Acos(wt)输出QPSK言号 图12QPSK调制器框图QPS

6、K军调QPSKJ号的解调原理如图3-5的方框图所示。解调是从已调信号中提取信号的过程，在某种意义上解调是调制的逆过程。由于QPSC号可以看作是两正交2PSK信号的叠加，故用两路正交的相干载波去解调，这样能够很容易地分离出这两路正交的2PSK信号。相干解调后的两路并行码元a和b经过“并/用”转换后成为串行数据输出。图3-5 QPSK信号解调原理方框图BPS缄制解调程序%勾造载波，产生8个码元，生成已调信号a=randsrc(1,8,0:1);% l=linspace(0,2*pi,50);% f=sin(2*l);%生成载波t=linspace(0,10*pi,400);% 点数out=1:40

7、0;%规定已调信号产生8个随机的二进制数利用linspace函数创建数组,2pi长度取点50个模拟一个码元定义时轴length为10pi,取点400个，代表8个码元的总取样b=1:400;%规定基带信号 w=1:400;%规定载波 length %fc成PSK信号 for i=1:8if a(i)=0for j=1:50out(j+50*(i-1)=f(j); % endelselength length若码元为0则将载波输出for j=1:50out(j+50*(i-1)=-f(j); % end若码元为1则将载波反相输出endend%俞出载波和基带信号 %for i=1:8for j=1:

8、50b(j+50*(i-1)=a(i); %b作为调制信号输出w(j+50*(i-1)=f(j); %w end作为载波输出endsubplot(3,3,1),plot(t,b),axis(0 10*pi 基带信号 ');grid on;subplot(3,3,2),plot(t,w),axis(0 10*pi 载波 '); grid on;subplot(3,3,3),plot(t,out),axis(0 '),title('PSK 波形 ');grid on;%已调信号加入高斯白噪声%noise=awgn(out,80,'measured&#

9、39;)比 80), xlabel('t'),ylabel('), xlabel('t'),ylabel('10*pi),xlabel('t'),ylabel('幅度 '),title('幅度 '),title('幅度产生噪音并加入到已调信号 out 中, 信噪subplot(334);plot(t,noise);ylabel('幅度');title('噪音+信号');xlabel('t');axis(0 10*pi ); grid on;%信

10、号通过BPF%Fs=400;%t=(1:400)*10*40/Fs;%b,a=ellip(4,40,10,25*2/Fs);%sf=filter(b,a,noise);%subplot(335); plot(t,sf);%xlabel('t'); ylabel(' 幅度 ');title(' axis(0 10*pi );grid on; %信号经过相乘器%抽样频率 400HZ时轴步进设计 IIR-BPF信号通过该滤波器画出信号通过该 BPF的波形通过BPF后的波形)；f=f f f f f f f f;%调整载波函数的长度，与 BPF输出函数统leng

11、ths=sf.*f;%信号与载波相乘s=(-1).*s;subplot(336);plot(t,s);%画出信号通过该相乘器的波形xlabel('t');ylabel('幅度');title('通过相乘器后波形');axis(010*pi-11);gridon;%信号通过LPF%Fs=400;%t=(1:400)*10*pi/Fs;%b,a=ellip(4,40,10*2/Fs);%sf=filter(b,a,s);%subplot(337);plot(t,sf);%xlabel('t'); ylabel(' 幅度 

12、9;);title('抽样频率 400HZ时轴步进设计 IIR-LPF信号通过该滤波器画出信号通过该低通滤波器的波形通过LPF后的波形)；axis(010*pi-11);gridon;%抽样判决%b=;%fori=1:8forj=1:50ifsf(j+50*(i-1)>bsf(j+50*(i-1)=1;%elsesf(j+50*(i-1)=0;%endendendsubplot(338);plot(t,sf);%xlabel('t');ylabel('幅度');title('抽样判决后波形');axis(310*pi);gridon

13、;QP颁制解调程序(1)调制%调相法clearallcloseallt=-1:;tt=length(t);x1=ones(1,800);fori=1:ttif(t(i)>=-1&t(i)<=1)|(t(i)>=5&t(i)<=7);x1(i)=1;elsex1(i)=-1;endendt1=0:;t2=0:;t3=-1:i=1:tt1if(t1(i)>=0&t1(i)<=2)|(t1(i)>=4&t1(i)<=8);x2(i)=1;elsex2(i)=-1;end设置判决门限若sf>判决门限，说明此时码元为1

14、若sf<判决门限，说明此时码元为0画出信号通过抽样判决器的波形endf=0:1;xrc=+*cos(pi*f);y1=conv(x1,xrc)/;y2=conv(x2,xrc)/;n0=randn(size(t2);f1=1;i=x1.*cos(2*pi*f1*t);q=x2.*sin(2*pi*f1*t1);I=i(101:800);Q=q(1:700);QPSK=sqrt(1/2).*I+sqrt(1/2).*Q;QPSK_n=(sqrt(1/2).*I+sqrt(1/2).*Q)+n0;n1=randn(size(t2);i_rc=y1.*cos(2*pi*f1*t3);q_rc=

15、y2.*sin(2*pi*f1*t4);I_rc=i_rc(101:800);Q_rc=q_rc(1:700);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);QPSK_rc_n1=QPSK_rc+n1;figure(1)subplot(4,1,1);plot(t3,i_rc);axis(-18-11);ylabel('a序列');subplot(4,1,2);plot(t4,q_rc);axis(-18-11);ylabel('b序列');subplot(4,1,3);plot(t2,QPSK_rc);axis(-18-11

16、);ylabel('合成序列');subplot(4,1,4);plot(t2,QPSK_rc_n1);axis(-18-11);ylabel('加入噪声');(2)解调%设定T=1,不加噪声clearallcloseall%调制bit_in=randint(1e3,1,01);bit_I=bit_in(1:2:1e3);bit_Q=bit_in(2:2:1e3);data_I=-2*bit_I+1;data_Q=-2*bit_Q+1;data_I1=repmat(data_I',20,1);data_Q1=repmat(data_Q',20,1

17、);fori=1:1e4data_I2(i)=data_I1(i);data_Q2(i)=data_Q1(i);end;t=0:;f=0:1;xrc=+*cos(pi*f);data_I2_rc=conv(data_I2,xrc)/;data_Q2_rc=conv(data_Q2,xrc)/;f1=1;t1=0:1e3+;I_rc=data_I2_rc.*cos(2*pi*f1*t1);Q_rc=data_Q2_rc.*sin(2*pi*f1*t1);QPSK_rc=(sqrt(1/2).*I_rc+sqrt(1/2).*Q_rc);%解调I_demo=QPSK_rc.*cos(2*pi*f1*t1);Q_demo=QPSK_rc.*sin(2*pi*f1*t1);I_recover=conv(I_demo,xrc);Q_recover=conv(Q_demo,xrc);I=I_recover(11:10010);Q=Q_recover(11:10010);t2=0:;t3=0:;data_recover=;fori=1:20:10000data_recover=data_recoverI(i:1:i+19)Q(i:1:i+19);end;ddd=-2*bit_in+1;ddd1=repmat(ddd',10,1);fori=1:1e4ddd2(