目录
0.OFDM子载波互拍效应
1.half-cycled技术原理
2.half-cycled算法实现的步骤
2.1 发送端的频域交叉OFDM信号生成
2.2 发送端的half-cycled截断操作
2.3 发送端的符号合并——频谱效率恢复
2.4 光纤传输与接收端的半周期恢复
2.5 接收端的FFT解调与互拍噪声隔离
3.Matlab程序和仿真
在直接检测正交频分复用(DDO-OFDM)系统中，子载波互拍噪声(由平方检波过程中不同子载波的非线性混频产生)是限制系统性能的核心瓶颈之一。传统抑制方案(如预留保护间隔)往往以牺牲频谱效率为代价，而half-cycled技术通过利用OFDM符号的时域对称性，在不降低频谱效率的前提下实现了子载波互拍噪声的有效隔离。
0.OFDM子载波互拍效应
OFDM子载波互拍效应，是直接检测正交频分复用(DDO-OFDM)系统中，由光电探测器的平方律检波特性引发的一种非线性噪声干扰现象，也是制约DDO-OFDM系统传输性能的核心瓶颈之一。在DDO-OFDM系统中，接收端的光电探测器(PD)遵循平方律检波特性，即探测器输出的光电流与入射光功率的平方成正比X，将OFDM时域信号代入平方运算后可得：
这个表达式中包含两类成分：
互拍噪声的频率为fk−fm，由于子载波数量N通常很大，不同子载波对的频率差会覆盖整个OFDM 信号的频带范围。这意味着互拍噪声会弥散在所有子载波的传输带宽内，无法通过简单的滤波手段消除。
1.half-cycled技术原理
定义DDO-OFDM系统中，OFDM符号的快速傅里叶逆变换(IFFT)尺寸为N ，单个符号的时间长度为T ，则经过IFFT后的时域信号为：
将该OFDM符号分为前半段（1≤t1≤2T）和后半段（2T+1≤t2≤T），分别表示为：
在频域交叉OFDM方案中，奇数子载波的调制系数ck=0 ，仅偶数子载波携带数据（记偶数子载波索引为n）。因此上式可简化为：
由此可得频域交叉OFDM符号的时域奇对称关系：
这一关系表明：频域交叉OFDM符号的后半段是前半段的 “反相复制”，即符号关于时间中点呈奇对称。这是half-cycled技术的核心理论基础。
2.half-cycled算法实现的步骤
half-cycled技术的实现分为发送端处理和接收端处理两个部分，通过“截断-合并-恢复”的流程，实现子载波互拍噪声的隔离与频谱效率的保持。
2.1 发送端的频域交叉OFDM信号生成
首先在发送端生成频域交叉OFDM信号：仅在偶数子载波上调制数据，奇数子载波全部置零。
假设IFFT尺寸为N，则频域数据可表示为：
其中ck为QPSK、16QAM等调制格式对应的复数符号。
对该频域数据执行IFFT，得到时域OFDM符号：
根据前文推导，该时域符号满足s(t+T/2)=−s(t)的奇对称关系。
2.2 发送端的half-cycled截断操作
由于频域交叉OFDM符号的后半段是前半段的反相复制，因此发送端仅需传输符号的前半段，直接丢弃后半段。设原OFDM符号的时间长度为T ，截断后符号的时间长度变为T/2，对应的时域信号为：
2.3 发送端的符号合并——频谱效率恢复
为了恢复频谱效率，发送端将两个连续的half-cycled符号合并为一个常规长度的OFDM符号：
第一个half-cycled符号的前半段作为合并后符号的前半段；
第二个half-cycled符号的前半段作为合并后符号的后半段。
合并后的符号时间长度为T，与常规OFDM符号长度一致，因此频谱效率与常规OFDM相同(解决了频域交叉 OFDM 频谱效率低的问题)。
2.4 光纤传输与接收端的半周期恢复
合并后的half-cycled符号经电光调制、光纤传输后，被接收端的光电探测器(PD)接收并转换为电信号。由于DDO-OFDM系统在短距离传输中可近似为时不变系统，接收端通过半周期复制操作恢复完整的OFDM符号：
设接收端接收到的half-cycled符号前半段为rhalf (t)(对应发送端的shalf(t) )，则通过奇对称关系恢复后半段：
恢复后的完整时域符号为：
2.5 接收端的FFT解调与互拍噪声隔离
对恢复后的完整时域符号执行快速傅里叶变换(FFT)，转换回频域：
由于发送端仅在偶数子载波上调制数据，接收端经FFT后，子载波互拍噪声会集中在奇数子载波上(偶数子载波为有效数据)。此时只需提取偶数子载波上的信号，即可避开互拍噪声的干扰，实现误码性能的改善。
3.Matlab程序和仿真
以下是half-cycled算法的核心流程：
根据这个流程图，我们编写如下的matlab程序：
clear; clc; close all; rng('default'); % 固定随机种子，保证结果可复现 rng(1); %% ====================== 1. 系统参数配置 ====================== N_fft = 256; % IFFT/FFT尺寸（与前文一致） N_subcarrier = 100; % 有效子载波数（偶数子载波） mod_order = 4; % 调制阶数：4=QPSK，16=16QAM T_symbol = 1e-6; % 单个OFDM符号时长（1μs） t = linspace(0, T_symbol, N_fft); % 时域采样点 CP_len = 4; % 循环前缀长度（half-cycled缩短为4点） SNR = 25; % 接收端信噪比(dB) fiber_length = 10; % 光纤传输距离(km)，仅仿真损耗，短距离忽略色散 %% ====================== 2. 发送端处理 ====================== % 2.1 生成频域数据（仅偶数子载波调制，奇数子载波置零） c_k = zeros(1, N_fft); % 频域系数初始化 even_idx = 2:2:N_subcarrier; % 偶数子载波索引（有效子载波） % 生成调制符号（QPSK/16QAM） if mod_order == 4 mod_sym = pskmod(randi([0,3],1,N_subcarrier/2),4); % QPSK elseif mod_order == 16 mod_sym = qammod(randi([0,15],1,N_subcarrier/2),16); % 16QAM end c_k(even_idx) = mod_sym; % 偶数子载波赋值，奇数子载波保持0 % 2.2 IFFT生成时域OFDM符号（频域交叉OFDM） s_time = ifft(c_k, N_fft); % IFFT变换 % 验证时域奇对称特性（前半段与后半段反相） s_half1 = s_time(1:N_fft/2); % 前半段 s_half2 = s_time(N_fft/2+1:end); % 后半段 fprintf('时域奇对称验证：后半段与前半段反相差值=%.4f\n', norm(s_half2 + s_half1)); % 2.3 half-cycled截断：仅保留前半段，丢弃后半段 s_half = s_half1; % 截断后的half-cycled符号 % 2.4 符号合并：2个half-cycled符号拼接（恢复频谱效率） % 生成第二个half-cycled符号（模拟连续符号传输） c_k2 = zeros(1, N_fft); mod_sym2 = pskmod(randi([0,3],1,N_subcarrier/2),4); % 第二个符号调制数据 c_k2(even_idx) = mod_sym2; s_time2 = ifft(c_k2, N_fft); s_half2 = s_time2(1:N_fft/2); % 合并为完整T长度符号（前半=符号1，后半=符号2） s_tx = [s_half, s_half2]; % 2.5 添加循环前缀（CP） s_tx_cp = [s_tx(end-CP_len+1:end), s_tx]; %% ====================== 3. 信道传输（简化DDO-OFDM信道） ====================== % 3.1 电光转换+光纤传输（仅考虑损耗和加性高斯噪声） fiber_loss = 0.2*fiber_length; % 光纤损耗：0.2dB/km s_tx_power = 10^(-fiber_loss/20) * s_tx_cp; % 功率衰减 % 加性高斯白噪声（模拟接收机噪声） s_rx_cp = awgn(s_tx_power, SNR, 'measured'); % 3.2 平方律检波（DDO-OFDM核心：光电探测器平方特性，产生互拍噪声） s_rx_cp_sq = abs(s_rx_cp).^2; %% ====================== 4. 接收端处理 ====================== % 4.1 去除循环前缀 s_rx = s_rx_cp_sq(CP_len+1:end); % 4.2 拆分合并后的符号：恢复两个half-cycled符号 s_rx_half1 = s_rx(1:N_fft/2); % 第一个符号的前半段 s_rx_half2 = s_rx(N_fft/2+1:end); % 第二个符号的前半段 % 4.3 半周期恢复：利用奇对称恢复完整符号 % 恢复第一个符号的后半段（反相复制） s_rx_full1 = [s_rx_half1, -s_rx_half1]; % 恢复第二个符号的后半段 s_rx_full2 = [s_rx_half2, -s_rx_half2]; % 4.4 FFT解调（频域恢复） r_k1 = fft(s_rx_full1, N_fft); % 第一个符号FFT r_k2 = fft(s_rx_full2, N_fft); % 第二个符号FFT % 4.5 提取有效子载波（偶数子载波，隔离奇数子载波互拍噪声） r_sym1 = r_k1(even_idx); % 第一个符号解调结果 r_sym2 = r_k2(even_idx); % 第二个符号解调结果 % 4.6 解调判决（恢复原始比特） if mod_order == 4 demod_sym1 = pskdemod(r_sym1,4); demod_sym2 = pskdemod(r_sym2,4); elseif mod_order == 16 demod_sym1 = qamdemod(r_sym1,16); demod_sym2 = qamdemod(r_sym2,16); end %% ====================== 5. 性能分析与结果展示 ====================== % 5.1 原始调制符号与解调符号对比（验证正确性） orig_sym1 = pskdemod(mod_sym,4); % 原始符号 orig_sym2 = pskdemod(mod_sym2,4); ber1 = sum(orig_sym1 ~= demod_sym1)/length(orig_sym1); % 误码率 ber2 = sum(orig_sym2 ~= demod_sym2)/length(orig_sym2); % 5.2 结果输出 fprintf('===== 系统仿真结果 =====\n'); fprintf('调制格式：%dQAM\n', mod_order); fprintf('第一个符号误码率：%.4f\n', ber1); fprintf('第二个符号误码率：%.4f\n', ber2); fprintf('接收端信噪比：%ddB\n', SNR); % 5.3 时域/频域波形可视化 figure('Color','w'); subplot(2,2,1); plot(t, real(s_time)); title('原始频域交叉OFDM时域信号'); xlabel('时间(s)'); ylabel('幅度'); grid on; subplot(2,2,2); plot(t(1:N_fft/2), real(s_half)); title('half-cycled截断后信号（前半段）'); xlabel('时间(s)'); ylabel('幅度'); grid on; subplot(2,2,3); plot(abs(c_k)); title('发送端频域数据（偶数子载波）'); xlabel('子载波索引'); ylabel('幅度'); grid on; subplot(2,2,4); plot(abs(r_k1)); title('接收端FFT后频域数据（互拍噪声在奇数子载波）'); xlabel('子载波索引'); ylabel('幅度'); grid on;仿真结构如下图所示：
