QPSK调制后进行载波同步实现

QPSK（Quadrature Phase Shift Keying，正交相移键控）是一种常见的数字调制方式，它通过在正交的两个载波上分别调制两个比特的信息来实现高效的信号传输。在接收端，为了正确解调QPSK信号，需要进行载波同步，以确保接收机能够准确地恢复出发送信号的载波相位和频率。

QPSK调制和载波同步的Matlab代码

%% QPSK调制和载波同步
% 清空环境变量
clear all
close all
clc

%% 参数设置
fs = 1e6; % 采样频率
fc = 100e3; % 载波频率
Tb = 1e-6; % 比特周期
N = 1000; % 总比特数
Eb = 1; % 每比特能量
noise_power = 0.1; % 噪声功率

%% 生成随机比特序列
data = randi([0, 1], N, 2); % 每个符号2比特
I_data = 2 * data(:, 1) - 1; % 映射到±1
Q_data = 2 * data(:, 2) - 1; % 映射到±1

%% QPSK调制
t = (0:N-1) * Tb; % 时间序列
I_signal = I_data' .* cos(2 * pi * fc * t); % I路信号
Q_signal = Q_data' .* sin(2 * pi * fc * t); % Q路信号
QPSK_signal = I_signal + 1j * Q_signal; % QPSK信号

%% 添加噪声
noise = sqrt(noise_power / 2) * (randn(1, N) + 1j * randn(1, N)); % 复高斯噪声
received_signal = QPSK_signal + noise; % 接收到的信号

%% 载波同步 - 锁相环（PLL）
% PLL参数
Kp = 0.01; % 比例增益
Ki = 0.001; % 积分增益
theta = 0; % 初始相位估计
omega = fc; % 初始频率估计
phase_error = 0; % 相位误差
phase_integral = 0; % 相位积分

% PLL处理
synchronized_signal = zeros(1, N);
for k = 1:N
    % 相位检测
    phase_error = angle(received_signal(k) * exp(-1j * theta));
    
    % 环路滤波
    phase_integral = phase_integral + Ki * phase_error;
    theta = theta + omega * Tb + Kp * phase_error + phase_integral;
    
    % 限制相位在[-pi, pi]范围内
    theta = mod(theta + pi, 2 * pi) - pi;
    
    % 解调信号
    synchronized_signal(k) = received_signal(k) * exp(-1j * theta);
end

%% 解调
I_demodulated = real(synchronized_signal);
Q_demodulated = imag(synchronized_signal);
I_demodulated = I_demodulated > 0;
Q_demodulated = Q_demodulated > 0;
demodulated_data = [I_demodulated; Q_demodulated]';

%% 计算误码率
errors = sum(abs(demodulated_data - data));
BER = errors / (2 * N);
fprintf('误码率（BER）: %.4f\n', BER);

%% 绘图
figure;
subplot(3, 1, 1);
plot(real(QPSK_signal), imag(QPSK_signal), 'o');
title('QPSK调制信号星座图');
xlabel('I');
ylabel('Q');

subplot(3, 1, 2);
plot(real(received_signal), imag(received_signal), 'o');
title('接收到的信号星座图');
xlabel('I');
ylabel('Q');

subplot(3, 1, 3);
plot(real(synchronized_signal), imag(synchronized_signal), 'o');
title('载波同步后的信号星座图');
xlabel('I');
ylabel('Q');

代码说明

1. 参数设置：设置采样频率、载波频率、比特周期、总比特数、每比特能量和噪声功率。
2. 生成随机比特序列：随机生成比特序列，并将其映射到QPSK符号上。
3. QPSK调制：将比特序列调制成QPSK信号。
4. 添加噪声：在调制后的信号中添加复高斯噪声，模拟信道噪声。
5. 载波同步 - 锁相环（PLL）：使用锁相环进行载波同步，恢复出信号的载波相位和频率。
6. 解调：对同步后的信号进行解调，恢复出原始比特序列。
7. 计算误码率：计算解调后的误码率。
8. 绘图：绘制调制信号、接收到的信号和载波同步后的信号的星座图，以便直观观察同步效果。

注意事项

• 在实际应用中，PLL的参数（如比例增益Kp和积分增益Ki）需要根据具体情况进行调整，以确保良好的同步性能。
• 代码中的噪声功率可以根据信噪比（SNR）的要求进行调整。
• 该示例假设信号的采样率足够高，能够满足奈奎斯特采样定理的要求。在实际应用中，需要根据具体信号和信道条件进行采样率的设计。
• QPSK调制后进行载波同步 代码