qpsk_all1.m

QPSK调制解调

%***********QPSK调制****************************
clear all
Rb=2400;   %信息速率
fc=19200;  %载波速率
fs=8*19200;%采样率
T=0.005;   %时间
s=randsrc(1,100*T*Rb,[0,1,2,3]);%信源序列
for k=1:length(s)
  st(((k-1)*fs/Rb+1):(k*fs/Rb))=sin(2*pi*(0.05*(1:fs/Rb)*fc/Rb+0.25*s(k)));
end
N=length(s);
N1=length(st);
%figure(1);
%plot(st);ylim([-2,2]);% 加未加噪声
%title('未加噪声QPSK信号');
%***********************************************
for k=1:61
SNR=61-k;
sout=AWGN(st,SNR);      %加9db的噪声
%figure(2);
%plot(sout);ylim([-2,2]);% 加噪声
%title('加噪声后的QPSK信号');
%***************加带通滤波器*********************
%*******************相干解调*********************
t=0:pi*Rb/fs:((N1-1)*Rb/fs)*pi;
y1=sout.*cos(2*t);
y2=sout.*sin(2*t);
scadd=zeros(1,N1);
ssadd=zeros(1,N1);
for i=1:N
    for j=2:fs/Rb
        scadd((i-1)*fs/Rb+j)=scadd((i-1)*fs/Rb+j-1)+y1((i-1)*fs/Rb+j-1);
        ssadd((i-1)*fs/Rb+j)=ssadd((i-1)*fs/Rb+j-1)+y2((i-1)*fs/Rb+j-1);
    end
end
%figure(3);
%subplot(3,1,1);plot(scadd);
%subplot(3,1,2);plot(ssadd);
%subplot(3,1,3);stem(s);
%*****************加低通滤波器******************
%*******************判决************************ 
for i=1:N
    if scadd(i*fs/Rb)>0
          if ssadd(i*fs/Rb)>0
          sget(i)=3;
          else sget(i)=0;
          end
    else
        if ssadd(i*fs/Rb)>0
        sget(i)=2;
        else 
        sget(i)=1;
        end
    end
end
sum=0;
for i=1:N
    if sget(i)==s(i);
        sum=sum;
    else sum=sum+1;
    end
end
Pe(k)=sum/N;
snr(k)=SNR;
disp(k);
end
figure(4);
plot(snr,Pe);grid on;ylim([0,1]);
title('QPSK的误码率与信噪比的关系图');