用MATLAB写一个QPSK调制解调的代码

以下是一个简单的QPSK调制解调的MATLAB代码示例：

% QPSK调制函数
function signal = qpsk_modulation(bits)
    % 将输入比特流按2个比特一组进行分组
    bit_pairs = reshape(bits, 2, []);
    % 将每组2个比特转换为一个QPSK符号，使用调制映射表
    symbol_indices = bi2de(bit_pairs', 'left-msb') + 1;
    modulation_table = [1+1i, -1+1i, -1-1i, 1-1i] / sqrt(2);
    symbols = modulation_table(symbol_indices);
    % 将符号序列串联为信号向量
    signal = symbols(:);
end

% QPSK解调函数
function bits = qpsk_demodulation(signal)
    % 将信号向量按4个符号一组进行分组
    symbols = reshape(signal, 4, []);
    % 将每组4个QPSK符号解调为2个比特，使用解调映射表
    demodulation_table = [0 0; 0 1; 1 1; 1 0];
    symbol_indices = zeros(1, size(symbols, 2));
    for i = 1:size(symbols, 2)
        distances = abs(symbols(:,i) - modulation_table);
        [~, index] = min(distances);
        symbol_indices(i) = index - 1;
    end
    bit_pairs = demodulation_table(symbol_indices+1, :);
    % 将比特对转换为比特流
    bits = bit_pairs(:)';
end

% 示例
bits = [0 1 1 0 1 0 0 1];
signal = qpsk_modulation(bits);
demodulated_bits = qpsk_demodulation(signal);
disp("原始比特流："); disp(bits);
disp("调制后信号："); disp(signal);
disp("解调后比特流："); disp(demodulated_bits);

该代码中，`qpsk_modulation`函数将输入的比特流进行QPSK调制，即将每2个二进制比特转换为一个QPSK符号。`qpsk_demodulation`函数实现了QPSK解调，即将接收到的QPSK符号解调为二进制比特。

示例中，我们将比特流`[0 1 1 0 1 0 0 1]`进行QPSK调制，得到调制信号`[0.7071+0.7071i -0.7071+0.7071i -0.7071-0.7071i 0.7071-0.7071i]`，然后进行解调操作，得到`[0 1 1 0 1 0 0 1]`，与原始比特流一致。