1. Introduction : « Calculer » signifie-t-il nécessairement compter des chiffres ?

Lorsque nous entendons le mot « ordinateur » aujourd’hui, nous imaginons presque exclusivement un calculateur numérique (Digital Computer). Dans le paradigme numérique, chaque phénomène continu de l’univers physique — ondes sonores, turbulences fluides, orbites planétaires ou inférence neuronale — doit d’abord être discrétisé par échantillonnage et quantification en suites finies de chiffres binaires ($0$ et $1$). Des milliards de portes logiques commutent ensuite au rythme d’une horloge globale pour approximer les mathématiques continues par des milliards d’opérations arithmétiques élémentaires.

Pourtant, à l’échelle de l’histoire des sciences et de la physique, réduire tout calcul au dénombrement de chiffres discrets (du latin digitus, « doigt ») est un détour relativement récent. Les lois fondamentales de la nature — mécanique newtonienne, électrodynamique de Maxwell, équations de Navier-Stokes, équation de Schrödinger — s’expriment nativement dans le langage des variables réelles continues gouvernées par des équations différentielles.

Bien avant l’hégémonie du numérique au milieu du XXe siècle, les scientifiques et ingénieurs ont développé une philosophie computationnelle radicalement différente : le calculateur analogique (Analog Computer).

Un calculateur analogique ne compte aucun chiffre. Il ne possède ni adresses mémoire, ni compteur ordinal, ni cycles d’horloge discrets, ni boucles d’intégration numérique itératives (telles que les méthodes d’Euler ou de Runge-Kutta). Au lieu de cela, il résout un problème en construisant un système physique — composé d’engrenages, d’intégrateurs mécaniques à disque, d’amplificateurs opérationnels (AOP), de condensateurs ou d’interféromètres optiques — dont les équations physiques gouvernantes sont mathématiquement isomorphes aux équations différentielles du problème étudié. Dès la mise sous tension, les lois de Kirchhoff et la conservation de la charge forcent les tensions continues $V(t)$ et les courants $I(t)$ du circuit à évoluer en temps réel selon la solution exacte de l’équation différentielle, sans erreur de troncature liée à la discrétisation et sans le goulot d’étranglement de von Neumann lié au transfert de données entre mémoire et processeur.

Cet article retrace l’épopée complète du calcul analogique : des intégrateurs mécaniques de James Thomson, Lord Kelvin et Vannevar Bush à l’invention de l’amplificateur opérationnel et de la méthode de bouclage (Bootstrap) de Kelvin, jusqu’aux limites thermodynamiques et théoriques (bruit thermique de Johnson-Nyquist) et à la spectaculaire renaissance du calcul analogique au XXIe siècle à travers le calcul en mémoire (memristors), la photonique sur silicium et l’ingénierie neuromorphique.

2. La lignée mécanique : exécuter le calcul infinitésimal avec des disques et des engrenages

2.1 De la machine d’Anticythère à la règle à calcul

Le plus ancien calculateur analogique complexe connu est la machine d’Anticythère, conçue en Grèce antique vers le IIe siècle av. J.-C. Grâce à plus de trente engrenages de bronze (dont un train épicycloïdal différentiel), elle calculait les positions écliptiques continues du Soleil, de la Lune et des planètes ainsi que les cycles d’éclipses de Saros. Vers 1622, à la suite de l’invention des logarithmes par John Napier, William Oughtred inventa la règle à calcul (Slide Rule) : en associant à tout réel positif $x$ une longueur physique proportionnelle à $\log_{10} x$, elle transforme la multiplication $x \cdot y$ en l’addition mécanique continue de deux déplacements linéaires $\log_{10} x + \log_{10} y = \log_{10}(xy)$.

2.2 James Thomson, Lord Kelvin et l’intégrateur à disque et galet

Le saut décisif permettant de résoudre mécaniquement des équations différentielles eut lieu en 1876 lorsque James Thomson, frère aîné de William Thomson (Lord Kelvin), inventa l’intégrateur à disque et galet (Wheel-and-Disk Integrator) :

