Concept

Marching — where it appears

Solving a network forward in time by replacing each reactive element with a discrete companion and solving the resulting resistive network at every step. It answers questions no steady state contains, and its own error is a property of the step size that has to be measured rather than assumed.

Named by 17 essays across 6 fields — each of them below, with the objects they name alongside it.

Three capacitances, all correct: 2.000 µF, 5.814 µF and 2.105 µF at 5 V. computed by solving, not by drawing. The charge is C∞·v + Q_s·tanh(v/V_k) — a linear backbone and a polarisation that saturates — with both parameters pinned by the capacitance at zero volts and at the rated voltage, so there is no third degree of freedom to tune the answer with. The three curves are three questions. The small-signal value is the slope at the bias, which is what a ripple sees. The charge-average is the total charge moved from zero divided by the voltage, which is what a reservoir or a hold capacitor obeys. What a bridge reads is neither: it is the fundamental of the charge waveform under a one-volt test, which is a measurement condition. At 5 V they are 2.000 µF, 5.814 µF and 2.105 µF — a factor of 2.91 between the extremes, and every one of them is the capacitance.

The capacitance that is not one number

A ten-microfarad ceramic at its rated five volts is 2.000 µF as a slope, 5.814 µF as a charge average and 2.105 µF as a bridge reads it — three answers to three different questions, all correct, all called the capacitance. At zero bias the standard one-volt test alone reads 5.2 per cent low. Marched in a circuit the same part distorts as the square of the drive with no bias and in proportion to it with a bias, because the bias is what puts a second harmonic there.

frequency · Real capacitor
Twenty nanoseconds of skew is 0.28% of duty, not 0.40. computed by solving, not by drawing. The duty cycle of a relaxation oscillator whose comparator takes longer to go one way than the other, marched, against the expression that lengthens each half cycle by its own delay. The two part company immediately and by a constant factor of about 1.42: during a delay the capacitor keeps charging past the threshold it already crossed, so the next half cycle starts further out and takes longer, and the two halves partly cancel. At 200 ns of skew on a 2585 ns period the duty is 52.805 per cent where the expression says 54.004. The open circles are the same quantity in closed form — the overshoot is V(1 − (1 − β)e^(−d/τ)) and the next half starts from it — which the march reproduces to parts in ten thousand.

The delay that is two delays

Modelling one comparator delay applied to both transitions makes the two half cycles equal by construction, and the fix is a change of one line. Made, the duty cycle moves by 0.28 per cent for twenty nanoseconds of skew rather than the 0.40 the obvious expression gives, more hysteresis improves the duty cycle and worsens the volt-seconds at the same time, and ten nanoseconds of skew saturates a hundred-turn core in 231 cycles.

applied · Hysteresis
The frequency at which a pulse train becomes an average. computed by solving, not by drawing. The same 5 watts of average dissipation at every frequency, delivered 2 per cent at a time. The flat line is the steady-state answer, which does not know about the frequency. The falling curve is the marched peak junction temperature, which does. They meet at 308 Hz, and that frequency is not a property of the converter: it is a fraction of one junction time constant per period — f·τ = 0.738 at this duty, with τ = 2.40 ms, and between 0.78 and 0.56 across the duties on the slider. A hundred-kilohertz converter fits 240 periods inside that time constant, and at the top of the sweep — 10.0 kHz — the steady state is already exact to 0.46 per cent, so the averaged-power fixed point is right and this is the measurement that says why. The march puts 48 steps inside each pulse, which is what the answer is sensitive to: at six it put the boundary 19 per cent too high.

The pulse the heatsink does not feel

A thermal resistance iterated to a fixed point with a diode or a switch is a statement about a power — so it assumes that a hundred and fifty watts for two per cent of the time is three watts. The die's own heat capacity decides whether that is true, and it decides it at a frequency: above 308 hertz the junction integrates, by a hundred kilohertz the fixed point is exact to five parts in ten thousand, and at one hertz the same average power on the same heatsink puts the junction three hundred kelvin hotter.

transients · Reverse-recovery
The charge that comes back: a 0.2% dielectric, 10 s shorted, read at 900 s. computed by solving, not by drawing. The capacitor is charged to 10 V until every relaxation is complete, shorted for 10 seconds, then opened and watched. It climbs back to 20.00 millivolts — 0.2000 per cent of where it was — and the shape is the finding: it is a straight line on a logarithmic time axis, gaining 0.097 per cent of the charging voltage per decade. There is no time constant after which it is over, because there is no single time constant: one branch of the model comes to equilibrium per decade, for as many decades as the dielectric has. A decade before the reading it was at 0.1033 per cent.

