四分之piQPSK解調部分濾波

1、% pi/4-QPSK系統(tǒng)仿真clear;% -初始化數(shù)據(jù)-T=1; % 碼元周期fc=1/T; % 載波頻率ml=2; % 調制階數(shù) nb=100; % 傳輸比特數(shù)samp_T=T/100; % 抽樣間隔fs=1/samp_T; % 抽樣頻率SNR=8; % 信噪比t=0:samp_T:nb*T-samp_T;N=length(t); % 碼元個數(shù) % -產生隨機信號-data=randint(1,nb,0 1); datanrz=data.*2-1; % 轉換成雙極性非歸零碼data1=zeros(1,nb/samp_T); for q=1:nb data1(q-1)/samp_T/2+1

2、:q/samp_T/2)=datanrz(q);end % 傳輸信號 % -pi/4-QPSK調制-idata=datanrz(1:ml: (nb-1); % 串/并轉換qdata=datanrz(2:ml:nb);idata1=zeros(1,nb/2/samp_T);qdata1=zeros(1,nb/2/samp_T);for q=1:nb/2 idata1(q-1)/samp_T+1:q/samp_T)=idata(q); qdata1(q-1)/samp_T+1:q/samp_T)=qdata(q);endidata2(1)=-1;qdata2(1)=0;for i=1:nb/2 %

3、 映射邏輯 if (idata(i)=0) & (qdata(i)=0) angl=-3*pi/4; elseif (idata(i)=0) & (qdata(i)=1) angl=3*pi/4; elseif (idata(i)=1) & (qdata(i)=0) angl=-pi/4; else angl=pi/4; endidata2(i+1)=idata2(i)*cos(angl)-qdata2(i)*sin(angl);qdata2(i+1)=idata2(i)*sin(angl)+qdata2(i)*cos(angl);endich=zeros(1,nb/2

4、/samp_T);for i=1:nb/2 ich(i-1)/samp_T+1:i/samp_T)=idata2(i);endfor ii=1:N/2 a(ii)=cos(2*pi*fc*t(ii);endidata3=ich.*a; % I支路調制qch=zeros(1,nb/2/samp_T);for j=1:nb/2 qch(j-1)/samp_T+1:j/samp_T)=qdata2(j);endfor jj=1:N/2 b(jj)=sin(2*pi*fc*t(jj);endqdata3=qch.*b; % Q支路調制s=idata3+qdata3; % pi/4-QPSK已調信號%-

5、pi/4-QPSK解調-ich1=s.*a; % 相乘載波qch1=s.*b;b=fir1(300,0.5); % 低通濾波ich2=filter(b,1,ich1);qch2=filter(b,1,qch1);for i=1:nb/2/samp_T-1 % 基帶差分變換ich3(i+1)=ich2(i+1)*ich2(i)+qch2(i+1)*qch2(i);qch3(i+1)=qch2(i+1)*ich2(i)-ich2(i+1)*qch2(i);ich3(i)=ich3(i+1);qch3(i)=qch3(i+1);endich4=zeros(1,nb/2);qch4=zeros(1,n

6、b/2);for n=1:nb/2 % 判決輸出 ichsum(n)=sum(ich3(n-1)/samp_T+1:n/samp_T)*samp_T; if ichsum(n)>=0 ich4(n)=1; else ich4(n)=0; end qchsum(n)=sum(qch3(n-1)/samp_T+1:n/samp_T)*samp_T; if qchsum(n)>=0 qch4(n)=1; else qch4(n)=0; endendich4=ich4.*2-1;qch4=qch4.*2-1;idata4=zeros(1,nb/samp_T/2); qdata4=zeros

7、(1,nb/samp_T/2);for q=1:nb/2 idata4(q-1)/samp_T+1:q/samp_T)=ich4(q); qdata4(q-1)/samp_T+1:q/samp_T)=qch4(q);end demodata=zeros(1,nb); % 并/串轉換demodata(1:ml: (nb-1)=ich4;demodata(2:ml:nb)=qch4;demodata1=zeros(1,nb/samp_T); for q=1:nb demodata1(q-1)/samp_T/2+1:q/samp_T/2)=demodata(q);end畫圖部分figure(1)su

8、bplot(311);plot(data1);axis(0 1e3 -2 2);title('基帶信號');grid onsubplot(312);plot(idata1);axis(0 1e3 -2 2);title('I支路信號');grid onsubplot(313);plot(qdata1);axis(0 1e3 -2 2);title('Q支路信號');grid onfigure(2)subplot(211);plot(ich);axis(0 3e3 -2 2);title('I支路映射后信號');grid onsub

9、plot(212);plot(qch);axis(0 3e3 -2 2);title('Q支路映射后信號');grid onfigure(3)subplot(311);plot(idata3);axis(0 3e3 -2 2);title('I支路載波調制后信號');grid onsubplot(312);plot(qdata3);axis(0 3e3 -2 2);title('Q支路載波調制后信號');grid onsubplot(313);plot(s);axis(0 3e3 -2 2);title('pi/4-QPSK已調信號');grid onfigure(4)subplot(211);plot(ich1);axis(0 3e3 -2 2);title('相乘載波后I支路信號');grid onsubplot(212);plot(qch1);axis(0 3e3 -2 2);title('相乘載波后Q支路信號');grid onfigure(5)s