Перейти к основному содержимому

Реконструкция и идентифицируемость Γ

Область действия и исправление (2026-10-03)

Эта страница задаёт основание обратной задачи для измерительного и экспериментального протоколов. Она заменяет снятый вывод «ригидность пространства состояний обеспечивает единственность нейронного кодировщика». Теоремы ниже относятся к заданной модели наблюдения. Соответствие биологических или ИИ-измерений этой модели остаётся эмпирической гипотезой [Г].

1. Четыре разных отображения​

Пусть D7={Γ=Γ†⪰0:Tr⁡Γ=1}\mathsf D_7=\{\Gamma=\Gamma^\dagger\succeq0:\operatorname{Tr}\Gamma=1\}. Различаем:

  1. Модель наблюдения Oθ:D7→∏a∈AP(Ya)\mathsf O_\theta:\mathsf D_7\to\prod_{a\in\mathsf A}\mathcal P(\mathsf Y_a): вероятностный закон данных для каждой заранее объявленной настройки измерения или вмешательства aa. Параметры калибровки θ\theta и функциональный репер фиксируются до подтверждающего опыта.
  2. Отображение статистик Fθ:D7→RmF_\theta:\mathsf D_7\to\mathbb R^m, например средних измеренных признаков. Равенство статистик обеспечивает равенство законов лишь при соответствующем условии модели, например в гауссовской модели положения с известной фиксированной ковариацией.
  3. Оценщик π^θ:y↦Γ^\hat\pi_\theta:y\mapsto\hat\Gamma либо многозначная реконструкция y↦C(y)⊆D7y\mapsto\mathcal C(y)\subseteq\mathsf D_7 по конечным шумным данным. Линейность и обратимость не предполагаются.
  4. Физический канал E:B(H)→B(C7)\mathcal E:\mathcal B(\mathcal H)\to\mathcal B(\mathbb C^7): если он называется CPTP, он должен быть линейным, вполне положительным и сохраняющим след.

Softmax, нормированное разложение Холецкого и оптимизация правдоподобия дают допустимые состояния; их композиция от этого не становится CPTP-каналом. Правило подготовки отдельно может задать классико-квантовый канал p↦∑xpxπ^θ(x)p\mapsto\sum_xp_x\hat\pi_\theta(x) на указанной классической входной алгебре. Этот канал не доказывает, что оценщик восстанавливает онтологическое состояние системы. См. Watrous, глава 2, определение 2.13.

Если θ\theta неизвестно, идентифицируемость проверяется для совместной модели (Γ,θ)↦Oθ(Γ)(\Gamma,\theta)\mapsto\mathsf O_\theta(\Gamma). Свободная калибровка и свободное скрытое состояние обычно допускают компенсирующие изменения; фиксация подогнанного θ\theta является правилом протокола, а не доказательством единственности скрытой интерпретации.

2. Сильнейшее общее утверждение о единственности​

Определим наблюдательную эквивалентность:

Γ∼θΣ⟺Oθ(Γ)=Oθ(Σ)при всех a∈A.\Gamma\sim_\theta\Sigma\quad\Longleftrightarrow\quad \mathsf O_\theta(\Gamma)=\mathsf O_\theta(\Sigma) \quad\text{при всех }a\in\mathsf A.
Теорема ID-1 (идентификация по слоям наблюдения) [Т]

При фиксированной модели точные популяционные данные определяют ровно класс эквивалентности [Γ]θ[\Gamma]_\theta. Полное состояние идентифицируемо на допустимом множестве S⊆D7\mathsf S\subseteq\mathsf D_7 тогда и только тогда, когда Oθ∣S\mathsf O_\theta|_{\mathsf S} инъективно. Величина qq идентифицируема тогда и только тогда, когда она постоянна на каждом слое наблюдения в S\mathsf S.

Доказательство. Одинаковые законы дают одинаковое распределение результатов всякого разрешённого планом эксперимента, поэтому состояния неразличимы. Разные законы различают состояния на популяционном уровне по определению. Величина спускается на Oθ(S)\mathsf O_\theta(\mathsf S) ровно тогда, когда принимает одинаковое значение на каждом слое. □\square

Это статистическая идентификация; она отличается от точности и успеха на конечной выборке. Классическое различение развито в Rothenberg, Identification in Parametric Models (1971).

Почему ригидность кодировщика не следует​

Прежнее доказательство образовывало π2∘π1−1\pi_2\circ\pi_1^{-1} из непрерывных отображений пространства признаков. Непрерывность не обеспечивает ни обратимости, ни сюръективности. Даже при обратимости π1\pi_1 на образе композиция сначала определена лишь на этом образе и не обязана продолжаться до автоморфизма всего D7\mathsf D_7, сохраняющего структуру.

