学习指南 /

生物医学信号:滤波不能撤销含混采样

平滑可抑制某成分,也可衰减和延迟所需成分;不同模拟分量甚至会在采样后相消。信任曲线之前先明确输入历史和采集模型。

课程专题: 生物医学信号处理

学习计划与例题由 University Reading List 编辑整理。 编辑与选书方法 · 来源确认日期:

本页目录

这份指南适合谁

虚构离散信号与理想模拟采集比较,不是患者波形或临床检测器。Smith为免费可读通用DSP,版权许可有限;Rangayyan与Krishnan加入生医推理,Semmlow与Griffel连接实现及影像。

选择你的起点

需要生医分析判断

信号、概率和基础生理后读Rangayyan与Krishnan。

需要数值实现

具微积分、信号与MATLAB或等效编程后读Semmlow与Griffel。

三本书如何取舍

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采样。

检查标准: 在同相原点t=n/200,两者为相反序列并相消。

书目信息与获取方式 →

开始前需要会什么

序列有限还是仅显示有限?

x[n]对所有整数定义,只显示n=0,…,3,滤波可用x[−1]。

二抽头平均如何作用正弦?

H(ω)=e^(−iω/2)cos(ω/2),π/2处幅值1/√2、相移−π/4。

抗混叠滤波放哪里?

采样之前的模拟段;含混采样后数字滤波不能识别已丢的模拟分量。

六步学习计划

勾选用于自我记录,不是成绩评定;标记不会发送给本站。

  1. 声明历史

    建议阅读范围: 序列和卷积。

    要做的练习: 写全整数公式并列x[−1]至x[3]。

    检查标准: 0、3、2、3、0。

  2. 直接卷积

    建议阅读范围: 移动平均。

    要做的练习: 手算或精确算四输出。

    检查标准: 1.5、2.5、2.5、1.5。

  3. 分解成分

    建议阅读范围: 频率响应。

    要做的练习: 分别求直流、交替与所需正弦。

    检查标准: 直流保留,交替消失,正弦衰减移相。

  4. 检查输出公式

    建议阅读范围: 线性系统。

    要做的练习: 代四索引核闭式。

    检查标准: 含n=0都吻合直接卷积。

  5. 改变采集

    建议阅读范围: 采样与后续题。

    要做的练习: 采指定两模拟正弦并推相消。

    检查标准: 无模拟滤波序列恒为2。

  6. 放预防滤波

    建议阅读范围: 采集推理。

    要做的练习: 先理想模拟低通,再用相同数字平均。

    检查标准: 预防混叠相消不同于事后恢复。

原创例题与解析

有明确无限历史的移动平均

全整数n定义x[n]=2+cos(πn)+sin(πn/2),用y[n]=(x[n]+x[n−1])/2。求闭式、n=0,…,3及各成分效果。

前提假设

精确算术,线性时不变二抽头及指定全整数历史,无启动重置、补零、有限截断或患者测量。

  1. 常数2不变,相邻cos(πn)互反,交替项消失。
  2. 正弦平均为(1/√2)sin(πn/2−π/4)。
  3. 所以y=2+(1/√2)sin(πn/2−π/4),四输出1.5、2.5、2.5、1.5;x[−1]=0,故y[0]=(3+0)/2。

答案与解释

直流保留、交替去除,π/2正弦幅值降至1/√2,相移−π/4,相当于半样本延迟。

检验结果

H(0)=1,H(π)=0,H(π/2)=(1−i)/2;去掉干扰不等于精确保住所需成分。

改变假设,再做一题

改变模拟采集暴露混叠相消

改采xa(t)=2+sin(2π50t)+sin(2π150t),fs=200 Hz,同步同相原点t=n/fs。比较无模拟滤波与ADC前理想单位增益、通过|f|<90 Hz的砖墙低通,之后均用相同二抽头数字滤波。

  1. 无前滤波时两项成为sin(πn/2)与sin(3πn/2)=−sin(πn/2),全整数相消,采样和输出均恒为2。
  2. 指定模拟低通保留直流与50 Hz,采前移除150 Hz,序列2+sin(πn/2)。
  3. 相同数字平均输出2+(1/√2)sin(πn/2−π/4),四样本又为1.5、2.5、2.5、1.5,但输入及采集模型不同于主例。

答案与解释

无前滤波输出常数;理想ADC前滤波保留50 Hz采样分量,数字平均仍衰减延迟它。

检验结果

纯常数模拟输入也对应同一未滤样本,事后算法无法仅凭这些样本区分两来源;模拟滤波预防相消,不是事后恢复丢失信息。

常见误区

不说明启动历史

补零会改变已规定持续信号的初始输出。

平滑就保留波形

必须同时查幅值与相位。

数字滤波修复所有混叠

不同模拟输入可有完全相同采样。

进入下一阶段之前

  • 用真实题设历史卷积。
  • 逐成分说明幅度和相位变化。
  • 展示采样含混及采前滤波作用。

参考来源与版本

来源确认日期:

这些是本站建议的学习活动,并非大学官方作业,也不代表对每本书逐页通读。

编辑与选书方法 →