不同多普勒频移下的QPSK误码率代码

以下是一个MATLAB实现的QPSK调制解调器的误码率仿真代码，考虑了不同多普勒频移的影响：

clear all; close all; clc;

% QPSK参数定义
mod_order = 4; % 调制阶数

% 信道参数定义
fd = [0 20 100]; % 多普勒频移
tau = [0 3 6] * 1e-6; % 时延
pdb = [0 -3 -6]; % 衰减

% 生成随机数据
tx_data = randi([0 mod_order-1], 1, 1e6);

% QPSK调制
tx_mod = qammod(tx_data, mod_order);

% 信道模型
for i = 1:length(fd)
    h(:,i) = rayleigh_chan(1, fd(i), tau(i), pdb(i), length(tx_mod));
end

% 加噪声并接收
rx_sig = zeros(size(tx_mod));
for i = 1:length(fd)
    rx_sig = rx_sig + awgn(conv(h(:,i), tx_mod), 10, 'measured');
end

% 信道估计
h_est = zeros(size(h));
for i = 1:length(fd)
    h_est(:,i) = ls_est(rx_sig, tx_mod, fd(i), tau(i), pdb(i));
end

% QPSK解调
rx_demod = qamdemod(rx_sig, mod_order);

% 解调得到的符号
rx_data = reshape(rx_demod, 1, length(tx_data));

% 误码率计算
num_err = sum(abs(tx_data-rx_data) > 0);
ber = num_err / length(tx_data);
fprintf('误码率为 %f\n', ber);

其中，`rayleigh_chan` 函数用于生成瑞利衰落信道，`ls_est` 函数用于进行信道估计，这里只是一个简单的实现，可以根据需要进行修改。