1. Un disque horizontal tourne d’un angle infinitésimal $dx$ proportionnel à la variable indépendante $x$.
2. Un petit galet de friction vertical de rayon $a$ repose sur la surface du disque à une distance radiale $y(x)$ du centre (représentant la fonction à intégrer).
3. Lorsque le disque tourne de $dx$, l’arc parcouru sous le point de contact vaut $y(x)\,dx$. En supposant un roulement sans glissement, la rotation $dz$ du galet satisfait $a\,dz = y(x)\,dx$, de sorte que l’angle total $z(x)$ du galet évalue physiquement l’intégrale définie :

Lord Kelvin utilisa cet intégrateur pour construire son célèbre prédicteur de marées (Tide-Predicting Machine), sommant les composantes harmoniques lunaires et solaires $\sum_k A_k \cos(\omega_k t + \phi_k)$ pour tracer en quatre heures les courbes de marées continues d’une année entière.

2.3 L’analyseur différentiel de Vannevar Bush et l’amplificateur de couple

Le couplage de plusieurs intégrateurs mécaniques resta longtemps bloqué par la faiblesse du couple de sortie du galet de friction, qui glissait dès qu’on lui imposait une charge mécanique. En 1931 au MIT, Vannevar Bush résolut ce problème grâce à l’amplificateur de couple mécanique (Torque Amplifier) fondé sur le principe du cabestan, amplifiant le couple des milliers de fois via un moteur externe sans charger le galet. En 1941, Claude Shannon démontra dans sa Théorie mathématique de l’analyseur différentiel que les fonctions générables par un nombre fini d’intégrateurs coïncident exactement avec les solutions des équations différentielles algébriques (ADE).

  flowchart TD
    subgraph Mechanical ["Calculateurs analogiques mécaniques (1876–1940 : Kelvin & Bush)"]
        M1["Grandeur physique : Angles de rotation d'arbre θ(t) & déplacements y(t)"]
        M2["Opérateur clé : Intégrateur à disque et galet (dz = (1/a) y dx)"]
        M3["Amplification : Amplificateur de couple mécanique à cabestan"]
        M4["Limites : Jeu d'engrenage, glissement, inertie mécanique (minutes par calcul)"]
        M1 --> M2 --> M3 --> M4
    end

    subgraph Electronic ["Calculateurs analogiques électroniques (1940–1970 : Bell Labs & AOP)"]
        E1["Grandeur physique : Tensions de nœud V(t) & courants de branche I(t)"]
        E2["Opérateur clé : Amplificateur opérationnel (AOP) + Résistances R + Condensateurs C"]
        E3["Circuits clés : Sommateur inverseur & Intégrateur de Miller (V_out = -(1/RC) ∫ V_in dt)"]
        E4["Avantages : Calcul en millisecondes/microsecondes & câblage sur panneau de brassage"]
        E1 --> E2 --> E3 --> E4
    end

    Mechanical ==>|"Invention de la contre-réaction & des AOP électroniques"| Electronic

3. Le cœur du calcul analogique électronique : l’amplificateur opérationnel (AOP)

Pendant la Seconde Guerre mondiale, Karl D. Swartzel Jr. (Bell Labs), Helmut Hoelzer (concepteur du calculateur analogique de guidage de la fusée V2) et John R. Ragazzini (Université Columbia) remplacèrent les arbres métalliques par des circuits électroniques à tubes. Le terme même d’« amplificateur opérationnel (Operational Amplifier, Op-Amp) » fut forgé par Ragazzini en 1947 pour désigner un amplificateur à courant continu à grand gain conçu pour exécuter des opérations mathématiques (addition, intégration, mise à l’échelle).

3.1 L’AOP idéal et la magie de la « masse virtuelle (Virtual Ground) »

