用SP061A实现心电数据的FFT与压缩

式中m为样本系列的数据个数，x[n]、x分别为原始数据系列及其平均值，y[n]、y分别为重建后的系列及其平均值。CR为压缩比，CR=1（m×10）/(压缩后的字节数×8)

处理时间为SP061A完成FFT与压缩花费的时间，CPU时钟设置为49．152MHz。

表1 实验结果数据

样本编号 代表波形 CR PRD（%） CC 原始数据个数 压缩后字节数 处理时间（ms） 1 图4（a） 3.93 23.132 0.97103 17926 5694 254 3 图4（b） 3.88 14.058 0.98999 19290 6086 263 6 图4（c） 4.08 3.0731 0.99953 65126 20884 857 9 图4（d） 3.92 7.4203 0.99592 12804 3978 179

实验表明，本方法用价格低廉的单片机实现了复杂的FFY与数据压缩，计算耗时少，所得结果满足实用要求。由图4可见，重建后的波形在段间结合点无畸变。噪声较弱时PRD和CC参数较为理想；而当噪声很强时，如图4(a)、4(b)，因滤除了高频噪声而使得重建波形与原始波形差距较大，PRD和CC参数已不能说明问题。压缩算法简便，CR约为4。顺便说明，本方法未实现50Hz干扰滤波、肌电干扰滤波和基线漂移，这些处理可在头端PC上进行。

《用SP061A实现心电数据的FFT与压缩(第3页)》