CN106093885A - 线性调频信号调频率的估计方法及装置 - Google Patents

Abstract

本发明公开了一种线性调频信号调频率的估计方法及装置，其中，方法包括以下步骤：获取待分析信号的调频率最大值；获取调频率误差，以得到分数阶傅里叶变换的转角及其个数；将信号均分，以统计信号对应的第一过零数与第二过零数；判断第一过零数是否小于第二过零数；如果是，则根据转角及其个数得到待分析信号的分数阶傅里叶变换的前半段，否则得到后半段；获取转角及其个数下变换结果的方差及对应的方差值；获取最大方差值，并得到对应的旋转角度；根据调频率最大值、旋转角度得到调频率。本发明实施例的估算方法不但计算速度快，而且可以根据绝对误差要求调整计算速度，提高调频率的估计精度。

Description

线性调频信号调频率的估计方法及装置

技术领域

本发明涉及信号处理技术领域，特别涉及一种线性调频信号调频率的估计方法及装置。

背景技术

线性调频信号具有低峰值功率、大时宽带宽积的特性，有较好的可压缩性，从而使得线性调频信号在声呐和雷达等领域得到广泛应用。其中，调频率是表征线性调频信号频率特性的基本参数，其估计问题是一直是信号处理领域的重要研究内容。

相关技术中，常用的线性调频信号调频率的估计方法有子孔径相关法、PGA(PhaseGradient Autofocus，相位梯度自聚焦法)法、最小熵方法和Wigner-Hough变换等等。然而，子孔径相关法要求子孔径间的相关性较强，才能实现较高精度的参数估计，一旦场景内分辨单元反射率低、无明显特征或强散射点，则计算误差较大；PGA法虽然精度较高，但计算量较大；最小熵值方法可以利用图像熵和相位误差的函数关系估算调频率，但估计相位误差时所需迭代次数较多，计算量较大；Wigner-Hough变换可以利用基于回波信号的二维时频平面进行Hough变换来估算调频率，该算法的优点是在高噪声条件下仍能得到高精度的结果，但是同样也有计算量大的缺点。

发明内容

本发明旨在至少在一定程度上解决相关技术中的技术问题之一。

为此，本发明的一个目的在于提出一种线性调频信号调频率的估计方法，该方法可以调整计算速度，并且提高调频率的估计精度。

本发明的另一个目的在于提出一种线性调频信号调频率的估计装置。

为达到上述目的，本发明一方面实施例提出了一种线性调频信号调频率的估计方法，包括以下步骤：根据信号采样频率与信号采样时长获取待分析信号的调频率最大值；获取调频率误差，以得到分数阶傅里叶变换的转角及其个数；将所述待分析信号均分为第一分析信号与第二分析信号；统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数；判断所述第一过零数是否小于所述第二过零数；如果所述第一过零数小于所述第二过零，则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段，否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段；获取所述转角及其个数下变换结果的方差及对应的方差值；获取最大方差值，并得到对应的旋转角度；根据所述调频率最大值、所述旋转角度得到调频率。

本发明实施例的线性调频信号调频率的估计方法，可以基于分数阶傅里叶变换，以信号的部分变换系数的方差追踪最佳旋转角，来实现对调频率的估算，具有计算量小的优点，且可以根据调频率的精度要求，调节计算量，使计算速度具备适应性，提高了线性调频信号调频率的计算速度。

另外，根据本发明上述实施例的线性调频信号调频率的估计方法还可以具有以下附加的技术特征：

进一步地，在本发明的一个实施例中，所述调频率最大值的计算公式为：

$k_{max}=\frac{f_s}{\mathrm{\Δ}t},$

其中，k_max为所述调频率最大值，f_s为所述信号采样频率，Δt为所述信号采样时长。

进一步地，在本发明的一个实施例中，所述转角及其个数的计算公式为：

其中，c为常系数，α和n分别为所述转角及其个数。

进一步地，在本发明的一个实施例中，所述获取所述转角及其个数下变换结果的方差及对应的方差值进一步包括：计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)，其中，方差计算公式为：

D＝E{X(u)-E[X(u)]}²，

其中，X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段，N为信号的采样点数。

另外，在本发明的一个实施例中，如果所述第一过零数小于所述第二过零，则所述调频率的计算公式为：

K_p＝k_max·cot(α_o)，

如果所述第一过零数大于所述第二过零，则所述调频率的计算公式为：

