Home/Projects/zprofile
D

zprofile

Project D · An inverse-problem algorithm on the foundation layer
04 / 07

Reconstructs a cable's impedance profile along its physical length from frequency-domain measurements.

+4 to +9 Ωbias of the plain transform±0.3 Ωafter reconstruction5reference cables14tests17 pptechnical report
Problem

The problem

A frequency-domain measurement says how much is wrong with a cable; it cannot say where. For a manufacturer, the difference matters enormously: "return loss fails at 391 MHz" starts an investigation, while "there is a defect 6.0 m from the near end, 0.25 Ω ripple every 26 cm on top" ends one.

The obvious tool — an inverse FFT into the time domain, as a VNA's "time-domain option" does — is quietly biased: on a lossy line the step response keeps rising as √t, so a 15 m automotive pair reads 4 to 9 Ω high by the far end. Real TDR instruments show the same slope. The instrument is not lying; the plain transform is just not the right inverse.

System

What was built

A reconstruction pipeline in which every processing choice is explicit, and the one that matters most — undoing loss — is an algorithm, not a correction factor.

S_dd11(f), and S_dd21(f) when availablethe differential reflection of the cable, fixture already removed
DC anchoringfour estimators; the default ties Γ(0) to the cable's own loop resistance
Windowing + rise-time matchinga catalogue of windows with numerically computed rise time, side lobes and overshoot at the actual bandwidth
Loss modelα = c/√f + a√f + b·f fitted from transmission, with a consistent complex Z₀(f) and γ(f)
Loss-aware layer peelingeach reflection referred back through the attenuation and delay of everything the wave already crossed — with regularised inverse propagation
Connector masking + featuresfixed masks and a settling detector; ripple period and amplitude; defect position, depth, extent
Z(x): impedance against positionwith deviation statistics against the nominal impedance
Takes
  • Differential S-parameters from cablecheck's pipeline
  • The cable's physical length and, optionally, its loop resistance
  • A window and bandwidth choice — with the trade-off quantified
Produces
  • Characteristic impedance metre by metre along the cable
  • Ripple period and amplitude — the stranding lay, made visible
  • Defect position, depth and extent
  • A sensitivity budget: how each choice moves each statistic
Engineering

The hard part

Layer peeling walks the cable cell by cell: read the first reflection, deduce the first impedance step, remove that interface's effect from everything that follows, repeat. On a lossless line the recursion is textbook. On a real line, every later echo has been attenuated and dispersed by everything before it, so the algorithm must amplify each buried reflection by exactly the loss it suffered — and amplification of a measured signal amplifies noise with it. The reconstruction therefore regularises the inverse propagation (Wiener-style, with an adaptive read time) and corrects each cell's read against a lossy remainder — the difference between an algorithm that works on paper and one that survives −85 dB of instrument noise.

Just as important is what surrounds the algorithm: a sensitivity budget. Every processing choice — window, bandwidth, DC estimator, mask — is varied, and the shift of every reported statistic is tabulated. The result does not just say "the mean is 100.2 Ω"; it says how much of that could be the processing.

Truth, plain transform and loss-aware reconstruction of a 15 m cable with ripple and a defect at 6 m.
FIG. 01  The whole argument in one figure: the plain transform (biased, rising) against the loss-aware reconstruction (flat, on truth), with the planted defect recovered at 6.0 m.
Evidence

Checked against ground truth

Five synthetic reference cables with exactly known profiles, plus a simulated TDR instrument as an independent path.

What was checkedResult
Plain transform, 15 m automotive pairbiased +4 to +9 Ω by the far end
Loss-aware reconstruction, same cableswithin a few tenths of an ohm
Ripple period and amplituderecovered once bandwidth resolves them
Defect position, depth and extentrecovered — position to the resolution cell
Simulated 35 ps TDR instrument (independent path)same story: instrument shows the bias, reconstruction removes it
Deviation statistics vs nominalimplemented and reported per profile
Reconstruction of a thin, lossy 15 m cable where the plain transform fails badly.
FIG. 02  The hard case: a thin, lossy pair — the plain transform is unusable; the reconstruction holds.
Sensitivity of the reported statistics to each processing choice at 600 MHz.
FIG. 03  The honest table: what each processing choice does to each number. Few tools publish this.
Declared limits

What it does not claim

From the report's own limitations section:

  • Resolution is bandwidth-bound: features shorter than the resolved rise distance blur, stated per sweep.
  • The loss model is smooth (√f + f); strongly resonant structures would need more.
  • Validation is against synthetic truth and a simulated instrument — the same honest boundary as the whole toolchain.
Context

Where it sits in the toolchain

Built directly on cablecheck's network layer — and it repaid the debt: building this project exposed a half-cell TDR bug in the foundation, worth up to 1 Ω, now fixed and pinned by test. labauto runs it as an analyser on every archived measurement, and cablecheck can use it as its impedance-profile engine (--profile-engine zprofile) when installed.

← PreviousC · shieldeval