The capacitor that remembers

Charge a capacitor, short it for ten seconds, open it, and it climbs back to a fifth of a per cent of where it was. Nothing leaked and nothing was gained: some of the dielectric had not finished discharging. The same defect measured as an admittance says the part is 0.593 per cent more capacitance at a tenth of a millihertz than at a kilohertz, and measured in a sample-and-hold it says a millisecond of hold costs a hundred parts per million — thirteen bits, on a part specified at nothing.

transients · Dielectric absorption
One part, three saturation currents: 1.90 A, 2.40 A, 2.69 A. computed by solving, not by drawing. B(H) is μ₀H plus a saturating magnetisation, written as a flux linkage, and the inductance drawn here is dλ/di — the slope of that flux, which is what a small signal on a direct current actually meets. It is 18.92 µH at no current, 16.73 µH at two amps and 10.43 µH at three, and it never reaches zero: the vacuum is still there, so the part falls to its air-core 0.05 µH and stays. The three marks are the ten, twenty and thirty per cent drops different manufacturers print as the saturation current — 1.896, 2.395, 2.694 amps, a spread of 42 per cent on one part.

The inductance the current decides

The rung below bounded the flux and the boundary is exact: the volt-seconds decide the flux swing whatever the material does, and a cycle whose current ripple runs from 519 mA to 75 A has the same flux excursion to better than two per cent. What saturation breaks is the relationship between that flux and the current — so a ripple the design expression puts at 514 mA is 1,022 mA at 2.56 A of load, the peak reaches 3.57 A where the part is at a fifth of its nameplate inductance, and the same part has three saturation currents depending on which per cent it was quoted at.

magnetics · Saturation
A major hysteresis loop at 9.0 A/m of coercivity, and the anhysteretic curve it closes onto. computed by solving, not by drawing. The B–H loop of a core driven sinusoidally to ±400 A/m, marched through a superposition of twenty-four play operators and drawn over the single-valued curve the two rungs below this one measured. The loop encloses 12.481 joules per cubic metre per cycle, which is the core loss and which no single-valued model can produce, because a curve has no area. The coercivity is 8.96 amperes per metre and the remanence 22.3 millitesla; both are read off the marched descending branch rather than handed in. The slider takes the threshold spread down to zero, where the two branches become one, the area falls to 9.8e-15 J/m³, and the object is exactly the core the field already had.

The area a curve cannot have

Every magnetic model in this collection is a single-valued B(H), and a single-valued B(H) cannot dissipate: ∮H dB around a curve is zero, so the transformers and inductors here have all run cold. Giving the material a second branch costs one object — a play operator, which lags the field by a threshold and is otherwise nothing — and a superposition of twenty-four of them produces a loop with 12.481 joules per cubic metre in it, a coercivity of 8.958 amperes per metre and a remanence of 22.3 millitesla, none of which was handed in. At zero threshold the whole thing collapses onto the curve the field already had, to fifteen figures.

magnetics · Magnetic loss
Which limit binds is a property of the load, and they change places near 22 nF. computed by solving, not by drawing. Each limit measured on its own, as the departure of its march from the linear one, at a load step of half the output stage's rating. The input pair's departure falls with load capacitance — a bigger reservoir holds the node while the loop responds, which is the sixth rung's own result — and the output stage's does not fall nearly as fast, because what it has to supply is the charge the capacitor wants. Below about 22 nanofarads the thermal voltage decides the answer and above it the output stage does, and nothing about the amplifier changed.

The current above which there is no impedance

The sixth rung found the impedance leaving at 10.6 mA, where the input pair's own tanh takes over and the slew rate is twice the thermal voltage times the gain-bandwidth in radians, with no design choice in it. A real output stage has a second limit that is nothing but design choice, and the two do not bind at the same load: at 0.47 nF the input pair's departure is 19.4 per cent against the output stage's 4.2, at 22 nF it is 0.9 against 2.3, and above the output stage's rating the excursion does not come back at all — 2,254 Ω for a quantity that was 37.