K_p＝-k_max·cot(α_o)，

其中，K_p为所述调频率，α_o为分数阶傅里叶变换方差取得最大值对应的旋转角度。

为达到上述目的，本发明另一方面实施例提出了一种线性调频信号调频率的估计装置，包括：第一获取模块，用于根据信号采样频率与信号采样时长获取待分析信号的调频率最大值，并且获取调频率误差，以得到分数阶傅里叶变换的转角及其个数；处理模块，用于将所述待分析信号均分为第一分析信号与第二分析信号；，并且统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数；判断模块，用于判断所述第一过零数是否小于所述第二过零数；计算模块，当所述第一过零数小于所述第二过零时，用于根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段，否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段；第二获取模块，用于获取所述转角及其个数下变换结果的方差及对应的方差值，并且获取最大方差值，并得到对应的旋转角度，以根据所述调频率最大值、所述旋转角度得到调频率。

本发明实施例的线性调频信号调频率的估计装置，可以基于分数阶傅里叶变换，以信号的部分变换系数的方差追踪最佳旋转角，来实现对调频率的估算，具有计算量小的优点，且可以根据调频率的精度要求，调节计算量，使计算速度具备适应性，提高了线性调频信号调频率的计算速度。

另外，根据本发明上述实施例的线性调频信号调频率的估计装置还可以具有以下附加的技术特征：

进一步地，在本发明的一个实施例中，所述调频率最大值的计算公式为：

$k_{max}=\frac{f_s}{\mathrm{\Δ}t},$

其中，k_max为所述调频率最大值，f_s为所述信号采样频率，Δt为所述信号采样时长。

进一步地，在本发明的一个实施例中，所述转角及其个数的计算公式为：

其中，c为常系数，α和n分别为所述转角及其个数。

进一步地，在本发明的一个实施例中，所述第二获取模块还用于计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)，其中，方差计算公式为：

D＝E{X(u)-E[X(u)]}²，

其中，X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段，N为信号的采样点数。

另外，在本发明的一个实施例中，当所述第一过零数小于所述第二过零时，所述调频率的计算公式为：

K_p＝k_max·cot(α_o)，

当所述第一过零数大于所述第二过零时，所述调频率的计算公式为：

K_p＝-k_max·cot(α_o)，

其中，K_p为所述调频率，α_o为分数阶傅里叶变换方差取得最大值对应的旋转角度。

本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：

图1为根据本发明实施例的线性调频信号调频率的估计方法的流程图；

图2为根据本发明一个实施例的线性调频信号调频率的估计方法的流程图；

图3为根据本发明一个实施例的不同转角下的信号s的前半段分数阶傅里叶变换的方差示意图；

图4为根据本发明一个实施例的不同c值下的调频率估计值示意图；

图5为根据本发明一个实施例的不同c值下的调频率估计值的误差百分比示意图；以及

图6为根据本发明实施例的线性调频信号调频率的估计装置的结构示意图。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。

下面参照附图描述根据本发明实施例提出的线性调频信号调频率的估计方法及装置，首先将参照附图描述根据本发明实施例提出的线性调频信号调频率的估计方法。

图1是本发明实施例的线性调频信号调频率的估计方法的流程图

如图1所示，该线性调频信号调频率的估计方法包括以下步骤：

在步骤S101中，根据信号采样频率与信号采样时长获取待分析信号的调频率最大值。

其中，在本发明的一个实施例中，调频率最大值的计算公式为：

$k_{max}=\frac{f_s}{\mathrm{\Δ}t},$

其中，k_max为调频率最大值，f_s为信号采样频率，Δt为信号采样时长。

也就是说，如图2所示，第一步：计算待分析信号s调频率最大值。其中，第一步包括信号调频率最大值k_max的计算过程。最大调频率k_max计算公式为：

在步骤S102中，获取调频率误差，以得到分数阶傅里叶变换的转角及其个数。

其中，在本发明的一个实施例中，转角及其个数的计算公式为：

其中，c为常系数，α和n分别为转角及其个数。

进一步地，第二步：输入调频率误差ε，计算分数阶傅里叶变换转角α及其个数n，其中，第二步包括根据最大调频率绝对误差计算分数阶傅里叶变换旋转角度α个数n的计算过程。

旋转角度α及其个数n的计算公式为：

其中，c为常系数，为了增加分辨率，可适当增加n值，c一般大于1，为2的整数幂较好，α单位为°。

