[发明专利]一种基于短时傅里叶变换的脉内频率编码信号参数测量方法

权利要求书

1.一种基于短时傅里叶变换的脉内频率编码信号参数测量方法，其特征在于，包括：

对待检信号做短时傅里叶变换STFT，根据STFT的结果计算一个子码的长度，并将所计算的子码的长度设置为后续STFT处理的窗长；

对STFT的结果做二维形态学滤波平滑噪声基底；

对平滑噪声基底后的STFT的变换结果做二维过门限检测；所述二维过门限检测为：以恒虚警检测门限作为自适应门限值，以所述待检信号的幅度最大值乘以预先设定的衰减系数作为大信号门限，将自适应门限值与大信号门限进行比较，取二者中较大值作为当前检测门限，对平滑噪声基底后的STFT的变换结果进行过门限检测；

在二维过门限检测结果中检测子码；

根据检测得到的子码的起止位置在所述待检信号的原始数据中截取每个子码对应的数据，并做傅里叶变换FFT计算子码频率，并根据所述子码频率进一步计算所述待检信号的其他参数。

2.如权利要求1所述的测量方法，其特征在于，所述以恒虚警检测门限作为自适应门限值，具体为：

所计算的自适应门限值为

上式中，为获得噪声样本的均值，k为调节因子，δ为偏置补偿因子。

3.如权利要求1所述的测量方法，其特征在于，所述二维过门限检测结果为二维结果图；

所述在二维过门限检测结果中检测子码，具体为：

保留所述二维结果图中每一列的最大值，则每一个子码在二维时频图上的形状转换为一条线，即获得二维线形时频图；

在所述二维线形时频图上，搜索并确定子码的数量及起止位置。

4.如权利要求3所述的测量方法，其特征在于，所述在所述二维线形时频图上，搜索并确定子码的数量及起止位置，采用如下具体步骤：

在所述二维线形时频图上，检测非零点，当前非零点为(i,j)，若对于当前非零点(i,j)，其在时间方向的前方设定范围内的点均为0，则当前非零点(i,j)为当前子码的起始位置，以当前当前非零点(i,j)为起始点，在时间方向的后方设定范围内检测获得的第一个非零点作为当前子码的结束位置，由此获得当前子码的起止位置；

以此类推，检测得到所有子码及其起止位置。

5.如权利要求1～4任一所述的测量方法，其特征在于，所述在二维过门限检测结果中检测子码，之后，还包括噪声子码滤除，具体为：

通过子码的起止位置计算子码长度；

设定噪声子码长度门限值为所述根据STFT的结果计算一个子码的长度的设定比例，将子码长度小于所述噪声子码长度门限值的子码滤除。

6.如权利要求5所述的测量方法，其特征在于，所述噪声子码滤除还包括：

对所有子码的起始位置进行检测，若两子码的起始位置之间的距离小于设定的距离阈值，则滤除两子码中子码长度较短的一个。

7.如权利要求1、2、3、4或6所述的测量方法，其特征在于，所述在二维过门限检测结果中检测子码，之后还包括：

将所有子码的长度值从小到大排列得到子码长度序列，在所述子码长度序列中选取设定点开始依次计算前设定数量个点的方差，若方差发生突变时，则将所述设定点对应的子码长度数值修改为所述前设定数量个点的均值。

8.如权利要求1、2、3、4或6所述的测量方法，其特征在于，所述在二维过门限检测结果中检测子码，之后还包括：

以所有子码长度的平均值作为子码时宽；

对两两相邻子码的起始位置进行比较，若存在两相邻子码起始位置相差子码时宽的N倍，则当前两相邻子码间缺失了N个子码，根据当前两相邻子码位置以及子码时宽计算缺失子码的位置；

对缺失子码进行补全。