8qam调制解调matlab

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

% 8QAM 调制和解调 MATLAB 代码
clc; clear all; close all;

% 参数设置
M = 8; % 8QAM
k = log2(M); % 每个符号的比特数
n = 30000; % 符号数
numSamplesPerSymbol = 16; % 每个符号的采样数

% 发送端
dataIn = randi([0 1], n*k, 1); % 随机生成二进制数据
dataInMatrix = reshape(dataIn, k, length(dataIn)/k)';
% 将二进制数据按照 k 比特进行分组，并转换为十进制数字
dataSymbolsIn = bi2de(dataInMatrix);

% 8QAM 调制
dataMod = qammod(dataSymbolsIn, M);

% 信道
SNR_dB = 20; % 信噪比
dataModNoisy = awgn(dataMod, SNR_dB, 'measured');

% 接收端
dataSymbolsOut = qamdemod(dataModNoisy, M);

% 将十进制数字转换为二进制数据
dataOutMatrix = de2bi(dataSymbolsOut, k);
dataOut = dataOutMatrix(:);

% BER 计算
[numErrors, ber] = biterr(dataIn, dataOut);

% 结果显示
disp(['符号数 = ' num2str(n)]);
disp(['符号每个比特数 = ' num2str(k)]);
disp(['总比特数 = ' num2str(n*k)]);
disp(['采样数 = ' num2str(numSamplesPerSymbol)]);
disp(['信噪比 = ' num2str(SNR_dB) ' dB']);
disp(['误比特率 = ' num2str(ber)]);
disp(['错误比特数 = ' num2str(numErrors)]);

% 信号显示
figure;
subplot(2,1,1);
plot(dataModNoisy(1:100*numSamplesPerSymbol),'b.');
title('8QAM 调制信号');
xlabel('实部');
ylabel('虚部');
subplot(2,1,2);
plot(dataSymbolsOut(1:100),'ro');
title('8QAM 解调信号');
xlabel('实部');
ylabel('虚部');