请写一个符合该文献的图6的matlab代码

为了生成类似于图6所示的64-GBaud DP-16-QAM系统设置的MATLAB仿真代码，可以参考以下示例代码框架：

% 参数初始化
symbol_rate = 64e9; % 符号率（每秒64千兆符号）
num_subcarriers = 8; % 子载波数量
prbs_length = 2^15 - 1; % PRBS长度
roll_off_factor = 0.1; % 滚降系数

% 生成伪随机比特序列
data_bits = prbs15_generator(prbs_length * log2(num_subcarriers));

% 映射为16-QAM符号
qam_symbols = qam_modulate(data_bits);

% 实现I-DFT-S
qam_symbols_parallel = reshape(qam_symbols, [], num_subcarriers);
dft_symbols = idft(qam_symbols_parallel);

% 将IDFT后的符号串行化并分配给子载波
tx_symbols = dft_symbols(:);

% 上采样和根升余弦脉冲成形
tx_symbols_upsampled = upsample(tx_symbols, oversampling_factor);
tx_waveform = rcosfilter(tx_symbols_upsampled, 'sqrt', roll_off_factor, symbol_period);

% 频率偏移补偿
tx_waveform_shifted = freq_offset_compensation(tx_waveform, freq_offset);

% 复用子载波信号
tx_signal = mux_subcarriers(tx_waveform_shifted);

% 调制到光载波上
optical_signal = iq_mzm_modulation(tx_signal, bias_voltage, pi_voltage);

% 发送端处理结束，后续部分涉及传输信道模拟等步骤...

function bits = prbs15_generator(length)
    % PRBS15生成器实现...
end

function qam_symbols = qam_modulate(bits)
    % 16-QAM调制函数实现...
end

function dft_symbols = idft(symbols_parallel)
    % 计算N点IDFT...
end

function waveform = rcosfilter(signal, type, alpha, Ts)
    % 根升余弦滤波器实现...
end

function compensated_signal = freq_offset_compensation(signal, freq_offset)
    % 频率偏移补偿算法实现...
end

function muxed_signal = mux_subcarriers(signals)
    % 子载波复用函数实现...
end

function optical_signal = iq_mzm_modulation(electrical_signal, bias, Vpi)
    % 双极性相位调制器IQ-MZM模型...
end

请注意，上述代码仅为示意性的，并未包含所有细节及实际参数值。需要具体实现时应参照论文中的详细说明进行调整和完善。此外，`prbs15_generator`, `qam_modulate`, `idft`, `rcosfilter`, `freq_offset_compensation`, `mux_subcarriers`, 和 `iq_mzm_modulation` 函数都需要具体的实现逻辑来完成各自的功能。