在步骤S103中，将待分析信号均分为第一分析信号与第二分析信号。

进一步地，第三步：将信号s均分为两段，且为s₁和s₂。其中，第三步包括了判断调频率符号的过程，需要先将信号均分为两部份s₁和s₂。

在步骤S104中，统计第一分析信号与第二分析信号对应的第一过零数与第二过零数。

进一步地，第四步：统计s₁和s₂的过零数，且对应分别为num₁和num₂。其中，第四步需要计算s₁和s₂过零次数，分别为num₁和num₂，计算过程为：

待处理信号为x(j)(j＝1、2、3、…、J，J为信号长度)，构建两个信号x₁(j₁)(j₁＝2、3、4、…、J_-1，J为信号x(j)长度)和x₂(j₂)(j₂＝2、3、4、…、J，J为信号x(j)长度)，生成新的信号x₀(j₀)：

x₀(j₀)＝x₁(j₁)·x₂(j₂)， 式4

统计x₀(j₀)中值小于0元素的个数。分别对s₁和s₂进行上述过程处理，得到num₁和num₂。

在步骤S105中，判断第一过零数是否小于第二过零数。

进一步地，第五步：对比num₁是否小于num₂。即言，第五步包括了对比num₁是否小于num₂大小的过程。

在步骤S106中，如果第一过零数小于第二过零，则根据转角及其个数得到待分析信号的分数阶傅里叶变换的前半段，否则根据转角及其个数得到待分析信号的分数阶傅里叶变换的后半段。

进一步地，第六步：当num₁小于num₂时，则计算n个转角α下，信号s分数阶傅里叶变换的前半段，即，X(u)，(为信号的采样点数)；否则，计算信号s分数阶傅里叶变换的后半段，即，X(u)，( 为信号的采样点数)。

具体地，第六步包括了根据num₁和num₂大小，分别计算信号s对应范围的分数阶傅里叶变换的过程，且为：

当num₁小于num₂时，则计算n个转角α下，信号s分数阶傅里叶变换的前半段，即，X(u)，(N为信号的采样点数)；否则，计算信号s分数阶傅里叶变换的后半段，即，X(u)，(为信号的采样点数)。

以此减少计算量，加快计算速度。

在步骤S107中，获取转角及其个数下变换结果的方差及对应的方差值。

进一步地，在本发明的一个实施例中，获取转角及其个数下变换结果的方差及对应的方差值进一步包括：计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)，其中，方差计算公式为：

D＝E{X(u)-E[X(u)]}²，

其中，X(u)为分数阶傅里叶变换的前半段或者分数阶傅里叶变换的后半段， N为信号的采样点数。

进一步地，第七步：计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)。其中，第七步包括计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)；方差计算公式为：

D＝E{X(u)-E[X(u)]}²。 式5

在步骤S108中，获取最大方差值，并得到对应的旋转角度。

进一步地，第八步：找到最大方差值D_max，并找到对应的旋转角度α_o。即言，第八步包括根据最大方差值D_max，并找到对应的旋转角度α_o。

在步骤S109中，根据调频率最大值、旋转角度得到调频率。

具体地，在本发明的一个实施例中，如果第一过零数小于第二过零，则调频率的计算公式为：

K_p＝k_max·cot(α_o)，

如果第一过零数大于第二过零，则调频率的计算公式为：

K_p＝-k_max·cot(α_o)，

其中，K_p为调频率，α_o为分数阶傅里叶变换方差取得最大值对应的旋转角度。。

进一步地，第九步：计算调频率K_p。即言，第九步包括线性调频信号调频率K_p的计算过程：

调频率K_p的计算公式为，

当num₁小于num₂时，

K_p＝k_max·cot(α_o)， 式6

当num₁大于等于num₂时，

K_p＝-k_max·cot(α_o)。 式7

为了便于本领域技术人员理解，下面以一个具体实施例进行详细描述。

该方法可以在matlab软件平台下进行实现。示例信号表达式可以为：

s＝exp[200·j·2π·t²]，0≤t≤1， 式8

即中心频率为0，调频率K_p为200Hz/s，t∈[0,1]的一组扫频信号，采样频率f_s为2000Hz。为验证方法的适用性，加上强度为0.2的白噪声信号，信噪比SNR为26.95db。

根据第一步中式1，最大调频率为2000Hz/s。

