尧图网站建设 
简介面向光纤传感与OFDR技术研究用的MATLAB、LabVIEW仿真源码包聚焦光学频率域反射OFDR分布式测量方法适合科研人员、工程师及高校学生学习OFDR信号处理与传感解算。压缩包共5个文件均为m脚本整体仅2KB代码精炼涵盖OFDR基本仿真、温度响应分析、空间分辨率计算与测量距离解算等核心模块便于按需调用与二次开发。已有1877人学习下载具备一定参考热度。通过对照源码可以理解啁啾脉冲生成、傅里叶变换解调、频域到空间域映射等关键算法掌握从回波频率变化反演温度、应变及位置信息的完整流程同时可结合温度、分辨率、测距不同功能脚本快速搭建自己的OFDR数据处理原型为电力电缆热监测、桥梁结构健康检测等实际分布式光纤传感系统设计提供理论支持。1. 从拍频到距离OFDR 信号处理为什么离不开 MATLAB 和 LabVIEW光频域反射OFDR和传统 OTDR 最大的区别在于它不靠时间飞行而靠拍频频率来定位光纤上的反射点和应变点。扫频激光器发出的光经过辅助干涉仪做时钟触发主干涉仪的信号被等光频间隔采样后做一次 FFT 就能把频域映射到空间域。这个思路听起来简单但真正在 MATLAB 里把干涉信号变成距离曲线时会遇到频谱泄漏、色散失配、扫频非线性残留等一系列问题。而 LabVIEW 的角色恰恰相反它更适合在采集端做实时显示和硬件控制把 FPGA 采集卡的数据边采边画同时把原始干涉数据落盘留给 MATLAB 做离线精算。这篇文章围绕 OFDR 这个主题把 MATLAB 和 LabVIEW 两条处理链路各自该干什么、参数怎么设、坑在哪讲清楚面向的是已经接手 OFDR 系统、手里有采集数据但还没跑通完整处理流程的工程师。2. OFDR 信号模型与拍频解调的几个关键参数2.1 拍频频率与空间分辨率的映射关系OFDR 的原始信号是扫频光源经过马赫-曾德尔干涉仪后产生的拍频信号。设扫频范围为 ΔF扫频周期为 T则扫频速率 γ ΔF / T。光纤中位置 z 处的反射点与参考臂之间的延时为 τ 2nz/c对应的拍频频率为 f_b γτ。对时域信号做 FFT 后频率轴乘以系数 c/(2nγ) 就得到距离轴。空间分辨率 δz c/(2nΔF)这意味着扫频范围直接决定分辨率而不是采样率。比如 ΔF 10 nm 在 1550 nm 波段对应约 1.25 THz理论上分辨率能到 80 μm 左右但实际还要受窗函数和色散补偿精度的影响。% 参数定义扫频范围、中心波长、光纤折射率 c 3e8; % 光速 n 1.4682; % 单模光纤群折射率 lambda_c 1550e-9; % 中心波长 delta_lambda 10e-9; % 扫频范围 10 nm delta_f c / lambda_c^2 * delta_lambda; % 换算成频率扫频范围 dz_theory c / (2 * n * delta_f); % 理论空间分辨率 fprintf(理论分辨率: %.2f um\n, dz_theory * 1e6);这段代码先把波长扫频范围换算成频率扫频范围再算理论分辨率。注意群折射率 n 不是相折射率OFDR 定位用的是光在光纤中的群速度取值通常在 1.468 附近不同厂家光纤有细微差别这个值不准会直接导致距离轴整体偏移。分辨率算出来只是理论值实际 FFT 加窗后主瓣展宽相邻反射点要分开通常需要 23 倍理论分辨率的间隔。2.2 辅助干涉仪时钟触发与重采样扫频激光器不可能是理想线性的频率随时间的变化总有非线性。如果不做校正FFT 之后的峰值会展宽甚至劈裂直接毁掉分辨率。最常见的做法是用一个辅助干涉仪生成等光频间隔的时钟信号用它作为采集卡的采样时钟或触发信号让主干涉仪信号在等频率间隔上被采样。这种硬件重采样方案在 LabVIEW 里实现时通常是把辅助干涉仪的过零点检测信号接到采集卡的 PFI 端口作为外部时钟。// LabVIEW 伪代码外部时钟采样配置 // 使用 NI 采集卡的 External Clock 模式 // 时钟源: PFI0 (辅助干涉仪过零脉冲) // 采样模式: Continuous, 采样率由外部时钟决定 // 主干涉仪信号接入 AI0, 采集数据写入循环缓冲区外部时钟采样的关键点在于辅助干涉仪的延时越长每个扫频周期内的时钟脉冲数越多重采样点数越密但过长的延时会增加光源相位噪声的影响。常见的辅助干涉仪延时设在 50200 米光纤对应长度折中之后每扫频周期能采到几万到几十万个点。如果采集卡不支持外部时钟也可以在内时钟模式下采集后在 MATLAB 里用辅助干涉仪信号做插值重采样效果稍差但可行。% 基于辅助干涉仪的软件重采样 % t: 等时间间隔时间轴, ai_sig: 辅助干涉仪信号 % mz_sig: 主干涉仪信号 [~, locs] findpeaks(ai_sig, MinPeakDistance, 10); % 相邻过零点对应等频率间隔 phase unwrap(angle(hilbert(ai_sig))); % 提取瞬时相位 freq_idx round(phase / (2*pi) * max_phase); % 映射到等频率索引 resampled interp1(t, mz_sig, t(locs), spline);这里用希尔伯特变换提取辅助干涉仪的瞬时相位相位每增加 2π 对应光源频率走过一个固定间隔用这些过零点位置对主干涉仪信号做插值。MinPeakDistance要按采样率设置太小会检测到噪声毛刺太大丢失有效过零点。插值方法建议用spline而不是linear因为相位映射后相邻点间距不均匀线性插值会引入额外的高频噪声。2.3 色散失配补偿的必要性参考臂和测试臂的光纤长度不一致时色散会导致不同频率成分经历不同的群延时表现为 FFT 后峰值展宽和位置偏移。特别是在测试距离超过几百米时这个效应不可忽略。补偿方法是在频域乘上一个二次相位因子参数可以通过对已知反射点的峰值宽度做优化来标定。% 色散补偿频域二次相位因子 N length(resampled); f_axis (0:N-1) / N * fs; D 17e-6; % 色散系数, 单位 s/m^2, 典型值需要标定 phase_comp exp(1i * pi * D * (f_axis - fs/2).^2); compensated ifft(fft(resampled) .* phase_comp);色散系数 D 不是光纤 datasheet 上的色散参数而是根据系统光路算出来的等效值最靠谱的标定方式是放一个反射峰扫 D 的值让峰最窄。实践中可以把 D 的搜索范围设定在理论值的 ±50%步长逐步减半两三轮就能收敛。3. OFDR MATLAB 处理链路从原始干涉数据到应变曲线3.1 数据预处理去直流、加窗与零填充采集卡拿到的原始干涉信号带有直流偏置来自光电探测器的平均光功率不先去掉会占据 FFT 动态范围。直流分量可以直接减去信号均值但如果扫频过程中光源功率有波动均值法不够准需要做高速滤波或用多项式拟合基线。% 去直流与加窗 mz_ac mz_sig - mean(mz_sig); % 去直流 win hanning(N, periodic); % 周期汉宁窗 mz_win mz_ac .* win; % 加窗 % 零填充到 2 的幂次方便 FFT Nfft 2^nextpow2(N * 4); spectrum fft(mz_win, Nfft);汉宁窗是 OFDR 里的默认选择主瓣比矩形窗宽但旁瓣抑制好能压住距离轴上反射峰旁边的假目标。periodic选项比symmetric更适合频谱分析因为它在离散傅里叶变换下的泄漏特性更好。零填充不会提高真实分辨率但让频谱看起来更平滑峰值定位精度能提升到亚像素级别。3.2 频域转距离域与峰值提取做完 FFT 之后频谱的横轴频率通过 f γτ 映射到距离。这里有一个容易被忽略的问题FFT 结果包含正负频率OFDR 信号经过希尔伯特变换或 IQ 解调后通常是单边带的但如果是直接用光电探测器采集的实数信号频谱会左右对称需要取单边谱并把幅值乘以 2。% 频域转距离域 f_res fs / Nfft; % 频率分辨率 f_axis (0:Nfft/2-1) * f_res; % 单边频率轴 dist_axis f_axis * c / (2 * n * gamma_sweep); % 距离轴 amp_spectrum abs(spectrum(1:Nfft/2)) * 2 / N; % 单边幅值谱 % 峰值检测 [pks, locs] findpeaks(amp_spectrum, MinPeakHeight, threshold, ... MinPeakDistance, round(dz_theory / (dist_axis(2)-dist_axis(1))));MinPeakDistance取值要大于理论分辨率对应的像素数否则同一个反射峰可能被检测成多个峰。阈值threshold一般设为噪声底以上 10 dB噪声底可以用距离轴上无反射区段的幅值统计出来。峰值检测出来的位置配合相位信息就能算出亚像素精度的距离值——OFDR 的定位精度可以做到远小于空间分辨率靠的就是相位。3.3 滑动窗 FFT 与分布式应变解调分布式应变测量的核心思路是把采集到的干涉数据分成若干段每段做一次 FFT得到不同位置处的反射光谱。当光纤某一段发生应变时对应位置的反射光谱会发生频移通过互相关计算这个频移量就能还原出应变分布。窗长决定了空间分辨率和应变分辨率之间的取舍。% 滑动窗 FFT 提取局部反射光谱 win_len 512; % 窗长, 决定空间分辨率 step 64; % 滑动步长, 决定采样间距 num_windows floor((N - win_len) / step); for k 1:num_windows idx (k-1)*step 1 : (k-1)*step win_len; seg mz_ac(idx) .* hanning(win_len, periodic); spec fft(seg); % 记录局部光谱, 用于后续互相关计算频移 spectra(:, k) spec(1:win_len/2); end窗长选 512 还是 2048取决于你要的定位粒度。假设每段对应的距离长度是 ΔL应变分辨率大约正比于 ΔL 的平方根所以追求高空间分辨率时应变灵敏度会下降。互相关计算频移时频移量与应变的换算系数大约为 0.78 GHz/με1550 nm 波段这个系数受光纤光弹系数影响不同类型的传感光纤有区别最好先在拉伸台上标定。3.4 频谱拼接与长距离测量优化测试距离越长单次扫频的数据量越大FFT 需要的点数也越多。距离超过 1 km 时一次扫频往往包含几百万个采样点直接 FFT 虽然能做但速度慢。常见做法是把数据分段做处理取平均或者在频域做子带滤波后降采样只保留感兴趣的距离范围。% 频域子带提取只关心特定距离范围 dist_min 100; dist_max 200; % 关心的距离范围, 米 freq_min 2 * n * dist_min * gamma_sweep / c; freq_max 2 * n * dist_max * gamma_sweep / c; idx_range find(f_axis freq_min f_axis freq_max); % 对子带做 IFFT 得到时间域信号, 降低后续处理数据量 sub_band spectrum(idx_range); sub_time ifft(sub_band);这种做法的好处是后续的应变计算只需要处理这一个子带的信号计算量能降一个数量级。注意 IFFT 之后距离轴的起点不再从零开始需要把偏移量加回去。4. 用 LabVIEW 搭 OFDR 实时采集与显示系统4.1 采集卡选型与采样参数配置LabVIEW 在 OFDR 系统里的定位是实时采集平台。采集卡的选型要看三个指标采样率至少 100 MS/s 以上取决于扫频速率和最大拍频频率、分辨率 14 bit 或更高、支持外部时钟输入。拍频频率的上限由最远测量距离决定如果最远测 100 米2nLγ/c 算出来的频率要在采集卡带宽范围内同时要留 20% 余量给滤波器滚降。// 关键参数配置示意 // DAQmx Create Channel AI Voltage // 采样率: 200 MS/s (内部时钟) // 触发源: PFI0 (辅助干涉仪时钟) // 采样模式: Continuous // 每通道采样数: 0 (连续模式不设单次采样数)外部时钟模式下AI 通道的采样率参数变成参考时钟的最大允许速率实际采样率由外部时钟频率决定。数据读取用循环缓冲方式每次读取最新的一帧数据做 FFT 显示同时把原始数据写入 TDMS 文件。连续采集时最怕的是读写速度不匹配缓冲区溢出会导致数据丢帧。读取速率要略快于产生速率并用队列机制做生产者消费者模式。4.2 生产者-消费者架构与实时频谱显示实时显示 OFDR 距离曲线核心是一个标准的生产者消费者模式采集循环往队列里写入原始数据块处理循环从队列取出数据做 FFT 和坐标映射再把结果送到波形图控件。队列深度要按采样率和每块数据的处理时间来定处理时间超过获取时间时队列就会持续增长最终耗尽内存。// 生产者循环: 读取 DAQ 数据 // DAQmx Read (NChan NSamp) - 队列插入 // 队列大小: 100 个数据块, 超出则丢弃最旧块 // 消费者循环: // 队列取出 - FFT - 距离轴映射 - XY 图显示 // FFT 使用 VI: Spectrum Measurements.vi // 窗口类型: Hanning, 输出单位: 幅值谱XY 图显示距离曲线时要注意FFT 输出的频率轴是均匀分布的但映射到距离轴后仍然是均匀的直接作为 X 轴数组传入就行。如果显示 10000 个点以上的曲线用Waveform Chart刷新率会下降建议用XY Graph一次性刷新整帧数据或者用Image控件做强度图横轴距离、纵轴扫频次数能直接看到整个测量过程中的变化。4.3 TDMS 文件落盘与 MATLAB 离线数据对接LabVIEW 采集的原始数据最终要交给 MATLAB 做精细处理。TDMS 格式是 NI 的推荐格式MATLAB 原生不支持直接读取但可以通过tdmsread函数MATLAB R2019a 之后内置读取或者用第三方工具包。更稳妥的方案是在 LabVIEW 里直接把数据写成二进制文件同时附一个包含采样率、扫频参数、折射率等元数据的文本文件。// TDMS 写入配置 // TDMS Open - TDMS Write - TDMS Close // 分组名: OFDR_Data, 通道名: MainInterferometer, AuxiliaryInterferometer // 属性: SweepRate, CenterWavelength, RefractiveIndex% MATLAB 读取 LabVIEW 写入的 TDMS 文件 data tdmsread(scan_001.tdms); mz data.MainInterferometer.Data; ai data.AuxiliaryInterferometer.Data; info data.MainInterferometer.Property; % 读取元数据落盘时最容易被忽略的是数据类型和字节序。LabVIEW 默认小端序MATLAB 的fread也要对应指定ieee-le。如果用 TDMS 格式加上属性通道后文件会包含完整的记录信息但文件体积会比纯二进制大 10%20%。对长时间连续采集的场景建议 LabVIEW 每扫频周期写入一个文件避免单个文件过大导致读取时内存不足。4.4 扫频光源触发的同步策略OFDR 系统需要保证数据采集和扫频光源严格同步。通常的做法是让光源输出一个扫频开始触发电平接到采集卡的 PFI 触发端口采集卡从第一个触发沿开始记录一帧完整的数据。连续测量时每次触发代表一个扫频周期开始帧间通过触发沿对齐数据处理时每帧独立做 FFT。5. 光纤瑞利散射信号的衰落噪声抑制与图像化显示5.1 偏振衰落与频谱平均的取舍OFDR 信号来自瑞利散射散射光的偏振态随机变化如果探测光的偏振方向与参考光不一致拍频效率会降低表现为局部信号衰落甚至完全消失。最常见的处理办法是偏振分集接收把信号分成两个正交偏振态分别探测再在数据处理时合成。偏振分集接收后信号强度会起伏但不会完全消失。5.2 频域降噪的数值实现频谱平滑是另一个抑制衰落噪声的手段。对反射光谱做移动平均能有效压低噪声底但会牺牲空间分辨率。更推荐的做法是对幅度谱做中值滤波去除孤立的噪声尖峰同时保持反射峰的形状。% 幅度谱中值滤波 smoothed_amp medfilt1(amp_spectrum, med_filt_len); % med_filt_len 设为 3 或 5, 过大则弱反射峰会被滤掉中值滤波的窗口长度是 3 或 5窗口太大弱反射峰会被当噪声滤掉。这里和窗函数的选择互为犄角窗函数抑制频谱泄漏中值滤波抑制随机噪声两者作用不同不要混为一谈。5.3 二维强度图距离-时间与异常定位把多次扫频的距离曲线堆叠成二维矩阵用imagesc画出来横轴距离纵轴扫频次数颜色代表反射强度。这种可视化方式能一眼看出来哪个位置的反射强度在随时间变化对应光纤上发生应变或温度变化的位置。% 距离-时间强度图 % dist_matrix: 每次扫频的距离曲线, 行扫频序号, 列距离 figure; imagesc(dist_axis, scan_index, 20*log10(dist_matrix)); xlabel(距离 (m)); ylabel(扫频次数); colorbar; caxis([-80 -20]); % 动态范围设置, 根据噪声底调整caxis的范围要把噪声底压到色标之外否则微弱的散射信号变化会被噪声的颜色范围淹没。动态范围 60 dB 是 OFDR 系统比较典型的水平如果噪声底太高先回头检查辅助干涉仪的重采样效果。本文还有配套的精品资源点击获取 OFDR信号处理实战：MATLAB与LabVIEW协同实现高精度分布式光纤传感
简介面向光纤传感与OFDR技术研究用的MATLAB、LabVIEW仿真源码包聚焦光学频率域反射OFDR分布式测量方法适合科研人员、工程师及高校学生学习OFDR信号处理与传感解算。压缩包共5个文件均为m脚本整体仅2KB代码精炼涵盖OFDR基本仿真、温度响应分析、空间分辨率计算与测量距离解算等核心模块便于按需调用与二次开发。已有1877人学习下载具备一定参考热度。通过对照源码可以理解啁啾脉冲生成、傅里叶变换解调、频域到空间域映射等关键算法掌握从回波频率变化反演温度、应变及位置信息的完整流程同时可结合温度、分辨率、测距不同功能脚本快速搭建自己的OFDR数据处理原型为电力电缆热监测、桥梁结构健康检测等实际分布式光纤传感系统设计提供理论支持。1. 从拍频到距离OFDR 信号处理为什么离不开 MATLAB 和 LabVIEW光频域反射OFDR和传统 OTDR 最大的区别在于它不靠时间飞行而靠拍频频率来定位光纤上的反射点和应变点。扫频激光器发出的光经过辅助干涉仪做时钟触发主干涉仪的信号被等光频间隔采样后做一次 FFT 就能把频域映射到空间域。这个思路听起来简单但真正在 MATLAB 里把干涉信号变成距离曲线时会遇到频谱泄漏、色散失配、扫频非线性残留等一系列问题。而 LabVIEW 的角色恰恰相反它更适合在采集端做实时显示和硬件控制把 FPGA 采集卡的数据边采边画同时把原始干涉数据落盘留给 MATLAB 做离线精算。这篇文章围绕 OFDR 这个主题把 MATLAB 和 LabVIEW 两条处理链路各自该干什么、参数怎么设、坑在哪讲清楚面向的是已经接手 OFDR 系统、手里有采集数据但还没跑通完整处理流程的工程师。2. OFDR 信号模型与拍频解调的几个关键参数2.1 拍频频率与空间分辨率的映射关系OFDR 的原始信号是扫频光源经过马赫-曾德尔干涉仪后产生的拍频信号。设扫频范围为 ΔF扫频周期为 T则扫频速率 γ ΔF / T。光纤中位置 z 处的反射点与参考臂之间的延时为 τ 2nz/c对应的拍频频率为 f_b γτ。对时域信号做 FFT 后频率轴乘以系数 c/(2nγ) 就得到距离轴。空间分辨率 δz c/(2nΔF)这意味着扫频范围直接决定分辨率而不是采样率。比如 ΔF 10 nm 在 1550 nm 波段对应约 1.25 THz理论上分辨率能到 80 μm 左右但实际还要受窗函数和色散补偿精度的影响。% 参数定义扫频范围、中心波长、光纤折射率 c 3e8; % 光速 n 1.4682; % 单模光纤群折射率 lambda_c 1550e-9; % 中心波长 delta_lambda 10e-9; % 扫频范围 10 nm delta_f c / lambda_c^2 * delta_lambda; % 换算成频率扫频范围 dz_theory c / (2 * n * delta_f); % 理论空间分辨率 fprintf(理论分辨率: %.2f um\n, dz_theory * 1e6);这段代码先把波长扫频范围换算成频率扫频范围再算理论分辨率。注意群折射率 n 不是相折射率OFDR 定位用的是光在光纤中的群速度取值通常在 1.468 附近不同厂家光纤有细微差别这个值不准会直接导致距离轴整体偏移。分辨率算出来只是理论值实际 FFT 加窗后主瓣展宽相邻反射点要分开通常需要 23 倍理论分辨率的间隔。2.2 辅助干涉仪时钟触发与重采样扫频激光器不可能是理想线性的频率随时间的变化总有非线性。如果不做校正FFT 之后的峰值会展宽甚至劈裂直接毁掉分辨率。最常见的做法是用一个辅助干涉仪生成等光频间隔的时钟信号用它作为采集卡的采样时钟或触发信号让主干涉仪信号在等频率间隔上被采样。这种硬件重采样方案在 LabVIEW 里实现时通常是把辅助干涉仪的过零点检测信号接到采集卡的 PFI 端口作为外部时钟。// LabVIEW 伪代码外部时钟采样配置 // 使用 NI 采集卡的 External Clock 模式 // 时钟源: PFI0 (辅助干涉仪过零脉冲) // 采样模式: Continuous, 采样率由外部时钟决定 // 主干涉仪信号接入 AI0, 采集数据写入循环缓冲区外部时钟采样的关键点在于辅助干涉仪的延时越长每个扫频周期内的时钟脉冲数越多重采样点数越密但过长的延时会增加光源相位噪声的影响。常见的辅助干涉仪延时设在 50200 米光纤对应长度折中之后每扫频周期能采到几万到几十万个点。如果采集卡不支持外部时钟也可以在内时钟模式下采集后在 MATLAB 里用辅助干涉仪信号做插值重采样效果稍差但可行。% 基于辅助干涉仪的软件重采样 % t: 等时间间隔时间轴, ai_sig: 辅助干涉仪信号 % mz_sig: 主干涉仪信号 [~, locs] findpeaks(ai_sig, MinPeakDistance, 10); % 相邻过零点对应等频率间隔 phase unwrap(angle(hilbert(ai_sig))); % 提取瞬时相位 freq_idx round(phase / (2*pi) * max_phase); % 映射到等频率索引 resampled interp1(t, mz_sig, t(locs), spline);这里用希尔伯特变换提取辅助干涉仪的瞬时相位相位每增加 2π 对应光源频率走过一个固定间隔用这些过零点位置对主干涉仪信号做插值。MinPeakDistance要按采样率设置太小会检测到噪声毛刺太大丢失有效过零点。插值方法建议用spline而不是linear因为相位映射后相邻点间距不均匀线性插值会引入额外的高频噪声。2.3 色散失配补偿的必要性参考臂和测试臂的光纤长度不一致时色散会导致不同频率成分经历不同的群延时表现为 FFT 后峰值展宽和位置偏移。特别是在测试距离超过几百米时这个效应不可忽略。补偿方法是在频域乘上一个二次相位因子参数可以通过对已知反射点的峰值宽度做优化来标定。% 色散补偿频域二次相位因子 N length(resampled); f_axis (0:N-1) / N * fs; D 17e-6; % 色散系数, 单位 s/m^2, 典型值需要标定 phase_comp exp(1i * pi * D * (f_axis - fs/2).^2); compensated ifft(fft(resampled) .* phase_comp);色散系数 D 不是光纤 datasheet 上的色散参数而是根据系统光路算出来的等效值最靠谱的标定方式是放一个反射峰扫 D 的值让峰最窄。实践中可以把 D 的搜索范围设定在理论值的 ±50%步长逐步减半两三轮就能收敛。3. OFDR MATLAB 处理链路从原始干涉数据到应变曲线3.1 数据预处理去直流、加窗与零填充采集卡拿到的原始干涉信号带有直流偏置来自光电探测器的平均光功率不先去掉会占据 FFT 动态范围。直流分量可以直接减去信号均值但如果扫频过程中光源功率有波动均值法不够准需要做高速滤波或用多项式拟合基线。% 去直流与加窗 mz_ac mz_sig - mean(mz_sig); % 去直流 win hanning(N, periodic); % 周期汉宁窗 mz_win mz_ac .* win; % 加窗 % 零填充到 2 的幂次方便 FFT Nfft 2^nextpow2(N * 4); spectrum fft(mz_win, Nfft);汉宁窗是 OFDR 里的默认选择主瓣比矩形窗宽但旁瓣抑制好能压住距离轴上反射峰旁边的假目标。periodic选项比symmetric更适合频谱分析因为它在离散傅里叶变换下的泄漏特性更好。零填充不会提高真实分辨率但让频谱看起来更平滑峰值定位精度能提升到亚像素级别。3.2 频域转距离域与峰值提取做完 FFT 之后频谱的横轴频率通过 f γτ 映射到距离。这里有一个容易被忽略的问题FFT 结果包含正负频率OFDR 信号经过希尔伯特变换或 IQ 解调后通常是单边带的但如果是直接用光电探测器采集的实数信号频谱会左右对称需要取单边谱并把幅值乘以 2。% 频域转距离域 f_res fs / Nfft; % 频率分辨率 f_axis (0:Nfft/2-1) * f_res; % 单边频率轴 dist_axis f_axis * c / (2 * n * gamma_sweep); % 距离轴 amp_spectrum abs(spectrum(1:Nfft/2)) * 2 / N; % 单边幅值谱 % 峰值检测 [pks, locs] findpeaks(amp_spectrum, MinPeakHeight, threshold, ... MinPeakDistance, round(dz_theory / (dist_axis(2)-dist_axis(1))));MinPeakDistance取值要大于理论分辨率对应的像素数否则同一个反射峰可能被检测成多个峰。阈值threshold一般设为噪声底以上 10 dB噪声底可以用距离轴上无反射区段的幅值统计出来。峰值检测出来的位置配合相位信息就能算出亚像素精度的距离值——OFDR 的定位精度可以做到远小于空间分辨率靠的就是相位。3.3 滑动窗 FFT 与分布式应变解调分布式应变测量的核心思路是把采集到的干涉数据分成若干段每段做一次 FFT得到不同位置处的反射光谱。当光纤某一段发生应变时对应位置的反射光谱会发生频移通过互相关计算这个频移量就能还原出应变分布。窗长决定了空间分辨率和应变分辨率之间的取舍。% 滑动窗 FFT 提取局部反射光谱 win_len 512; % 窗长, 决定空间分辨率 step 64; % 滑动步长, 决定采样间距 num_windows floor((N - win_len) / step); for k 1:num_windows idx (k-1)*step 1 : (k-1)*step win_len; seg mz_ac(idx) .* hanning(win_len, periodic); spec fft(seg); % 记录局部光谱, 用于后续互相关计算频移 spectra(:, k) spec(1:win_len/2); end窗长选 512 还是 2048取决于你要的定位粒度。假设每段对应的距离长度是 ΔL应变分辨率大约正比于 ΔL 的平方根所以追求高空间分辨率时应变灵敏度会下降。互相关计算频移时频移量与应变的换算系数大约为 0.78 GHz/με1550 nm 波段这个系数受光纤光弹系数影响不同类型的传感光纤有区别最好先在拉伸台上标定。3.4 频谱拼接与长距离测量优化测试距离越长单次扫频的数据量越大FFT 需要的点数也越多。距离超过 1 km 时一次扫频往往包含几百万个采样点直接 FFT 虽然能做但速度慢。常见做法是把数据分段做处理取平均或者在频域做子带滤波后降采样只保留感兴趣的距离范围。% 频域子带提取只关心特定距离范围 dist_min 100; dist_max 200; % 关心的距离范围, 米 freq_min 2 * n * dist_min * gamma_sweep / c; freq_max 2 * n * dist_max * gamma_sweep / c; idx_range find(f_axis freq_min f_axis freq_max); % 对子带做 IFFT 得到时间域信号, 降低后续处理数据量 sub_band spectrum(idx_range); sub_time ifft(sub_band);这种做法的好处是后续的应变计算只需要处理这一个子带的信号计算量能降一个数量级。注意 IFFT 之后距离轴的起点不再从零开始需要把偏移量加回去。4. 用 LabVIEW 搭 OFDR 实时采集与显示系统4.1 采集卡选型与采样参数配置LabVIEW 在 OFDR 系统里的定位是实时采集平台。采集卡的选型要看三个指标采样率至少 100 MS/s 以上取决于扫频速率和最大拍频频率、分辨率 14 bit 或更高、支持外部时钟输入。拍频频率的上限由最远测量距离决定如果最远测 100 米2nLγ/c 算出来的频率要在采集卡带宽范围内同时要留 20% 余量给滤波器滚降。// 关键参数配置示意 // DAQmx Create Channel AI Voltage // 采样率: 200 MS/s (内部时钟) // 触发源: PFI0 (辅助干涉仪时钟) // 采样模式: Continuous // 每通道采样数: 0 (连续模式不设单次采样数)外部时钟模式下AI 通道的采样率参数变成参考时钟的最大允许速率实际采样率由外部时钟频率决定。数据读取用循环缓冲方式每次读取最新的一帧数据做 FFT 显示同时把原始数据写入 TDMS 文件。连续采集时最怕的是读写速度不匹配缓冲区溢出会导致数据丢帧。读取速率要略快于产生速率并用队列机制做生产者消费者模式。4.2 生产者-消费者架构与实时频谱显示实时显示 OFDR 距离曲线核心是一个标准的生产者消费者模式采集循环往队列里写入原始数据块处理循环从队列取出数据做 FFT 和坐标映射再把结果送到波形图控件。队列深度要按采样率和每块数据的处理时间来定处理时间超过获取时间时队列就会持续增长最终耗尽内存。// 生产者循环: 读取 DAQ 数据 // DAQmx Read (NChan NSamp) - 队列插入 // 队列大小: 100 个数据块, 超出则丢弃最旧块 // 消费者循环: // 队列取出 - FFT - 距离轴映射 - XY 图显示 // FFT 使用 VI: Spectrum Measurements.vi // 窗口类型: Hanning, 输出单位: 幅值谱XY 图显示距离曲线时要注意FFT 输出的频率轴是均匀分布的但映射到距离轴后仍然是均匀的直接作为 X 轴数组传入就行。如果显示 10000 个点以上的曲线用Waveform Chart刷新率会下降建议用XY Graph一次性刷新整帧数据或者用Image控件做强度图横轴距离、纵轴扫频次数能直接看到整个测量过程中的变化。4.3 TDMS 文件落盘与 MATLAB 离线数据对接LabVIEW 采集的原始数据最终要交给 MATLAB 做精细处理。TDMS 格式是 NI 的推荐格式MATLAB 原生不支持直接读取但可以通过tdmsread函数MATLAB R2019a 之后内置读取或者用第三方工具包。更稳妥的方案是在 LabVIEW 里直接把数据写成二进制文件同时附一个包含采样率、扫频参数、折射率等元数据的文本文件。// TDMS 写入配置 // TDMS Open - TDMS Write - TDMS Close // 分组名: OFDR_Data, 通道名: MainInterferometer, AuxiliaryInterferometer // 属性: SweepRate, CenterWavelength, RefractiveIndex% MATLAB 读取 LabVIEW 写入的 TDMS 文件 data tdmsread(scan_001.tdms); mz data.MainInterferometer.Data; ai data.AuxiliaryInterferometer.Data; info data.MainInterferometer.Property; % 读取元数据落盘时最容易被忽略的是数据类型和字节序。LabVIEW 默认小端序MATLAB 的fread也要对应指定ieee-le。如果用 TDMS 格式加上属性通道后文件会包含完整的记录信息但文件体积会比纯二进制大 10%20%。对长时间连续采集的场景建议 LabVIEW 每扫频周期写入一个文件避免单个文件过大导致读取时内存不足。4.4 扫频光源触发的同步策略OFDR 系统需要保证数据采集和扫频光源严格同步。通常的做法是让光源输出一个扫频开始触发电平接到采集卡的 PFI 触发端口采集卡从第一个触发沿开始记录一帧完整的数据。连续测量时每次触发代表一个扫频周期开始帧间通过触发沿对齐数据处理时每帧独立做 FFT。5. 光纤瑞利散射信号的衰落噪声抑制与图像化显示5.1 偏振衰落与频谱平均的取舍OFDR 信号来自瑞利散射散射光的偏振态随机变化如果探测光的偏振方向与参考光不一致拍频效率会降低表现为局部信号衰落甚至完全消失。最常见的处理办法是偏振分集接收把信号分成两个正交偏振态分别探测再在数据处理时合成。偏振分集接收后信号强度会起伏但不会完全消失。5.2 频域降噪的数值实现频谱平滑是另一个抑制衰落噪声的手段。对反射光谱做移动平均能有效压低噪声底但会牺牲空间分辨率。更推荐的做法是对幅度谱做中值滤波去除孤立的噪声尖峰同时保持反射峰的形状。% 幅度谱中值滤波 smoothed_amp medfilt1(amp_spectrum, med_filt_len); % med_filt_len 设为 3 或 5, 过大则弱反射峰会被滤掉中值滤波的窗口长度是 3 或 5窗口太大弱反射峰会被当噪声滤掉。这里和窗函数的选择互为犄角窗函数抑制频谱泄漏中值滤波抑制随机噪声两者作用不同不要混为一谈。5.3 二维强度图距离-时间与异常定位把多次扫频的距离曲线堆叠成二维矩阵用imagesc画出来横轴距离纵轴扫频次数颜色代表反射强度。这种可视化方式能一眼看出来哪个位置的反射强度在随时间变化对应光纤上发生应变或温度变化的位置。% 距离-时间强度图 % dist_matrix: 每次扫频的距离曲线, 行扫频序号, 列距离 figure; imagesc(dist_axis, scan_index, 20*log10(dist_matrix)); xlabel(距离 (m)); ylabel(扫频次数); colorbar; caxis([-80 -20]); % 动态范围设置, 根据噪声底调整caxis的范围要把噪声底压到色标之外否则微弱的散射信号变化会被噪声的颜色范围淹没。动态范围 60 dB 是 OFDR 系统比较典型的水平如果噪声底太高先回头检查辅助干涉仪的重采样效果。本文还有配套的精品资源点击获取