Элементарный контрпример: X=[−1,1]2\mathsf X=[-1,1]^2, A=diag⁡(1,−1,0,…,0)A=\operatorname{diag}(1,-1,0,\ldots,0) и

π1(x,z)=I/7+x2A/28,π2(x,z)=I/7+x2A/56.\pi_1(x,z)=I/7+x^2A/28,\qquad \pi_2(x,z)=I/7+x^2A/56.

Оба кодировщика непрерывны, положительны, имеют единичный след и необратимы. При x≠0x\ne0 их чистоты различны, поэтому они не связаны унитарным сопряжением. Это опровергает вывод из одной непрерывности и допустимости состояний; всякую более сильную функциональную или феноменологическую посылку нужно записать как явное ограничение и доказать её достаточность.

Ковариантность тоже не определяет кодировщик. На D7×[0,1]\mathsf D_7\times[0,1] отображения (Γ,s)↦Γ(\Gamma,s)\mapsto\Gamma и (Γ,s)↦(1−s)Γ+sI/7(\Gamma,s)\mapsto(1-s)\Gamma+sI/7 оба эквивариантны относительно одновременного унитарного сопряжения. Они в общем случае неэквивалентны и необратимы.

Остаточная симметрия​

Пусть объявленная группа KK действует на состояниях. Единственность с точностью до KK требует отдельного утверждения

Oθ(Γ)=Oθ(Σ)⟺Σ∈K⋅Γ.\mathsf O_\theta(\Gamma)=\mathsf O_\theta(\Sigma) \quad\Longleftrightarrow\quad \Sigma\in K\cdot\Gamma.

Инвариантность наблюдения даёт только включение K⋅Γ⊆[Γ]θK\cdot\Gamma\subseteq[\Gamma]_\theta, а не равенство. Применимая остаточная группа должна сохранять весь план наблюдения, метки, калибровку и заявленные величины. Она может быть меньше G2G_2, меньше Γ ⁣oct\Gamma_{\!\mathrm{oct}} или тривиальной. Ненаблюдаемая степень свободы, меняющая физическую величину, — неидентифицируемость, а не автоматически калибровочная свобода. Выравнивание репера не должно вращать уже фиксированный функциональный репер или скрывать расхождения Φ\Phi и CohE\mathrm{Coh}_E.

3. Полная и устойчивая линейная модель наблюдения​

Пусть H1,…,HmH_1,\ldots,H_m — калиброванные эрмитовы операторы наблюдения и

M(Γ)a=Tr⁡(HaΓ),H0={Δ=Δ†:Tr⁡Δ=0}.M(\Gamma)_a=\operatorname{Tr}(H_a\Gamma),\qquad \mathsf H_0=\{\Delta=\Delta^\dagger:\operatorname{Tr}\Delta=0\}.

Вещественная размерность H0\mathsf H_0 равна 72−1=487^2-1=48.

Теорема ID-2 (информационная полнота и устойчивость) [Т]

Эквивалентны следующие условия:

  1. MM инъективно на всём D7\mathsf D_7.
  2. ker⁡M∩H0={0}\ker M\cap\mathsf H_0=\{0\}.
  3. Проекции Ha−Tr⁡(Ha)I/7H_a-\operatorname{Tr}(H_a)I/7 линейно порождают H0\mathsf H_0.

Если они выполнены, положим α=inf⁡Δ∈H0,∥Δ∥HS=1∥M(Δ)∥2>0\alpha=\inf_{\Delta\in\mathsf H_0,\|\Delta\|_{\mathrm{HS}}=1}\|M(\Delta)\|_2>0. Тогда

∥Γ−Σ∥HS≤α−1∥M(Γ)−M(Σ)∥2.\|\Gamma-\Sigma\|_{\mathrm{HS}} \leq\alpha^{-1}\|M(\Gamma)-M(\Sigma)\|_2.

Доказательство. Разность эрмитовых матриц с единичным следом бесследова. Если ненулевое бесследовое Δ\Delta принадлежит ядру, I/7±εΔI/7\pm\varepsilon\Delta при достаточно малом ε>0\varepsilon>0 — различные матрицы плотности с одинаковыми наблюдениями. Это доказывает (1) iff (2); ортогональность ко всем проекциям доказывает (2) iff (3). Инъективность и компактность единичной сферы H0\mathsf H_0 дают α>0\alpha>0 и оценку. □\square

Конечномерное доказательство приведено здесь полностью. Операторные фреймы и информационно полные квантовые измерения систематически рассматриваются у Scott (2006). Для физических квантовых измерений HaH_a должны быть реализуемыми наблюдаемыми либо соответствующие эффекты должны образовывать указанные POVM; произвольный нейронный признак не становится известным эффектом POVM от присвоения ему метки оси.

Явный координатный план использует Eij=∣i⟩⟨j∣E_{ij}=|i\rangle\langle j| и измеряет

