QPSK调制与解调（matlab实现）

一、QPSK调制与解调的原理

二、仿真

画出QPSK信号产生的信息、I路、Q路信号对比图；显示接收端解调的I路、Q 路，信息恢复的对比图。

三、matlab代码

clear all;clc;
%参数设定
N=20;%比特数
T=1;%比特周期
fc=2;%载波频率，载波频率是比特频率的两倍
Fs=100;%抽样频率
bitstream=randi([0,1],1,N);%产生随机的比特流。返回区间[0,1]的离散分布整数，组成1×N的数组
bitstream=2*bitstream-1;%单极性码变双极性码(0，1变为-1,1)
I=[];Q=[];

%奇数进I路，偶数进Q路。串并变换
for i=1:N
    if mod(i,2)~=0
        I=[I,bitstream(i)];
    else
        Q=[Q,bitstream(i)];
    end
end
%绘图比较I、Q比特流
bit_data=[];
for i=1:N
    bit_data=[bit_data,bitstream(i)*ones(1,T*Fs)];%一个比特周期里面有T*Fs个采样点
end
I_data=[];Q_data=[];
%I路和Q路是原来比特周期的两倍，因此采样点数为2*比特流采样点数
for i=1:N/2
    I_data=[I_data,I(i)*ones(1,T*Fs*2)];
    Q_data=[Q_data,Q(i)*ones(1,T*Fs*2)];
end
%绘图
figure();
%时间轴
t=0:1/Fs:T*N-1/Fs;
subplot(3,1,1)
plot(t,bit_data);legend('Bitstream')%比特信息
subplot(3,1,2)
plot(t,I_data);legend('I Bitstream')
subplot(3,1,3)
plot(t,Q_data);legend('Q Bitstream')

%载波信号
bit_t=0:1/Fs:2*T-1/Fs;%定义时间轴，I、Q两路码元长度为2T，为了好显示
%定义I路和Q路载波
I_carrier=[];Q_carrier=[];
for i=1:N/2
    I_carrier=[I_carrier,I(i)*cos(2*pi*fc*bit_t)];
    Q_carrier=[Q_carrier,Q(i)*cos(2*pi*fc*bit_t+pi/2)];
end
%传输信号
QPSK_signal=I_carrier+Q_carrier;
%绘图
figure();%产生一个新图
subplot(3,1,1)
plot(t,I_carrier);legend('I signal')%I路信号
subplot(3,1,2)
plot(t,Q_carrier);legend('Q signal')%Q路信号
subplot(3,1,3)
plot(t,QPSK_signal);legend('QPSK signal')%I路、Q路的和信号
snr=1;%信噪比
%接收信号
QPSK_receive=awgn(QPSK_signal,snr);%awgn()添加噪声

%解调
for i=1:N/2
    I_output=QPSK_receive((i-1)*length(bit_t)+1:i*length(bit_t)).*cos(2*pi*fc*bit_t);
    %当i=1时，QPSK_receive(1,200),每bit_t这个时间长度解调一次。matlab中数组下标从1开始，区别与C语言
    if sum(I_output)>0
        I_recover(i)=1;
    else
        I_recover(i)=-1;
    end
    Q_output=QPSK_receive(1,(i-1)*length(bit_t)+1:i*length(bit_t)).*cos(2*pi*fc*bit_t+pi/2)
    if sum(Q_output)>0
        Q_recover(i)=1;
    else 
        Q_recover(i)=-1;
    end
end
%并串变换
bit_recover=[];
for i=1:N
    if mod(i,2)~=0
        bit_recover=[bit_recover,I_recover((i-1)/2+1)];%奇数取I路信息
    else
        bit_recover=[bit_recover,Q_recover(i/2)];%偶数取Q路信息
    end
end
%绘图比较I、Q比特流
recover_data=[];
for i=1:N
    recover_data=[recover_data,bit_recover(i)*ones(1,T*Fs)];
end
I_recover_data=[];Q_recover_data=[];
for i=1:N/2
    I_recover_data=[I_recover_data,I_recover(i)*ones(1,2*T*Fs)];
    Q_recover_data=[Q_recover_data,Q_recover(i)*ones(1,2*T*Fs)];
end
%绘图
figure();
t=0:1/Fs:N*T-1/Fs;
subplot(3,1,1)
plot(t,recover_data);legend('Bitstream');
subplot(3,1,2)
plot(t,I_recover_data);legend('I Bitstream');
subplot(3,1,3)
plot(t,Q_recover_data);legend('Q Bitstream');

四、思考与总结

1.有关bit_t的定义。解调时，分I、Q两路，此时两路信号码元周期为2T，因此定义时间轴bit_t=2T。根据定义length(bit_t)=2*T*Fs=2×1×100=200，I、Q两路信号的每个0-1序列对应长度200，那么每200就解调一次。

2.在matlab中，数组索引从1开始，区别于C语言。