生成正交调制解调QPSK信号MATLAB程序

以下是一个简单的 MATLAB 程序，用于生成正交调制解调 QPSK 信号：

% 定义常量和参数
fc = 1000; % 载波频率
fs = 10000; % 采样率
T = 1/fs; % 采样间隔
Ts = 1/10; % 符号间隔
Ns = fs * Ts; % 每个符号的采样数
L = 1000; % 信号长度
SNR = 10; % 信噪比

% 定义调制信号序列
data = randi([0 1], 1, L/2); % 生成随机的二进制数据
s = qammod(data, 4); % 使用 QAM 调制将二进制数据转换为 QPSK 符号

% 将 QPSK 符号解调为二进制数据
received_signal = awgn(s, SNR); % 添加高斯白噪声
data_hat = qamdemod(received_signal, 4); % 使用 QAM 解调将接收信号转换回二进制数据

% 将二进制数据转换为 QPSK 符号
s_hat = qammod(data_hat, 4); % 使用 QAM 调制将二进制数据转换为 QPSK 符号

% 生成正交载波
t = 0:T:(L*Ns-1)*T;
carrier1 = cos(2*pi*fc*t);
carrier2 = sin(2*pi*fc*t);

% 将 QPSK 符号调制到正交载波上
x1 = real(s) .* carrier1;
x2 = imag(s) .* carrier2;
x = x1 - x2;

% 将调制信号解调回 QPSK 符号
y1 = x .* carrier1;
y2 = x .* carrier2;
y = conv(y1 - y2, ones(1, Ns)); % 低通滤波
y = y(1:Ns:end); % 降采样
y = y(1:length(s)); % 去除多余部分
s_hat2 = y./abs(y); % 能量归一化
data_hat2 = qamdemod(s_hat2, 4); % 将符号转换回二进制数据

% 输出结果
disp("原始数据：");
disp(data);
disp("解调后的数据：");
disp(data_hat);
disp("解调后的数据（第二种方法）：");
disp(data_hat2);

该程序生成一个长度为 L 的随机二进制数据序列，将其转换为 QPSK 符号并调制到正交载波上。然后，程序添加高斯白噪声并将接收信号解调回二进制数据。最后，程序将接收到的 QPSK 符号解调回二进制数据，以便与原始数据进行比较。