生体信号：曖昧な標本化は後のフィルタでは戻せない
平滑化は不要成分を抑える一方、必要な成分も減衰・遅延させる。異なるアナログ成分が標本化で相殺する場合もある。入力履歴と取得条件を先に定めよう。
学習順路・例題の編集責任：University Reading List編集プロジェクト。 編集と選定の方法 · 出典確認日:
このページの目次
対象と編集責任
架空離散信号と理想アナログ取得の比較。患者波形や臨床検出器ではない。無料閲覧のSmith、専門的RangayyanとKrishnan、実装と画像のSemmlowとGriffelを比較する。
出発点を選ぶ
無料閲覧でDSP基礎
Smithの標本化、畳み込み、移動平均を使う。
三冊を比較する
The Scientist and Engineer's Guide to Digital Signal Processing（科学者と技術者のためのディジタル信号処理ガイド）
Steven W. Smith · 1997
一般DSPの基礎
- 選ぶ場面
- 正弦波、和、添字付き列を扱える。
- 別の教材が適する場面
- 臨床検証手順や無制限再配布を求める。
- 読む範囲
- Sampling、Convolution、Moving Average Filtersの時間・周波数両視点。
- 限界と注意点
- 分かりやすい基礎だが臨床性能研究でもCC作品でもない。
選ぶ前に試す
主例の最初の出力前にx[−1]を求める。
確認すること: 全整数の式から0。ゼロ埋めは与えた履歴ではない。
Biomedical Signal Analysis（生体信号解析）
Rangaraj M. Rangayyan; Sridhar Krishnan · 版: 3 · 2024
生体解析主教材
- 選ぶ場面
- 線形系、確率、生理的信号起源を知る。
- 別の教材が適する場面
- 代数入門や玩具例からの診断が必要。
- 読む範囲
- 信号起源、同時・結合過程、アーチファクト除去から検出へ。
- 限界と注意点
- 適用例は判断を促すが数理応答だけで臨床指標は検証されない。
選ぶ前に試す
h=[1/2,1/2]を(−1)^nへ適用。
確認すること: 隣接項が相殺、H(π)=0。
Biosignal and Medical Image Processing（生体信号と医用画像の処理）
John L. Semmlow; Benjamin Griffel · 版: 3 · 2014
実装と画像
- 選ぶ場面
- 式をコード化し取得条件を調べられる。
- 別の教材が適する場面
- 概念だけ、またはソフトに取得設計を任せたい。
- 読む範囲
- 測定、標本化、離散系、フィルタの後に画像。
- 限界と注意点
- 実装は境界条件を露呈するが後半画像は負荷増。ソフト利用権は別。
選ぶ前に試す
50 Hzと150 Hzを200 Hzで標本化。
確認すること: 同じ位相原点では逆符号の列になり相殺する。
準備を確認する
有限列か表示だけ有限か
全整数で定義され、表示のみn=0,…,3。x[−1]も使う。
二タップ平均の応答は何か
H(ω)=e^(−iω/2)cos(ω/2)。π/2では振幅1/√2、位相−π/4。
アンチエイリアス処理はどこか
アナログ段、標本化前。後処理だけで失われた成分は識別できない。
六つのステップで学ぶ
完了したステップ: 0 / 6
自己記録用のチェックであり、評価ではありません。チェック情報は当サイトに送信されません。
-
履歴を宣言
読む範囲: 列と畳み込み。
取り組む課題: x[−1]からx[3]を書く。
確認すること: 0,3,2,3,0。
-
直接計算
読む範囲: 移動平均。
取り組む課題: 四出力を正確に計算。
確認すること: 1.5,2.5,2.5,1.5。
-
成分別に解く
読む範囲: 周波数応答。
取り組む課題: 直流、交互項、正弦項を分ける。
確認すること: 直流通過、交互項ゼロ、正弦項は変化。
-
閉形式を確認
読む範囲: 線形系。
取り組む課題: 四添字で評価。
確認すること: n=0も直接計算と一致。
-
取得を変える
読む範囲: 標本化と追試。
取り組む課題: 二アナログ正弦波の標本相殺を導く。
確認すること: 未処理列は定数2。
-
予防処理を置く
読む範囲: 取得設計。
取り組む課題: ADC前理想低域処理後に同じデジタル平均。
確認すること: 防止と事後回復は異なる。
解答付きの例題
全整数履歴がある移動平均
全整数nでx[n]=2+cos(πn)+sin(πn/2)。y[n]=(x[n]+x[n−1])/2の式、n=0,…,3、各成分への作用を求めよ。
前提
正確計算、線形時不変二タップ、指定履歴。初期リセット、ゼロ埋め、有限切出し、患者測定なし。
- 直流2は不変。隣接cos(πn)は逆符号で相殺。
- 正弦平均は(1/√2)sin(πn/2−π/4)。
- y=2+(1/√2)sin(πn/2−π/4)。四出力1.5,2.5,2.5,1.5。x[−1]=0で最初は(3+0)/2。
解答と解釈
交互項は消えるが正弦は1/√2倍、位相−π/4、半標本遅延。
結果を検証する
H(0)=1、H(π)=0、H(π/2)=(1−i)/2。不要成分除去は必要成分の完全保存ではない。
前提を変えて考える
取得条件を変えてエイリアス相殺を示す
xa(t)=2+sin(2π50t)+sin(2π150t)を共通位相原点t=n/200で200 Hz標本化する。ADC前処理なしと、|f|<90 Hzを利得1で通す理想矩形低域フィルタを比較し、後で同じ二タップ平均を適用。
- 処理なしでは二項sin(πn/2)とsin(3πn/2)=−sin(πn/2)が全整数で相殺。列も出力も2。
- アナログ低域は直流と50 Hzを通し150 Hzを除く。取得列は2+sin(πn/2)。
- 同じデジタル平均で2+(1/√2)sin(πn/2−π/4)。四出力は1.5,2.5,2.5,1.5だが入力モデルは主例と異なる。
解答と解釈
未処理取得は定数。ADC前理想処理なら50 Hzが残るがデジタル平均で減衰・遅延する。
結果を検証する
定数アナログ入力も同じ未処理標本を生む。標本だけでは二源を区別できず、前処理は相殺を防ぐもので事後回復ではない。
よくある誤り
履歴を省く
初期ゼロ埋めで出力が変わる。
平滑化は形を保存
振幅と位相の両方を確認。
デジタル処理で全エイリアスを修復
別アナログ入力が同じ列を持ちうる。
次に進む前に
- 実際の履歴で畳み込む。
- 成分ごとの振幅と位相を説明。
- 標本化の曖昧さと事前処理を示す。
出典と版の確認
出典確認日:
- Author: DSP guide and editions
著者は1997紙本と後の別形式を分け、無料閲覧の許諾範囲を示す。
- Wiley: third-edition Biomedical Signal Analysis
第3版は起源、結合過程、除去、検出を扱う。
- CRC: exact third-edition hardback
第3版の実装・画像課程は架空信号の臨床妥当性を保証しない。
提案する学習活動であり、公式の課題や全巻通読の主張ではありません。