Backpropagation | From a Dot to a Chatbot

The network after 28 updates

Train it

Start 1 step 2 steps 10 steps 100 steps 1,000 steps 10,000 steps One more step

0.107285total error E

0.4417, 0.8220outputs (want 0.01, 0.99)

1. Forward pass: compute the answer

net_h1 = w1·i1 + w2·i2 + b1 = 0.1564·0.05 + 0.2129·0.10 + 0.35 = 0.379112

h1 = σ(net_h1) = 0.593659 h2 = σ(0.393907) = 0.597223

net_o1 = w5·h1 + w6·h2 + b2 = −0.234299 o1 = 0.441692

net_o2 = w7·h1 + w8·h2 + b2 = 1.530176 o2 = 0.822032

E = ½(t1 − o1)² + ½(t2 − o2)² = 0.093179 + 0.014107 = 0.107285459

What to notice

Output weights are easy. w5's blame is just three numbers multiplied: how wrong o1 is, how steep o1's sigmoid is, and how big h1 was (a big input means the weight mattered more).

Hidden weights reuse that work. w1 can't see the error directly. Its blame arrives through both outputs: output 1's blame times w5, plus output 2's blame times w7. Then the same local slopes. That reuse is what makes backprop fast enough for billions of weights.

The check column agrees. "Nudge and measure" gives the same numbers without any calculus, but it needs two extra forward passes per weight. For an LLM that would be billions of passes per step; backprop needs one.

2. Backward pass: blame for every weight

Weight	Local slopes multiplied together (chain rule)	∂E/∂w	Check (nudge)	New value = w − 0.50 × ∂E/∂w
w5	(o₁ − t₁) 0.43169 × o₁(1 − o₁) 0.24660 × h₁ 0.59366	0.0631981	0.0631981	−0.72224 → −0.75384104
w6	(o₁ − t₁) 0.43169 × o₁(1 − o₁) 0.24660 × h₂ 0.59722	0.0635775	0.0635775	−0.67903 → −0.71082163
w7	(o₂ − t₂) −0.16797 × o₂(1 − o₂) 0.14630 × h₁ 0.59366	−0.0145879	−0.0145879	0.75523 → 0.76252607
w8	(o₂ − t₂) −0.16797 × o₂(1 − o₂) 0.14630 × h₂ 0.59722	−0.0146755	−0.0146755	0.80678 → 0.81411469
w1	δo₁·w5 + δo₂·w7 −0.09544 × h₁(1 − h₁) 0.24123 × i₁ 0.05000	−0.0011512	−0.0011512	0.15645 → 0.15702302
w2	δo₁·w5 + δo₂·w7 −0.09544 × h₁(1 − h₁) 0.24123 × i₂ 0.10000	−0.0023024	−0.0023024	0.21289 → 0.21404604
w3	δo₁·w6 + δo₂·w8 −0.09211 × h₂(1 − h₂) 0.24055 × i₁ 0.05000	−0.0011079	−0.0011079	0.25563 → 0.25618271
w4	δo₁·w6 + δo₂·w8 −0.09211 × h₂(1 − h₂) 0.24055 × i₂ 0.10000	−0.0022157	−0.0022157	0.31126 → 0.31236542

At the start with learning rate 0.5 you should see w5 → 0.35891648 and w1 → 0.14978072, the same as the original worked example. δo means (o − t)·o(1 − o), the output neuron's blame.

The code that just ran

Everything on this page came from calling this code, which lives at the top of this same file.

Show or hide the source

final class TinyNet
{
    public array $w = ['w1' => 0.15, 'w2' => 0.20, 'w3' => 0.25, 'w4' => 0.30,
                       'w5' => 0.40, 'w6' => 0.45, 'w7' => 0.50, 'w8' => 0.55,
                       'b1' => 0.35, 'b2' => 0.60];   // biases stay fixed, as in the original example
    public array $in = [0.05, 0.10];
    public array $target = [0.01, 0.99];

