企业官网建设流程全解析

最近有朋友问起用MATLAB/Simulink做行波故障测距仿真到底该怎么上手。我前前后后把模型、算法、误差来源都捋过一遍，踩了不少坑，也积累了一些可以复用经验。这篇文章就把整个探索过程摊开来讲，从原理、模型搭建、波头识别到故障场景验证，按我实际操作的顺序走一遍。想用Simulink复现行波测距算法的朋友，或者正在做输电线路故障定位相关课题的电力方向学生，应该都能从中找到直接用得上的东西。

先说清楚一句话：行波测距不是“测电压算距离”，而是把故障产生的电磁波传播过程本身当作测量工具。它比传统阻抗法精度高得多，但对采样率、线路模型、波头标定的要求也苛刻很多。仿真里大量细节没做对，结果就完全不对。下面这些内容基本覆盖了一条完整可行的实践路线。

1. 行波测距的原理与方案选型

1.1 输电线路故障定位为什么这么难

输电线路动辄几十上百公里，故障发生之后，保护装置能很快判断“有故障”，但判断“故障到底在哪”是另一回事。早期用得最多的阻抗法，本质是靠测量故障回路阻抗来估算距离。这个思路在理论上是通的，但一到现场就有很多干扰因素：过渡电阻会抬升测量阻抗，对侧电源注入电流会改变阻抗角，线路负荷也会影响计算精度。实际运行中，阻抗法定位误差到公里级甚至十几公里都是常态，巡线人员拿着仪器在几十公里范围内找故障点，效率非常低。

行波法切入的角度完全不同。它是这样想的：故障发生的瞬间，故障点会出现一个电压突降，这个突降会以电磁波的形式沿线路向两端传播。这个电磁波的传播速度接近光速，在架空线中大概每微秒跑0.3公里左右。如果我们能精确测出这个波到达线路两端的时间差，或者测出它从故障点反射回来的时间，就可以用非常简单的几何关系算出故障距离。

这里面的核心逻辑可以类比成一维空间里的声波测距：你在隧道里喊一声，听到回声的时间差乘以声速，就能估算离墙壁多远。行波测距也是如此，只不过把声波换成了电磁波，把声速换成了波速，而且波速比声速快一百万倍。

1.2 行波测距的A/B/C/D型方法差在哪

行波测距在教科书里通常分成四种典型方式，工程上叫A型、B型、C型、D型。我把它整理成一张表，方便对照理解：

我在仿真里主要验证的是A型和D型。A型能完整展示行波的反射过程，适合搞清原理；D型精度高、实现简单，是研究算法的主线。B型从算法本质上看和D型是一样的，差别只在工程实现上如何同步两端时钟，仿真里不需要刻意分开处理。C型涉及高压脉冲注入设备，仿真时需要额外建模激励源，后续有机会再单独展开。

先说D型的测距公式，它比A型简单，也更容易理解。设线路长度为L，波速为v，故障点距离A端为d。故障发生后，行波从故障点分别向A端和B端传播，到达时间分别为tA和tB。显然，如果故障点离A端近，tA就小。两者的时间差Δt = |tA - tB|，对应的路程差就是v×Δt。因为tA和tB对应的路程分别是d和L-d，所以路程差等于|d - (L-d)| = |2d - L|。把这个式子反解出来就是：

• 若A端先收到行波，则故障点离A端更近，d = (L - v×Δt) / 2
• 若B端先收到行波，则故障点离B端更近，d = (L + v×Δt) / 2

这个推导看似简单，实际隐藏了一个关键点：只需要识别第一个到达的初始行波，不需要管后续的反射波。所以D型的可靠性天然比A型高，因为初始行波是故障发生后最先到达测量点的信号，幅值最大、特征最明显，不容易被其他折反射波干扰。

A型的公式也简单，但理解上要多一层。单端检测到初始行波后，这个波会继续向故障点方向传，到达故障点后反射回来，再次到达测量端。设初始行波到达时刻为t0，故障点反射波到达时刻为t1，那么行波从测量端到故障点再返回的往返时间就是t1 - t0。单程距离就是d = v×(t1 - t0) / 2。注意这里有个因子2，因为波走了两趟。

1.3 为什么选MATLAB/Simulink而不是其他工具