Avec un gain en boucle ouverte $A_{\text{OL}} \sim 10^5 \text{–} 10^8$, la tension de sortie vaut $V_{\text{out}} = A_{\text{OL}}(V_+ - V_-)$. En reliant l’entrée non inverseuse à la masse ($V_+ = 0\,\text{V}$) et en appliquant une contre-réaction (Negative Feedback) vers l’entrée inverseuse $V_-$, la limite $A_{\text{OL}} \to \infty$ (et l’impédance d’entrée infinie $Z_{\text{in}} \to \infty$) impose :

Le nœud $V_-$ est maintenu à $0\,\text{V}$ sans être connecté physiquement à la terre : c’est la masse virtuelle (Virtual Ground), qui rend la fonction de transfert indépendante des fluctuations internes de l’amplificateur.

3.2 Dérivation exacte du sommateur et de l’intégrateur de Miller

1. Sommateur inverseur pondéré (Summing Amplifier) : En appliquant la loi des nœuds de Kirchhoff (KCL) à la masse virtuelle ($V_- = 0$) pour $N$ tensions d’entrée $V_k(t)$ à travers les résistances $R_k$ et une résistance de contre-réaction $R_f$ : $$\sum_{k=1}^N \frac{V_k(t)}{R_k} + \frac{V_{\text{out}}(t)}{R_f} = 0 \implies V_{\text{out}}(t) = -\sum_{k=1}^N \left(\frac{R_f}{R_k}\right) V_k(t)$$
2. Intégrateur continu de Miller (Miller Integrator) : En remplaçant $R_f$ par un condensateur de capacité $C$, le courant de charge est $I_C(t) = C \frac{d V_{\text{out}}(t)}{dt}$, ce qui donne par KCL : $$\frac{V_{\text{in}}(t)}{R} + C \frac{d V_{\text{out}}(t)}{dt} = 0 \implies V_{\text{out}}(t) = -\frac{1}{RC} \int_0^t V_{\text{in}}(\tau)\,d\tau + V_{\text{out}}(0)$$
3. Pourquoi les calculateurs analogiques proscrivent les dérivateurs au profit exclusif des intégrateurs : En fréquence ($s = j\omega$), l’intégrateur possède un gain $|H_{\text{int}}(\omega)| = \frac{1}{\omega RC}$ qui agit comme un filtre passe-bas lissant le bruit thermique haute fréquence. À l’inverse, un dérivateur a un gain $|H_{\text{diff}}(\omega)| = \omega RC$ (filtre passe-haut) qui amplifie tout bruit haute fréquence $\varepsilon \sin(\omega_{\text{high}} t)$ d’un facteur $\omega_{\text{high}}$, saturant immédiatement l’AOP.

4. Résoudre les équations différentielles par câblage : la méthode Bootstrap de Kelvin

Pour résoudre l’équation de l’oscillateur harmonique amorti $m \ddot{x} + c \dot{x} + k x = F(t)$ uniquement avec des intégrateurs, la méthode de bouclage (Bootstrap) de Kelvin procède en quatre étapes :

1. Isoler la dérivée d’ordre le plus élevé : $\ddot{x}(t) = -\frac{c}{m}\dot{x}(t) - \frac{k}{m}x(t) + \frac{1}{m}F(t)$.
2. Supposer que le fil portant la tension $\ddot{x}(t)$ est déjà disponible.
3. Intégrer en cascade : Un premier intégrateur inverseur (condition initiale $-v_0$) produit $-\dot{x}(t)$, et un second intégrateur inverseur (condition initiale $+x_0$) produit la position $+x(t)$.
4. Fermer la boucle de rétroaction : Combiner $-\dot{x}(t)$, $-x(t)$ et $F(t)$ via des potentiomètres et un sommateur pour former le membre de droite, puis relier la sortie de ce sommateur directement au nœud d’entrée $\ddot{x}(t)$ supposé à l’étape 2 !

  flowchart LR
    F["Tension de force +F(t)"] -->|"Potentiomètre 1/m"| SUM["Sommateur inverseur (Σ)"]
    SUM -->|"Nœud Bootstrap : d²x/dt²"| INT1["Intégrateur 1 (-∫dt)<br/>Cond. initiale : -v₀"]
    INT1 -->|"Sortie : -dx/dt"| INT2["Intégrateur 2 (-∫dt)<br/>Cond. initiale : +x₀"]
    INT2 -->|"Solution : +x(t)"| OUT["Oscilloscope / Traceur"]
    INT1 -->|"Amortissement : -c/m"| SUM
    INT2 -->|"Raideur : +k/m"| INV["Inverseur (-1)"]
    INV -->|"- (k/m) x(t)"| SUM

