[发明专利]基于特征谱的胃磁慢波信号频率检测方法

权利要求书

1.一种基于特征谱的胃磁慢波信号频率检测方法，其特征在于，包括以下步骤：

数据预处理，读取巨磁阻抗传感器在特表测量到的原始胃磁数据，并对原始胃磁数据进行重新采样和滤波处理；

特征谱计算，计算所述滤波处理后的数据的特征谱；

特征谱主峰识别，对上述特征谱进行主峰识别。

2.根据权利要求1所述的基于特征谱的胃磁慢波信号频率检测方法，其特征在于，所述数据预处理包括以下步骤：

重采样，将原始胃磁数据重新采样为1Hz进行数据采样；

数字滤波，对重采样后的数据进行0.01Hz-0.5Hz数字带通滤波处理。

3.根据权利要求1所述的基于特征谱的胃磁慢波信号频率检测方法，其特征在于，所述特征谱计算包括以下步骤：

自相关序列估计，对所述滤波处理后的数据进行数据加窗处理，利用快速傅里叶变换的方法计算出功率谱密度函数的估值，然后再计算它的傅里叶反变换，即得自相关估计函数；

自相关矩阵估计，利用上面得到的自相关估计函数，定义一个N×N的方阵，其主对角线为自相关函数的常数项，主对角线旁边的第一对角线为自相关函数的一次项，然后为自相关函数的二次项，最右上角和最左下角为自相关函数的N次项，从而得到胃磁数据的自相关矩阵；

特征向量计算，计算所述自相关矩阵的特征值，将计算处理的特征值按降序排列得到一组特征向量；

谱计算，将所述特征向量通过预设的参数计算出特征谱。

4.根据权利要求3所述的基于特征谱的胃磁慢波信号频率检测方法，其特征在于，所述谱计算采用采用MUSIC方法：

Pmusic(w)=1Σk=p+1M|eHvk|2]]>

e=[1ejwek2w...ej(M-1)w]T]]>

其中，e为自然对数的底数，j为虚数单位，w为角频率，M为特征向量的维数，e^H为e的共轭转置，p为信号子空间的维数，v_k为相关矩阵的第k个特征向量，这些特征向量构成噪声子空间。

5.根据权利要求1所述的基于特征谱的胃磁慢波信号频率检测方法，其特征在于，所述特征谱主峰识别具体为：获取特征谱的谱函数，基于函数的最大值计算方法，计算谱函数最大值所对应的频率位置。

6.根据权利要求1所述的基于特征谱的胃磁慢波信号频率检测方法，其特征在于，如果特征谱主峰识别不成功，则修改特征谱计算过程中的参数，重新生成特征谱。

7.根据权利要求1所述的基于特征谱的胃磁慢波信号频率检测方法，其特征在于，如果特征谱主峰识别成功，则将所述特征谱输出或存储。