用Matlab写出QPSK误比特率性能仿真的程序

以下是一个简单的QPSK误比特率性能仿真程序：

clear all;
close all;

% 设置参数
M = 4; % 调制阶数
k = log2(M); % 每个符号的比特数
n = 1e6; % 模拟的符号数
SNRdB = 0:2:16; % 信噪比范围（dB）

% 生成QPSK调制表
qpskMod = comm.QPSKModulator();

% 生成AWGN信道
awgnChan = comm.AWGNChannel('NoiseMethod', 'Signal to noise ratio (SNR)', 'SNR', SNRdB(1), 'SignalPower', 1);

% 误码率变量
ber = zeros(1, length(SNRdB));

% 循环信噪比
for i = 1:length(SNRdB)
    % 显示当前信噪比
    disp(['SNR: ' num2str(SNRdB(i)) ' dB']);
    
    % 设置当前信噪比
    awgnChan.SNR = SNRdB(i);
    
    % 生成随机比特流
    bits = randi([0 1], n*k, 1);
    
    % 将比特流变成符号
    tx = qpskMod(bits);
    
    % 通过AWGN信道发送符号
    rx = awgnChan(tx);
    
    % 将接收的符号变成比特流
    rxBits = qpskDemod(rx);
    
    % 计算误码率
    ber(i) = sum(bits ~= rxBits)/length(bits);
end

% 绘制误码率性能曲线
semilogy(SNRdB, ber, 'b.-');
grid on;
xlabel('SNR (dB)');
ylabel('BER');
title('QPSK误比特率性能仿真');

注：其中 `qpskDemod` 函数是 `comm.QPSKDemodulator()` 的简单封装，代码如下：

function bits = qpskDemod(symbols)
% QPSK调制解调器
qpskDemod = comm.QPSKDemodulator();
bits = qpskDemod(symbols);
bits = reshape(de2bi(bits, 2, 'left-msb').', [], 1);
end