16QAM调制解调 matlab代码

以下是一个简单的16QAM调制解调的matlab代码示例：

% 16QAM调制解调
clc;
clear all;
close all;

%% 参数设置
M = 16; % QAM调制阶数
k = log2(M); % 每个符号的比特数

%% 生成随机数据
data = randi([0 1], 1, 1000*k); % 产生随机比特流

%% 16QAM调制
dataInMatrix = reshape(data, k, length(data)/k)';
dataSymbolsIn = bi2de(dataInMatrix); % 二进制转十进制
dataMod = qammod(dataSymbolsIn, M); % 调制

%% 16QAM解调
dataSymbolsOut = qamdemod(dataMod, M); % 解调
dataOutMatrix = de2bi(dataSymbolsOut, k); % 十进制转二进制
dataOut = dataOutMatrix(:)'; % 展开

%% 计算误码率
[numErrors,ber] = biterr(data,dataOut); % 计算误码率
disp(['误码率 = ',num2str(ber)]);

%% 绘制调制前后的星座图
scatterplot(dataSymbolsIn); % 调制前的星座图
title('16QAM调制前星座图');
scatterplot(dataMod); % 调制后的星座图
title('16QAM调制后星座图');

%% 绘制调制前后的眼图
eyediagram(data,5*k); % 调制前的眼图
title('16QAM调制前眼图');
eyediagram(dataMod,5*k); % 调制后的眼图
title('16QAM调制后眼图');