5. Physique de la précision et du bruit : pourquoi le numérique a détrôné l’analogique

La limite fondamentale du calcul analogique découle de la thermodynamique statistique et de la théorie de l’information de Shannon. Toute résistance $R$ à température absolue $T > 0\,\text{K}$ génère un bruit thermique de Johnson-Nyquist sur une bande passante $\Delta f$ :

Pour représenter une précision de $b$ bits ($2^b = V_{\max}/v_{\text{noise}}$) sur un fil analogique continu, le rapport signal sur bruit en puissance doit atteindre $\text{SNR}_{\text{power}} = 2^{2b}$. Chaque bit de précision supplémentaire multiplie par quatre la puissance énergétique ou la capacité requise : le coût physique du calcul analogique croît exponentiellement en $O(2^{2b})$ ! En pratique, la précision analogique plafonne à 3 ou 4 chiffres décimaux ($10\text{–}14\text{ bits}$, erreur de $0{,}01\% \text{ à } 0{,}1\%$). Dans un calculateur numérique, chaque transistor ne distingue que deux états ($0$ et $1$) avec une large marge de bruit, de sorte que le coût matériel pour $b$ bits ne croît que linéairement ou polynomialement en $O(b \text{ à } b^2)$.

Lors du programme spatial Apollo, les calculateurs hybrides (Hybrid Computers) associaient la vitesse temps réel de l’analogique (simulation du vol à 6 degrés de liberté avec matériel dans la boucle) à la précision de 64 bits du numérique pour la navigation translunaire.

6. La renaissance au XXIe siècle : calcul en mémoire (AIMC), photonique et puces neuromorphiques

Aujourd’hui, face à la fin de la loi de Moore et au mur énergétique du transfert de données mémoire-processeur (goulot d’étranglement de von Neumann) dans l’apprentissage profond et les grands modèles de langage (LLM), le calcul analogique connaît une renaissance spectaculaire :

• Calcul analogique en mémoire (AIMC) à matrices de memristors : Plus de $99\%$ des calculs d’IA consistent en multiplications matrice-vecteur $y_i = \sum_{j=1}^N W_{ij} x_j$ qui ne requièrent qu’une précision de $4\text{ à }8\text{ bits}$. En stockant les poids $W_{ij}$ sous forme de conductances $G_{ij} = 1/R_{ij}$ dans une grille de memristors (ReRAM/PCM), la loi d’Ohm ($I_{ij} = G_{ij} V_j$) et la loi des nœuds de Kirchhoff ($I_i = \sum_j G_{ij} V_j$) exécutent la multiplication matricielle entière en quelques nanosecondes en temps $O(1)$ directement dans la mémoire !
• Photonique sur silicium, puces neuromorphiques et Neural ODEs : Les processeurs photoniques utilisent l’interférence laser dans des réseaux d’interféromètres de Mach-Zehnder (MZI) pour calculer à la vitesse de la lumière ; les circuits neuromorphiques sous le seuil (Subthreshold CMOS) de Carver Mead exploitent la physique exponentielle des transistors pour émuler directement les équations différentielles des neurones biologiques, tandis que les Neural ODEs et réseaux liquides réconcilient équations différentielles continues et intelligence artificielle moderne.