根据第二步，输入调频率误差ε为4Hz/s，根据式2可算得转角个数n为500(c＝1)。

根据第三步，将信号均分为两段，且为s₁和s₂。

根据第四步，统计s₁和s₂的过零数，且对应分别为num₁＝56和num₂＝150。

根据第五步，比较num₁和num₂，则可知num₁小于num₂。

根据第六步，因为num₁小于num₂，需要计算n个转角α下，信号s分数阶傅里叶变换的前半段，即，X(u)，(N为信号的采样点数)。

根据第七步，计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、500)。图3为不同转角下的信号s的前半段分数阶傅里叶变换的方差。

根据第八步，找到最大方差值D_max＝0.58，并找到对应的旋转角度α_o＝84.24°。

根据第九步，因为num₁小于num₂，需根据式5计算调频率，K_p＝201.74Hz/s。绝对误差为1.74Hz/s，误差百分比为0.87％，为最大调频率的0.087％。由此可知，该方法具备较高精度。

根据第二步可知，增大c值可以提高精度，为此，设定c＝2、4、6、8、16、32、64，分别对扫频率K_p进行估算，其估测值和误差值分别如图4和图5所示。由图4可知，随着c值的增大，调频率估计值越来越趋近于200Hz/s，误差百分比也越来越小，当c＝32时，误差百分比为0.02％。由此可知，本发明实施例的方法可以根据精度要求调节计算量，更具灵活性。

根据本发明实施例的线性调频信号调频率的估计方法，可以基于分数阶傅里叶变换，以信号的部分变换系数的方差追踪最佳旋转角，来实现对调频率的估算，具有计算量小的优点，且可以根据调频率的精度要求，调节计算量，使计算速度具备适应性，提高了线性调频信号调频率的计算速度。

其次参照附图描述根据本发明实施例提出的线性调频信号调频率的估计装置。

图6是本发明实施例的线性调频信号调频率的估计装置的结构示意图。

如图6所示，该线性调频信号调频率的估计装置10包括：第一获取模块100、处理模块200、判断模块300、计算模块400和第二获取模块500。

其中，第一获取模块100用于根据信号采样频率与信号采样时长获取待分析信号的调频率最大值，并且获取调频率误差，以得到分数阶傅里叶变换的转角及其个数。处理模块200用于将待分析信号均分为第一分析信号与第二分析信号；，并且统计第一分析信号与第二分析信号对应的第一过零数与第二过零数。判断模块300用于判断第一过零数是否小于第二过零数。当第一过零数小于第二过零时，计算模块400用于根据转角及其个数得到待分析信号的分数阶傅里叶变换的前半段，否则根据转角及其个数得到待分析信号的分数阶傅里叶变换的后半段。第二获取模块500用于获取转角及其个数下变换结果的方差及对应的方差值，并且获取最大方差值，并得到对应的旋转角度，以根据调频率最大值、旋转角度得到调频率。本发明实施例的估计装置10不但计算速度快，而且可以根据绝对误差要求调整计算速度，提高调频率的估计精度。

进一步地，在本发明的一个实施例中，调频率最大值的计算公式为：

$k_{max}=\frac{f_s}{\mathrm{\Δ}t},$

其中，k_max为调频率最大值，f_s为信号采样频率，Δt为信号采样时长。

进一步地，在本发明的一个实施例中，转角及其个数的计算公式为：

其中，c为常系数，α和n分别为转角及其个数。

进一步地，在本发明的一个实施例中，第二获取模块500还用于计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)，其中，方差计算公式为：

D＝E{X(u)-E[X(u)]}²，

其中，X(u)为分数阶傅里叶变换的前半段或者分数阶傅里叶变换的后半段， N为信号的采样点数。

另外，在本发明的一个实施例中，当第一过零数小于第二过零时，调频率的计算公式为：

K_p＝k_max·cot(α_o)，

当第一过零数大于第二过零时，调频率的计算公式为：

K_p＝-k_max·cot(α_o)，

其中，K_p为调频率，α_o为分数阶傅里叶变换方差取得最大值对应的旋转角度。

需要说明的是，前述对线性调频信号调频率的估计方法实施例的解释说明也适用于该实施例的线性调频信号调频率的估计装置，此处不再赘述。

根据本发明实施例的线性调频信号调频率的估计装置，可以基于分数阶傅里叶变换，以信号的部分变换系数的方差追踪最佳旋转角，来实现对调频率的估算，具有计算量小的优点，且可以根据调频率的精度要求，调节计算量，使计算速度具备适应性，提高了线性调频信号调频率的计算速度。

