feat(ops): 大気圧・ブレーキ・勾配を調べる3本のスクリプト

[RyoheiHashimoto]

狙い

いずれも「車が既に持っているはずだが、読めていない量」を対象にする。センサーを買い足す前に、まず車が流しているものを確認する。

1. 大気圧 — 一度も問い合わせていなかった

probe-standard-pids.py（30秒・停車で可）

生ログに現れた Mode01 の応答は 05/0B/0F/10/42 の5つだけ。対応PIDビットマップ (0x00/0x20/0x40/0x60) も取っていない。0x33 (大気圧) は「非対応」ではなく未確認。

対応していれば燃料モデルの実バグが直る:

atmosphericKPa = 101.3                  // ハードコード
ve := intakeMAP / atmosphericKPa        // 体積効率の近似

比叡山 (標高848m) の実大気圧は約 91.6kPa。山では体積効率を1割過大に見積もっている。

外気温 (0x46)、エンジン油温 (0x5C)、絶対負荷 (0x43) も同時に確認する。

2. ブレーキ — 対照群を含む設計にした

capture-brake.py（15分・要走行） → find-brake-bit.py（Mac で解析）

操作ラベル付きで candump を録り、踏む/離すに連動するビットを探す。

要点は「ブレーキを踏まずにアクセルオフだけで減速する」区間を必ず含めること。 これが対照群になる。ブレーキ信号ならその区間で立たず、単に減速と相関するだけのビットはそこでも立つ。

手持ちの生ログで 0x200 を調べて決着しなかったのは、走行データが最高17.2km/hの駐車場走行しか無く、対照が取れなかったため（相関は r=+0.20/-0.30/+0.08 と全部弱い）。

3. 勾配 — 17個から2個に絞れた

find-grade-pid.py（Mac で解析、追加操作不要）

FN4A-EL は road load / grade logic で変速線を動かす。登坂で不要なシフトアップを避け、降坂でエンジンブレーキを残すには、PCM が勾配相当の量を持っていなければ制御できない。

未同定 PID 17個 (#150) から「停車中は固定・走行中に動く・車速と無相関」を満たすものを探す。

手持ちの平地データ (1,269,416サンプル) で実際に走らせた結果:

PID	範囲	車速相関	判定
098E	0〜40	+0.243	★停車中は常に0、走行中のみ動く
1746	0〜5376	+0.200	★停車中は固定、走行中に動く

除外されたもの: 16E8 (r=-0.830)、1782 (+0.655)、16CD (+0.651)、17C1 (-0.645) は車速の写し。096E/1101/1103/1104/16E9/3201 は停車中も動くため別の量。

平地データなので「性質が合う」までで同定ではない。 登坂で値が跳ねるかを見る必要がある。スクリプト自身にもその旨を出力させている。

検証

• 3本とも ast.parse 構文チェック通過
• find-grade-pid.py は手持ちの実データ62ファイルで実行し、上記の絞り込み結果を得ている（書きっぱなしにしていない）
• probe-standard-pids.py と capture-brake.py は車両が要るため未実行

使い方

scripts/ops/README.md に順序と所要時間、安全上の注意を書いた。

関連

• GPS が一度も測位できていない（480km・全12ファイルで捕捉衛星0個） #164 GPS が一度も測位できていない
• IMU を追加して勾配・横G・加速度を直接測る #165 IMU を追加して勾配・横G・加速度を直接測る
• ブレーキ信号が CAN にあるか未確認 — 0x200 も未同定のまま #166 ブレーキ信号が CAN にあるか未確認
• chore(obd): 未同定の Mode 22 PID を追跡する（87個中10個が未同定） #150 未同定の Mode 22 PID 17個

🤖 Generated with Claude Code

https://claude.ai/code/session_01JeHX6eyTnRPgubr7XLL8ZB