CN107561357A - 一种基于信道化的高精度瞬时测频方法和装置 - Google Patents
一种基于信道化的高精度瞬时测频方法和装置 Download PDFInfo
- Publication number
- CN107561357A CN107561357A CN201710543287.6A CN201710543287A CN107561357A CN 107561357 A CN107561357 A CN 107561357A CN 201710543287 A CN201710543287 A CN 201710543287A CN 107561357 A CN107561357 A CN 107561357A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- frequency
- signal
- channelizing
- 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.)
- Granted
Links
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出了一种基于信道化的高精度瞬时测频方法和装置,包括:对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理;分别采用相位差法和插值法对数字信道化处理后的多路信号进行瞬时频率测量;根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率。本发明在数字信道化不过采时仍然具有较高的频率测量精度。
Description
技术领域
本发明涉及电子侦察中的雷达信号侦察和干扰技术领域,特别涉及一种基于信道化的高精度瞬时测频方法和装置。
背景技术
现代电子战需要截获敌方的电磁频谱信息,其中最重要的是测量敌方电磁信号的频率参数,根据此频率参数,可以得到频率捷变雷达的捷变频范围和频谱宽度等重要指标,它能反映雷达的功能和用途。在现代复杂的电磁环境下,为了实施有效的干扰,必须首先对信号进行分选和威胁识别,雷达的频率信息是进行分选和威胁识别的重要依据。
比较常见的测频方法有纯信道化测频和比相法瞬时测频。纯信道化测频具有灵敏度高的优点,但是受信道个数不能太多的限制,测频精度较低;比相法瞬时测频具有很高的测频精度,但是对信号信噪比要求较高,灵敏度较差。两者相结合的方法既能提高信号检测灵敏度,也能提高测频精度。即先使用数字信道化对信号进行信道化分,然后再使用相位差法测量信号瞬时频率。
但是当数字信道化不过采时,受噪声影响,相位差瞬时测频会存在相位测量模糊问题,极大地影响测频精度。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种基于信道化的高精度瞬时测频方法和装置。
为了实现上述目的,本发明一方面的实施例提供一种基于信道化的高精度瞬时测频方法,包括如下步骤:
步骤S1,对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理;
步骤S2,分别采用相位差法和插值法对步骤S1中数字信道化处理后的多路信号进行瞬时频率测量;
步骤S3,根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率。
进一步,在所述步骤S1中,采用偶排列、临界抽取的多相相带通滤波器组实现数字信道化处理。
进一步,在所述步骤S2中,
假设所述步骤S1数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
采用所述相位差法进行瞬时频率测量的计算步骤如下:
步骤1:对数字信道化之后第i个信道相邻两个采样点xi(n)与xi(n-1)进行共轭相乘:
yi(n)=xi(n)·conj(xi(n-1))=Ii(n)+jQi(n)
步骤2:利用正交信号Ii(n)与Qi(n)的关系得到相角φi(n),即
步骤3:通过相位与频率之间的转换运算得到每个信道中对应信号的瞬时频率:
其中,Δt是采样间隔,即当前信道两采样点之间的时间间隔,fi(n)是第i个信道相位差法得到的信号频率。
进一步,在所述步骤S2中,
假设所述步骤S1数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
采用所述插值法,进行瞬时频率测量的计算步骤如下:
步骤1:计算第i个信道以及第i-1个信道和第i+1个信道信号幅度:
其中,Ii(n)与Qi(n)是第i个信道第n个点的正交信号,Ai(n)是第i个信道第n个点的信号幅度;
步骤2:使用插值法根据信号幅度计算第i个信道信号频率:
其中,Ai(n)是第i个信道信号幅度,Aneighbour(n)为Ai-1(n)和Ai+1(n)中的较大信号幅度,Δt是采样间隔。当Ai-1(n)大于Ai+1(n)时,λ=-1;当Ai-1(n)小于Ai+1(n)时,λ=1。
进一步,在所述步骤S3中,
假设相位差法得到的频率为f1(n),插值法得到的频率为f2(n),最终解模糊之后的频率为f(n),解模糊过程如下:
如果则
如果则
如果则f(n)=f1(n),
其中,0≤α≤1,α比较合适的取值是0.75。
本发明实施例还提出一种基于信道化的高精度瞬时测频装置,包括:数字信道化模块、瞬时频率测量模块和频率解模糊模块,其中,
所述数字信道化模块用于对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理;
所述瞬时频率测量模块用于分别采用相位差法和插值法对步骤S1中数字信道化处理后的多路信号进行瞬时频率测量;
所述频率解模糊模块用于根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率。
进一步,所述数字信道化模块采用偶排列、临界抽取的多相相带通滤波器组实现数字信道化处理。
进一步,假设数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
所述瞬时频率测量模块采用相位差法进行瞬时频率测量的计算步骤如下:
步骤1:对数字信道化之后第i个信道相邻两个采样点xi(n)与xi(n-1)进行共轭相乘:
yi(n)=xi(n)·conj(xi(n-1))=Ii(n)+jQi(n)
步骤2:利用正交信号Ii(n)与Qi(n)的关系得到相角φi(n),即
步骤3:通过相位与频率之间的转换运算得到每个信道中对应信号的瞬时频率:
其中,Δt是采样间隔,即当前信道两采样点之间的时间间隔,fi(n)是第i个信道相位差法得到的信号频率。
进一步,假设数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
所述瞬时频率测量模块采用所述插值法进行瞬时频率测量的计算步骤如下:
步骤1:计算第i个信道以及第i-1个信道和第i+1个信道信号幅度:
其中,Ii(n)与Qi(n)是第i个信道第n个点的正交信号,Ai(n)是第i个信道第n个点的信号幅度;
步骤2:使用插值法根据信号幅度计算第i个信道信号频率:
其中,Ai(n)是第i个信道信号幅度,Aneighbour(n)为Ai-1(n)和Ai+1(n)中的较大信号幅度,Δt是采样间隔。当Ai-1(n)大于Ai+1(n)时,λ=-1;当Ai-1(n)小于Ai+1(n)时,λ=1。
进一步,所述频率解模糊模块假设相位差法得到的频率为f1(n),插值法得到的频率为f2(n),最终解模糊之后的频率为f(n),解模糊过程如下:
如果则
如果则
如果则
其中,0≤α≤1,α比较合适的取值是0.75。
根据本发明实施例的基于信道化的高精度瞬时测频方法和装置,其目的是为了能够在数字信道化不过采时仍然具有较高的频率测量精度。本发明对接收到的外部信号进行信道化处理,能够滤除带外噪声,提高信噪比,并提高信号检测灵敏度和测频精度;通过对信道化之后的信号进行基于相位差和插值法的瞬时频率测量,能够较准确地得到信号的频率信息;通过频率解模糊处理,即根据插值法得到的频率对相位差法测得的频率进行解模糊,可以确保得到无模糊的信号频率,极大地提高信号频率的测量精度。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的基于信道化的高精度瞬时测频方法的流程图;
图2为根据本发明实施例的基于信道化的高精度瞬时测频装置的结构图;
图3为根据本发明实施例的基于多相滤波结构的数字信道化原理图;
图4为根据本发明实施例的相位差法瞬时测频原理图;
图5为根据本发明实施例的频率解模糊示例图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
如图1所示,本发明实施例的基于信道化的高精度瞬时测频方法,包括如下步骤:
步骤S1,对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理,以滤除信道之外的噪声。
具体地,采用偶排列、临界抽取的多相相带通滤波器组实现数字信道化处理。对输入雷达信号使用一组性能相同的带通滤波器组,将原始全频带信号分割为多个子频带信号。数字信道化有一种高效的实现方法,即基于多相滤波器组的信道化结构。
如图3所示,图中所示为复信号的偶排列、临界抽取多相滤波结构,主要过程是先对输入信号进行D倍抽取,然后对抽取信号进行滤波,最后经过D点FFT输出D个信道的样本点。本实施例采用偶排列、临界抽取的多相滤波结构实现数字信道化。
步骤S2,分别采用相位差法和插值法对步骤S1中数字信道化处理后的多路信号进行瞬时频率测量。
(1)相位差法
假设步骤S1数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
如图4所示,采用相位差法进行瞬时频率测量的计算步骤如下:
步骤1:对数字信道化之后第i个信道相邻两个采样点xi(n)与xi(n-1)进行共轭相乘:
yi(n)=xi(n)·conj(xi(n-1))=Ii(n)+jQi(n)
步骤2:利用正交信号Ii(n)与Qi(n)的关系得到相角φi(n),即
步骤3:通过相位与频率之间的转换运算得到每个信道中对应信号的瞬时频率:
其中,Δt是采样间隔,即当前信道两采样点之间的时间间隔,fi(n)是第i个信道相位差法得到的信号频率。
(2)插值法
假设步骤S1数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
采用插值法,进行瞬时频率测量的计算步骤如下:
步骤1:计算第i个信道以及第i-1个信道和第i+1个信道信号幅度:
其中,Ii(n)与Qi(n)是第i个信道第n个点的正交信号,Ai(n)是第i个信道第n个点的信号幅度;
步骤2:使用插值法根据信号幅度计算第i个信道信号频率:
其中,Ai(n)是第i个信道信号幅度,Aneighbour(n)为Ai-1(n)和Ai+1(n)中的较大信号幅度,Δt是采样间隔。当Ai-1(n)大于Ai+1(n)时,λ=-1;当Ai-1(n)小于Ai+1(n)时,λ=1。
由于相位差法瞬时测频得到的频率精度较高,但是存在相位模糊现象,导致测频模糊;插值法测频得到的频率精度较低,但是不会存在测频模糊现象。
步骤S3,根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率,确保频率测量的准确性。
具体地,假设相位差法得到的频率为f1(n),插值法得到的频率为f2(n),最终解模糊之后的频率为f(n),解模糊过程如下:
如果则
如果则
如果则f(n)=f1(n),
其中,0≤α≤1,α比较合适的取值是0.75。
经过解模糊处理之后,原始可能存在模糊的频率均被校正回正确的频率,从而使得测量的频率更准确。图5显示的是采用频率解模糊算法前后的测频效果图。图5中,信道带宽为18.75MHz,即频率范围是-9.375MHz到+9.375之间。
采用相位差法瞬时测频时,测频模糊现象很严重,即-9MHz附近的频率被测成了+9MHz左右的频率,经过解模糊处理之后,+9MHz附近的频率被校正回正确的频率。
如图2所示,本发明实施例还提出一种基于信道化的高精度瞬时测频装置,包括:数字信道化模块100、瞬时频率测量模块200和频率解模糊模块300。
具体地,数字信道化模块100用于对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理。
在本发明的一个实施例中,数字信道化模块100采用偶排列、临界抽取的多相相带通滤波器组实现数字信道化处理。对输入雷达信号使用一组性能相同的带通滤波器组,将原始全频带信号分割为多个子频带信号。数字信道化有一种高效的实现方法,即基于多相滤波器组的信道化结构。
如图3所示,图中所示为复信号的偶排列、临界抽取多相滤波结构,主要过程是先对输入信号进行D倍抽取,然后对抽取信号进行滤波,最后经过D点FFT输出D个信道的样本点。本实施例采用偶排列、临界抽取的多相滤波结构实现数字信道化。
瞬时频率测量模块200用于分别采用相位差法和插值法对步骤S1中数字信道化处理后的多路信号进行瞬时频率测量。
(1)相位差法
假设数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
瞬时频率测量模块200采用相位差法进行瞬时频率测量的计算步骤如下:
步骤1:对数字信道化之后第i个信道相邻两个采样点xi(n)与xi(n-1)进行共轭相乘:
yi(n)=xi(n)·conj(xi(n-1))=Ii(n)+jQi(n)
步骤2:利用正交信号Ii(n)与Qi(n)的关系得到相角φi(n),即
步骤3:通过相位与频率之间的转换运算得到每个信道中对应信号的瞬时频率:
其中,Δt是采样间隔,即当前信道两采样点之间的时间间隔,fi(n)是第i个信道相位差法得到的信号频率。
(2)插值法
假设数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
瞬时频率测量模块200采用插值法进行瞬时频率测量的计算步骤如下:
步骤1:计算第i个信道以及第i-1个信道和第i+1个信道信号幅度:
其中,Ii(n)与Qi(n)是第i个信道第n个点的正交信号,Ai(n)是第i个信道第n个点的信号幅度;
步骤2:使用插值法根据信号幅度计算第i个信道信号频率:
其中,Ai(n)是第i个信道信号幅度,Aneighbour(n)为Ai-1(n)和Ai+1(n)中的较大信号幅度,Δt是采样间隔。当Ai-1(n)大于Ai+1(n)时,λ=-1;当Ai-1(n)小于Ai+1(n)时,λ=1。
频率解模糊模块300用于根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率,确保频率测量的准确性。
具体地,频率解模糊模块300假设相位差法得到的频率为f1(n),插值法得到的频率为f2(n),最终解模糊之后的频率为f(n),解模糊过程如下:
如果则
如果则
如果则f(n)=f1(n),
其中,0≤α≤1,α比较合适的取值是0.75。
经过解模糊处理之后,原始可能存在模糊的频率均被校正回正确的频率,从而使得测量的频率更准确。图5显示的是采用频率解模糊算法前后的测频效果图。图5中,信道带宽为18.75MHz,即频率范围是-9.375MHz到+9.375之间。
采用相位差法瞬时测频时,测频模糊现象很严重,即-9MHz附近的频率被测成了+9MHz左右的频率,经过解模糊处理之后,+9MHz附近的频率被校正回正确的频率。
根据本发明实施例的基于信道化的高精度瞬时测频方法和装置,其目的是为了能够在数字信道化不过采时仍然具有较高的频率测量精度。本发明对接收到的外部信号进行信道化处理,能够滤除带外噪声,提高信噪比,并提高信号检测灵敏度和测频精度;通过对信道化之后的信号进行基于相位差和插值法的瞬时频率测量,能够较准确地得到信号的频率信息;通过频率解模糊处理,即根据插值法得到的频率对相位差法测得的频率进行解模糊,可以确保得到无模糊的信号频率,极大地提高信号频率的测量精度。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (10)
1.一种基于信道化的高精度瞬时测频方法,其特征在于,包括如下步骤:
步骤S1,对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理;
步骤S2,分别采用相位差法和插值法对步骤S1中数字信道化处理后的多路信号进行瞬时频率测量;
步骤S3,根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率。
2.如权利要求1所述的基于信道化的高精度瞬时测频方法,其特征在于,在所述步骤S1中,采用偶排列、临界抽取的多相相带通滤波器组实现数字信道化处理。
3.如权利要求1所述的基于信道化的高精度瞬时测频方法,其特征在于,在所述步骤S2中,
假设所述步骤S1数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
采用所述相位差法进行瞬时频率测量的计算步骤如下:
步骤1:对数字信道化之后第i个信道相邻两个采样点xi(n)与xi(n-1)进行共轭相乘:
yi(n)=xi(n)·conj(xi(n-1))=Ii(n)+jQi(n)
步骤2:利用正交信号Ii(n)与Qi(n)的关系得到相角φi(n),即
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>a</mi>
<mi>r</mi>
<mi>c</mi>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
步骤3:通过相位与频率之间的转换运算得到每个信道中对应信号的瞬时频率:
<mrow>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,Δt是采样间隔,即当前信道两采样点之间的时间间隔,fi(n)是第i个信道相位差法得到的信号频率。
4.如权利要求1所述的基于信道化的高精度瞬时测频方法,其特征在于,在所述步骤S2中,
假设所述步骤S1数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
采用所述插值法,进行瞬时频率测量的计算步骤如下:
步骤1:计算第i个信道以及第i-1个信道和第i+1个信道信号幅度:
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
其中,Ii(n)与Qi(n)是第i个信道第n个点的正交信号,Ai(n)是第i个信道第n个点的信号幅度;
步骤2:使用插值法根据信号幅度计算第i个信道信号频率:
<mrow>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>i</mi>
<mi>g</mi>
<mi>h</mi>
<mi>b</mi>
<mi>o</mi>
<mi>u</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>i</mi>
<mi>g</mi>
<mi>h</mi>
<mi>b</mi>
<mi>o</mi>
<mi>u</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
<mo>)</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,Ai(n)是第i个信道信号幅度,Aneighbour(n)为Ai-1(n)和Ai+1(n)中的较大信号幅度,Δt是采样间隔。当Ai-1(n)大于Ai+1(n)时,λ=-1;当Ai-1(n)小于Ai+1(n)时,λ=1。
5.如权利要求1所述的基于信道化的高精度瞬时测频方法,其特征在于,在所述步骤S3中,
假设相位差法得到的频率为f1(n),插值法得到的频率为f2(n),最终解模糊之后的频率为f(n),解模糊过程如下:
如果则
如果则
如果则f(n)=f1(n),
其中,0≤α≤1,α比较合适的取值是0.75。
6.一种基于信道化的高精度瞬时测频装置,其特征在于,包括:数字信道化模块、瞬时频率测量模块和频率解模糊模块,其中,
所述数字信道化模块用于对输入的雷达信号进行数字信道化处理,包括:采用多相带通滤波器组对输入的雷达信号进行抽取和滤波,并进行傅里叶FFT变换后输出多个信道的样本点,以实现数字信道化处理;
所述瞬时频率测量模块用于分别采用相位差法和插值法对步骤S1中数字信道化处理后的多路信号进行瞬时频率测量;
所述频率解模糊模块用于根据插值法得到的频率对相位差法得到的频率进行解模糊处理,以将原始存在模糊的频率校正回正确的频率。
7.如权利要求6所述的基于信道化的高精度瞬时测频装置,其特征在于,所述数字信道化模块采用偶排列、临界抽取的多相相带通滤波器组实现数字信道化处理。
8.如权利要求6所述的基于信道化的高精度瞬时测频装置,其特征在于,
假设数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
所述瞬时频率测量模块采用相位差法进行瞬时频率测量的计算步骤如下:
步骤1:对数字信道化之后第i个信道相邻两个采样点xi(n)与xi(n-1)进行共轭相乘:
yi(n)=xi(n)·conj(xi(n-1))=Ii(n)+jQi(n)
步骤2:利用正交信号Ii(n)与Qi(n)的关系得到相角φi(n),即
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>a</mi>
<mi>r</mi>
<mi>c</mi>
<mi>t</mi>
<mi>a</mi>
<mi>n</mi>
<mfrac>
<mrow>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
步骤3:通过相位与频率之间的转换运算得到每个信道中对应信号的瞬时频率:
<mrow>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&phi;</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,Δt是采样间隔,即当前信道两采样点之间的时间间隔,fi(n)是第i个信道相位差法得到的信号频率。
9.如权利要求6所述的基于信道化的高精度瞬时测频装置,其特征在于,
假设数字信道化之后的信号为基带单一频率信号,其表达式为:
其中,是真实信号,A是信号幅度,φ是初始相位,f是信号频率,Δt是采样时间间隔,ω(n)是零均值高斯白噪声;
所述瞬时频率测量模块采用所述插值法进行瞬时频率测量的计算步骤如下:
步骤1:计算第i个信道以及第i-1个信道和第i+1个信道信号幅度:
<mrow>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>I</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
其中,Ii(n)与Qi(n)是第i个信道第n个点的正交信号,Ai(n)是第i个信道第n个点的信号幅度;
步骤2:使用插值法根据信号幅度计算第i个信道信号频率:
<mrow>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>i</mi>
<mi>g</mi>
<mi>h</mi>
<mi>b</mi>
<mi>o</mi>
<mi>u</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>n</mi>
<mi>e</mi>
<mi>i</mi>
<mi>g</mi>
<mi>h</mi>
<mi>b</mi>
<mi>o</mi>
<mi>u</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
<mo>)</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中,Ai(n)是第i个信道信号幅度,Aneighbour(n)为Ai-1(n)和Ai+1(n)中的较大信号幅度,Δt是采样间隔。当Ai-1(n)大于Ai+1(n)时,λ=-1;当Ai-1(n)小于Ai+1(n)时,λ=1。
10.如权利要求6所述的基于信道化的高精度瞬时测频装置,其特征在于,
所述频率解模糊模块假设相位差法得到的频率为f1(n),插值法得到的频率为f2(n),最终解模糊之后的频率为f(n),解模糊过程如下:
如果则
如果则
如果则f(n)=f1(n),
其中,0≤α≤1,α比较合适的取值是0.75。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710543287.6A CN107561357B (zh) | 2017-07-05 | 2017-07-05 | 一种基于信道化的高精度瞬时测频方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710543287.6A CN107561357B (zh) | 2017-07-05 | 2017-07-05 | 一种基于信道化的高精度瞬时测频方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107561357A true CN107561357A (zh) | 2018-01-09 |
CN107561357B CN107561357B (zh) | 2020-08-14 |
Family
ID=60973462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710543287.6A Active CN107561357B (zh) | 2017-07-05 | 2017-07-05 | 一种基于信道化的高精度瞬时测频方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107561357B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061296A (zh) * | 2018-07-17 | 2018-12-21 | 南京恒电电子有限公司 | 一种射频脉冲信号的高精度载频估计方法 |
CN109270344A (zh) * | 2018-10-07 | 2019-01-25 | 扬州大学 | 脉冲丢失下的相参脉冲信号频率估计方法 |
CN109387813A (zh) * | 2018-08-31 | 2019-02-26 | 中国电子科技集团公司第二十九研究所 | 一种基于宽带雷达信号接收的高精度频率测量方法 |
CN109617631A (zh) * | 2018-12-28 | 2019-04-12 | 华航高科(北京)技术有限公司 | 基于数字信道化瞬时参数测量的侦察***自适应接收方法 |
CN113219424A (zh) * | 2021-04-29 | 2021-08-06 | 中国船舶重工集团公司第七二三研究所 | 一种并行输出的多通道瞬时测频*** |
CN115694529A (zh) * | 2023-01-05 | 2023-02-03 | 扬州宇安电子科技有限公司 | 一种基于数字信道化接收机的测频精度优化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5235287A (en) * | 1991-07-05 | 1993-08-10 | Sanderson Richard B | Frequency measurement receiver with bandwidth improvement through phase shifted sampling of real signals using sampling rate selection |
CN202334540U (zh) * | 2011-11-17 | 2012-07-11 | 四川蓝讯宝迩电子科技有限公司 | 一种用于数字信道化接收机的信道化测量*** |
CN104459319A (zh) * | 2014-11-27 | 2015-03-25 | 中国船舶重工集团公司第七二四研究所 | 一种短基线干涉仪矢量叠加鉴相方法 |
CN104735007A (zh) * | 2015-03-12 | 2015-06-24 | 电子科技大学 | 一种对中心频率无约束的直接中频信道化方法 |
CN105486920A (zh) * | 2015-11-27 | 2016-04-13 | 北京振兴计量测试研究所 | 一种基于数字信道化的窄脉冲频率测量方法 |
CN105510706A (zh) * | 2015-12-30 | 2016-04-20 | 中国航天时代电子公司 | 一种高精度欠采样测频方法 |
-
2017
- 2017-07-05 CN CN201710543287.6A patent/CN107561357B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5235287A (en) * | 1991-07-05 | 1993-08-10 | Sanderson Richard B | Frequency measurement receiver with bandwidth improvement through phase shifted sampling of real signals using sampling rate selection |
CN202334540U (zh) * | 2011-11-17 | 2012-07-11 | 四川蓝讯宝迩电子科技有限公司 | 一种用于数字信道化接收机的信道化测量*** |
CN104459319A (zh) * | 2014-11-27 | 2015-03-25 | 中国船舶重工集团公司第七二四研究所 | 一种短基线干涉仪矢量叠加鉴相方法 |
CN104735007A (zh) * | 2015-03-12 | 2015-06-24 | 电子科技大学 | 一种对中心频率无约束的直接中频信道化方法 |
CN105486920A (zh) * | 2015-11-27 | 2016-04-13 | 北京振兴计量测试研究所 | 一种基于数字信道化的窄脉冲频率测量方法 |
CN105510706A (zh) * | 2015-12-30 | 2016-04-20 | 中国航天时代电子公司 | 一种高精度欠采样测频方法 |
Non-Patent Citations (2)
Title |
---|
周志军: "电子侦察仿真***的建模与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
曾德国 等: "基于相位差分的脉内调制信号类型识别", 《电子测量与仪器学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109061296A (zh) * | 2018-07-17 | 2018-12-21 | 南京恒电电子有限公司 | 一种射频脉冲信号的高精度载频估计方法 |
CN109061296B (zh) * | 2018-07-17 | 2020-11-27 | 南京恒电电子有限公司 | 一种射频脉冲信号的高精度载频估计方法 |
CN109387813A (zh) * | 2018-08-31 | 2019-02-26 | 中国电子科技集团公司第二十九研究所 | 一种基于宽带雷达信号接收的高精度频率测量方法 |
CN109270344A (zh) * | 2018-10-07 | 2019-01-25 | 扬州大学 | 脉冲丢失下的相参脉冲信号频率估计方法 |
CN109270344B (zh) * | 2018-10-07 | 2021-01-08 | 扬州大学 | 脉冲丢失下的相参脉冲信号频率估计方法 |
CN109617631A (zh) * | 2018-12-28 | 2019-04-12 | 华航高科(北京)技术有限公司 | 基于数字信道化瞬时参数测量的侦察***自适应接收方法 |
CN109617631B (zh) * | 2018-12-28 | 2021-09-14 | 华航高科(北京)技术有限公司 | 基于数字信道化瞬时参数测量的侦察***自适应接收方法 |
CN113219424A (zh) * | 2021-04-29 | 2021-08-06 | 中国船舶重工集团公司第七二三研究所 | 一种并行输出的多通道瞬时测频*** |
CN115694529A (zh) * | 2023-01-05 | 2023-02-03 | 扬州宇安电子科技有限公司 | 一种基于数字信道化接收机的测频精度优化方法 |
CN115694529B (zh) * | 2023-01-05 | 2023-08-11 | 扬州宇安电子科技有限公司 | 一种基于数字信道化接收机的测频精度优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107561357B (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107561357A (zh) | 一种基于信道化的高精度瞬时测频方法和装置 | |
CN111693944B (zh) | 雷达有源干扰信号参数提取与干扰样式识别方法及装置 | |
CN109459732B (zh) | 一种宽带雷达干扰信号的识别分类方法 | |
CN103454495B (zh) | 自适应高精度快速频谱分析方法 | |
CN101135726B (zh) | 一种星载sar内定标信号处理平台***及实现方法 | |
CN102043091B (zh) | 数字化高精度相位检测器 | |
CN103929387B (zh) | 基于fpga的大动态范围数字信道化接收机及工作方法 | |
CN106385288B (zh) | 一种雷达脉冲跨信道融合检测装置及其检测方法 | |
CN105067650A (zh) | 一种利用小波计算导数检测谱特征峰的方法 | |
CN104375132A (zh) | 数字阵列雷达多路模拟通道相对延时测量装置及方法 | |
CN106680585B (zh) | 谐波/间谐波的检测方法 | |
CN107528646B (zh) | 一种基于宽带频谱的干扰信号识别及提取方法 | |
CN105629224A (zh) | 调频连续波雷达高精度测距方法 | |
CN110868432B (zh) | 一种射频指纹提取方法 | |
CN110208601B (zh) | 基于fpga的瞬时测频方法及其数字接收机 | |
CN107272026B (zh) | 一种导航信号分量相位偏差测试方法 | |
CN104991171B (zh) | 基于超高频信号的gis局部放电分频故障图谱绘制方法 | |
CN107132511B (zh) | 一种精确的雷达线性调频源预失真方法 | |
CN106353776A (zh) | 一种非理想信道下测量零值无偏的频域抗干扰***及其方法 | |
CN105807204B (zh) | 基于频谱细化的硬件木马检测方法 | |
CN106198008B (zh) | 一种电机轴承外圈故障辨识方法 | |
CN111030959A (zh) | 一种nbiot的频域时频同步方法 | |
CN106353594A (zh) | 一种快速多分辨率频谱分析***及方法 | |
CN101527584A (zh) | 基于频域幅值处理的扩频伪码检测方法 | |
CN107329123A (zh) | 一种弱雷达脉冲包络信号检测方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |