MATLAB实现QPSK调制解调及频谱分析

QPSK调制解调：

QPSK调制：

% QPSK调制
M = 4; % QPSK调制方式
N = 10000; % 产生的码元数
data = randi([0 M-1],1,N); % 随机产生0~M-1之间的整数
% 码元映射为QPSK符号
s = pskmod(data,M,pi/M);

QPSK解调：

% QPSK解调
% 产生QPSK调制信号
M = 4;
N = 10000;
data = randi([0 M-1],1,N);
s = pskmod(data,M,pi/M);
% 加入AWGN噪声
EbNo = 10; % 信噪比，单位dB
snr = EbNo + 10*log10(log2(M)); % SNR，单位dB
r = awgn(s,snr,'measured');
% QPSK解调
r_demod = pskdemod(r,M,pi/M);
% 计算误码率
error = sum(data~=r_demod);
ber = error/N;

QPSK频谱分析：

% QPSK频谱分析
M = 4;
N = 10000;
data = randi([0 M-1],1,N);
s = pskmod(data,M,pi/M);
Fs = 10000; % 采样率
t = 0:1/Fs:(N-1)/Fs; % 时域
f = (-N/2:N/2-1)/N*Fs; % 频域
% 计算QPSK调制信号的频谱
S = fftshift(fft(s))/N;
% 绘制频谱图
plot(f,abs(S));
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('QPSK Modulation Spectrum');