feedback · Capacitive load
A switched-capacitor low-pass driven past half its own clock. computed by solving, not by drawing. The clock is 1.00 MHz and the corner the rung below fitted is 1.59 kHz. The falling dashed curve is that continuous model, which knows nothing about a clock and goes on falling. The circles are the marched circuit, read at the frequency the output actually appears at. They part company past half the clock and by 992 kHz the model is 42.1 decibels wrong — an input just below the clock arrives just above direct current, in the middle of the passband, with the passband's own gain. The third curve is the exact discrete transfer function evaluated at the folded frequency, and it agrees with the march to 0.26 decibels, which is what says the march is measuring the folding rather than an artefact of itself.

The filter that samples

The rung below built a resistor out of a clock and measured two ways it is not one: a settling time, and a corner that is a capacitor ratio rather than an R–C product. Both are errors in a value and both get smaller as the design gets better. This is an error of a different kind — the arrangement is not a continuous system at all, and nothing below half the clock shows it. An input at 992 kilohertz arrives at 7.8 kilohertz with the passband's own gain, where the continuous model the rung below fitted says it is 56 decibels down.

filters · Switched capacitor
A staircase costs one Nth, computed rather than quoted. computed by solving, not by drawing, marched, with the energies rebuilt from the element laws. The charge is broken into N equal risers, each held for 16 time constants so that it completes. The measured losses are 1.00000, 0.500001, 0.250000, 0.125000, 0.0625000, 0.0312500 of ½CV² — which is 1.000004, 1.000002, 1.000001, 1.000001, 1.000000, 1.000000 times 1/N, so the law is exact to four parts in a million at the worst rather than approximately true. The fitted exponent is -1.00000 and the energy account closes to 2.17e-6 at the worst.

The half a switch keeps

The rung below found that charging a capacitor from a step loses half the delivered energy whatever the resistance, and that a ramp takes the loss down as 2τ/T with no floor. A staircase of N settled risers costs one Nth of the step, exact to four parts in a million, and the law ends at a dwell of 5.272 time constants. A switch is the other half of the same product and buys nothing at all: with the supply held at five volts and the channel conductance ramped over a thousand time constants, the loss is 1.00000000 of ½CV².

transients · Switching energy
The reading is a count of decades: 1.0288 parts per thousand of them. computed by solving, not by drawing. 17 marched tests, three families of absolute time — a tenth of a second, one second and ten seconds of short — plotted against the number of decades between the short and the reading. The families lie on one another, which is the finding: the answer is not a property of the part alone and not a property of either duration, it is a count of the decades of relaxation time the test leaves in. The line is a least-squares fit through the origin at 1.0288e-3 per decade; the model's own capacitance per decade of relaxation time is α = 1.0343e-3, which nothing in the fit was told — the fit sits 0.53 per cent under it, because the charge that comes back is shared with the slow branches it came off. The worst residual is 4.80 per cent, at the narrowest ratio drawn, and 1.23 per cent over the 8 tests that are two decades wide and read before the slowest relaxation the model has; the 3 read after it fall away to 4.49 per cent, which is where the law ends. Families a hundred times apart in absolute time differ by at most 0.84 per cent, which is the whole of the collapse.

Ten seconds, and fifteen minutes

A data sheet's dielectric absorption is quoted as a property of the part. It is not: the same modelled capacitor reads 0.4050 per cent with a tenth-of-a-second short and 0.0305 per cent with a thousand-second one, and 0.0047 against 0.2948 depending on when the reading is taken. Seventeen marched tests collapse onto one line — the recovery is 1.0288 parts per thousand for every decade between the two durations — and four dielectrics the specified test declares identical read a factor of 3.31 apart one decade away from it.

transients · Dielectric absorption
Switched on at 0.95× resonance, a Q 50 capacitor reaches 1.550 times its settled voltage. computed by solving, not by drawing. The capacitor voltage of the series circuit, switched on from rest at the crest of the drive, drawn as the tip of its arrow in the frame that turns with the drive, so that the settled state is a fixed point — the arrow from the centre, 10.067 V long, with the circle of that radius around the centre. The path is the settled arrow plus a second one turning at the circuit's own frequency and shrinking, drawn until the second is a hundredth of its first length. The voltage reaches 15.603 V in cycle 9, 1.5499 times the settled amplitude, and the tip's farthest point is 1.5499 times it. The approximation 1 + exp(−π/2Q|δ|) gives 1.5335. Marched in time by the trapezoidal rule, the network agrees with the exact solution to 1.2e-3 of the settled amplitude over its first 11 cycles.