Tr⁡(EiiΓ)=γii,Tr⁡(Eij+Eji2Γ)=Re⁡γij,Tr⁡(i(Eij−Eji)2Γ)=Im⁡γij.\operatorname{Tr}(E_{ii}\Gamma)=\gamma_{ii},\quad \operatorname{Tr}\left(\frac{E_{ij}+E_{ji}}2\Gamma\right)=\operatorname{Re}\gamma_{ij},\quad \operatorname{Tr}\left(\frac{i(E_{ij}-E_{ji})}2\Gamma\right)=\operatorname{Im}\gamma_{ij}.

Семь диагональных величин связаны условием следа; 21+2121+21 вещественных и мнимых частей когерентностей дают в сумме 48 независимых вещественных координат. Для произвольного полного состояния в линейном плане необходимо m≥48m\geq48 независимых статистик. Одна конечная информационно полная POVM требует как минимум 49 исходов, поскольку их вероятности суммируются в единицу. Ни один подсчёт не означает достаточности 48 образцов или глобальной идентификации любыми 48 нелинейными признаками.

4. Фазы нельзя восстановить из модулей​

Если модель зависит лишь от γii\gamma_{ii}, ∣γij∣|\gamma_{ij}| и Gapij=∣sin⁡arg⁡γij∣\mathrm{Gap}_{ij}=|\sin\arg\gamma_{ij}| (нуль при нулевой когерентности), она информационно неполна.

Возьмём равномерную диагональ 1/71/7, a=1/28a=1/28 и ненулевые когерентности только на ориентированном треугольнике (1,2),(2,3),(3,1)(1,2),(2,3),(3,1): все три равны aeiπ/6ae^{i\pi/6}, на обратных рёбрах — эрмитово сопряжённые значения. Строгое диагональное преобладание (1/7>2a1/7>2a) обеспечивает положительную определённость. У этой матрицы и её комплексного сопряжения совпадают диагонали, модули и Gap, но противоположны голономии треугольника:

arg⁡(γ12γ23γ31)=±π/2.\arg(\gamma_{12}\gamma_{23}\gamma_{31})=\pm\pi/2.

Их спектры равны; первый пример показывает потерю ориентированной фазовой информации, а не различие энтропии. Чтобы увидеть потерю спектральной информации, поменяем только ориентированную фазу (3,1)(3,1) на −π/6-\pi/6. Те же диагональ, модули и все попарные Gap сохраняются, а голономия становится π/6\pi/6. Для этого семейства

Tr⁡(Γ3)=149+18a27+6a3cos⁡(θ12+θ23+θ31),\operatorname{Tr}(\Gamma^3)=\frac1{49}+\frac{18a^2}{7}+6a^3\cos(\theta_{12}+\theta_{23}+\theta_{31}),

поэтому у последних состояний спектры различны. Правило arcsin⁡(Gap)\arcsin(\mathrm{Gap}) выбирает ветвь, а не восстанавливает фазу. Величина фазовой синхронизации ∣⟨ei(ϕi−ϕj)⟩∣|\langle e^{i(\phi_i-\phi_j)}\rangle| вещественна и неотрицательна. Знаковая фаза требует комплексного среднего ⟨ei(ϕi−ϕj)⟩\langle e^{i(\phi_i-\phi_j)}\rangle, объявленного эталона и проверенной связи с γij\gamma_{ij}. Отсутствующие знаковые измерения должны оставаться отсутствующими.

PP, RR и каноническая Φ\Phi могут быть идентифицируемы по полной диагональной и модульной информации в фиксированном репере, даже если полное состояние не идентифицируемо. Спектральные и фазочувствительные величины проверяются отдельно; положительность ограничивает допустимые фазы, но не устраняет приведённые контрпримеры.

5. Нелинейная, граничная и динамическая идентификация​

Для C1C^1-отображения статистик FθF_\theta вблизи полнорангового состояния используем локальные координаты q∈R48q\in\mathbb R^{48}. Если DqFθD_qF_\theta имеет столбцовый ранг 48, у некоторого набора из 48 выходных координат якобиан обратим. Теорема об обратной функции даёт локальную идентификацию и локальную оценку устойчивости. Это достаточное, не необходимое условие в особых точках: x↦x3x\mapsto x^3 инъективно при нулевой производной в нуле. Полный ранг не исключает удалённых состояний с одинаковыми данными. Дефект ранга во всей окрестности постоянного ранга даёт локальные слои; один вырожденный якобиан сам по себе не доказывает неидентифицируемость.

