技術 / テクニカルガイド

数値微分:信号の傾きとノイズの増幅

微分は、ある量がどれだけ速く変化するかを推定します。また、測定誤差の寄与も変化させます。時間的に近い測定値の差を取り、その値を時間間隔で割るためです。したがって、微分値の変化が鋭く見えても、元の変化をより適切に表しているとは限りません。ノイズ、相関、近似誤差のすべてが関与します。

傾きには、差分と時間スケールを用います。

連続する測定誤差の関係を説明しなければ、一回の測定値の不確かさから微分値の不確かさを推定することはできません。

01 / 差分

離散的な傾きは、ある区間にわたる推定値です。

Ts秒間隔の測定値y[n]について、前進差分はd̂[n] = (y[n+1] − y[n])/Tsです。

Ts秒間隔の測定値y[n]について、前進差分はd̂[n] = (y[n+1] − y[n])/Tsです。yをメートルで測定する場合、結果の単位はm/sです。測定間隔は式の一部であり、任意の表示換算ではありません。

元の信号が滑らかであれば、間隔が小さくなるにつれて、この商は微分値に近づきます。間隔が有限の場合は、その区間の平均勾配を表します。その値を区間の始点、中点、終点のどこに対応付けるかによって、時間上の解釈が変わります。

02 / エラーモデル

減算では、二つの誤差がその共分散とともに組み合わされます。

各読み取り値を、真値と誤差εの和として表してください。

各測定値を、真値と誤差 ε の和として表します。微分値の誤差は (ε[n+1] − ε[n])/Ts です。誤差の分散が等しく σ² で、隣接する誤差の共分散が c の場合、その分散は (2σ² − 2c)/Ts² です。結果には、これらの仮定を併記する必要があります。

隣接する誤差が互いに独立している場合、c はゼロとなり、標準偏差は √2 σ/Ts になります。正の相関は、この独立モデルと比べて差分誤差を減らし、負の相関は増やします。したがって、単一チャネルのノイズ振幅だけでは微分の不確かさは決まりません。

03 / 簡略化した計算例

より短い間隔は、傾きのノイズを増やす可能性があります。

説明のためだけの仮定として、どちらの間隔でも、互いに独立した読み取り誤差の標準偏差が0.1単位であるとします。

あくまで例として、どちらのサンプリング間隔でも、独立した読み取り誤差の標準偏差が0.1単位であると仮定します。Ts = 0.2 sでは、傾き誤差の標準偏差は約0.707単位/sです。Ts = 0.1 sでは、約1.414単位/sになります。

この比較では、一回の読み取り当たりの不確かさを一定としています。実際のセンサーの不確かさがサンプリング、平均化、帯域幅によってどう変化するかを示すものではありません。データ収集モデルを変更すると、σと相関の両方が変わる可能性があるため、この計算は普遍的なサンプリング規則ではありません。

04 / 周波数による重み付け

微分は、急速な変化により大きな重みを与えます。

MathWorksは、差分に基づく微分が、ノイズの多い信号の要素をより際立たせる理由について説明しています。

MathWorksは、差分に基づく微分が信号中のノイズ成分をより強調する理由を説明しています。選択した離散差分の振幅応答は2|sin(Ω/2)|/Tsであり、Ωはサンプルあたりのラジアンで表した角周波数です。周波数がゼロに近い領域では、連続時間の微分に近似します。

有界な離散応答を、連続周波数に対して無制限に拡張したものと混同してはいけません。サンプリング後にすでに区別できなくなった周波数は、微分によって復元できません。微分演算は、情報の欠落や誤った表現も含め、利用可能な表現に対して作用します。

05 / 競合する誤差

ノイズと有限間隔による近似は、別々の項です。

曲線状に変化する信号の傾きは、区間全体で一定ではありません。

曲線状の信号は、区間全体にわたって一定の傾きを持ちません。そのため、測定値にノイズがなくても、有限の間隔によって近似誤差が生じることがあります。元の関数が十分に滑らかであれば、間隔を小さくすることで、この誤差を減らせます。

一方、独立した読み取りノイズは、その小さくなった間隔で割られます。これらの効果は、互いに逆方向に作用することがあります。一定の共通オフセットは差を取ると相殺されますが、ドリフトするオフセットは必ずしも相殺されません。ある種類の誤差を除去しても、微分値が不偏であるとは確認できません。

06 / 対象範囲

微分量は制御器の設定値ではありません。

高階微分では、さらに差分を取る必要があり、時間間隔や滑らかさへの依存性が強まります。

高次微分では差分をさらに重ねるため、タイミングと滑らかさへの依存性が強くなります。そのため、独自の誤差モデルが必要です。また、数値的な傾きは、追加の関係や制約を伴う完全なフィードバック構成内の微分項とは異なります。

MARSTONの微分方法、ノイズレベル、微分制御の挙動は指定されていません。このページでは、実測データの変換、フィルタの調整、センサーの調整は行いません。明示された差分商の数学的な帰結を説明します。

関連する技術ガイド: 推定量のバイアス、分散、平均二乗誤差 · 最小二乗法:残差、重み、および入力誤差 · 振動の変位・速度・加速度.

一般的な背景情報: MathWorks: 信号の微分と計測ノイズ。これらの参考資料は業界の考え方を説明するものであり、適用するMARSTONの製品資料やプロジェクト資料とは別のものです。

主な概念と条件。

各量を、その基準、単位、明示された条件と一緒に読んでください。一般的な関係から、特定モデルの値が確定するわけではありません。

概念 / 量意味 / 関係条件 / 限界
間隔Ts対象とする測定値間の秒数想定した表示レートではなく、実際のモデルのタイミングを使用してください
傾きの単位対象量の単位を秒で割ったもの微分には物理次元の変化が伴う
測定値の分散 σ²入力単位の二乗記載された誤差過程について定義される
近傍共分散 c誤差間で共通する変動独立は、別々のタイムスタンプから推測することはできません。
近似誤差区間の傾きと微分の違い滑らかさと区間の配置に依存する

技術検討上の質問。

サンプリング周波数を高くすれば、常に微分の精度は向上しますか?

いいえ。結果は、間隔、各測定値の不確かさ、および相関に依存します。

固定のオフセットがなぜ消えるのですか?

同じオフセットが、反対の符号で差に現れます。

ドリフトするオフセットも打ち消されますか?

たいていはそうではありません。変化は、推定される傾斜に寄与します。

微分によって、エイリアシングが生じた情報を復元できますか?

いいえ、すでにサンプリングされたデータに存在する情報損失を回復することはできません。

微分の誤差は常に √2 σ/Ts ですか?

記載された、隣接する誤差が等分散かつ互いに独立である場合に限ります。

PIDの微分設定は推奨されていますか?

いいえ。このページで扱うのは有限差分の数学であり、現場での制御変更ではありません。

技術的なお問い合わせ

技術上の疑問について相談する。

微分は、ある量がどれだけ速く変化するかを推定します。技術的な協議のために、モデルと適用される運転情報または設計情報を共有してください。