做行波测距仿真，可选工具其实不少。PSCAD是电力系统电磁暂态仿真的老牌工具，建模能力强，很多论文里的行波测距仿真都是用PSCAD做的；EMTP系列也一样。但我的实际体会是，如果研究重点在“测距算法”而不在“电磁暂态建模”本身，MATLAB/Simulink的组合有它独特的优势。

行波测距的完整流程是“建模仿真 → 数据导出 → 波头标定 → 测距计算”。Simulink负责线路建模和波形生成，MATLAB脚本负责小波变换、模极大值搜索、公式计算这些算法工作。这两个环境之间是无缝衔接的，Simulink的数据可以直接通过To Workspace模块送到MATLAB工作区，不需要像PSCAD那样额外导出数据文件再处理。

另外，MATLAB自带Wavelet Toolbox，而小波变换恰恰是行波波头标定最常用的工具。如果不用MATLAB，还得另找小波分析库，或者手写Mallat算法，工作量会明显增加。Simulink里的Simscape Electrical专门有分布参数线路模型，可以模拟行波的传播和折反射过程，这正好是行波测距仿真最核心的建模需求。

做这个方向，我建议先接受一个观点：仿真不是越复杂越好。初期用无损分布参数线路模型把整个流程跑通，后期再逐渐加入线路损耗、频变参数、噪声等因素，这是比较稳妥的路线。上来就堆复杂模型，波形乱七八糟，反而很难验证算法对不对。

2. Simulink行波测距模型搭建过程

2.1 线路模型选择：分布参数是底线

行波测距仿真模型里最关键的模块就是输电线路模型。很多人第一次搭模型时习惯从Simulink的元件库里拖一个“Three-Phase PI Section Line”出来，毕竟这个模块参数简单，看起来也直观。但我要特别提醒：做行波测距仿真，千万别用PI分段集中参数模型。集中参数模型把线路的分布电感、电容合并成几个集中元件，它只能反映线路在工频附近的电气特性，无法正确模拟行波在导线上的传播过程。用PI模型跑出来的故障波形，要么看不到清晰的行波波头，要么反射波完全对不上，后面所有算法都没法验证。

需要用的是Simscape Electrical里的分布参数线路模型。这个模型内部用特征线法（Bergeron方法）求解波动方程，能够正确描述行波沿线传播、反射、折射的物理过程。在Simscape > Electrical > Specialized Power Systems > Power Grid Elements里可以找到“Distributed Parameters Line”模块，三相和单相都有。模块参数里主要设置每公里电阻、电感、电容、电导，以及线路总长度。

对于行波测距的原理验证，我建议先用无损模型，也就是把R和G设成0。这样做的原因是波速只由L和C决定，而且波形中没有衰减和色散，波头非常清晰，适合先把测距算法跑通。等算法没问题了，再把R和G加回去，看算法在有损耗情况下的鲁棒性。如果你一上来就用带损耗的频变参数模型，波形中高频分量被衰减得很厉害，波头识别难度会大很多，不利于排查到底是算法问题还是模型问题。

2.2 电源、故障与测量模块的搭建细节

标准的三相行波测距仿真模型，按信号流动方向大致是：三相电源 → 线路A段 → 故障点模块 → 线路B段 → 线路末端负载/电源。故障点放在两条线路段之间，通过改变两条线路段的长度比例，就可以控制故障距离。

三相电源用Specialized Power Systems里的“Three-Phase Source”模块。这里有一个比较容易忽略的细节：电源内阻抗会影响行波在母线端的反射特性。实际系统中，母线侧有变压器、断路器、电压互感器等设备，这些设备对高频行波呈现的等效阻抗很复杂。仿真里至少要设置一个合理的内部阻抗，不能把电源理想化到内阻抗为零。我用的是变压器等效阻抗折算后的数值，具体大小对行波测距仿真结果有影响，但不影响算法的验证逻辑。

故障模块用“Three-Phase Fault”模块。故障类型、故障相、故障起始时间、过渡电阻都要在这里设置。行波测距对故障起始时间很敏感，因为我们要用它作为时间零点来测量行波到达时刻。在仿真里，故障触发时刻是精确已知的，但在实际系统中这个时刻未知，需要通过波形特征来推断。这一点在做结果分析时要区分清楚。

