QPSK调制与解调在MATLAB平台上的实现

1、下载可编辑.专业 .整理 .qpsk 调制与解调在matlabf 台上的实现李悦qpsi 即四进制移向键控 (quaternary phase shift keying), 它利用载波的四种不同相位来表示数字信息，由于每一种载波相位代表两个比特信息，因此每个四进制码元可以用两个二进制码元的组合来表示。两个二进制码元中的前一个码元用 a 表示，后一个码元用b 表示。qpsk 言号可以看作两个载波正交2psk 信号的合成，下图表示qpske交调制器。由 qpsk! 号的调制可知，对它的解调可以采用与2psk 信号类似的解调方法进行解调。解调原理图如下所示，同相支路和正交支路分别采用相干解调方式解调

2、，得到l(t)和q(t)，经过抽样判决和并 /串交换器，将上下支路得到的并行数据恢复成串行数据输入磯波振沏输出输入抽样刿决输出滤波器变换低通. 抽样判决载波低通滤液器位定时下载可编辑.专业 .整理 .%调相法clear all close all t=-1:0.01:7-0.01; tt=le ngth(t); x1=o nes(1,800); for i=1:tt if (t(i)=-1 & t(i)=5 & t(i)=0 & t1(i)= 4& t1(i)0 data_recover_a(i:i+19)=1; bit_recover=bit_recover

3、1; else data_recover_a(i:i+19)=-1; bit_recover=bit_recover -1; end end error=0; dd = -2*bit_i n+1; ddd=dd; ddd 仁 repmat(ddd,20,1); for i=1:2e4 ddd2(i)=ddd1(i); end for i=1:1e3 if bit_recover(i)=ddd(i) error=error+1; end end p=error/1000; figure(1) subplot(2,1,1);plot(t2,ddd2);axis(0 100 -2 2);title(

4、 subplot(2,1,2);plot(t2,data_recover_a);axis(0 100 -2 2);title(原序列 ); 解调后序列 );下载可编辑.专业 .整理 .效果图 : 下载可编辑.专业 .整理 .%设定 t=1,不加噪声clear all close all %调制bit_in = randin t(1e3, 1, 0 1); bit=bit_i n(1:2:1e3); bit_q = bit_i n(2:2:1e3); data_i = -2*bit_l+1; data_q = -2*bit_q+1; data_i1=repmat(data_i,20,1); da

5、ta_q 1=repmat(data_q,20,1); for i=1:1e4 data_i2(i)=data_i1(i); data_q2(i)=data_q1(i); en d; t=0:0.1:1e3-0.1; f=0:0.1:1; xrc=0.5+0.5*cos(pi*f); data_i2_rc=c on v(data_i2,xrc)/5.5; data_q2_rc=co nv(data_q2,xrc)/5.5; f1=1; t1=0:0.1:1e3+0.9; i_rc=data_i2_rc.*cos(2*pi*f1*t1); q_rc=data_q2_rc.*si n(2*pi*f

6、1*t1); qpsk_rc=(sqrt(1/2).*l_rc+sqrt(1/2).*q_rc); %解调i_demo=qpsk_rc.*cos(2*pi*f1*t1); q_demo=qpsk_rc.*si n(2*pi*f1*t1); i_recover=c on v(i_demo,xrc); q_recover=c on v(q_demo,xrc); l=l_recover(11:10010); q=q_recover(11:10010); t2=0:0.05:1e3-0.05;dafalrecove卫一八for耳20=0000 daqrlrecove卫daqrlrecover-(z+l

7、9)q(z+19)kdddh2*brn+l八ddchhrepmasdd二03-forlll二e4figures subpof(411)八p-ofsdaxisao2006一)八subpof(4l2)八pof(?q)oi.s(o20下载可编辑.专业 .整理 .% qpsk 误码率分析snri ndb 仁 0:2:10; snri ndb2=0:0.1:10; for i=1:le ngth(snri ndb1) pb,ps=cm_sm32(snr in db1(i); smld_bit_err_prb(i)=pb; smld_symbol_err_prb(i)=ps; en d; for i=1:

8、le ngth(snr in db2) snr=exp(snri ndb2(i)*log(10)/10); theo_err_prb(i)=qfu nct(sqrt(2*snr); en d; title(qpsk 误码率分析 ); semilogy(snri ndb1,smld_bit_err_prb,*); axis(0 10 10e-8 1); hold on; % semilogy(snri ndb1,smld_symbol_err_prb,o); semilogy(snr in db2,theo_err_prb); legend( 仿真比特误码率 ,理论比特误码率 ); hold o

9、ff; fun ctio ny=qfu nct(x) y=(1/2)*erfc(x/sqrt(2); fun ctio n pb,ps=cm_sm32(snri ndb) n=10000; e=1; snr=10a(snri ndb/10); sgma=sqrt(e/snr)/2; s00=1 0; s0 仁0 1; si 仁-1 0; s10=0 -1; for i=1:n temp=ra nd; if (temp0.25) dsource1(i)=0; dsource2(i)=0; elseif (temp0.5) dsource1(i)=0; 下载可编辑.专业 .整理 .dsource2

10、(i)=1; elseif (temp0.75) dsource1(i)=1; dsource2(i)=0; else dsource1(i)=1; dsource2(i)=1; en d; en d; numo fsymbolerror=0; num ofbiterror=0; for i=1:n n=sgma*ra ndn( size(soo); if(dsource1(i)=0)&( dsource2(i)=0) r=soo+n; elseif(dsource1(i)=0)&(dsource2(i)=1) r=s01+ n; elseif(dsource1(i)=1)&a

11、mp;(dsource2(i)=0) r=s10+n; else r=s11+ n; en d; coo=dot(r,soo); c01= dot(r,s01); c10=dot(r,s10); c11=dot(r,s11); c_max=max(coo c01 c10 c11); if (c00=c_max) decis1=0;decis2=0; elseif(c01=c_max) decis1=0;decis2=1; elseif(c10=c_max) decis1=1;decis2=0; else decis1=1;decis2=1; en d; symbolerror=0; if(decis1=dsource1(i) num ofbiterror= num ofbiterror+1; symbolerror=1; en d; if(decis2=dsource2(i) num ofbiterro