В гауссовской модели среднего y∼N(Fθ(q),Σ)y\sim\mathcal N(F_\theta(q),\Sigma) при известной Σ≻0\Sigma\succ0 матрица Фишера равна JTΣ−1JJ^T\Sigma^{-1}J. Невырожденность даёт предыдущий локальный критерий. Общие эквивалентности через информацию Фишера требуют регулярности; их нельзя без изменений переносить на границу неполного ранга или зависящий от состояния носитель распределения. Состояния ранга rr образуют страт вещественной размерности 14r−r2−114r-r^2-1; идентификация при ограничении ранга должна отличать его от соседних состояний других рангов. См. Heinosaari, Mazzarella и Wolf (2013), Quantum Tomography under Prior Information.

Реконструкция состояния, идентификация скоростей динамики и валидация кодировщика — разные задачи. При известной линейной динамике Γ(t)=Tt(Γ0)\Gamma(t)=\mathcal T_t(\Gamma_0) наблюдения в моменты tkt_k определяют начальное состояние тогда и только тогда, когда

Δ⟼(MTtk(Δ))k\Delta\longmapsto\big(M\mathcal T_{t_k}(\Delta)\big)_k

инъективно на H0\mathsf H_0. Неизвестные скорости требуют совместной идентификации с Γ0\Gamma_0. Достаточно богатый временной ряд может дать недостающую информацию, но ни существование ОДУ, ни единственность его решений не доказывают наблюдаемость. Диссипативная обратная задача может быть плохо обусловлена даже при единственности.

6. Реконструкция должна содержать неопределённость​

В линейной гауссовской модели решаем выпуклую задачу

Γ^∈arg⁡min⁡Γ∈D7∥Σ−1/2(M(Γ)−y)∥22.\hat\Gamma\in\arg\min_{\Gamma\in\mathsf D_7} \|\Sigma^{-1/2}(M(\Gamma)-y)\|_2^2.

Непрерывность и компактность обеспечивают существование минимума; ID-2 и Σ≻0\Sigma\succ0 дают строго выпуклую цель на аффинном пространстве единичного следа, а значит единственность минимума. При неполном плане минимумов может быть несколько. Априорное ограничение или правило выбора может оставить один; это не информация, полученная из данных.

Публикуем калиброванное доверительное множество, например

Cε(y)={Γ∈D7:∥Σ−1/2(M(Γ)−y)∥2≤ε},\mathcal C_\varepsilon(y)=\{\Gamma\in\mathsf D_7: \|\Sigma^{-1/2}(M(\Gamma)-y)\|_2\leq\varepsilon\},

с заранее объявленным правилом покрытия для ε\varepsilon, проверенным в модели шума и калибровки. Если множество пусто, сообщаем несовместимость модели наблюдения; не принуждаем состояние попасть в проверяемое окно жизнеспособности. Если ID-2 выполнена и α\alpha вычислена для отбелённого отображения, диаметр множества не превосходит 2ε/α2\varepsilon/\alpha. Также на D7\mathsf D_7 выполнены ∣P(Γ)−P(Σ)∣≤2∥Γ−Σ∥HS|P(\Gamma)-P(\Sigma)|\leq2\|\Gamma-\Sigma\|_{\mathrm{HS}} и ∣R(Γ)−R(Σ)∣≤7∣P(Γ)−P(Σ)∣|R(\Gamma)-R(\Sigma)|\leq7|P(\Gamma)-P(\Sigma)|.

Для каждой величины публикуем диапазон на Cε\mathcal C_\varepsilon. Вердикт определён только при постоянстве предиката на всём множестве; иначе сообщаем не определено данным измерением. Представитель, выбранный оптимизатором, не заменяет эту проверку.

7. Требования подтверждающего протокола​

Защиты SUB-1 … SUB-6 сохраняются: заморозить калибровку на бодрствовании, исключить PCI и поведение из подтверждающего предиктора, положить λ1=λ2=0\lambda_1=\lambda_2=0. Дополняем:

  • ID-A: предрегистрировать закон наблюдения, измерения знаковых фаз, функциональный репер, обработку пропусков, модель шума и независимые настройки; доказать глобальную идентификацию на заявленном допустимом множестве либо явно сообщать частичную идентификацию.
  • ID-B: публиковать ранг и наименьшее сингулярное число калиброванного плана/якобиана, неопределённость калибровки, доверительные множества состояний и тесты восстановления синтетических состояний, включая фазовые контрпримеры.
  • ID-C: предрегистрировать диапазоны величин и неопределённые вердикты; не заменять отсутствующие знаковые фазы нулями или arcsin⁡(Gap)\arcsin(\mathrm{Gap}).
  • ID-D: сравнивать предсказания вне выборки с равным образом заданными альтернативными кодировщиками и моделями сопоставимой сложности. Различать качество замороженного классификатора, адекватность его модели наблюдения и онтологическое утверждение, что его состояние есть Γ\Gamma системы.

Эти требования усиливают эксперимент, делая его неудачу содержательной. Они не задают феноменологический или онтологический мост по определению.