测量模块用“Three-Phase V-I Measurement”来测量故障点两侧的电压和电流。行波测距通常以电压行波为主，因为电压行波幅值大、波头陡，更容易标定。电流行波在仿真里也能看到，但实际工程中电流互感器对高频分量的传变特性远不如电压互感器，所以现场装置普遍优先处理电压行波。仿真阶段可以同时记录电压和电流，但算法验证时优先用电压数据。

数据导出用“To Workspace”模块。这个模块要把存储格式设置为Timeseries，变量名根据后续MATLAB脚本的读取方便来命名。我习惯用vA、vB、vC这样带相别和后缀的名字，后面写脚本时一眼就能看出对应关系。

2.3 采样步长与求解器设置的数值依据

行波测距对采样率的要求非常高。现场测距装置的采样率通常在1MHz到10MHz这个量级，也就是每微秒采1到10个点，原因很简单：波速约为0.3 km/μs，如果采样间隔是1微秒，那么时间测量误差带来的距离误差最大就是0.3公里。换个说法，采样率每提高一个数量级，测距误差上限就降一个数量级。

Simulink里固定步长求解器的步长就相当于采样间隔。我把仿真步长设为1e-7秒，也就是100纳秒，等效采样率10MHz。这时候时间量化误差对应的距离误差大约是30米左右，对于验证行波测距算法来说已经够了。如果步长设成1e-6秒（1MHz采样率），时间量化误差是300米，作为对比实验观察采样率对精度的影响也很有价值。

固定步长求解器我选择ode4或ode5。选择固定步长而不是变步长，一是保证数据点均匀，便于后续小波变换处理；二是避免求解器自适应调整步长时对高频行波产生数值畸变。行波仿真本身就是模拟高频过程，必须用固定步长才能得到可靠结果。

仿真时长的设置也有讲究。以100公里线路为例，行波从一端传到另一端大约需要337微秒。如果我们要观察初始行波、故障点反射波、对端反射波等多个波头，仿真时间至少要有几毫秒。我一般设置成0.02秒，也就是20毫秒。这个时长足够观察到多次折反射，又不会因为仿真时间过长导致数据文件过大。需要说明的是，故障前还需要留一段预故障时间，让模型进入稳态。故障起始时间一般设在5毫秒左右，故障前的5毫秒数据可以作为正常波形参照。

2.4 三相线路的线模提取与数据导出

三相线路中，故障行波在A、B、C三相之间是相互耦合的，直接拿某一相的电压波形去分析，波头混叠比较严重。工程上行波测距装置会利用相模变换（Karenbauer变换或Clarke变换）把三相电压解耦成线模分量和零模分量，然后选取线模分量进行测距。

为什么选线模而不是零模？这里有个物理原因。零模分量是通过大地回路传播的，受大地电阻率影响很大，传播速度不稳定，高频分量衰减快；线模分量通过导线之间传播，特性更接近理想的无损传输线，波速稳定、衰减小，所以测距精度更高。在仿真中做线模提取很简单，一个近似的线模分量可以用两相电压之差来近似，比如vMode = vA - vB。更严格的做法是用Karenbauer变换矩阵，把三相量变换成0模、1模、2模分量。实际算法验证时，我推荐用完整的相模变换，这样更贴近现场装置的处理流程。

模型仿真结束后，电压行波数据通过To Workspace模块导出到MATLAB工作区。导出的数据是Simulink的Timeseries格式，包含时间向量和数值向量。在MATLAB脚本里，用vA.Data获取电压数值，用vA.Time获取时间值。这里要提醒一句：如果Timeseries变量的维度不对，常见的错误是每个通道的数据被存成了多维数组，需要先用squeeze函数把多余的维度去掉，否则后面小波变换会报维度不匹配的错误。

3. 波头识别与测距算法的MATLAB实现

3.1 小波变换为什么适合标定行波波头

波形拿到手之后，下一个问题就是：故障发生后，行波到底在哪个时刻到达测量点。听起来简单，实际做起来并不容易。故障后的原始波形里同时包含工频分量、高频行波分量、可能的噪声和多次折反射波。直接从原始波形找幅值最大的突变点，往往会找错，因为第一个到达的初始行波不一定幅值最大，后续反射波叠加后可能更大。

