# 初回 clock/event 分散 encoder：保存済み科学結果の独立解析

2026-10-19 UTC。新しい合成・fit・科学条件・閾値調整・テスト・実装変更・Git操作なし。旧結果は保持。

## 結論

モデル内C3では4/4 seedが正しい周期AMと96 ms・8 ms・c=4の単一eventを選択した。一方、最大の失敗はC5の非定常meanである。真のq=0・eventなしに対し全4例がq index 10（半径0.6、角度π/4）を選び、背景分散を7.13–7.18倍へ増大させた。eventの偽検出ではなく、mean変化を広域の周期noiseへ誤帰属した。次の単一変更はmeanの時間変化への対応を優先する。C4相互作用は局所的に明瞭な残課題だが、C4だけを理由に予定順を機械的に進めない。

## 対象・再計算の範囲

- 実行条件と全保存数値の公開表は index.html と compact-results.json を参照。元の高精度配列・実行ソースは別途保全し、計測の識別情報を measurement-provenance.json に記録した。
- 受領archiveの545 memberのSHA/bytesと10 source entryを確認済み。本解析は保存配列に対する独立した算術再集計。codec・truth実装をimportせず、推定も復号も再実行していない。
- 6条件×4 input seed＝24 encode。native120/target140の48個の解析分布。candidate96波形/PML96波形/WORLD48波形は実現数であり、解析的KLの反復数ではない。
- 独立再計算したKL、mean/variance/local誤差、積分power、重心、支持、q誤差、event数、C4局所誤差など計2376 scalar照合の最大絶対差は8.88e-16。各target内の2 decode seedで解析配列は完全一致。
- 真のevent支持内は96–104 msの閉区間129/4096標本（3.1494%）。非ゼロ推定支持は127/4096標本（3.1006%）。endpoint 0を含む定義差であり、支持違反ではない。C0/C1/C5のinsideは存在せずnull。

## 条件平均：4 input seedのみを平均

KL単位はnats/sample。各セルはnative120 / target140。

|条件|全区間KL|真event内KL|event外KL|mean由来KL|variance由来KL|
|---|---:|---:|---:|---:|---:|
|C0|0.000848355 / 0.000848918|— / —|0.000848355 / 0.000848918|0.000727869 / 0.000728432|0.000120486 / 0.000120486|
|C1|0.00167486 / 0.00168082|— / —|0.00167486 / 0.00168082|0.00158664 / 0.0015926|8.82144e-05 / 8.82144e-05|
|C2|0.00139774 / 0.00139686|0.021635 / 0.0216353|0.000739663 / 0.00073874|0.00066165 / 0.000660765|0.000736095 / 0.000736095|
|C3|0.00112682 / 0.0011298|0.000674641 / 0.00064228|0.00114152 / 0.00114565|0.00107359 / 0.00107657|5.32305e-05 / 5.32305e-05|
|C4|0.00192408 / 0.00191427|0.031201 / 0.0310446|0.000972042 / 0.000967003|0.000934634 / 0.000931794|0.000989441 / 0.000982477|
|C5|0.930304 / 0.928339|— / —|0.930304 / 0.928339|0.410867 / 0.40938|0.519437 / 0.518959|

分解はKL_mean=0.5(m_true−m_hat)²/v_hat、KL_variance=0.5[log(v_hat/v_true)+v_true/v_hat−1]。加算すると全KLとなる。これは誤差の因果分離ではなく、推定分散を共通分母にした分解である。分散を膨らませればmean項は小さくなり得るため、mean MSEも併記する。

## 24 encodeの選択とseed別KL

強度index 0/1/2はc=1/4/16。群193は96 ms・8 ms、196は97 ms・8 ms。

