1. Введение: означает ли «вычисление» обязательный счёт цифр?

Когда сегодня мы слышим слово «компьютер», мы почти всегда представляем себе цифровой компьютер (Digital Computer). В цифровой парадигме любое непрерывное явление физического мира — звуковые волны, вихри жидкости, орбиты планет или работа нейронных сетей — сначала должно быть дискретизировано путём квантования по времени и уровню в конечные последовательности двоичных разрядов (\(0\) и \(1\)). Затем миллиарды логических вентилей переключаются по тактам глобального генератора тактовой частоты, аппроксимируя непрерывную математику миллиардами элементарных арифметических шагов.

Однако с точки зрения физики и истории науки сведение всех вычислений к счёту дискретных цифр (от лат. digitus — «палец») — это относительно недавний обходной путь. Фундаментальные законы природы — уравнения механики Ньютона, электродинамики Максвелла, гидродинамики Навье — Стокса и уравнение Шрёдингера — изначально записаны на языке непрерывных вещественных величин и дифференциальных уравнений.

Задолго до того, как цифровые ЭВМ завоевали господство в середине XX века, учёные и инженеры развивали принципиально иную вычислительную философию — аналоговый компьютер (Analog Computer).

Аналоговый компьютер вообще не считает цифры. У него нет адресов памяти, счётчика команд, дискретных тактов и итерационных циклов численного интегрирования (таких как методы Эйлера или Рунге — Кутты). Вместо этого аналоговый компьютер решает задачу путём создания физической системы — из шестерён, дисковых механических интеграторов, операционных усилителей (ОУ), конденсаторов или оптических интерферометров, — чьи физические уравнения движения математически изоморфны дифференциальным уравнениям исследуемой задачи. В момент замыкания ключа законы Кирхгофа и закон сохранения заряда заставляют непрерывные узловые напряжения \(V(t)\) и токи ветвей \(I(t)\) изменяться в реальном времени строго по точному решению дифференциального уравнения — без ошибок дискретизации и без «бутылочного горлышка фон Неймана» при пересылке данных между памятью и процессором.

В этой статье мы подробно разберём историю и физику аналоговых вычислений: от механических интеграторов Джеймса Томсона, лорда Кельвина и Ванневара Буша до изобретения операционного усилителя, метода обратной связи (Bootstrap) Кельвина, термодинамических пределов теплового шума Джонсона — Найквиста и триумфального возрождения аналоговых вычислений в XXI веке в мемристорных кроссбар-матрицах для ИИ, кремниевой фотонике и нейроморфных чипах.

2. Механическая эпоха: математический анализ с помощью шестерён и дисков

2.1 От Антикитерского механизма до логарифмической линейки

Древнейший известный сложный аналоговый вычислитель — Антикитерский механизм, созданный в эллинистической Греции во II веке до н. э. Содержа более тридцати бронзовых шестерён (включая планетарный дифференциал), он вычислял непрерывные эклиптические координаты Солнца, Луны и планет, а также циклы затмений Сароса. В 1620-х годах на основе логарифмов Джона Непера Уильям Отред изобрёл логарифмическую линейку (Slide Rule): отображая положительное число \(x\) в физический отрезок длины \(\log_{10} x\), она превратила умножение \(x \cdot y\) в непрерывное механическое сложение двух длин \(\log_{10} x + \log_{10} y = \log_{10}(xy)\).

2.2 Джеймс Томсон, лорд Кельвин и диско-роликовый интегратор

Решающий прорыв к механическому решению дифференциальных уравнений произошёл в 1876 году, когда Джеймс Томсон, старший брат Уильяма Томсона (лорда Кельвина), изобрёл диско-роликовый фрикционный интегратор (Wheel-and-Disk Integrator):

1. Горизонтальный диск поворачивается на малый угол \(dx\), пропорциональный независимой переменной \(x\).
2. Маленький вертикальный фрикционный ролик радиуса \(a\) прижат к поверхности диска на радиальном расстоянии \(y(x)\) от центра (представляющем подынтегральную функцию).
3. При повороте диска на \(dx\) дуга пути под точкой контакта равна \(y(x)\,dx\). При качении без проскальзывания угол поворота ролика \(dz\) удовлетворяет соотношению \(a\,dz = y(x)\,dx\), так что суммарный угол поворота ролика \(z(x)\) физически вычисляет определённый интеграл:

Лорд Кельвин использовал этот интегратор в своей знаменитой машине для предсказания приливов (Tide-Predicting Machine), суммировавшей десятки лунных и солнечных гармоник \(\sum_k A_k \cos(\omega_k t + \phi_k)\).

2.3 Дифференциальный анализатор Ванневара Буша и усилитель крутящего момента

