CN103441975B - 一种基于功率谱的二相编码信号参数估值方法 - Google Patents
一种基于功率谱的二相编码信号参数估值方法 Download PDFInfo
- Publication number
- CN103441975B CN103441975B CN201310390121.7A CN201310390121A CN103441975B CN 103441975 B CN103441975 B CN 103441975B CN 201310390121 A CN201310390121 A CN 201310390121A CN 103441975 B CN103441975 B CN 103441975B
- Authority
- CN
- China
- Prior art keywords
- tau
- signal
- power spectrum
- bandwidth
- pulse width
- 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.)
- Expired - Fee Related
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000005070 sampling Methods 0.000 claims abstract description 6
- 238000010606 normalization Methods 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims description 8
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 4
- GNFTZDOKVXKIBK-UHFFFAOYSA-N 3-(2-methoxyethoxy)benzohydrazide Chemical compound COCCOC1=CC=CC(C(=O)NN)=C1 GNFTZDOKVXKIBK-UHFFFAOYSA-N 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 125000004122 cyclic group Chemical group 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于功率谱的二相编码信号参数估值方法,该方法包括:求出二相编码信号不同码长时对应的子脉冲宽度τ和3dB带宽B3dB的拟合关系;检测待估值的二相编码信号的功率谱,对所述功率谱进行幅值归一化得到P(f),对所述P(f)均匀采样M个点,变成长度为M的序列P(m),对所述P(m)做离散余弦变换得到Ypre(P(m)),对所述Ypre(P(m))进行阈值处理得到Y(P(m));求所述Y(P(m))的上包络,计算所述上包络的峰值个数n,根据信号码长N与所述峰值个数n的关系求出所述信号码长N;对所述Y(P(m))进行离散余弦逆变换,得到Pi(m),计算所述Pi(m)的载频fc、3dB带宽B3dB和子脉冲宽度τ。本发明克服现有技术中低信噪比环境下二相编码信号估值结果受噪声影响大而且计算量大的缺陷,实现未知先验信息条件下的二相编码信号参数的估值。
Description
技术领域
本发明涉及一种基于功率谱离散余弦变换的二相编码信号参数估值方法。
背景技术
相位编码信号,其相位调制函数是离散的有限的状态,属于离散编码脉冲压缩信号,由于相位编码采用伪随机序列,故又称为伪随机编码信号。相位编码信号常用于多普勒变化较小的场合。伪随机编码信号可以按相位取值数目分类,如果相位只取两个值,称为二相编码信号。
二相编码信号是指信号相位调制函数是离散的两个值,相位编码一般采用伪随机序列,时域表达式为:
其中,τ是子脉冲宽度,N是码长,子脉冲函数
二相编码信号功率谱表达式为:
二相编码信号功率谱特点如下:
①功率谱关于中心频率对称;
②带宽内呈现多峰、码长和带宽内的峰或谷个数相等。
由于二相编码信号的参数估计问题较一般的信号复杂,因此成为研究的难点。二相编码信号常用的参数如下所示:
①载频fc;
②3dB带宽B3dB;
③子脉冲宽度τ,也就是时宽。
现有的二相编码信号参数估值方法大都基于时域的信号分析,常用的参数估计方法有时域相位差分法,时频分析和循环谱法等。2011年3月9日公开的专利CN101984613A中提出了一种低速率BPSK突发信号码速率估计方法,利用该方法可获得二相编码信号的码速率。2011年,崔伟亮等人发表在《电子与信息学报》的文献“改进的循环谱估计快速算法与性能分析”。2012年,徐会法等人发表在《兵工学报》的文献“基于FRFT的几种典型相位编码信号检测与参数估计”。基于时频分析和循环谱的二相编码信号参数估计的计算量很大,无法解决估计精度与运算量之间的矛盾,而且不适用于低信噪比环境。
发明内容
本发明目的在于:实现未知先验信息条件下的二相编码信号参数的估值,克服现有技术中低信噪比环境下二相编码信号估值结果受噪声影响大而且计算量大的缺陷。
本发明的技术方案是:本发明提供了一种基于功率谱离散余弦变换的二相编码信号参数估值方法,其特征在于,该方法包括:
步骤1、求出二相编码信号不同码长N时对应的子脉冲宽度τ和3dB带宽B3dB的拟合关系,其具体包括:
步骤a、设二相编码信号码长N取值范围为{N1,N2,…,Nl,…,Nr},l∈{1,2,…,r},设定l=1;
步骤b、设定所述二相编码信号码长为Nl,子脉冲宽度τ变化区间为[τ1,τ2]、变化步长为Δτ;
步骤c、在所述τ变化区间内,计算不同τ值对应的所述二相编码信号功率谱,设所述功率谱最大值为Pmax,搜索所述功率谱第一个值为Pmax/2处对应的频率f1和最后一个值为Pmax/2处对应的频率f2,根据定义B3dB=f2-f1计算出每个所述功率谱对应的3dB带宽B3dB,则获得一一对应的τ-B3dB数据点;
步骤d、选择指数函数对所述τ-B3dB数据点进行数据拟合;
步骤e、求出拟合函数分别对应的系数a、b,得到所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式,其中,拟合函数指的是步骤d中数据拟合得到的拟合函数;
步骤f、令l加1;
步骤g,判断l是否等于r,如果l等于r,执行步骤h,如果l不等于r,返回步骤b;
步骤h,得出不同码长N分别对应的所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式。
步骤2、检测待估值的二相编码信号的功率谱,对所述功率谱进行幅值归一化得到P(f),对所述P(f)均匀采样M个点,变成长度为M的序列P(m),对所述P(m)做离散余弦变换得到Ypre(P(m)),对所述Ypre(P(m))进行阈值处理得到Y(P(m));
步骤3、求所述Y(P(m))的上包络,计算所述上包络的峰值个数n,根据信号码长N与所述峰值个数n的关系求出所述信号码长N;
步骤4、对所述Y(P(m))进行离散余弦逆变换,得到Pi(m),计算所述Pi(m)的载频fc、3dB带宽B3dB和子脉冲宽度τ,其具体包括:
a、设所述Pi(m)的最大值为Pmax’,搜索所述Pi(m)第一个值为Pmax’/2处对应的频率f’1和最后一个值为Pmax’/2处对应的频率f2’,根据定义计算出所述Pi(m)的载频fc;
b、根据定义B3dB=f2’-f1’计算出所述Pi(m)的3dB带宽B3dB;
c、根据所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式求出所述子脉冲宽度τ。
本发明的有益效果是:本发明公开了一种基于功率谱离散余弦变换的二相编码信号参数估值方法,该方法包括:求出二相编码信号不同码长时对应的子脉冲宽度τ和3dB带宽B3dB的拟合关系;检测待估值的二相编码信号的功率谱,对所述功率谱进行幅值归一化得到P(f),对所述P(f)均匀采样M个点,变成长度为M的序列P(m),对所述P(m)做离散余弦变换得到Ypre(P(m)),对所述Ypre(P(m))进行阈值处理得到Y(P(m));求所述Y(P(m))的上包络,计算所述上包络的峰值个数n,根据信号码长N与所述峰值个数n的关系求出所述信号码长N;对所述Y(P(m))进行离散余弦逆变换,得到Pi(m),计算所述Pi(m)的载频fc、3dB带宽B3dB和子脉冲宽度τ。本发明方法的计算量主要为一次FFT变换的复数乘法量和两次实数DCT变换的计算量,克服了目前常用的时频分析方法计算量大的缺陷,在低信噪比环境下,实现未知先验信息条件的二相编码信号参数的准确估值。
附图说明
图1本发明公开的一种基于功率谱的二相编码信号参数估值方法流程图;
图2本发明公开的确定二相编码信号子脉冲宽度τ和3dB带宽B3dB关系流程图;
图3本发明公开的选取阈值流程图;
图4本发明公开的二相编码信号码长N={5,7,11,13}时子脉冲宽度τ和3dB带宽B3dB关系的拟合曲线;
其中,左上为码长N=5时;右上为码长N=7时;左下为码长N=11时;右下为码长N=13时;
图5本发明实施例中不同阈值时降噪后信号均方根误差曲线图;
其中,·代表阈值T=0.01·max(Y(m))时;*代表阈值T=0.02·max(Y(m))时;△代表阈值T=0.03·max(Y(m))时;○代表阈值T=0.04·max(Y(m))时;虚线代表不取阈值时;
图6本发明实施例中不同参数时降噪后信号的均方根误差曲线图;
其中,代表参数N=5,τ=0.01μs时;代表参数N=5,τ=0.1μs时;○代表参数N=7,τ=0.01μs时;代表参数N=7,τ=0.1μs时;□代表参数N=11,τ=0.01μs时;△代表参数N=11,τ=0.1μs时;*代表参数N=13,τ=0.01μs时;·代表参数N=13,τ=0.1μs时;
图7本发明实施例中二相编码信号功率谱经离散余弦变换后的信号图;
图8本发明实施例中二相编码信号功率谱经离散余弦变换后信号的包络图;
图9本发明实施例中二相编码信号码长估计的准确率图;
其中,*代表码长N=5时;△代表码长N=7时;◇代表码长N=11时;○代表码长N=13时。
图10本发明实施例中二相编码信号经实验处理前后功率谱图与不含噪声功率谱图;
其中,上图为不含噪声时的功率谱图;中图为SNR=0dB时的功率谱图;下图为上图经过本发明方法处理后的功率谱图。
图11本发明实施例中二相编码信号子脉冲宽度估值的均方根误差图;
其中,*代表脉宽τ=0.01μs时;·代表τ=0.05μs时;○代表脉宽τ=0.1μs时;
图12本发明实施例中二相编码信号载频估值的均方根误差图。
其中,*代表脉宽τ=0.01μs时;·代表τ=0.05μs时;○代表脉宽τ=0.1μs时。
具体实施方式
以下将参照图1-12对本发明的实施方式进行说明。
如图1所示,本发明实施例进行基于功率谱的二相编码信号参数估值方法包括下列步骤:
步骤1、求出二相编码信号不同码长N时对应的子脉冲宽度τ和3dB带宽B3dB的拟合关系,如图2所示,其具体包括:
步骤a、设二相编码信号码长N取值范围为{N1,N2,…,Nl,…,Nr},l∈{1,2,…,r},设定l=1;
步骤b、设定所述二相编码信号码长为Nl,子脉冲宽度τ变化区间为[τ1,τ2]、变化步长为Δτ;
步骤c、在所述τ变化区间内,计算不同τ值对应的所述二相编码信号功率谱,设所述功率谱最大值为Pmax,搜索所述功率谱第一个值为Pmax/2处对应的频率f1和最后一个值为Pmax/2处对应的频率f2,根据定义B3dB=f2-f1计算出每个所述功率谱对应的3dB带宽B3dB,则获得一一对应的τ-B3dB数据点;
步骤d、选择指数函数对所述τ-B3dB数据点进行数据拟合;
步骤e、求出拟合函数分别对应的系数a、b,得到所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式;
步骤f、令l加1;
步骤g,判断l是否等于r,如果l等于r,执行步骤h,如果l不等于r,返回步骤b;
步骤h,得出不同码长N分别对应的所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式。
步骤2、检测待估值的二相编码信号的功率谱,对所述功率谱进行幅值归一化得到P(f),对所述P(f)均匀采样M个点,变成长度为M的序列P(m),对所述P(m)做离散余弦变换得到Ypre(P(m)),对所述Ypre(P(m))进行阈值处理得到Y(P(m));
其中,步骤2中,通过以下步骤获得所述Y(P(m)):
在二相编码信号功率谱表达式中,令τ为子脉冲宽度,N为码长,fc为载频,二相编码信号功率谱表达式为:
对P(f)均匀采样M个点,变成长度为M的序列P(m),其中所述离散余弦变换式中
所述离散余弦变换式为:
而实际收到的信号是带有噪声的,经过离散余弦变换后的信号也会受到噪声影响,因此设定一个阈值T,将所述Ypre(P(m))值为(-T,T)的部分置零就得到所述Y(P(m))。
需要说明的是,如图3所示,阈值的选取步骤如下:
a、对所述Y(P(m))做离散余弦逆变换得到信号Pi(m),设定所述Pi(m)与模板信号Pnn(m)的均方根误差阈值为σr;
b、设定所述阈值T的粗范围(0,Δ1),在所述粗范围内以10-nΔ1为步长,其中n=1;
c、分别计算所述均方根误差σi,i=1,2,3,…;
d、选取出min{σi}所对应的阈值Ta;
e、判断所述min{σi}是否小于所述σr,是,则执行步骤h,否,则执行步骤f;
f、重新设置阈值范围(Ta-10-nΔ1,Ta+10-nΔ1),步长为10-(n+1)Δ1;
g、令n加1,返回执行步骤c;
h、选取Ta为最优阈值。
其中,在步骤a中,本发明专利中所述Pi(m)是通过离散余弦逆变换计算得到的,其表达式如下:
需要说明的是,本发明专利中所述模板信号Pnn(m)定义为不含噪声时的信号;
需要说明的是,本发明专利中所述Pi(m)与所述模板信号Pnn(m)的均方根误差定义如下:
均方根误差即所述Pi(m)与所述模板信号Pnn(m)之间的方差的平方根,其定义式为:
步骤3、求所述Y(P(m))的上包络,计算所述上包络的峰值个数n,根据信号码长N与所述峰值个数n的关系求出所述信号码长N;
其中,步骤3中,信号码长N与所述Y(P(m))上包络峰值个数n的关系为:
N=2n+1
需要说明的是,本发明中上包络的定义如下:
将所述Y(P(m))中的峰值点连成一条曲线,即为所述Y(P(m))的上包络;
需要说明的是,上包络的峰值个数判断方法如下:
对所述上包络的每个点进行判断,当被判断的点同时大于它左侧5个连续的点和它右侧5个连续的点时,便可判定这个点是峰值点;
步骤4、对所述Y(P(m))进行离散余弦逆变换,得到Pi(m),计算所述Pi(m)的载频fc、3dB带宽B3dB和子脉冲宽度τ,其具体包括:
步骤4、对所述Y(P(m))进行离散余弦逆变换,得到Pi(m),计算所述Pi(m)的载频fc、3dB带宽B3dB和子脉冲宽度τ,其具体包括:
a、设所述Pi(m)的最大值为Pmax’,搜索所述Pi(m)第一个值为Pmax’/2处对应的频率f’1和最后一个值为Pmax’/2处对应的频率f2’,根据定义计算出所述Pi(m)的载频fc;
b、根据定义B3dB=f2’-f1’计算出所述Pi(m)的3dB带宽B3dB;
c、根据所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式求出所述子脉冲宽度τ。
本发明以二相编码信号为实施例对本发明进行详细说明;本发明具体实施例用缩写BPSK表示二相编码信号,用缩写DCT表示离散余弦变换,用缩写IDCT表示离散余弦逆变换,用缩写SNR表示信噪比。仿真参数设置:码长N={5,7,11,13},子脉冲宽度τ=0.05μs,载频为500MHz。
本发明实施例可以应用到军用电子对抗以及民用频谱监督领域,在低信噪比环境下准确地估计出信号参数,可以反映雷达的功能和用途,提高声纳跟踪***的定位精度与改善跟踪效果具有重要意义。
一、计算不同码长对应的τ-B3dB拟合关系式;
拟合关系曲线如图4所示,且拟合优度均大于99.5%,拟合公式为:
①5位BPSK信号τ-B3dB拟合公式:
τ=0.18604·B3dB -0.9922
②7位BPSK信号τ-B3dB拟合公式:
τ=0.82959·B3dB -0.99148
③11位BPSK信号τ-B3dB拟合公式:
τ=0.86469·B3dB -0.98822
④13位BPSK信号τ-B3dB拟合公式:
τ=0.21173·B3dB -0.98466
二、阈值T的选取;
设DCT后的信号最大值为max{Y(P(m}),设定所述σr=0.1,Δ1=0.1·max{Y(P(m))},根据阈值选取原则,当阈值取为0.02·max{Y(P(m))}时即可满足阈值选取条件,设定BPSK信号参数为码长N=13,SNR在-5dB~20dB变化时,所述均方根误差曲线如图5所示,当参数变化时,所述均方根误差曲线如图6所示。
由图5-6可知,当采取阈值处理时可以使信号具有降噪效果,而且当阈值取为0.02·max{Y(P(m))}时,所述均方根误差均小于0.1,因此本发明专利采取阈值T=0.02·max{Y(P(m))}。
三、检测待估值的BPSK信号功率谱P(f),并对所述BPSK信号进行码长估计;
对所述BPSK信号功率谱进行DCT变换,变换后信号如图7所示,并做阈值处理得到Y(P(m)),求所述经过阈值处理后信号的上包络,如图8所示,求出峰值个数n,于是得到所述BPSK信号的码长N的估计,当SNR变化范围为[-5dB~20dB]时,进行100次Monte Carlo试验。
如图9所示,当SNR大于-5dB时,BPSK信号码长估计准确率均为100%,因此可以进一步实现对所述BPSK信号的子脉冲宽度和载频的准确估值。
四、对所述BPSK信号进行载频和子脉冲宽度的估值;
对所述Y(P(m))进行IDCT变换,得到Pi(m),设定SNR=0dB,如图10所示,经过处理后的信号更接近于不含噪声的信号。
对SNR在区间为[-5dB,20dB]时的BPSK信号进行100次MonteCarlo试验,计算所述BPSK信号载频和子脉冲宽度的均方根误差(RMSE),结果如图11和图12所示。
本发明方法在SNR为[-5dB,20dB]时载频的估计误差如图11所示,而且在SNR>-5dB时载频的估计误差均小于0.5MHz。
本发明方法在SNR为[-5dB,20dB]时子脉冲宽度的估计误差如图12所示,而且在SNR>-5dB时子脉冲宽度的估计误差均小于0.032μs。
Claims (1)
1.一种基于功率谱的二相编码信号参数估值方法,其特征在于,该方法包括:
步骤1、求出二相编码信号不同码长N时对应的子脉冲宽度τ和3dB带宽B3dB的拟合关系,其具体包括:
步骤a、设二相编码信号码长N取值范围为{N1,N2,…,Nl,…,Nr},l∈{1,2,…,r},设定l=1;
步骤b、设定所述二相编码信号码长为Nl,子脉冲宽度τ变化区间为[τ1,τ2]、变化步长为Δτ;
步骤c、在所述τ变化区间内,计算不同τ值对应的所述二相编码信号功率谱,设所述二相编码信号功率谱最大值为Pmax,搜索所述二相编码信号功率谱第一个值为Pmax/2处对应的频率f1和最后一个值为Pmax/2处对应的频率f2,根据定义B3dB=f2-f1计算出每个所述二相编码信号功率谱对应的3dB带宽B3dB,则获得一一对应的τ-B3dB数据点;
步骤d、选择指数函数对所述τ-B3dB数据点进行数据拟合;
步骤e、求出拟合函数分别对应的系数a、b,得到所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式,其中,拟合函数指的是步骤d中数据拟合得到的拟合函数;
步骤f、令l加1;
步骤g,判断l是否等于r,如果l等于r,执行步骤h,如果l不等于r,返回步骤b;
步骤h,得出不同码长N分别对应的所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式;
步骤2、检测待估值的二相编码信号的功率谱,对所述功率谱进行幅值归一化得到P(f),对所述P(f)均匀采样M个点,变成长度为M的序列P(m),对所述P(m)做离散余弦变换得到Ypre(P(m)),对所述Ypre(P(m))进行阈值处理得到Y(P(m));
步骤3、求所述Y(P(m))的上包络,计算所述上包络的峰值个数n,根据信号码长N与所述峰值个数n的关系求出所述信号码长N,其中:N=2n+1;
步骤4、对所述Y(P(m))进行离散余弦逆变换,得到Pi(m),计算所述Pi(m)的载频fc、3dB带宽B3dB和子脉冲宽度τ,其具体包括:
a、设所述Pi(m)的最大值为Pmax’,搜索所述Pi(m)第一个值为Pmax’/2处对应的频率f’1和最后一个值为Pmax’/2处对应的频率f2’,根据定义计算出所述Pi(m)的载频fc;
b、根据定义B3dB=f2’-f1’计算出所述Pi(m)的3dB带宽B3dB;
c、根据所述子脉冲宽度τ与所述3dB带宽B3dB的拟合关系式求出所述子脉冲宽度τ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310390121.7A CN103441975B (zh) | 2013-08-30 | 2013-08-30 | 一种基于功率谱的二相编码信号参数估值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310390121.7A CN103441975B (zh) | 2013-08-30 | 2013-08-30 | 一种基于功率谱的二相编码信号参数估值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103441975A CN103441975A (zh) | 2013-12-11 |
CN103441975B true CN103441975B (zh) | 2016-12-28 |
Family
ID=49695641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310390121.7A Expired - Fee Related CN103441975B (zh) | 2013-08-30 | 2013-08-30 | 一种基于功率谱的二相编码信号参数估值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103441975B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106501787B (zh) * | 2016-11-02 | 2019-03-08 | 西安电子科技大学 | 基于平滑伪魏格纳分布的二相编码信号参数估计方法 |
CN109490862B (zh) * | 2018-11-23 | 2023-08-29 | 山东航天电子技术研究所 | 一种基于相位差分统计谱的载频估计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101984613A (zh) * | 2010-11-30 | 2011-03-09 | 中国工程物理研究院电子工程研究所 | 一种低速率bpsk突发信号码速率估计方法 |
CN102707157A (zh) * | 2012-05-18 | 2012-10-03 | 天津理工大学 | 一种基于功率谱的单频脉冲信号参数估值方法 |
CN102710560A (zh) * | 2012-05-18 | 2012-10-03 | 天津理工大学 | 一种td-scdma***的信道估计方法 |
CN103063909A (zh) * | 2012-12-18 | 2013-04-24 | 天津理工大学 | 一种基于功率谱的线性调频信号参数估值方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7778785B2 (en) * | 2008-02-14 | 2010-08-17 | Advantest Corporation | Signal-to-noise ratio measurement for discrete waveform |
-
2013
- 2013-08-30 CN CN201310390121.7A patent/CN103441975B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101984613A (zh) * | 2010-11-30 | 2011-03-09 | 中国工程物理研究院电子工程研究所 | 一种低速率bpsk突发信号码速率估计方法 |
CN102707157A (zh) * | 2012-05-18 | 2012-10-03 | 天津理工大学 | 一种基于功率谱的单频脉冲信号参数估值方法 |
CN102710560A (zh) * | 2012-05-18 | 2012-10-03 | 天津理工大学 | 一种td-scdma***的信道估计方法 |
CN103063909A (zh) * | 2012-12-18 | 2013-04-24 | 天津理工大学 | 一种基于功率谱的线性调频信号参数估值方法 |
Non-Patent Citations (2)
Title |
---|
automatic modulation recognition based on morphological operations;Yuan Zhang;《circuits,systems,and signal processing》;20130327;第32卷(第5期);2517-2525 * |
基于功率谱FFT的BPSK信号参数估计;马秀荣;《电子与信息学报》;20130515;第35卷(第5期);1252-1256 * |
Also Published As
Publication number | Publication date |
---|---|
CN103441975A (zh) | 2013-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105785324B (zh) | 基于mgcstft的线性调频信号参数估计方法 | |
CN106597408B (zh) | 基于时频分析和瞬时频率曲线拟合的高阶pps信号参数估计方法 | |
CN102680948B (zh) | 一种线性调频信号调频率和起始频率估计方法 | |
CN108469602B (zh) | 一种基于谱特征提取的脉冲信号类型自动判别方法 | |
CN101334469B (zh) | 基于分数阶傅立叶变换的风廓线雷达杂波抑制方法 | |
CN102508206B (zh) | 基于小波包去噪和功率谱熵的线性调频信号参数估计方法 | |
CN110865357A (zh) | 一种基于参数优化vmd的激光雷达回波信号降噪方法 | |
CN113378661A (zh) | 一种基于改进小波阈值和相关检测的直流电能信号去噪方法 | |
CN101984360B (zh) | 基于frft的归一化泄露lms自适应动目标检测方法 | |
CN103063909B (zh) | 一种基于功率谱的线性调频信号参数估值方法 | |
CN103675758B (zh) | 一种双曲调频信号周期斜率和起始频率估计方法 | |
CN103454622A (zh) | 基于稀疏约束的宽带雷达目标复回波去噪方法 | |
CN103995950A (zh) | 基于空域相关修正阈值的变小波系数局部放电信号消噪方法 | |
CN104901909B (zh) | 一种α非高斯噪声下chirp信号的参数估计方法 | |
WO2016004687A1 (zh) | 超高频局放信号初始时刻判别方法 | |
CN113156206B (zh) | 时频结合的含噪信号参数估计新算法 | |
CN110988834A (zh) | 一种基于自适应阈值双门限的脉冲到达时间测量方法 | |
CN111159888A (zh) | 一种基于互相关函数的协方差矩阵稀疏迭代时延估计方法 | |
CN106772271A (zh) | 一种基于匹配自相关算法的雷达信号toa估计方法 | |
CN103441975B (zh) | 一种基于功率谱的二相编码信号参数估值方法 | |
CN107315714B (zh) | 一种去卷积功率谱估计方法 | |
CN110632563B (zh) | 一种基于短时傅里叶变换的脉内频率编码信号参数测量方法 | |
CN106501787B (zh) | 基于平滑伪魏格纳分布的二相编码信号参数估计方法 | |
CN108957416A (zh) | 脉冲噪声环境下基于分数阶功率谱密度的线性调频信号参数估计方法 | |
CN116643250A (zh) | 基于稀疏线性调频信号的宽带雷达高速目标相参积累方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161228 Termination date: 20200830 |
|
CF01 | Termination of patent right due to non-payment of annual fee |