小波变换能派上用场，是因为它具有“局部化”分析能力。行波到达测量端时，波形上会出现一个很陡峭的突变点，这个突变点在各个频率尺度上都会产生较大的小波系数。小波变换可以把信号分解成不同频率成分，并在每个频率分量上定位突变发生的时刻。特别地，行波波头在本质上是一个阶跃或冲击信号，它的奇异性很强，小波变换的模极大值正好对应这个突变点的位置。

用生活化的方式理解：小波变换就像一个“多重分辨率的放大镜”。低频分量看波形的整体趋势，高频分量看波形中突然变化的地方。故障行波到达的瞬间，高频分量会突然增强，小波系数在这个时刻出现一个峰值。通过查找这个峰值对应的时间，就能确定行波到达时刻。

这里要说一个常用的判断原则：用小波变换标定波头时，通常看的是初始行波的第一个模极大值，而不是全周期里的最大值。后续反射波、对端折射波都可能产生很大的小波系数，如果直接在全时间范围内找最大值，很容易标定到后续的反射波上。正确做法是把分析窗口截断在故障发生后的短时间内，比如1到2毫秒内，在这个窗口内找第一个显著模极大值。

3.2 用MATLAB实现波头标定的参考代码

下面这段代码是我在探索过程中整理出来的核心标定流程，以较新的MATLAB版本为例。假设已经通过Simulink的To Workspace模块导出了线模电压信号vMode，采样率为2MHz，故障发生时刻已知为tFault。

% 读取Simulink导出的行波数据 % vMode: 线模电压信号，Timeseries格式 % fs: 采样率，单位Hz fs = 2e6; tAll = vMode.Time; vAll = vMode.Data; vAll = squeeze(vAll); % 去掉多余的维度 % 截取故障后1ms窗口内的数据 tFault = 0.005; % 故障发生时刻，单位s idxWindow = tAll >= tFault & tAll <= tFault + 1e-3; tWin = tAll(idxWindow); vWin = vAll(idxWindow); % 连续小波变换，关注高频段 [coef, freqs] = cwt(vWin, fs); % 选择50kHz以上的尺度，取对应小波系数的绝对值 scaleIdx = find(freqs > 50e3, 1, 'last'); absCoef = abs(coef(scaleIdx, :)); % 在窗口内找第一个局部模极大值，作为波头时刻 [pks, locs] = findpeaks(absCoef); tHead = tWin(locs(1)); % 第一个峰值对应初始行波到达时刻

这里有几个细节值得说明。findpeaks函数在默认情况下会返回所有局部峰值，如果波形中有其他干扰，可以在函数参数里指定MinPeakHeight或MinPeakDistance来过滤掉一些不明显的峰值。波头标定的稳定性和窗口长度关系很大，窗口太短可能把多次反射波截掉，窗口太长又会引入后续大波峰干扰判断，需要根据线路长度来调整。

如果你不想用小波变换，MATLAB里还有一个简单实用的函数findchangepts，可以检测信号中的突变点。它不需要选择小波基函数和尺度，上手快，但精度和对参数的敏感性控制不如小波变换细腻。我的体验是，先用findchangepts粗标定，再用小波变换细验证，两者结果一致时基本可以确定波头位置没问题。

3.3 单端测距与双端测距的计算实现

波头标定完成之后，测距计算就非常简单了，核心就是套公式。双端测距需要两个端点的波头到达时间，单端测距需要一个端点上的初始波头和反射波波头。

双端测距的MATLAB实现步骤大致如下。假设A端波头标定时刻为tA，B端波头标定时刻为tB，线路总长为L，波速为v，采样间隔为dt。先判断哪端先收到行波：

% 判断先到端 if tA < tB dtArrive = tB - tA; dA = (L - v * dtArrive) / 2; dB = L - dA; else dtArrive = tA - tB; dA = (L + v * dtArrive) / 2; dB = L - dA; end

需要注意，上面的tA和tB都必须是绝对时间，而不是以故障时刻为基准的相对时间。如果两个端点的数据采集起始时刻不一致，会直接导致测距错误。在仿真里，两端数据往往来自同一个模型，时间基准天然一致，但在实际系统中，两端必须通过GPS或北斗授时来同步时钟。

