帮我写一个用matlab实现4qam调制的通信系统

好的，下面是一个简单的用 Matlab 实现 4QAM 调制的通信系统的示例代码：

% 信号生成
M = 4; % 调制阶数
k = log2(M); % 每个符号的比特数
num_bits = 10000; % 数据位数
data = randi([0 1], num_bits, 1); % 生成随机比特流
data_in = reshape(data, k, length(data)/k)'; % 将比特流重塑为 k 位比特块
symbols = bi2de(data_in, 'left-msb')'; % 将二进制数据块转换为符号

% 符号映射
constellation = qammod(0:M-1, M); % 生成调制符号映射
mapped_symbols = constellation(symbols+1); % 符号映射

% 添加高斯噪声
EbNo = 10; % 信噪比（以 dB 为单位）
snr = EbNo + 10*log10(k) - 10*log10(M); % 将信噪比转换为信道信噪比
noisy_symbols = awgn(mapped_symbols, snr, 'measured'); % 添加高斯噪声

% 解调和译码
noisy_data = qamdemod(noisy_symbols, constellation); % 解调调制符号
noisy_data_in = de2bi(noisy_data, k, 'left-msb')'; % 将解调后的数据转换为二进制
received_data = noisy_data_in(:); % 将解调后的数据重塑为一维比特流

% 比特误码率计算
[num_err, ber] = biterr(data, received_data); % 计算比特误码率

% 结果输出
fprintf('比特误码率：%d/%d = %g\n', num_err, num_bits, ber);

以上代码实现了一个简单的 4QAM 调制通信系统，包括数据生成、符号映射、高斯噪声添加、解调译码和比特误码率计算等步骤。请根据自己的需求进行修改和调整。