|条件・input seed|q index|event [群,強度index]|KL120|KL140|次点ΔL|
|---|---:|---|---:|---:|---:|
|C0_seed0|0|[]|0.0005960204|0.0005944138|8.196|
|C0_seed20261015|0|[]|0.0009204877|0.000919588|7.271|
|C0_seed20261019|0|[]|0.0009156317|0.0009159748|6.465|
|C0_seed104729|0|[]|0.0009612796|0.0009656951|8.337|
|C1_seed0|10|[]|0.001646122|0.001666417|9.858|
|C1_seed20261015|10|[]|0.001589286|0.001584837|6.452|
|C1_seed20261019|10|[]|0.002108548|0.002114111|7.402|
|C1_seed104729|10|[]|0.001355465|0.001357898|6.109|
|C2_seed0|0|[[193, 1]]|0.0003765515|0.0003768968|7.926|
|C2_seed20261015|0|[[193, 1]]|0.00107562|0.001071518|2.821|
|C2_seed20261019|0|[[196, 1]]|0.003568069|0.003570937|0.2001|
|C2_seed104729|0|[[193, 1]]|0.0005707364|0.0005680863|6.838|
|C3_seed0|10|[[193, 1]]|0.001572949|0.001597771|5.023|
|C3_seed20261015|10|[[193, 1]]|0.001474332|0.001469382|5.163|
|C3_seed20261019|10|[[193, 1]]|0.0007613254|0.0007563211|5.542|
|C3_seed104729|10|[[193, 1]]|0.0006986744|0.0006957222|0.7803|
|C4_seed0|10|[[193, 0]]|0.003311813|0.003085312|1.32|
|C4_seed20261015|10|[[193, 1]]|0.001209608|0.001282341|7.542|
|C4_seed20261019|10|[[193, 1]]|0.001109497|0.001173675|6.866|
|C4_seed104729|10|[[193, 1]]|0.002065382|0.002115755|7.884|
|C5_seed0|10|[]|0.9304496|0.9285745|6.509|
|C5_seed20261015|10|[]|0.9301917|0.9281312|5.676|
|C5_seed20261019|10|[]|0.9303443|0.9284163|6.959|
|C5_seed104729|10|[]|0.9302298|0.9282357|6.428|

## 局所誤差と帰属

- C0：q=0、eventなしが4/4。C1：真q=10、eventなしが4/4。周期分散をevent枝へ退避する現象はこの小標本では見られない。C2もq=0が4/4で、局在増分の周期化は選択結果に見られない。
- C2：3/4は正しいonset、seed20261019だけ1 ms遅れ。全区間KL約0.003568に対しinside KL約0.085226。次点ΔL約0.2001と近く、単一実現の位置不確実性がある。局在積分powerは全4例−1.70%〜+1.98%でも、位置誤差は残る。遅れ例の支持外power比は約1.09e−5で、真onset前powerは0。これは小さい後尾漏れであり完全支持一致ではない。
- C3：q/位置/幅/強度/個数が4/4一致。局在積分power誤差−1.21%〜+1.71%。全区間KL約0.00113、inside約0.00064–0.00067。最小次点ΔLは0.7803（seed104729）なので、正しい選択を高確信posteriorと読み替えない。
- C4：全4例q=10、onset96 ms・幅8 ms・1 event。ただしseed0はc=1、他3例はc=4。局所在KL平均0.0312は区間全体0.001924より16倍程度大きい。外側KL約0.00097でqベクトル誤差は丸め精度。平均全体KLだけなら局所相互作用不足を見落とす。
- C4の局在power相対誤差平均はnative+0.46%に見えるが、実際の範囲は−68.59%〜+24.39%。正負の相殺である。targetでは平均+6.95%、範囲−66.56%〜+32.43%。推定局在分散形は絶対時間に固定、真値はtarget位相で変形する。重心はnative−0.484 ms、target−0.501 ms。周期qが正しくても局所分布は正しくない。
- C5：qのL2誤差0.6が4/4、phase誤差は真q=0なので採点不能。event数は正しい0でも帰属成功ではない。native KLのmean項0.4109、variance項0.5194、mean MSE0.002483（他条件約0.0000005–0.000001）。nativeと編集の両方で大きい誤差が再現し、単なるtarget時計移動の失敗ではない。native C3に対して全KLは約826倍。
- 全caseに過剰event・欠落eventはない。推定onsetより前および真onsetより前の局在増分は0。背景noiseをevent先行漏れに数えない。2 eventは今回一度も選ばれず、この結果は2 event容量の妥当性を実証しない。

## 通常API参照：F0・レベルのみ

各cellは出力Harvestの正F0中央値の平均Hz / 波形RMSの平均。candidate/PMLは条件・targetごと8波形、WORLDは4波形。seedを均等重みで扱う。これは波形診断であり、KL順位・知覚品質・自然音声優越ではない。本案はoracle位相・Hあり、参照は通常F0推定ありという非対称がある。

|条件|本案120|本案140|PML120|PML140|WORLD120|WORLD140|
|---|---:|---:|---:|---:|---:|---:|
|C0|119.941 / 0.1517|140.024 / 0.1511|119.163 / 0.1851|123.433 / 0.1893|119.912 / 0.1566|139.883 / 0.1507|
|C1|119.930 / 0.1516|140.033 / 0.1515|120.171 / 0.1483|139.966 / 0.1412|119.960 / 0.1567|139.953 / 0.1509|
|C2|119.956 / 0.1518|140.018 / 0.1514|119.526 / 0.1671|129.410 / 0.1725|119.906 / 0.1568|139.894 / 0.1510|
|C3|119.966 / 0.1513|140.023 / 0.1511|120.141 / 0.1487|139.972 / 0.1410|119.939 / 0.1562|139.929 / 0.1504|
|C4|119.951 / 0.1514|140.034 / 0.1512|120.066 / 0.1474|139.875 / 0.1399|119.927 / 0.1564|139.940 / 0.1506|
|C5|119.952 / 0.1612|139.948 / 0.1609|119.946 / 0.1610|139.776 / 0.1522|119.938 / 0.1729|139.905 / 0.1668|