The arrow that goes past where it settles

A phasor is where a driven resonator ends up, and counting the cycles it takes to get there says nothing about the path. Switched on from rest a little away from resonance, the capacitor's arrow circles its settled tip instead of approaching it: a Q of 50 driven at 0.8 of resonance reaches 1.854 times its settled voltage, and a Q of 200 at 1.25 switched on through zero reaches 2.182. At resonance exactly, where the settled voltage is largest, the arrow never passes its mark. So a resonance read by the largest voltage after switch-on is 1.20 times wider than its phasor says.

frequency · Phasors
The winding current a hysteretic core asks for, and the current the same circuit's single-valued core asks for. computed by solving, not by drawing. One netlist — a 22.6 volt peak sinusoid at 5.00 kHz through 4 ohms into a sixty-turn winding — marched twice, once with the core as a superposition of play operators and once with the single-valued saturating curve every earlier figure in this field used. The two currents differ in shape and not only in size: the hysteretic one leads the flux by an angle that is not ninety degrees, which is the whole of the core loss, and its peak is 59.45 milliamperes against 50.12. The energy the source delivers over a complete cycle is 17.19 microjoules for the loop and 2.5e-5 for the curve, which is zero to the resolution of the march. The slider takes the flux the drive demands from a twentieth of the material's saturation to most of it.

The core the solver has to remember

A saturating inductor's state is one number, because the current is a function of the flux linkage. A hysteretic core's is not: half an amp is one flux on the way up and a different flux on the way down, and there is no function of the current that returns the flux. So the march's state vector gains twenty-four more numbers, one per play operator, and the Newton loop is forbidden to touch them. The circuit and the bench then agree on the loss to three parts in ten million — one integral taken at two wires, the other inside the material.

magnetics · Magnetic loss
The ladder's step response, and the sum of its own stages — 0.95 per cent apart at worst. computed by solving, not by drawing. A step of power into a three-stage thermal ladder, and the junction's rise divided by it. The solid curve is exact: the impedance is a continued fraction in s, its denominator has 3 real negative roots, and the partial-fraction expansion of Z(s)/s is a sum of that many ordinary exponentials — no march, no step size. The dashed curve is the sum every account of a thermal path writes, each stage's own resistance times 1 − exp(−t/RC) with its own local time constant, and it is an approximation because the stages load each other. What that costs is 0.950 per cent, once, at 12.9 ms — between the fastest stage's 2.4 ms and the next one's 200 ms, which is the only place two stages are moving together. It is one-sided: the sum never reads low.

Two ladders the terminals cannot tell apart

A thermal path drawn as a ladder and the same path drawn as a sum of exponentials are called different models of one object, and the difference between them has never been priced because pricing it needs an exact answer. Solved in closed form, the sum is 0.950 per cent high at worst and never low; the marched netlist is right to a part in 21,169; and the largest disagreement in the picture was 2.919 per cent that has nothing to do with heat at all, which reading the curve one sample differently removes.

transients · Reverse-recovery
Two large-signal limits, each alone and then both, at 20 mA and half of it. computed by solving, not by drawing. Four marches of one netlist at each load: neither limit, the input pair's tanh alone, the output stage's 20 mA alone, and both, driven by a 10 mA step. The three curves are each limit's departure from the linear march and the departure with both present; the faint line is the two singles added. Both lies on the sum and a little above it — 1.112 times it at 0.47 nF and 1.022 at 47 nF — so the limits are present together rather than taking turns. Where the two singles cross, near 10 nanofarads, the pair costs 1.88 times what the worse of them costs alone.

The load that neither limit owns

Nine rungs of this argument asked which of an amplifier's two large-signal limits binds, and drew the load capacitance where the answer changes hands. Both are present at every load: the excursion with both in the netlist is the two departures added and between 2 and 12 per cent more, never the larger of them. So the crossing is not a handover but a maximum — at 12 nanofarads the pair costs 2.084 times what the worse of them costs alone, against 1.35 at 2.2 nanofarads and 1.07 at 47 — and the same peak sits on the resistance axis at 20 ohms and the gain-bandwidth axis at 50 megahertz.

