QAM调制：使用随机0/1序列生成数据并进行调制

本文将介绍如何将原代码修改为先产生随机0/1序列data，再进行调用MOD_QAM函数进行QAM调制，而不改变代码的功能。

修改后的代码

data_generated = randi([0 1], 1, Nfft-Cpilot);  // 生成随机0/1序列数据
data_modulated = MOD_QAM(data_generated, M);  // 调用MOD_QAM函数进行QAM调制

MOD_QAM函数代码

function data_modulated = MOD_QAM(data_generated, M)
% data_generated: 待调制数据
% M: 调制阶数（2/4/16/64/256 ...）
% data_modulated：已调制数据

% 生成星座图
[gray_map, gray_idx] = generate_gray_map(M);

% 计算需要调制的符号数量
num_symbols = ceil(length(data_generated)/log2(M));

% 初始化调制输出数组
data_modulated = zeros(1, num_symbols);

% 遍历每个符号，进行调制
for m = 1:num_symbols
    % 计算当前符号对应的二进制数据
    binary_data = dec2bin(data_generated(m), log2(M));
    
    % 查找对应的格雷码
    gray_code = '';
    for j = 1:log2(M)
        if j == 1
            gray_code(j) = binary_data(j);
        else
            gray_code(j) = xor(binary_data(j), binary_data(j-1)) + '0';
        end
    end
    
    % 查找对应的星座点索引
    idx = find(gray_idx==bin2dec(gray_code))-1;
    
    % 将当前符号对应的星座点赋值给调制输出数组
    if m <= length(data_generated)
        % 如果还有剩余符号没有调制，则正常赋值
        data_modulated(m) = gray_map(idx+1);
    else
        % 如果所有符号已经调制完毕，则只扩展数组大小，不赋值
        data_modulated = [data_modulated, gray_map(idx+1)];
    end
end
end

解释：

代码首先使用 randi([0 1], 1, Nfft-Cpilot) 生成随机的0/1序列，作为待调制数据。
然后调用MOD_QAM函数，将随机生成的0/1序列数据进行QAM调制，生成已调制数据。
MOD_QAM函数实现根据调制阶数M对输入的二进制数据进行格雷码映射，然后根据格雷码索引找到对应的星座点，并将星座点赋值给调制输出数组。

优点：

这种方法更灵活，可以方便地调整数据生成方式，例如使用其他随机数生成方法或自定义数据序列。
采用函数调用方式，代码结构更加清晰，便于维护和扩展。

注意：

需要确保Nfft、Cpilot、M等参数定义正确，以保证代码的功能和逻辑正确性。
MOD_QAM函数需根据实际需求进行实现，例如可以使用不同的格雷码映射方法或星座图定义。

通过以上修改，代码可以实现先产生随机0/1序列，再进行QAM调制的流程，并保持原有功能不变。这可以提高代码的灵活性，并使代码结构更加清晰。