基于SIMULINK的心电信号源系统设计分析

　　心电信号主要频率范围为0.05"100Hz，幅度约为0"4mV，信号十分微弱。由于心电信号中通常混杂有其它生物信号，加之体外以500Hz工频为主的电磁场干扰，使得心电噪声背景较强，测量条件比较复杂，为医学研究、救治和教学带来一定困难。本文通过对心电信号的特征分析，提出了基于MATLAB/SIMULINK的心电信号仿真方法，采用快速原型技术通过输入/输出卡，将虚拟仿真信号转化生成实际的物理电信号，并与实际硬件电路连接起来，构成心电信号的半实物仿真模型。

　　基本原理

　　心电信号

　　一个典型的心电信号如图1示，它由P波、Q波、R波、S波和T波组成，在诊断过程中可以根据这些特征波的间隔幅值以及P-R间隔、S-T间隔和Q-T间隔判断病人病情。

　　Fourier级数

　　若一个函数满足狄里赫利条件，则可以将它展开为傅立叶级数的形式。傅立叶级数的三角函数形式定义如下：

　　狄里赫利条件充分条件：

　　(1)函数在任意有限区间连续，或只有有限个第一类间断点(当t从左或右趋向于这个间断点时，函数存在有限的左极限和右极限);

　　(2)在一周期内，函数存在有限个极大值或极小值。

　　心电信号仿真及结果分析

　　心电信号仿真

　　本设计方案是将心电信号看作一个由各个三角波信号和正弦信号的组合，先经过计算得到各个特征波序列，再将各波合成为最终的仿真信号。具体实现方法如下：

　　(1)心电信号的QRS波、Q波和S波用三角波实现。

　　(2)心电信号的P波、T波和U波用正弦波实现。

　　三角波数学模型建立

　　三角波数学模型建立如下：

　　为方便计算，将此函数展开为傅立叶级数的形式：

　　式中，

　　正弦波数学模型建立

　　正弦波建立模型如下：

　　为方便计算，将此函数展开为傅立叶级数的形式：

　　仿真结果分析

　　通过上述仿真方案，仿真心电信号各特征波参数如表1示，心电信号仿真波形如图2示。

　　实时心电信号源设计

　　设计方案采用SIMULINK实现心电信号发生器的设计。SIMULINK可以方便地进行可视化动态建模，且仿真过程是交互的，可以随时修改参数，能够立即看到仿真结果，生成可以实际应用的实时信号。

　　心电信号发生器实现

　　设计方案中将心电信号看作是由各个特征波组合在一起的波形，这些特征波包括P波、Q波、QRS波、S波、T波和U波。仿真设计过程中，Q波、QRS波和S波用三角波实现，P波、T波和U波用正弦波实现。