单端测距实现的关键是正确识别故障点的反射波。设初始行波到达时刻为t0，反射波到达时刻为tRef，则测距距离为：

d = v * (tRef - t0) / 2;

这个公式看似简单，但tRef的识别如果选错了波头，结果会差得离谱。我后面专门总结了一些反射波识别的经验。

3.4 波速如何确定：理论计算与实际修正

波速是行波测距公式里另一个决定性参数。波速不准，再准的波头识别也白搭。理想情况下，无损线路上的行波波速等于：

v = 1 / sqrt(L×C)

其中L是线路单位长度的电感，C是单位长度电容。这个公式是从波动方程直接推出来的，和电磁波在介质中的传播速度公式形式上一致。以一条典型架空线路为例，取L=0.9716 mH/km，C=0.01216 μF/km，计算如下：

L×C = 0.9716×10⁻³ × 0.01216×10⁻⁶ ≈ 1.181×10⁻¹¹ (s²/km²)

开根号得到约3.437×10⁻⁶ s/km，取倒数就是约2.91×10⁵ km/s，换算成m/s就是2.91×10⁸，略小于真空光速。这个数值和仿真中观察到的波头传播时间基本吻合。

我实测过不同线路参数对测距结果的影响。如果把理论波速直接用进测距公式，在100公里线路上通常会有几百米量级的误差。这个误差不完全来自波速公式本身，而在于实际线路是有损线路，频变参数导致不同频率分量的传播速度略有差异，也就是色散效应。工程上行波测距装置一般会用两种方式处理：一是根据历史故障数据反推校准波速；二是用双端测距结果和已知故障距离进行比对，反算出等效波速。

在仿真阶段，我建议采取“理论计算波速 + 仿真结果校准”的方式。先用理论波速跑一次，看测距误差是正偏还是负偏，然后根据固定误差修正波速数值。这个修正后的波速在后续不同故障场景中应该保持一致，如果不同场景下修正系数差异很大，说明算法中还有其他系统性问题需要排查。

4. 仿真结果分析与关键影响规律

4.1 典型故障场景下的波形特征

以一条100公里、分布参数、无损三相线路为例，故障点设在距离A端30公里处，A相在0.005秒时刻发生经10欧姆过渡电阻的单相接地故障。仿真步长设为100纳秒，记录两个端点的线模电压信号。

故障发生后大约100微秒，A端就能看到第一个行波波头。为什么是100微秒而不是传播时间30公里除以0.3公里每微秒正好100微秒？因为这个波头就是故障行波沿着30公里线路直接传到A端的，传播时间为30/0.297约等于101微秒，和理论计算一致。紧接着，波头后面会跟着一串逐渐衰减的高频振荡，这些是行波在故障点、A端母线之间反复反射形成的。

在B端，初始行波大约在337微秒左右到达，因为故障点距离B端为70公里，传播时间约70/0.297约等于235微秒。等一下，这个计算似乎不对。故障距离A端30公里，那么距离B端是70公里。70/0.297约等于235微秒，不是337微秒。337微秒是行波跑完整条100公里线路的时间。所以B端初始行波到达时刻约为故障后235微秒左右。

从这个例子可以看到，行波测距的精度本质上取决于两个量：波头标定的时间精度和波速的准确性。两者换算成距离误差后叠加，就是最终的测距误差。

4.2 故障距离、过渡电阻和采样率的影响

我把一组典型实验结果整理在下面，目的是展示几个关键参数对测距精度的敏感程度。表中所用线路为100公里无损分布参数线路，双端法测距，理论波速校正后作为基准值。每个场景重复10次取平均值。

从数据里能看出几个规律。故障距离对误差没有明显的单调影响，远距离和近距离的误差都在同一量级，这说明只要波头标定正确，测距精度基本不随故障位置退化。过渡电阻从0.1欧姆增加到100欧姆，误差略有增加但不大，这是行波法相比阻抗法的核心优势。采样率的影响最直观：2MHz采样率下误差约为0.1到0.3公里，降到1MHz后误差明显变大，升到5MHz后误差进一步缩小，这符合时间量化误差与采样率的反比关系。

需要特别强调的是，这里的误差包括波头标定误差和波速近似误差的综合效果。如果只考虑采样率引起的量化误差，2MHz采样间隔0.5微秒，按波速0.297公里每微秒折算，最大量化误差约0.15公里，和表中观测到的误差量级匹配。