在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。

在本发明中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。

在本发明中，除非另有明确的规定和限定，第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触，或第一和第二特征通过中间媒介间接接触。而且，第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方，或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方，或仅仅表示第一特征水平高度小于第二特征。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (10)

1.一种线性调频信号调频率的估计方法，其特征在于，包括以下步骤：

根据信号采样频率与信号采样时长获取待分析信号的调频率最大值；

获取调频率误差，以得到分数阶傅里叶变换的转角及其个数；

将所述待分析信号均分为第一分析信号与第二分析信号；

统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数；

判断所述第一过零数是否小于所述第二过零数；

如果所述第一过零数小于所述第二过零，则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段，否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段；

获取所述转角及其个数下变换结果的方差及对应的方差值；

获取最大方差值，并得到对应的旋转角度；以及

根据所述调频率最大值、所述旋转角度得到调频率。

2.根据权利要求1所述的线性调频信号调频率的估计方法，其特征在于，所述调频率最大值的计算公式为：

$k_{max}=\frac{f_s}{\mathrm{\Δ}t},$

其中，k_max为所述调频率最大值，f_s为所述信号采样频率，Δt为所述信号采样时长。

3.根据权利要求2所述的线性调频信号调频率的估计方法，其特征在于，所述转角及其个数的计算公式为：

其中，(i＝1、2、…、n)，c为常系数，α和n分别为所述转角及其个数。

4.根据权利要求3所述的线性调频信号调频率的估计方法，其特征在于，所述获取所述转角及其个数下变换结果的方差及对应的方差值进一步包括：

计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)，其中，方差计算公式为：

D＝E{X(u)-E[X(u)]}²，

其中，X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段，u＝1、2、3、…、N为信号的采样点数。

5.根据权利要求4所述的线性调频信号调频率的估计方法，其特征在于，其中，

如果所述第一过零数小于所述第二过零，则所述调频率的计算公式为：

K_p＝k_max·cot(α_o)，

如果所述第一过零数大于所述第二过零，则所述调频率的计算公式为：

K_p＝-k_max·cot(α_o)，

其中，K_p为所述调频率，α_o为分数阶傅里叶变换方差取得最大值对应的旋转角度。

6.一种线性调频信号调频率的估计装置，其特征在于，包括：

第一获取模块，用于根据信号采样频率与信号采样时长获取待分析信号的调频率最大值，并且获取调频率误差，以得到分数阶傅里叶变换的转角及其个数；

处理模块，用于将所述待分析信号均分为第一分析信号与第二分析信号；，并且统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数；

判断模块，用于判断所述第一过零数是否小于所述第二过零数；

计算模块，当所述第一过零数小于所述第二过零时，用于根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段，否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段；以及

第二获取模块，用于获取所述转角及其个数下变换结果的方差及对应的方差值，并且获取最大方差值，并得到对应的旋转角度，以根据所述调频率最大值、所述旋转角度得到调频率。

7.根据权利要求6所述的线性调频信号调频率的估计装置，其特征在于，所述调频率最大值的计算公式为：

$k_{max}=\frac{f_s}{\mathrm{\Δ}t},$

其中，k_max为所述调频率最大值，f_s为所述信号采样频率，Δt为所述信号采样时长。

8.根据权利要求7所述的线性调频信号调频率的估计装置，其特征在于，所述转角及其个数的计算公式为：

其中，(i＝1、2、…、n)，c为常系数，α和n分别为所述转角及其个数。

9.根据权利要求8所述的线性调频信号调频率的估计装置，其特征在于，所述第二获取模块还用于计算n个转角α下变换结果的方差，得到对应方差值D(i)(i＝1、2、3、…、n)，其中，方差计算公式为：

D＝E{X(u)-E[X(u)]}²，

其中，X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段，u＝1、2、3、…、N为信号的采样点数。

10.根据权利要求9所述的线性调频信号调频率的估计装置，其特征在于，其中，

当所述第一过零数小于所述第二过零时，所述调频率的计算公式为：

K_p＝k_max·cot(α_o)，

当所述第一过零数大于所述第二过零时，所述调频率的计算公式为：

K_p＝-k_max·cot(α_o)，

其中，K_p为所述调频率，α_o为分数阶傅里叶变换方差取得最大值对应的旋转角度。