matlab实现QPSK

QPSK（Quadrature Phase Shift Keying）是一种数字调制技术，它可以将两个比特信息映射到一个符号上，从而实现信息的传输和接收。下面是MATLAB实现QPSK的示例代码：

% 生成随机比特序列
N = 1000; % 比特序列长度
bits = randi([0, 1], 1, N);

% 比特序列分组成2比特一组
symbols = zeros(1, N/2); % 初始化符号数组
for i = 1:N/2
    if bits(2*i-1) == 0 && bits(2*i) == 0
        symbols(i) = -1-1i; % 映射为QPSK符号
    elseif bits(2*i-1) == 0 && bits(2*i) == 1
        symbols(i) = -1+1i;
    elseif bits(2*i-1) == 1 && bits(2*i) == 0
        symbols(i) = 1-1i;
    elseif bits(2*i-1) == 1 && bits(2*i) == 1
        symbols(i) = 1+1i;
    end
end

% 显示QPSK星座图
scatterplot(symbols);
title('QPSK Constellation');

% 添加噪声
EbN0dB = 10; % 信噪比（dB）
EbN0 = 10^(EbN0dB/10); % 信噪比（线性值）
noiseVar = 1/(2*EbN0); % 噪声方差
noise = sqrt(noiseVar) * (randn(1, N/2) + 1i*randn(1, N/2)); % 生成噪声
rxSymbols = symbols + noise; % 接收到的符号

% 解调
rxBits = zeros(1, N); % 初始化接收到的比特序列
for i = 1:N/2
    if real(rxSymbols(i)) < 0
        rxBits(2*i-1) = 0;
    else
        rxBits(2*i-1) = 1;
    end
    if imag(rxSymbols(i)) < 0
        rxBits(2*i) = 0;
    else
        rxBits(2*i) = 1;
    end
end

% 计算误码率
numErrors = sum(rxBits ~= bits);
BER = numErrors/N;
fprintf('误码率为 %f\n', BER);

该代码首先生成随机比特序列，然后将比特序列分组成2比特一组，并映射为QPSK符号。接着，添加高斯噪声，并进行解调。最后，计算误码率并输出。运行该代码可以得到QPSK星座图，以及在10dB信噪比下的误码率。