    public function forward(): array
    {
        $w = $this->w; [$i1, $i2] = $this->in; [$t1, $t2] = $this->target;
        $nh1 = $w['w1'] * $i1 + $w['w2'] * $i2 + $w['b1'];  $h1 = sigmoid($nh1);
        $nh2 = $w['w3'] * $i1 + $w['w4'] * $i2 + $w['b1'];  $h2 = sigmoid($nh2);
        $no1 = $w['w5'] * $h1 + $w['w6'] * $h2 + $w['b2'];  $o1 = sigmoid($no1);
        $no2 = $w['w7'] * $h1 + $w['w8'] * $h2 + $w['b2'];  $o2 = sigmoid($no2);
        $E1 = 0.5 * ($t1 - $o1) ** 2;
        $E2 = 0.5 * ($t2 - $o2) ** 2;
        return compact('i1', 'i2', 't1', 't2', 'nh1', 'h1', 'nh2', 'h2', 'no1', 'o1', 'no2', 'o2', 'E1', 'E2') + ['E' => $E1 + $E2];
    }

    public function loss(): float { return $this->forward()['E']; }

    // BACKPROPAGATION. Each gradient is a product of simple local slopes:
    //   (o − t)      how the error changes with the output
    //   o·(1 − o)    how the sigmoid output changes with its input
    //   h or i       how that input changes with the weight
    // Hidden weights get the SUM of the blame flowing back from both outputs.
    public function gradients(): array
    {
        $f = $this->forward(); $w = $this->w;
        $do1 = ($f['o1'] - $f['t1']) * $f['o1'] * (1 - $f['o1']);   // δ of output 1
        $do2 = ($f['o2'] - $f['t2']) * $f['o2'] * (1 - $f['o2']);   // δ of output 2
        $back1 = $do1 * $w['w5'] + $do2 * $w['w7'];                 // blame reaching h1
        $back2 = $do1 * $w['w6'] + $do2 * $w['w8'];                 // blame reaching h2
        $s1 = $f['h1'] * (1 - $f['h1']);
        $s2 = $f['h2'] * (1 - $f['h2']);

        $g = fn($factors) => ['factors' => $factors, 'value' => array_product(array_column($factors, 1))];
        return [
            'w5' => $g([['(o₁ − t₁)', $f['o1'] - $f['t1']], ['o₁(1 − o₁)', $f['o1'] * (1 - $f['o1'])], ['h₁', $f['h1']]]),
            'w6' => $g([['(o₁ − t₁)', $f['o1'] - $f['t1']], ['o₁(1 − o₁)', $f['o1'] * (1 - $f['o1'])], ['h₂', $f['h2']]]),
            'w7' => $g([['(o₂ − t₂)', $f['o2'] - $f['t2']], ['o₂(1 − o₂)', $f['o2'] * (1 - $f['o2'])], ['h₁', $f['h1']]]),
            'w8' => $g([['(o₂ − t₂)', $f['o2'] - $f['t2']], ['o₂(1 − o₂)', $f['o2'] * (1 - $f['o2'])], ['h₂', $f['h2']]]),
            'w1' => $g([['δo₁·w5 + δo₂·w7', $back1], ['h₁(1 − h₁)', $s1], ['i₁', $f['i1']]]),
            'w2' => $g([['δo₁·w5 + δo₂·w7', $back1], ['h₁(1 − h₁)', $s1], ['i₂', $f['i2']]]),
            'w3' => $g([['δo₁·w6 + δo₂·w8', $back2], ['h₂(1 − h₂)', $s2], ['i₁', $f['i1']]]),
            'w4' => $g([['δo₁·w6 + δo₂·w8', $back2], ['h₂(1 − h₂)', $s2], ['i₂', $f['i2']]]),
        ];
    }

    // Sanity check without any calculus: nudge the weight up and down a hair
    // and see how much the loss moves. Should match backprop almost exactly.
    public function numericalGradient(string $name, float $eps = 1e-5): float
    {
        $keep = $this->w[$name];
        $this->w[$name] = $keep + $eps; $up = $this->loss();
        $this->w[$name] = $keep - $eps; $down = $this->loss();
        $this->w[$name] = $keep;
        return ($up - $down) / (2 * $eps);
    }

    public function step(float $lr): void
    {
        foreach ($this->gradients() as $name => $g) $this->w[$name] -= $lr * $g['value'];
    }
}