Каскадное соединение механических интеграторов долгое время было невозможным из-за малого выходного крутящего момента фрикционного ролика, который начинал проскальзывать под нагрузкой. В 1931 году в MIT Ванневар Буш (Vannevar Bush) решил эту проблему, создав механический усилитель крутящего момента (Torque Amplifier) по принципу кабестана (лебёдки), который увеличивал момент в тысячи раз за счёт внешнего электродвигателя без нагрузки на ролик. В 1941 году Клод Шеннон доказал в работе Mathematical Theory of the Differential Analyzer, что класс функций, воспроизводимых конечным числом интеграторов, в точности совпадает с решениями алгебраических дифференциальных уравнений (ADE).

  flowchart TD
    subgraph Mechanical ["Механические аналоговые вычислители (1876–1940-е: Кельвин и Буш)"]
        M1["Физическая величина: Углы поворота валов θ(t) и смещения y(t)"]
        M2["Базовый оператор: Диско-роликовый интегратор (dz = (1/a) y dx)"]
        M3["Усиление сигнала: Кабестановый механический усилитель крутящего момента"]
        M4["Ограничения: Люфт шестерён, проскальзывание и инерция (минуты на расчёт)"]
        M1 --> M2 --> M3 --> M4
    end

    subgraph Electronic ["Электронные аналоговые компьютеры (1940–1970-е: Bell Labs и ОУ)"]
        E1["Физическая величина: Узловые напряжения V(t) и токи ветвей I(t)"]
        E2["Базовый оператор: Операционный усилитель (ОУ) + Резисторы R + Конденсаторы C"]
        E3["Базовые схемы: Инвертирующий сумматор и интегратор Миллера (V_out = -(1/RC) ∫ V_in dt)"]
        E4["Преимущества: Миллисекундное быстродействие и наборные поля коммутации"]
        E1 --> E2 --> E3 --> E4
    end

    Mechanical ==>|"Изобретение отрицательной обратной связи и электронных ОУ"| Electronic

3. Сердце электронного аналогового компьютера: операционный усилитель (ОУ)

В годы Второй мировой войны Карл Сварцель (Karl D. Swartzel Jr.) в Bell Labs, Гельмут Хёльцер (создатель бортового аналогового автопилота ракеты Фау-2) и Джон Рагаццини (John R. Ragazzini) в Колумбийском университете заменили стальные валы ламповыми схемами. Сам термин «операционный усилитель (Operational Amplifier, Op-Amp)» был введён Рагаццини в 1947 году для обозначения усилителя постоянного тока с огромным коэффициентом усиления, предназначенного для выполнения математических операций (сложения, масштабирования и интегрирования).

3.1 Идеальный ОУ и магия «виртуальной земли (Virtual Ground)»

При коэффициенте усиления разомкнутой петли \(A_{\text{OL}} \sim 10^5 \text{–} 10^8\) выходное напряжение равно \(V_{\text{out}} = A_{\text{OL}}(V_+ - V_-)\). Если заземлить неинвертирующий вход (\(V_+ = 0\,\text{В}\)) и подать отрицательную обратную связь (Negative Feedback) на инвертирующий вход \(V_-\), то при \(A_{\text{OL}} \to \infty\) и бесконечном входном сопротивлении \(Z_{\text{in}} \to \infty\) получаем:

Узел \(V_-\) удерживается под потенциалом \(0\,\text{В}\) без физического соединения с землёй — это виртуальная земля (Virtual Ground), полностью исключающая внутренние нелинейности усилителя из передаточной функции!

3.2 Строгий вывод сумматора и интегратора Миллера