feedback · Capacitive load
Stepped at 20 of its time constant, a 1 µs pole rings between 1.818 and 0.331 V, and needs 23 steps to settle. Marched with the trapezoidal rule at a step of 20.0 µs. A 1 µs pole (1 kΩ, 1 nF) drives, through a unity buffer, a 1 ms pole (1 kΩ, 1 µF). The fast node's exact response reaches its final volt within a few microseconds; the march's first values are 1.8182, 0.3306, 1.5477, 0.5519, 1.3666 V. Its distance from its final volt is multiplied by (1 − h/2τ)/(1 + h/2τ) = −0.8182 every step, measured and checked against that form, so it changes sign every step and takes 23 steps to fall below 1% — 460 µs. The slow node it drives is 1.23e-5 V from exact at 1 ms, because a 1 ms pole averages an alternation at half the stepping rate to nothing.

The ringing that belongs to the rule

The trapezoidal rule is stable for every stable circuit and every step size, and it is not damping. March a one-microsecond pole with twenty-microsecond steps and its node reads 1.818, 0.331, 1.548, 0.552 volts — an oscillation at half the stepping rate, its distance from the final volt multiplied by exactly −0.8182 every step, taking twenty-three steps to fall below one per cent. The slow node that pole drives is right to 1.2 × 10⁻⁵ V at a millisecond. One backward-Euler step at the discontinuity cuts the first swing from 0.818 V to 0.048 and two to 0.0023, because backward Euler multiplies the same error by 1/(1 + h/τ) and the trapezoidal rule by (1 − h/2τ)/(1 + h/2τ), which approaches −1.

transients · Step response
At 20 steps a cycle, ten cycles of an undamped LC: the trapezoidal rule keeps the amplitude and falls 29.2° behind; backward Euler keeps 0.0082% of it. Marched, both rules, against 1 − cos ωt for a 1 kHz inductor–capacitor pair stepped with no resistance at all. At 20 steps a cycle the trapezoidal march's amplitude stays at 1.00000 a cycle and its frequency is slow: it loses 2.918° a cycle, measured from the march's own recurrence, against 2π − 2N·atan(π/N) = 2.918°, so after ten cycles it is 29.2° behind. Backward Euler keeps 0.3901 of its amplitude a cycle, against (1 + (2π/N)²)^(−N/2) = 0.3901, so 0.0082% is left after ten, and it loses 11.19° a cycle. No resistance is in the circuit; every loss is the rule's.

The phase the rule loses

An inductor and a capacitor with no resistance ring for ever, and two ways of marching them disagree about how. The trapezoidal rule keeps the amplitude exactly — its factor per step has a magnitude of one — and loses phase instead: 2π − 2N·atan(π/N) a cycle, 2.918° at twenty steps a cycle, so ten cycles later it is 29.2° behind the circuit. Backward Euler keeps 0.3901 of the amplitude a cycle at the same step, and after ten cycles 0.0082 per cent of the ringing is left, in a circuit that has no loss. The two errors fall at different rates: the trapezoidal rule's phase as the square of the steps a cycle, backward Euler's amplitude as the first power. A hundred cycles to within one per cent needs 182 steps a cycle of one and 196,404 of the other.

transients · Step response
The temperature a part cannot come back from, and how long it takes to leave. computed by solving, not by drawing. The same fixed-point equation as the rung below, marched in time with a thermal capacitance rather than solved for its steady states: C dT/dt = P(T) − (T − T_a)/R_th, stepped adaptively on the temperature change because dT/dt goes through zero at each fixed point. Every trajectory starting below 191.1 °C returns to 88.8, however far above the operating point it began; every one starting above it leaves the material's range entirely, the closest in 0.2 minutes. The two nearest starts are 3.0 kelvin apart. The ignition temperature is a boundary in the STARTING CONDITION, and no steady-state analysis contains one.

The boundary that is a starting point

A wound part with a stable operating point at 88.8 degrees and an ignition temperature at 191.1 will never reach the second, because nothing takes it there. Marched in time rather than solved for its steady states, the same equation says what does: a trajectory starting at 189.6 degrees settles back and one starting at 192.6 leaves the material's range in twelve seconds — two starts three kelvin apart. And an overload of four times the normal loss is survivable for ever, while seven times is survivable for seventeen minutes.

magnetics · Thermal feedback

Named alongside it

The objects these essays reach for when they reach for this one.

Model rangeDesign tradeoffLarge-signalSaturationVerificationCoreDielectric absorptionMagnetic lossSettling timeThermal resistanceTrapezoidal ruleAir gap

All concepts