- PMLはC0のseed0/20261015/20261019、C2のseed20261015/20261019の5/24 encodeでraw F0正値なし。fallback後パラメータ中央値は65 Hz、7/6編集後は75.833 Hz。C0/C2の出力が140へ届かない例は、この入口失敗と併記すべき。出力F0が約120でも正しい内部分析を意味しない。
- PMLの通常編集は7/6比であり絶対140保証ではない。fallbackのない条件では出力140付近。WORLDは全24通常分析で正F0あり、全48出力の正F0中央値は目標近傍。candidateもF0目標付近だが、C5の分布誤帰属を救わない。
- 全入力/出力の保存診断で±1超過なし。全出力F0測定は正F0あり。PMLの平均voiced frame率は条件・targetによって約0.875–0.964、candidate/WORLDは1.0。F0測定も推定器依存であり、波形の全品質判定にはしない。
- 共通固定尺度のまま比較。targetのH周波数応答、有限窓mean power、noise実現、参照再合成でRMSは変化し得る。方式別レベル正規化なし。C5では本案/PML native RMSが似ていてもmean時間構造の正しさは判定できない。

## 費用・数値・失敗

|方式|encode RTF 平均 [範囲]|decode RTF平均|保存payload bytes [範囲]|
|---|---:|---:|---:|
|candidate|4.5443 [3.6363, 6.2339]|0.00454|522.7 [503, 542]|
|pml|1.3549 [0.5622, 15.4409]|0.14353|22926.7 [21954, 23476]|
|world|0.2755 [0.2358, 0.4077]|0.01467|403644.3 [403447, 403861]|

- 本案encodeは1区間平均1.163秒、RTF4.54で非リアルタイム。24 encodeのcandidate総比較数2,004,093。有限探索・小payloadを高速encodeと同一視しない。PML初回を含む時間の外れ値があるため、この一回の測定から普遍的な速度比は主張しない。
- coreは実選択で38/41 bytes、source clock込み46/49 bytes。実JSONは503–542 bytesで、理論coreと混同しない。PMLは約22–23 KB、WORLD保存NPZは約404 KBだが、資産・形式・表現力・oracle条件が異なるため等品質圧縮比較ではない。
- 最大条件数12.336、最大abs RSS差9.95e−14、最大abs目的関数直接評価差1.24e−10、σ床近傍フラグ0。case/出力statusはすべてok。数値失敗の証拠は保存診断に見られない。
- RSSメモリはLinux同一processの生涯peakで、方式単独・call単独の割当ではない。最終最大約135,974,912 bytes（129.7 MiB）。方式メモリ優劣に使わない。

## 次の一変更：stationary mean仮定のみを緩める

推奨は、4調波meanに共通する、絶対時間上の小さなgain envelopeを混合から推定する次版。周期noise、非負局在分散、辞書、q格子、罰則を同時変更しない。C5の真gain/128 ms/遷移幅をencoderへ与えず、同じ推定を全条件へ適用する。保存波形や全残差を増やす救済はしない。

具体設計候補は「単一の滑らかなgain遷移」で、mean係数との尺度不定性を固定し、contrast・中心・幅を有限範囲で推定するもの。これならfloat32を3個追加する設計でcore +12 bytes（最大56/clock込み64 bytes）が一つの予算案となる。ただしこれは今回未採用の次版設計案であり、探索域・自由度・記述長罰則・追加演算を事前に決める必要がある。現在の結果から追加RTFや成功率は推定できない。既知C5形状への過適合の危険を明記する。

次版を実行する場合の問いは、C5でmean MSE・全KL・σ膨張・偽qが同時に改善し、C0–C3の帰属を壊さないか。旧版を保持し同じ固定24入力と両targetで比較する。今回ここでは実装もrunも行わない。C4の局所相互作用誤差は別の未解決項として残し、mean変更と同時に直さない。

自然/b/、未知H/位相、48 kHz、window接続、正式120録音/40話者の採用判断はこの実験の範囲外。4 input seedの有限人工条件の結果を母集団成功率や生理的同定へ一般化しない。