4.3 行波法的精度边界在哪里

行波法也不是万能的。它的精度边界主要受三个因素限制：采样率、波速准确性和波头标定能力。

采样率是硬件层面的天花板。现场测距装置用1MHz采样率，理论上时间量化误差是1微秒，对应约0.3公里的测距上限误差。如果装置用10MHz采样率，理论上限缩到约0.03公里。但采样率高了，数据量也大，对装置的存储和计算能力提出了更高要求。所以现场装置会做权衡，一般不是单纯追求最高采样率。

波速准确性的问题前面说过，理论计算波速和实际等效波速之间会有一个偏差。现场装置的惯用做法是定期用已知故障或者人工短路试验来校正波速。仿真里我们虽然可以用理想波速，但在做工程化评估时，一定要在校正后的波速基础上计算误差，不要把理想波速带来的“虚假高精度”当成算法性能。

波头标定则是算法层面的核心。行波在传播过程中会衰减，尤其是经过母线、拐点、T接点时会发生折射和反射，波头会变得不那么陡峭。这时候单纯靠小波模极大值容易产生时间偏移。我的经验是，不要只取一个小波尺度上的峰值，而是综合多个尺度的模极大值趋势来判断波头，这样对噪声和衰减的鲁棒性更好。

还有一个需要注意的点是故障初始角。故障发生在电压过零附近时，产生的行波幅值很小，波头识别难度剧增。这是一个在仿真中容易被忽略但其实非常现实的场景。我验证过电压过零附近故障和电压峰值附近故障两种情况，前者的波头标定时间明显偏移，测距误差大好几倍。做算法验证时，故障初始角至少要设为30度、60度、90度三组，看看算法在不同初始条件下的表现。

5. 仿真中的常见问题与排查实录

5.1 波形在故障后看不到明显行波波头

这是我被问得最多的一个问题。模型搭完，仿真跑起来，结果故障后的波形看起来就是一团乱麻，根本找不到那个“陡峭的波头”。排查顺序应该按下面这几步走。

第一，确认线路模型类型。如果用了PI Section Line模型，请换成分布参数线路模型，这是最常见的原因。PI模型无法正确模拟行波传播，波形当然不对。

第二，检查过渡电阻设置。如果过渡电阻设了很大，比如几千欧姆，故障产生的行波幅值会非常小，淹没在工频分量里。可以先设成0.01欧姆左右验证算法，再逐渐加大。

第三，检查故障初始角。故障点电压接近过零时，行波幅值很小，可以调整故障模块的触发时间，让故障发生在电压峰值附近。

第四，尝试提高采样率。如果采样率太低，比如只有10kHz，根本不可能观察到高频行波波头，因为每个行波脉冲只持续几百纳秒到几微秒，采样点太少就完全丢失了。行波测距至少需要1MHz以上采样率，这是硬性要求。

还有一个隐藏问题：如果观测窗口太长，MATLAB在绘制波形时会自动缩放，把初始行波那个窄脉冲压成几乎看不见的细线。处理方法是先用xlim限制时间轴到故障后几毫秒内，单独画一段窗口的波形，波头就显现出来了。

5.2 单端法总是找错反射波

单端测距的反射波识别是行波测距算法里最难的环节。故障点反射波和对端母线反射波在时间上、极性上都有差异，但实际波形里往往混在一起。我总结了几个有效的判据。

第一个判据是时间合理性。初始行波到达测量端后，故障点反射波返回测量端的往返时间不会超过2d/v，其中d是预估故障距离。如果用双端法先粗测出一个故障距离，可以用这个距离反推反射波应该出现的时刻范围，再在这个范围内找反射波。

第二个判据是幅值衰减规律。行波每经过一次反射或折射，幅值都会衰减。故障点反射波的幅值通常大于对端母线反射波在测量端观测到的幅值，因为后者经过了更长的传播路径和更多次折反射。如果看到的“候选反射波”幅值和初始行波差不多大，那很可能是对端反射波。

第三个判据是波形极性和形态。不同性质的波头在电压行波中的极性有一定规律。可以对比初始波头和候选波头的极性是否匹配预期，再用小波变换分别看两个波头在不同尺度下的系数分布，排除偶然的波形尖峰。

