查询网站空间商,建立公司网站的重点,网站主机域名,免费申请淘宝账号注册目录
0.OFDM子载波互拍效应
1.half-cycled技术原理
2.half-cycled算法实现的步骤
2.1 发送端的频域交叉OFDM信号生成
2.2 发送端的half-cycled截断操作
2.3 发送端的符号合并——频谱效率恢复
2.4 光纤传输与接收端的半周期恢复
2.5 接收端的FFT解调与互拍噪声隔离
3.…目录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和后半段2T1≤t2≤T分别表示为在频域交叉OFDM方案中奇数子载波的调制系数ck0 仅偶数子载波携带数据记偶数子载波索引为n。因此上式可简化为由此可得频域交叉OFDM符号的时域奇对称关系这一关系表明频域交叉OFDM符号的后半段是前半段的 “反相复制”即符号关于时间中点呈奇对称。这是half-cycled技术的核心理论基础。2.half-cycled算法实现的步骤half-cycled技术的实现分为发送端处理和接收端处理两个部分通过“截断-合并-恢复”的流程实现子载波互拍噪声的隔离与频谱效率的保持。2.1 发送端的频域交叉OFDM信号生成首先在发送端生成频域交叉OFDM信号仅在偶数子载波上调制数据奇数子载波全部置零。假设IFFT尺寸为N则频域数据可表示为其中ck为QPSK、16QAM等调制格式对应的复数符号。对该频域数据执行IFFT得到时域OFDM符号根据前文推导该时域符号满足s(tT/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; % 调制阶数4QPSK1616QAM 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/21: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_len1: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_len1:end); % 4.2 拆分合并后的符号恢复两个half-cycled符号 s_rx_half1 s_rx(1:N_fft/2); % 第一个符号的前半段 s_rx_half2 s_rx(N_fft/21: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;仿真结构如下图所示