CN106093885B - 线性调频信号调频率的估计方法及装置 - Google Patents
线性调频信号调频率的估计方法及装置 Download PDFInfo
- Publication number
- CN106093885B CN106093885B CN201610379315.0A CN201610379315A CN106093885B CN 106093885 B CN106093885 B CN 106093885B CN 201610379315 A CN201610379315 A CN 201610379315A CN 106093885 B CN106093885 B CN 106093885B
- Authority
- CN
- China
- Prior art keywords
- signal
- zero
- modulation rate
- frequency modulation
- corner
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/534—Details of non-pulse systems
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种线性调频信号调频率的估计方法及装置,其中,方法包括以下步骤:获取待分析信号的调频率最大值;获取调频率误差,以得到分数阶傅里叶变换的转角及其个数;将信号均分,以统计信号对应的第一过零数与第二过零数;判断第一过零数是否小于第二过零数;如果是,则根据转角及其个数得到待分析信号的分数阶傅里叶变换的前半段,否则得到后半段;获取转角及其个数下变换结果的方差及对应的方差值;获取最大方差值,并得到对应的旋转角度;根据调频率最大值、旋转角度得到调频率。本发明实施例的估算方法不但计算速度快,而且可以根据绝对误差要求调整计算速度,提高调频率的估计精度。
Description
技术领域
本发明涉及信号处理技术领域,特别涉及一种线性调频信号调频率的估计方法及装置。
背景技术
线性调频信号具有低峰值功率、大时宽带宽积的特性,有较好的可压缩性,从而使得线性调频信号在声呐和雷达等领域得到广泛应用。其中,调频率是表征线性调频信号频率特性的基本参数,其估计问题是一直是信号处理领域的重要研究内容。
相关技术中,常用的线性调频信号调频率的估计方法有子孔径相关法、PGA(PhaseGradient Autofocus,相位梯度自聚焦法)法、最小熵方法和Wigner-Hough变换等等。然而,子孔径相关法要求子孔径间的相关性较强,才能实现较高精度的参数估计,一旦场景内分辨单元反射率低、无明显特征或强散射点,则计算误差较大;PGA法虽然精度较高,但计算量较大;最小熵值方法可以利用图像熵和相位误差的函数关系估算调频率,但估计相位误差时所需迭代次数较多,计算量较大;Wigner-Hough变换可以利用基于回波信号的二维时频平面进行Hough变换来估算调频率,该算法的优点是在高噪声条件下仍能得到高精度的结果,但是同样也有计算量大的缺点。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种线性调频信号调频率的估计方法,该方法可以调整计算速度,并且提高调频率的估计精度。
本发明的另一个目的在于提出一种线性调频信号调频率的估计装置。
为达到上述目的,本发明一方面实施例提出了一种线性调频信号调频率的估计方法,包括以下步骤:根据信号采样频率与信号采样时长获取待分析信号的调频率最大值;获取调频率误差,以得到分数阶傅里叶变换的转角及其个数;将所述待分析信号均分为第一分析信号与第二分析信号;统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数;判断所述第一过零数是否小于所述第二过零数;如果所述第一过零数小于所述第二过零数,则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段,否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段;获取所述转角及其个数下变换结果的方差及对应的方差值;获取最大方差值,并得到对应的旋转角度;根据所述调频率最大值、所述旋转角度得到调频率。
本发明实施例的线性调频信号调频率的估计方法,可以基于分数阶傅里叶变换,以信号的部分变换系数的方差追踪最佳旋转角,来实现对调频率的估算,具有计算量小的优点,且可以根据调频率的精度要求,调节计算量,使计算速度具备适应性,提高了线性调频信号调频率的计算速度。
另外,根据本发明上述实施例的线性调频信号调频率的估计方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述调频率最大值的计算公式为:
其中,kmax为所述调频率最大值,fs为所述信号采样频率,Δt为所述信号采样时长。
进一步地,在本发明的一个实施例中,所述转角及其个数的计算公式为:
其中,c为常系数,α和n分别为所述转角及其个数。
进一步地,在本发明的一个实施例中,所述获取所述转角及其个数下变换结果的方差及对应的方差值进一步包括:计算n个转角α下变换结果的方差,得到对应方差值D(i)(i=1、2、3、…、n),其中,方差计算公式为:
D=E{X(u)-E[X(u)]}2,
其中,X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段,N为信号的采样点数。
另外,在本发明的一个实施例中,如果所述第一过零数小于所述第二过零数,则所述调频率的计算公式为:
Kp=kmax·cot(αo),
如果所述第一过零数大于所述第二过零数,则所述调频率的计算公式为:
Kp=-kmax·cot(αo),
其中,Kp为所述调频率,αo为分数阶傅里叶变换方差取得最大值对应的旋转角度。
为达到上述目的,本发明另一方面实施例提出了一种线性调频信号调频率的估计装置,包括:第一获取模块,用于根据信号采样频率与信号采样时长获取待分析信号的调频率最大值,并且获取调频率误差,以得到分数阶傅里叶变换的转角及其个数;处理模块,用于将所述待分析信号均分为第一分析信号与第二分析信号;,并且统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数;判断模块,用于判断所述第一过零数是否小于所述第二过零数;计算模块,当所述第一过零数小于所述第二过零数时,用于根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段,否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段;第二获取模块,用于获取所述转角及其个数下变换结果的方差及对应的方差值,并且获取最大方差值,并得到对应的旋转角度,以根据所述调频率最大值、所述旋转角度得到调频率。
本发明实施例的线性调频信号调频率的估计装置,可以基于分数阶傅里叶变换,以信号的部分变换系数的方差追踪最佳旋转角,来实现对调频率的估算,具有计算量小的优点,且可以根据调频率的精度要求,调节计算量,使计算速度具备适应性,提高了线性调频信号调频率的计算速度。
另外,根据本发明上述实施例的线性调频信号调频率的估计装置还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述调频率最大值的计算公式为:
其中,kmax为所述调频率最大值,fs为所述信号采样频率,Δt为所述信号采样时长。
进一步地,在本发明的一个实施例中,所述转角及其个数的计算公式为:
其中,c为常系数,α和n分别为所述转角及其个数。
进一步地,在本发明的一个实施例中,所述第二获取模块还用于计算n个转角α下变换结果的方差,得到对应方差值D(i)(i=1、2、3、…、n),其中,方差计算公式为:
D=E{X(u)-E[X(u)]}2,
其中,X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段,N为信号的采样点数。
另外,在本发明的一个实施例中,当所述第一过零数小于所述第二过零数时,所述调频率的计算公式为:
Kp=kmax·cot(αo),
当所述第一过零数大于所述第二过零数时,所述调频率的计算公式为:
Kp=-kmax·cot(αo),
其中,Kp为所述调频率,αo为分数阶傅里叶变换方差取得最大值对应的旋转角度。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的线性调频信号调频率的估计方法的流程图;
图2为根据本发明一个实施例的线性调频信号调频率的估计方法的流程图;
图3为根据本发明一个实施例的不同转角下的信号s的前半段分数阶傅里叶变换的方差示意图;
图4为根据本发明一个实施例的不同c值下的调频率估计值示意图;
图5为根据本发明一个实施例的不同c值下的调频率估计值的误差百分比示意图;以及
图6为根据本发明实施例的线性调频信号调频率的估计装置的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的线性调频信号调频率的估计方法及装置,首先将参照附图描述根据本发明实施例提出的线性调频信号调频率的估计方法。
图1是本发明实施例的线性调频信号调频率的估计方法的流程图
如图1所示,该线性调频信号调频率的估计方法包括以下步骤:
在步骤S101中,根据信号采样频率与信号采样时长获取待分析信号的调频率最大值。
其中,在本发明的一个实施例中,调频率最大值的计算公式为:
其中,kmax为调频率最大值,fs为信号采样频率,Δt为信号采样时长。
也就是说,如图2所示,第一步:计算待分析信号s调频率最大值。其中,第一步包括信号调频率最大值kmax的计算过程。最大调频率kmax计算公式为:
在步骤S102中,获取调频率误差,以得到分数阶傅里叶变换的转角及其个数。
其中,在本发明的一个实施例中,转角及其个数的计算公式为:
其中,c为常系数,α和n分别为转角及其个数。
进一步地,第二步:输入调频率误差ε,计算分数阶傅里叶变换转角α及其个数n,其中,第二步包括根据最大调频率绝对误差计算分数阶傅里叶变换旋转角度α个数n的计算过程。
旋转角度α及其个数n的计算公式为:
其中,c为常系数,为了增加分辨率,可适当增加n值,c一般大于1,为2的整数幂较好,α单位为°。
在步骤S103中,将待分析信号均分为第一分析信号与第二分析信号。
进一步地,第三步:将信号s均分为两段,且为s1和s2。其中,第三步包括了判断调频率符号的过程,需要先将信号均分为两部份s1和s2。
在步骤S104中,统计第一分析信号与第二分析信号对应的第一过零数与第二过零数。
进一步地,第四步:统计s1和s2的过零数,且对应分别为num1和num2。其中,第四步需要计算s1和s2过零次数,分别为num1和num2,计算过程为:
待处理信号为x(j)(j=1、2、3、…、J,J为信号长度),构建两个信号x1(j1)(j1=1、2、3、4、…、J-1,J为信号x(j)长度)和x2(j2)(j2=2、3、4、…、J,J为信号x(j)长度),生成新的信号x0(j0):
x0(j0)=x1(j1)·x2(j2), 式4
统计x0(j0)中值小于0元素的个数。分别对s1和s2进行上述过程处理,得到num1和num2。
在步骤S105中,判断第一过零数是否小于第二过零数。
进一步地,第五步:对比num1是否小于num2。即言,第五步包括了对比num1是否小于num2大小的过程。
在步骤S106中,如果第一过零数小于第二过零数,则根据转角及其个数得到待分析信号的分数阶傅里叶变换的前半段,否则根据转角及其个数得到待分析信号的分数阶傅里叶变换的后半段。
进一步地,第六步:当num1小于num2时,则计算n个转角α下,信号s分数阶傅里叶变换的前半段,即,X(u),(N为信号的采样点数);否则,计算信号s分数阶傅里叶变换的后半段,即,X(u),( N为信号的采样点数)。
具体地,第六步包括了根据num1和num2大小,分别计算信号s对应范围的分数阶傅里叶变换的过程,且为:
当num1小于num2时,则计算n个转角α下,信号s分数阶傅里叶变换的前半段,即,X(u),(N为信号的采样点数);否则,计算信号s分数阶傅里叶变换的后半段,即,X(u),(N为信号的采样点数)。
以此减少计算量,加快计算速度。
在步骤S107中,获取转角及其个数下变换结果的方差及对应的方差值。
进一步地,在本发明的一个实施例中,获取转角及其个数下变换结果的方差及对应的方差值进一步包括:计算n个转角α下变换结果的方差,得到对应方差值D(i)(i=1、2、3、…、n),其中,方差计算公式为:
D=E{X(u)-E[X(u)]}2,
其中,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)]}2。 式5
在步骤S108中,获取最大方差值,并得到对应的旋转角度。
进一步地,第八步:找到最大方差值Dmax,并找到对应的旋转角度αo。即言,第八步包括根据最大方差值Dmax,并找到对应的旋转角度αo。
在步骤S109中,根据调频率最大值、旋转角度得到调频率。
具体地,在本发明的一个实施例中,如果第一过零数小于第二过零数,则调频率的计算公式为:
Kp=kmax·cot(αo),
如果第一过零数大于第二过零数,则调频率的计算公式为:
Kp=-kmax·cot(αo),
其中,Kp为调频率,αo为分数阶傅里叶变换方差取得最大值对应的旋转角度。。
进一步地,第九步:计算调频率Kp。即言,第九步包括线性调频信号调频率Kp的计算过程:
调频率Kp的计算公式为,
当num1小于num2时,
Kp=kmax·cot(αo), 式6
当num1大于等于num2时,
Kp=-kmax·cot(αo)。 式7
为了便于本领域技术人员理解,下面以一个具体实施例进行详细描述。
该方法可以在matlab软件平台下进行实现。示例信号表达式可以为:
s=exp[200·j·2π·t2],0≤t≤1, 式8
即中心频率为0,调频率Kp为200Hz/s,t∈[0,1]的一组扫频信号,采样频率fs为2000Hz。为验证方法的适用性,加上强度为0.2的白噪声信号,信噪比SNR为26.95db。
根据第一步中式1,最大调频率为2000Hz/s。
根据第二步,输入调频率误差ε为4Hz/s,根据式2可算得转角个数n为500(c=1)。
根据第三步,将信号均分为两段,且为s1和s2。
根据第四步,统计s1和s2的过零数,且对应分别为num1=56和num2=150。
根据第五步,比较num1和num2,则可知num1小于num2。
根据第六步,因为num1小于num2,需要计算n个转角α下,信号s分数阶傅里叶变换的前半段,即,X(u),(N为信号的采样点数)。
根据第七步,计算n个转角α下变换结果的方差,得到对应方差值D(i)(i=1、2、3、…、500)。图3为不同转角下的信号s的前半段分数阶傅里叶变换的方差。
根据第八步,找到最大方差值Dmax=0.58,并找到对应的旋转角度αo=84.24°。
根据第九步,因为num1小于num2,需根据式5计算调频率,Kp=201.74Hz/s。绝对误差为1.74Hz/s,误差百分比为0.87%,为最大调频率的0.087%。由此可知,该方法具备较高精度。
根据第二步可知,增大c值可以提高精度,为此,设定c=2、4、6、8、16、32、64,分别对扫频率Kp进行估算,其估测值和误差值分别如图4和图5所示。由图4可知,随着c值的增大,调频率估计值越来越趋近于200Hz/s,误差百分比也越来越小,当c=32时,误差百分比为0.02%。由此可知,本发明实施例的方法可以根据精度要求调节计算量,更具灵活性。
根据本发明实施例的线性调频信号调频率的估计方法,可以基于分数阶傅里叶变换,以信号的部分变换系数的方差追踪最佳旋转角,来实现对调频率的估算,具有计算量小的优点,且可以根据调频率的精度要求,调节计算量,使计算速度具备适应性,提高了线性调频信号调频率的计算速度。
其次参照附图描述根据本发明实施例提出的线性调频信号调频率的估计装置。
图6是本发明实施例的线性调频信号调频率的估计装置的结构示意图。
如图6所示,该线性调频信号调频率的估计装置10包括:第一获取模块100、处理模块200、判断模块300、计算模块400和第二获取模块500。
其中,第一获取模块100用于根据信号采样频率与信号采样时长获取待分析信号的调频率最大值,并且获取调频率误差,以得到分数阶傅里叶变换的转角及其个数。处理模块200用于将待分析信号均分为第一分析信号与第二分析信号;,并且统计第一分析信号与第二分析信号对应的第一过零数与第二过零数。判断模块300用于判断第一过零数是否小于第二过零数。当第一过零数小于第二过零数时,计算模块400用于根据转角及其个数得到待分析信号的分数阶傅里叶变换的前半段,否则根据转角及其个数得到待分析信号的分数阶傅里叶变换的后半段。第二获取模块500用于获取转角及其个数下变换结果的方差及对应的方差值,并且获取最大方差值,并得到对应的旋转角度,以根据调频率最大值、旋转角度得到调频率。本发明实施例的估计装置10不但计算速度快,而且可以根据绝对误差要求调整计算速度,提高调频率的估计精度。
进一步地,在本发明的一个实施例中,调频率最大值的计算公式为:
其中,kmax为调频率最大值,fs为信号采样频率,Δt为信号采样时长。
进一步地,在本发明的一个实施例中,转角及其个数的计算公式为:
其中,c为常系数,α和n分别为转角及其个数。
进一步地,在本发明的一个实施例中,第二获取模块500还用于计算n个转角α下变换结果的方差,得到对应方差值D(i)(i=1、2、3、…、n),其中,方差计算公式为:
D=E{X(u)-E[X(u)]}2,
其中,X(u)为分数阶傅里叶变换的前半段或者分数阶傅里叶变换的后半段, N为信号的采样点数。
另外,在本发明的一个实施例中,当第一过零数小于第二过零数时,调频率的计算公式为:
Kp=kmax·cot(αo),
当第一过零数大于第二过零数时,调频率的计算公式为:
Kp=-kmax·cot(αo),
其中,Kp为调频率,αo为分数阶傅里叶变换方差取得最大值对应的旋转角度。
需要说明的是,前述对线性调频信号调频率的估计方法实施例的解释说明也适用于该实施例的线性调频信号调频率的估计装置,此处不再赘述。
根据本发明实施例的线性调频信号调频率的估计装置,可以基于分数阶傅里叶变换,以信号的部分变换系数的方差追踪最佳旋转角,来实现对调频率的估算,具有计算量小的优点,且可以根据调频率的精度要求,调节计算量,使计算速度具备适应性,提高了线性调频信号调频率的计算速度。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (6)
1.一种线性调频信号调频率的估计方法,其特征在于,包括以下步骤:
根据信号采样频率与信号采样时长获取待分析信号的调频率最大值;
获取调频率误差,以得到分数阶傅里叶变换的转角及其个数,其中,所述转角及其个数的计算公式为:
其中,i=1、2、…、n,c为常系数,α和n分别为所述转角及其个数,ε为调频率误差;
将所述待分析信号均分为两段,分别作为第一分析信号与第二分析信号;
统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数;
判断所述第一过零数是否小于所述第二过零数;
如果所述第一过零数小于所述第二过零数,则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段,否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段,其中,
如果所述第一过零数小于所述第二过零数,则所述调频率的计算公式为:
Kp=kmax·cot(αo),
如果所述第一过零数大于所述第二过零数,则所述调频率的计算公式为:
Kp=-kmax·cot(αo),
其中,Kp为所述调频率,kmax为所述调频率最大值,αo为分数阶傅里叶变换方差取得最大值对应的旋转角度;
获取所述转角及其个数下分数阶傅里叶变换结果的方差及对应的方差值;
获取最大方差值,并得到对应的旋转角度;以及
根据所述调频率最大值、所述旋转角度得到调频率。
2.根据权利要求1所述的线性调频信号调频率的估计方法,其特征在于,所述调频率最大值的计算公式为:
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,kmax为所述调频率最大值,fs为所述信号采样频率,Δt为所述信号采样时长。
3.根据权利要求1所述的线性调频信号调频率的估计方法,其特征在于,所述获取所述转角及其个数下变换结果的方差及对应的方差值进一步包括:
计算n个转角α下变换结果的方差,得到对应方差值D(i)i=1、2、3、…、n,其中,方差计算公式为:
D=E{X(u)-E[X(u)]}2,
其中,X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段,N为信号的采样点数。
4.一种线性调频信号调频率的估计装置,其特征在于,包括:
第一获取模块,用于根据信号采样频率与信号采样时长获取待分析信号的调频率最大值,并且获取调频率误差,以得到分数阶傅里叶变换的转角及其个数,其中,所述转角及其个数的计算公式为:
其中,i=1、2、…、n,c为常系数,α和n分别为所述转角及其个数,ε为调频率误差;
处理模块,用于将所述待分析信号均分为两段,分别作为第一分析信号与第二分析信号,并且统计所述第一分析信号与第二分析信号对应的第一过零数与第二过零数;
判断模块,用于判断所述第一过零数是否小于所述第二过零数;
计算模块,当所述第一过零数小于所述第二过零数时,用于根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的前半段,否则根据所述转角及其个数得到所述待分析信号的分数阶傅里叶变换的后半段,其中,
如果所述第一过零数小于所述第二过零数,则所述调频率的计算公式为:
Kp=kmax·cot(αo),
如果所述第一过零数大于所述第二过零数,则所述调频率的计算公式为:
Kp=-kmax·cot(αo),
其中,Kp为所述调频率,kmax为所述调频率最大值,αo为分数阶傅里叶变换方差取得最大值对应的旋转角度;以及
第二获取模块,用于获取所述转角及其个数下分数阶傅里叶变换结果的方差及对应的方差值,并且获取最大方差值,并得到对应的旋转角度,以根据所述调频率最大值、所述旋转角度得到调频率。
5.根据权利要求4所述的线性调频信号调频率的估计装置,其特征在于,所述调频率最大值的计算公式为:
<mrow>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,kmax为所述调频率最大值,fs为所述信号采样频率,Δt为所述信号采样时长。
6.根据权利要求4所述的线性调频信号调频率的估计装置,其特征在于,所述第二获取模块还用于计算n个转角α下变换结果的方差,得到对应方差值D(i)i=1、2、3、…、n,其中,方差计算公式为:
D=E{X(u)-E[X(u)]}2,
其中,X(u)为所述分数阶傅里叶变换的前半段或者所述分数阶傅里叶变换的后半段,N为信号的采样点数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610379315.0A CN106093885B (zh) | 2016-05-31 | 2016-05-31 | 线性调频信号调频率的估计方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610379315.0A CN106093885B (zh) | 2016-05-31 | 2016-05-31 | 线性调频信号调频率的估计方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106093885A CN106093885A (zh) | 2016-11-09 |
CN106093885B true CN106093885B (zh) | 2018-05-29 |
Family
ID=57229774
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610379315.0A Active CN106093885B (zh) | 2016-05-31 | 2016-05-31 | 线性调频信号调频率的估计方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106093885B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108594185B (zh) * | 2018-07-25 | 2021-06-04 | 电子科技大学 | 一种线性调频信号调频率的估计方法 |
CN117118536B (zh) * | 2023-10-25 | 2023-12-19 | 南京派格测控科技有限公司 | 调频稳定性的确定方法、装置、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833035A (zh) * | 2010-04-19 | 2010-09-15 | 天津大学 | 线性调频信号参数估计方法及其实施装置 |
CN102832942A (zh) * | 2011-06-16 | 2012-12-19 | 中国人民解放军海军航空工程学院 | 基于分数阶Fourier变换的三角线性调频连续波特征提取方法 |
CN102866391A (zh) * | 2012-09-05 | 2013-01-09 | 中北大学 | 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法 |
CN102999473A (zh) * | 2012-10-18 | 2013-03-27 | 中国人民解放军电子工程学院 | 一种线性调频信号的检测与参数估计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05297119A (ja) * | 1992-04-21 | 1993-11-12 | Japan Radio Co Ltd | ビートカウント型fm−cw測距装置 |
-
2016
- 2016-05-31 CN CN201610379315.0A patent/CN106093885B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101833035A (zh) * | 2010-04-19 | 2010-09-15 | 天津大学 | 线性调频信号参数估计方法及其实施装置 |
CN102832942A (zh) * | 2011-06-16 | 2012-12-19 | 中国人民解放军海军航空工程学院 | 基于分数阶Fourier变换的三角线性调频连续波特征提取方法 |
CN102866391A (zh) * | 2012-09-05 | 2013-01-09 | 中北大学 | 基于短时傅里叶变换和分数阶傅里叶变换的多目标检测方法 |
CN102999473A (zh) * | 2012-10-18 | 2013-03-27 | 中国人民解放军电子工程学院 | 一种线性调频信号的检测与参数估计方法 |
Non-Patent Citations (1)
Title |
---|
基于分数阶Fourier变换的多分量LFM信号的检测和参数估计;齐林等;《中国科学(E辑)》;20030831;第33卷(第8期);749-759 * |
Also Published As
Publication number | Publication date |
---|---|
CN106093885A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106093885B (zh) | 线性调频信号调频率的估计方法及装置 | |
CN107683423A (zh) | 用于调频连续波雷达中的对象检测的低复杂度超分辨率技术 | |
DE102005007057B4 (de) | Verfahren zur Regelung eines Fluidstroms sowie damit gesteuerte Verbrennungskraftmaschine | |
TW201513566A (zh) | 阻抗匹配方法與阻抗匹配系統 | |
DE102014114350A1 (de) | Radarvorrichtung | |
CN105137373B (zh) | 一种指数信号的去噪方法 | |
CN105652085B (zh) | 一种基于改进不完全s变换的基波和谐波参数估计方法 | |
JP2016148515A (ja) | レーダ装置 | |
DE19531540A1 (de) | Verfahren zur Messung des Füllstandes einer Flüssigkeit | |
CN110187320A (zh) | 一种改进雷达信号时频分析方法 | |
CN108646553A (zh) | 一种统计在线监控闭环控制系统模型质量的方法 | |
CN104407197B (zh) | 一种基于三角函数迭代的信号相量测量的方法 | |
Corsaro et al. | Contrasting microscopic interactions determine the properties of water/methanol solutions | |
CN113252881A (zh) | 一种油品含水分析系统和信息融合分析方法 | |
EP3729074B1 (de) | Verfahren und vorrichtung zur milchfettmessung | |
CN105911349B (zh) | 基于重排时频谱的线性扫频信号基本参数估算方法及装置 | |
KR101983873B1 (ko) | 임피던스 매칭 방법 및 임피던스 매칭 시스템 | |
CN106778001A (zh) | 基于改进时频单源区的欠定混合矩阵盲估计方法 | |
Meyne et al. | Quasi-lumped coplanar transmission-line sensors for broadband liquid characterization | |
CN106644050B (zh) | 在频谱质量好的情况下测量物体的振动基频的方法 | |
DE3732856C2 (zh) | ||
Nordon et al. | Determination of the ethylene oxide content of polyether polyols by low-field 1H nuclear magnetic resonance spectrometry | |
JP2010271088A (ja) | 位相情報を用いた高分解能距離測定方法及び距離測定装置 | |
Nian et al. | Anticipated rife interpolation algorithm for frequency estimation of sinusoid signal | |
Graener | Anharmonicity and overtone spectra of oxygen-hydrogen stretching vibrations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |