1.基于非延迟代价函数和PauTa检验的瞬时频率估计方法，其特征在于，包括以下步骤：步骤1：采用短时傅里叶变换算法将信号x（k）（k=1, 2, …,N）转换为时频谱图，N代表信号的长度；

步骤2：从时频谱图中选取一块具有较高信噪比的局部区域，采用Canny检测算法将该局部区域转换成二值图像，二值图像包含多条脊带；

步骤3：采用PauTa检验算法排除每条脊带上下边缘的异常值；

步骤4：将上述多条脊带按照相互之间的运动学比例关系叠加到其中一条轮廓最完整的脊带上，构建一条具有完整清晰边缘的合成脊带；

步骤5：采用PauTa检验算法排除上述合成脊带上下边缘的异常值；

步骤6：计算上述合成脊带的均值曲线，采用五点三次平滑算法对均值曲线进行平滑处理，得到平滑均值曲线，计算该平滑均值曲线在95%置信水平上的置信区间；

步骤7：将上述平滑均值曲线及其置信区间按照平滑均值曲线与待估计目标脊线之间的运动学比例关系映射到目标脊线上；

步骤8：将映射后的平滑均值曲线作为目标脊线的参考线，将映射后的置信区间作为目标脊线的局部搜索区间；

步骤9：采用非延迟代价函数在每个时刻所对应的局部搜索区间内搜索脊点，确定每个时刻所对应的瞬时频率，最后得到整个时间区间上的瞬时频率。

2.根据权利要求1所述的基于非延迟代价函数和PauTa检验的瞬时频率估计方法，其特征在于，所述步骤1中短时傅里叶变换算法包括以下步骤：

1）对信号x（k）进行短时傅里叶变换：，

TF(t, f)代表信号x（k）的短时傅里叶变换结果，t代表时间因子，f代表尺度因子，函数w(z)代表自变量为z的窗口函数；

2）计算信号x（k）的时频谱：

，

spectrogram(t, f)代表x（k）的时频谱。

3.根据权利要求1所述的基于非延迟代价函数和PauTa检验的瞬时频率估计方法，其特征在于，所述步骤2中Canny检测算法包括以下步骤：

1) 采用高斯滤波器对图像f(x, y)进行平滑处理，消除图像中的噪声和无关细节：，

，

g(x, y)代表平滑后的图像，G(x, y)代表二维高斯滤波器，x代表图像的时间点，y代表图像的频率点，符号*代表卷积计算，σ代表高斯标准差；

2) 计算g(x, y)强度梯度的幅值和方向，

，

，

M(x, y)代表强度梯度的幅值，θ(x, y)代表强度梯度的方向，gx(x, y)代表g(x, y)对x的偏导数，gy(x, y)代表g(x, y)对y的偏导数；

3) 采用非最大值镇压法消除虚假的边缘：在梯度方向θ(x, y)上，如果非零梯度值大于相邻的两个梯度值，则该非零梯度值保持不变，否则，该非零梯度值置零；

4) 采用大小不同的两个阈值对上述消除虚假边缘后的图像进行滤波，两个阈值记为T1和T2，T1

5) 从I2中剔除不与强边缘连接的弱边缘，然后连接I1和I2中的边缘形成连续边缘。

4.根据权利要求1所述的基于非延迟代价函数和PauTa检验的瞬时频率估计方法，其特征在于：所述步骤3中PauTa检验算法包括以下步骤：

1）估计信号x（n n=1, 2, …,N）的标准差，，

代表样本均值，σ代表样本标准差，N代表样本长度；

2）如果 ，则剔除xn。

5.根据权利要求1所述的基于非延迟代价函数和PauTa检验的瞬时频率估计方法，其特征在于：所述步骤9中非延迟代价函数包括以下步骤：

1）第k个时刻所对应的局部搜索区间FBk 定义为，

fk(pmc)代表映射后的平滑均值曲线在第k个时刻的值， 代表映射后的平滑均值曲线置信区间在第k个时刻宽度的一半，m代表目标脊线的长度；

2）第k个时刻所对应的非延迟代价函数CFk定义为：,

，

fk(i) 代表在FBk范围内所取的频率值，TF(tk, fk)代表TF(t, f)在第k个时刻的值，tk代表t在第k个时刻的值，fk代表f在第k个时刻的值，ek代表权重因子。

6.根据权利要求1所述的基于非延迟代价函数和PauTa检验的瞬时频率估计方法，其特征在于：相对误差≤0.621%，平均相对误差≤0.056%。