第四个判据也是最实用的：用双端法测距结果对单端法标定进行校验。仿真中两端数据都有，先算双端结果，再依据这个结果约束单端反射波的搜索区间。如果单端结果和双端结果差异超过一定范围，几乎可以肯定是反射波选错了。

5.3 Simulink模型运行与数据导出中的坑

建模执行过程中也有一些细节问题，单独列出来给大家避坑。

Bus Selector没有可选信号。这种情况通常是因为上游信号不是总线信号，或者总线中的信号没有被正确命名。Simulink里面，信号线必须显式命名，且上级模块的输出名称要和Bus Creator的输入名称一一对应，Bus Selector才能正确下拉出信号。如果发现下拉列表是空的，检查一下是不是连线没有对准，或者信号标签没命名。

To Workspace导出的变量维度不对。如果你记录的是三相电压，导出变量可能是一个三维数组，第几维对应哪一相需要捋清楚。我一般使用squeeze函数处理掉空维度，再按列提取各相信号。如果直接用原始数组做小波变换，很容易出现维度错误。

仿真速度太慢。分布参数线路模型加10MHz采样率，仿真0.02秒要跑很久。我的做法是先跑一个只覆盖故障后5毫秒的短仿真，用较短时长验证算法流程；等确认算法没问题后，再适当拉长仿真时间观察多次反射波。如果只是验证双端测距，故障后1毫秒的数据都足够了，没必要仿真全部折反射过程。

Simulink报“仿真停止”或“数值不稳定”的错误。这种情况多数是模型参数设置不合理，比如线路模型参数L和C数量级错误，或者是求解器步长设置得太小导致数值计算量过大。先检查线路参数的单位是否一致，再检查故障模块的过渡电阻是否设成了零（有时零电阻会让数值解出现奇异，设成0.001欧姆即可）。

5.4 从仿真到实际装置的距离有多远

最后聊一个很容易被忽略的问题：仿真跑通了，离真实系统还有多远。

真实行波测距装置的采样率通常在1MHz到10MHz之间，和仿真建议的设置一致。但有一个关键差异：仿真中的采样是理想等间隔的，实际装置的采样时钟会有抖动，这会直接增加波头标定的时间不确定性。我建议在仿真后期给波形加上随机噪声，并在波头标定算法里加入一些平滑处理，提高抗噪能力。

实际装置的传感器响应也是一个大问题。电压互感器和电流互感器对高频行波都有传变限制，实际波形里的波头比仿真波形圆滑得多。如果只用仿真中那种理想陡峭波头调算法，到现场很容易失灵。一种可行的方法是在仿真数据后端加一个低通滤波器，模拟互感器的带宽限制，再用滤波后数据测试标定算法。

双端测距还需要通信通道和时钟同步。现场两端之间的光纤通信延时、GPS对时误差都会影响测距精度。这些在Simulink里都可以通过人为加入延时和偏差来模拟。我的建议是至少加一个固定通信延时和一个随机对时误差，看看测距结果的变化范围，这样才能对算法的现场表现有个心理预期。

关于把Simulink模型生成C代码、在嵌入式装置上运行，确实是很多入门者想走通的路。但从行波测距的实际实现来看，装置核心是高速数据采集和FPGA上的并行处理逻辑，Simulink模型更多承担的是算法验证和离线诊断的角色。先把离线仿真做扎实，理解波头标定和反射波识别的每一个细节，再去碰装置实现，路线会顺畅很多。

整套仿真做下来，我最大的体会是：行波测距的数学原理极其简单，就是一个距离等于速度乘以时间的问题，但工程难度全部藏在“如何把时间测准”和“如何把波速定准”这两件事里。Simulink模型能帮你直观看到行波传播和反射的完整过程，但它不会告诉你哪个波头是故障点反射波，不会告诉你波速参数该调成多少。这些判断能力，只能在反复调整参数、反复对比误差的过程中慢慢积累。我给自己的习惯做法是每次只改一个变量，记录下测距误差的变化，而不是一次性把所有参数都调一遍，否则根本定位不到问题出在哪。这套方法也分享给你，做行波测距仿真时，保持这个耐心，结果会好很多。