1. Инвертирующий взвешенный сумматор (Summing Amplifier): По первому закону Кирхгофа (KCL) в узле виртуальной земли (\(V_- = 0\)) для \(N\) входных напряжений \(V_k(t)\) через резисторы \(R_k\) и резистор обратной связи \(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. Непрерывный интегратор Миллера (Miller Integrator): Заменяя \(R_f\) конденсатором ёмкости \(C\) с током \(I_C(t) = C \frac{d V_{\text{out}}(t)}{dt}\), по закону Кирхгофа получаем: \(\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. Почему в аналоговых компьютерах запрещены дифференциаторы и используются только интеграторы: В частотной области (\(s = j\omega\)) коэффициент передачи интегратора \(|H_{\text{int}}(\omega)| = \frac{1}{\omega RC}\) спадает с частотой, работая как фильтр нижних частот и подавляя высокочастотный шум. Дифференциатор же имеет передаточную функцию \(|H_{\text{diff}}(\omega)| = \omega RC\) (фильтр верхних частот), который усиливает любую микровольтовую высокочастотную помеху \(\varepsilon \sin(\omega_{\text{high}} t)\) в \(\omega_{\text{high}}\) раз, мгновенно вводя ОУ в насыщение и вызывая самовозбуждение.

4. Решение дифференциальных уравнений коммутацией проводов: метод Кельвина

Чтобы решить уравнение затухающего гармонического осциллятора \(m \ddot{x} + c \dot{x} + k x = F(t)\) только с помощью интеграторов, применяется метод обратной связи (Bootstrap) Кельвина:

1. Выразить старшую производную: \(\ddot{x}(t) = -\frac{c}{m}\dot{x}(t) - \frac{k}{m}x(t) + \frac{1}{m}F(t)\).
2. Предположить, что проводник с напряжением \(\ddot{x}(t)\) уже существует.
3. Каскадно проинтегрировать: Первый инвертирующий интегратор (начальное условие \(-v_0\)) выдаёт \(-\dot{x}(t)\), а второй инвертирующий интегратор (начальное условие \(+x_0\)) выдаёт искомое решение \(+x(t)\).
4. Замкнуть петлю обратной связи: Собрать правую часть из \(-\dot{x}(t)\), \(-x(t)\) и \(F(t)\) с помощью потенциометров и сумматора и соединить выход сумматора напрямую со входом первого интегратора \(\ddot{x}(t)\), предположенным на шаге 2!

  flowchart LR
    F["Внешняя сила +F(t)"] -->|"Потенциометр 1/m"| SUM["Инвертирующий сумматор (Σ)"]
    SUM -->|"Узел Кельвина: d²x/dt²"| INT1["Интегратор 1 (-∫dt)<br/>Нач. условие: -v₀"]
    INT1 -->|"Выход: -dx/dt"| INT2["Интегратор 2 (-∫dt)<br/>Нач. условие: +x₀"]
    INT2 -->|"Решение: +x(t)"| OUT["Осциллограф / Графопостроитель"]
    INT1 -->|"Демпфирование: -c/m"| SUM
    INT2 -->|"Жёсткость: +k/m"| INV["Инвертор (-1)"]
    INV -->|"- (k/m) x(t)"| SUM

5. Физика точности и шума: почему цифра победила аналог

Любой резистор \(R\) при абсолютной температуре \(T > 0\,\text{К}\) создаёт в полосе частот \(\Delta f\) неустранимый тепловой шум Джонсона — Найквиста:

Чтобы передать точность в \(b\) бит (\(2^b = V_{\max}/v_{\text{noise}}\)) по одному аналоговому проводу, отношение сигнал/шум по мощности должно составлять \(\text{SNR}_{\text{power}} = 2^{2b}\). Каждый дополнительный бит точности в аналоговых вычислениях требует четырёхкратного увеличения мощности сигнала или ёмкости конденсатора — физические затраты растут экспоненциально как \(O(2^{2b})\)! Поэтому практический потолок точности электронных АВМ составляет 3–4 десятичных знака (\(10\text{–}14\text{ бит}\), погрешность \(0{,}01\%\text{–}0{,}1\%\)). В цифровом же компьютере каждый транзистор различает лишь два состояния (\(0\) и \(1\)), и аппаратные затраты на \(b\) бит растут лишь линейно или полиномиально как \(O(b \sim b^2)\).

В лунной программе «Аполлон» широко применялись гибридные вычислительные машины (ГВМ), объединявшие мгновенное аналоговое интегрирование динамики полёта в реальном времени с 64-битной цифровой точностью баллистической навигации.

6. Возрождение в XXI веке: вычисления в памяти (AIMC), фотоника и нейроморфные чипы

Сегодня, когда закон Мура упёрся в атомные пределы, а глубокое обучение и большие языковые модели (LLM) столкнулись с колоссальными затратами энергии на пересылку данных между памятью и процессором (узкое место фон Неймана), аналоговые вычисления переживают бурное возрождение:

• Аналоговые вычисления в памяти (AIMC) на мемристорных кроссбарах: Более \(99\%\) вычислений в нейронных сетях составляют матрично-векторные умножения \(y_i = \sum_{j=1}^N W_{ij} x_j\), которым достаточно точности \(4\text{–}8\text{ бит}\). Если записать веса \(W_{ij}\) как проводимости \(G_{ij} = 1/R_{ij}\) мемристоров (ReRAM/PCM), то по закону Ома (\(I_{ij} = G_{ij} V_j\)) и закону Кирхгофа (\(I_i = \sum_j G_{ij} V_j\)) всё матричное умножение выполняется физически за наносекунды за время \(O(1)\) прямо внутри памяти!
• Кремниевая фотоника, нейроморфные чипы и Neural ODE: Когерентные фотонные процессоры перемножают матрицы со скоростью света за счёт интерференции лазеров в решётках интерферометров Маха — Цендера (MZI); подпороговые нейроморфные КМОП-схемы Карвера Мида напрямую моделируют дифференциальные уравнения ионных каналов биологических нейронов, а нейронные ОДУ (Neural ODEs) и жидкостные нейросети (Liquid Neural Networks) возвращают непрерывные дифференциальные уравнения в самый центр современного искусственного интеллекта.