目录
1.MZM调制器的原理
1.1 输入输出电场关系
1.2 不同检测方式下的传输函数
1.3 MZM的非线性与偏置点选择
2.MATLAB仿真
光通信系统的调制方式分为直接调制和外调制，两者的核心差异与适用场景如下：
直接调制：通过改变半导体激光器的注入电流调节光强，优点是结构简单、成本低，但存在带宽有限、消光比低、线宽较大的缺陷，仅适用于短距离(≤100km)直接检测光正交频分复用(DDO-OFDM)系统。
外调制：利用独立于激光源的外部调制器实现电光转换，具备带宽大、消光比高、激光器线宽可调的优势，适用于长距离(>100km)DDO-OFDM或大容量中长距离相干光正交频分复用(CO-OFDM)系统。目前主流外调制器包括马赫-曾德尔调制器(MZM)(基于电光效应)和电吸收调制器(EAM)(基于电吸收效应)。
1.MZM调制器的原理
MZM通过电光效应实现电信号对光信号的调制，其结构包含上下两个电极：通过调节加载在电极上的外加电压，改变光电材料的折射率，进而控制输出光信号的强度。
MZM参数定义：
1.1 输入输出电场关系
MZM输入光信号Ein(t)与输出光信号Eout(t)的电场满足：
将电信号的三角函数形式(ejθ=cosθ+jsinθ)代入，可分解为幅度调制分量和相位调制分量：
1.2 不同检测方式下的传输函数
根据接收端检测方式的不同，MZM的传输特性分为两类：
直接检测(DDO-OFDM系统)：传输函数为输出功率与输入功率的比值(功率与电场幅度的平方成正比)，令V(t)=V1(t)−V2(t)，则：
相干检测(CO-OFDM 系统)：传输函数为输出电场与输入电场的比值：
1.3 MZM的非线性与偏置点选择
为了使MZM高效高性能的工作，应该使调制信号尽量落在MZM线性度高的区域内，这可以通过调节调制信号的峰值和MZM的配置电压来控制。当信号的幅度位于MZM传输函数非线性高的区域时，信号将严重失真，从而降低系统的误码性能。为避免信号失真，需将调制信号限制在线性度高的区域，并选择合适的偏置点：
DDO-OFDM系统：最佳偏置点为正交点；
CO-OFDM系统：最佳偏置点为零点。
DDO-OFDM系统与CO-OFDM系统相比，由于在接收端不需要提供光混频器、平衡接收机等高成本的器件，成本将会得到很好的控制。但由于子载波互拍噪声和光纤传输中的色散会导致频率选择性衰落，DDO-OFDM 系统的光纤传输距离非常有限，因此主要应用于短距离的有线接入网和点对点的数据中心之间的传输。
2.MATLAB仿真
%% ==================== 1. 仿真参数配置 ==================== clear; clc; close all; % 系统基础参数 Fs = 10e9; % 采样频率 10GHz Ts = 1/Fs; % 采样周期 t_total = 1e-6; % 总仿真时间 1us t = 0:Ts:t_total-Ts;% 时间轴 % MZM核心参数 Vpi = 2; % 半波电压 2V delta = 20; % 消光比 20dB gamma = (sqrt(10^(delta/10)) - 1) / (sqrt(10^(delta/10)) + 1); % 对称因子 V_bias_ddo = Vpi/2; % DDO-OFDM最佳偏置（正交点 Vpi/2） V_bias_co = Vpi; % CO-OFDM最佳偏置（零点 Vpi） % OFDM基带信号参数 N_carrier = 64; % OFDM子载波数 N_sym = 16; % OFDM符号数 mod_order = 16; % 16QAM调制 % ==================== 2. 生成OFDM基带调制信号 ==================== % 生成16QAM星座图数据 data = randi([0 mod_order-1], N_carrier, N_sym); qam_mod = qammod(data, mod_order, 'UnitAveragePower', true); % IFFT变换生成OFDM时域信号 ofdm_sym = ifft(qam_mod, N_carrier, 1); ofdm_sym = [ofdm_sym(N_carrier/2+1:end, :); ofdm_sym(1:N_carrier/2, :)]; % 频谱搬移 ofdm_signal = reshape(ofdm_sym, [], 1); % 串行化 % 信号重采样匹配系统采样率 tx_signal = resample(ofdm_signal, Fs, N_carrier*Fs/(t_total*Fs)); tx_signal = tx_signal(1:length(t)); % 对齐时间轴 tx_signal = tx_signal / max(abs(tx_signal)); % 归一化幅度 % ==================== 3. MZM调制过程建模 ==================== % 输入光信号（假设为单频连续光） E_in = 1; % 输入光电场幅度归一化 % 调制电压 = 基带信号 + 直流偏置 V_ddo = tx_signal + V_bias_ddo; % DDO-OFDM调制电压 V_co = tx_signal + V_bias_co; % CO-OFDM调制电压 % 3.1 DDO-OFDM直接检测（功率检测） P_out_ddo = (abs(E_in))^2 .* cos(pi * V_ddo / (2*Vpi)).^2; P_in = (abs(E_in))^2; % 输入光功率 % 3.2 CO-OFDM相干检测（电场检测） E_out_co = E_in .* cos(pi * V_co / (2*Vpi)) .* exp(1j * pi * V_bias_co / Vpi); % 相干检测需本地振荡光（假设LO功率与输入光相同） LO_phase = 0; % 本地振荡光相位 E_lo = E_in * exp(1j * LO_phase); I_co = real(E_out_co .* conj(E_lo)); % 相干检测同相分量 % ==================== 4. 传输特性曲线绘制 ==================== % 生成扫描电压（用于绘制MZM传输曲线） V_scan = linspace(-2*Vpi, 2*Vpi, 1000); P_trans_ddo = cos(pi*(V_scan + V_bias_ddo)/(2*Vpi)).^2; % DDO传输曲线 E_trans_co = cos(pi*(V_scan + V_bias_co)/(2*Vpi)); % CO传输曲线 figure('Color','w'); subplot(2,1,1); plot(V_scan, P_trans_ddo, 'b-', 'LineWidth',1.5); hold on; grid on; xlabel('调制电压 V(t) [V]'); ylabel('归一化输出功率 P_{out}/P_{in}'); title('DDO-OFDM 直接检测 MZM传输特性'); xlim([-2*Vpi 2*Vpi]); ylim([0 1.1]); subplot(2,1,2); plot(V_scan, E_trans_co, 'r-', 'LineWidth',1.5); hold on; grid on; xlabel('调制电压 V(t) [V]'); ylabel('归一化输出电场 |E_{out}/E_{in}|'); title('CO-OFDM 相干检测 MZM传输特性'); xlim([-2*Vpi 2*Vpi]); ylim([-1.1 1.1]); % ==================== 5. 调制信号时域波形绘制 ==================== figure('Color','w'); subplot(2,1,1); plot(t*1e9, P_out_ddo, 'b-', 'LineWidth',1); xlabel('时间 [ns]'); ylabel('归一化功率'); title('DDO-OFDM 直接检测 输出功率波形'); xlim([0 0.2]); % 显示前0.2us便于观察 subplot(2,1,2); plot(t*1e9, I_co, 'r-', 'LineWidth',1); xlabel('时间 [ns]'); ylabel('相干检测同相分量'); title('CO-OFDM 相干检测 输出波形'); xlim([0 0.2]); % ==================== 6. 性能分析：计算信号失真度 ==================== % 计算基带信号与调制后信号的相关系数（衡量失真） corr_ddo = corrcoef(abs(tx_signal), sqrt(P_out_ddo)); corr_co = corrcoef(real(tx_signal), I_co); fprintf('DDO-OFDM调制信号相关系数：%.4f\n', corr_ddo(1,2)); fprintf('CO-OFDM调制信号相关系数：%.4f\n', corr_co(1,2));仿真结果如下图所示：
