CN105277801B - 一种测量调制脉冲参数的方法 - Google Patents
一种测量调制脉冲参数的方法 Download PDFInfo
- Publication number
- CN105277801B CN105277801B CN201510487559.6A CN201510487559A CN105277801B CN 105277801 B CN105277801 B CN 105277801B CN 201510487559 A CN201510487559 A CN 201510487559A CN 105277801 B CN105277801 B CN 105277801B
- Authority
- CN
- China
- Prior art keywords
- pulse
- peak
- sequence
- waveform
- method described
- 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
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明涉及一种测量调制脉冲参数的方法,采用如下步骤:测量系统的硬件射频采样速率大于被测信号的载波频率fc;对所述被测信号的波形序列数据求绝对值;设定所述被测信号的载波周期为Tc,采样器的采样周期为Ts,计算局部求峰值的点数Np,求取包络波形;包络波形作数字低通滤波;若经过低通滤波后的波形中含有多个脉冲,则进行以下处理:<1>找出平顶脉冲顶部幅度的中位值;<2>逐个获取单个脉冲;<3>从单个脉冲中切分上升沿部分和下降沿部分;对切分得到的所述上升沿部分和所述下降沿部分进行处理,求取上升时间、下降时间、脉冲宽度等参量。本发明的测量调制脉冲参数的方法能够降低信号中的噪声干扰,检出的波形轮廓清晰,没有谐波分量,从而能够准确地测量出脉冲的相关参量。
Description
技术领域
本发明涉及一种测量调制脉冲参数的方法,属于测量领域。
技术背景
传统上测量已调射频脉冲参数,如上升时间、下降时间等,主要依靠检波器将射频脉冲包络检出,进而使用示波器进行测量。但是检波器是一种非线性器件,在检波的同时引入了大量的谐波分量,导致检出的波形轮廓模糊,如图1所示,特别是如图2所示的局部放大图,对照图1可以看出其高次谐波分量的幅度可以达到基频分量的三分之一左右,这种现象将导致脉冲参数很难得到准确有效的测量。
实验证明,使用检波器测量调制脉冲信号的参数,在使用的检波器不同、采样率不同时,存在显著的测量结果不一致的现象。在实验中,我们使用射频调制信号发生器的基带波形设计功能,设计了具有不同上升沿和下降沿参数的射频脉冲,然后使用不同的检波器(简称为A型、B型),示波器设置不同的采样率进行测量,同时为了对照还测量了基带未调脉冲的参数,产生的部分测量结果如表1~表3所示。
表1载波频率2GHz,脉冲上升沿和下降时间设计值均为100ns
表2载波频率2GHz,脉冲上升沿和下降时间设计值均为500ns
表3载波频率13.6MHz,脉冲上升沿和下降时间设计值均为500ns
显然,采用检波器的测量结果与设计值(预期值)差别较大,且同样的信号采用不同检波器、甚至是采用同一个检波器但是采样率不同时,测量结果差别也很显著。这说明传统的采用检波器的脉冲测量方法是不完善的,测量结果也是缺乏可信度的,为此,提出了本发明的方法。
发明内容
本发明的目的是针对上述现有技术的不足,提供一种测量调制脉冲参数的方法。
一种测量调制脉冲参数的方法,包括以下步骤:
(1).测量系统的硬件射频采样速率大于被测信号的载波频率fc;
(2).对所述被测信号的波形序列数据求绝对值;
(3).设定所述被测信号的载波周期为Tc,采样器的采样周期为Ts,通过(1)式计算局部求峰值的点数Np,其中是向上取整数运算符:
从射频已调波形数据序列SPF的第一个点开始,逐次取Np个点的序列,然后求这Np个点的极大值,第一段Np个点的序列的极大值记做Mp1,下一段的极大值记做Mp2,再下一段的极大值记做Mp3,第n段的极大值记做Mpn,一直取到整个波形SPF的结束,形成了一个新的序列SP,从而得到包络波形;
(4).对所述包络波形作数字低通滤波,低通滤波器的带宽的初选值为Bp0,由
(2)式计算得到:
其中TErise和TEfall分别是上升时间和下降时间的预期值,min()为求最小值函数,低通滤波器最终选择的带宽为Bp,由(3)式计算得到:
Bp=min(Bp0,fc) (3);
(5).若经过低通滤波后的波形中含有多个脉冲,则进行以下处理:<1>找出平顶脉冲顶部幅度的中位值;<2>逐个获取单个脉冲;<3>从单个脉冲中切分上升沿部分和下降沿部分;
(6).对切分得到的所述上升沿部分和所述下降沿部分进行处理,求取上升时间、下降时间、脉冲宽度等参量;
进一步地,其中步骤(1)中的采样速率为所述被测信号的载波频率fc的8倍以上;
进一步地,所述步骤(5)中的步骤<1>的具体操作为:先求出整个脉冲波形的幅度极大值Amax,然后找出波形中所有幅度大于kmaxAmax的点的幅度序列然后求该幅度序列的中位值,记做Amaxmedian,其中,Kmax为小于1的系数;
进一步地,所述Kmax取0.95;
进一步地,所述步骤(5)中的步骤<2>的具体操作为:在所述被测信号的波形序列中找出所有幅值大于kmAmaxmedian的点,将这些点的位置形成各个群,群内各个位置数相邻,而后一个群的第一个位置数和前一个群的最后一个位置数是不连续的,存在差值,并且该差值等于脉冲宽度内的点数,利用该差值将各个群分离出来,从而得到单个脉冲,km为介于0到0.1之间的系数;
进一步地,所述km取0.05;
进一步地,所述步骤(6)中的所述上升时间取0.1Apeak到0.9Apeak,其中Apeak为脉冲峰值;
进一步地,所述采样器为数字示波器;
进一步地,当0.1Apeak和0.9Apeak这两个点不能准确定位时,先求kLApeak到kHApeak之间的时间trise0,其中kL是一个介于0.1到0.5之间的常数,kH是一个介于0.5到0.9之间的常数,然后利用式(4)估算0.1倍峰值到0.9倍峰值的上升时间trise:
进一步地,所述单个脉冲的脉冲参量包括下降时间、脉冲宽度。
采用上述技术方案,与传统的采用检波器检测相比,本发明具有以下的有益效果:本发明的测量调制脉冲参数的方法能够降低信号中的噪声干扰,检出的波形轮廓清晰,没有谐波分量,从而能够准确地测量出脉冲的相关参量。
附图说明
图1为现有技术已调射频脉冲经检波器检出信号波形;
图2为附图1椭圆圈出的局部放大图;
图3为本发明脉冲参量测量装置;
图4为本发明已调射频脉冲波形;
图5为本发明已调射频脉冲波形包络(基带信号)的检出;
图6为本发明直接数字检波后的波形图;
图7为图6椭圆标出的局部放大图;
图8为检波包络滤波后波形图;
图9为附图8的脉冲顶部局部放大图;
图10为切分后的单独脉冲图;
图11为切分后的上升沿部分和下降沿部分图。
具体实施方式
为了使本发明的内容更容易被清楚地理解,下面根据具体实施例并结合附图,对本发明作进一步详细的说明。
如图1-11所示,一种测量调制脉冲参数的方法,包括以下步骤:
(1).测量系统的硬件射频采样速率大于被测信号的载波频率fc;
(2).对所述被测信号的波形序列数据求绝对值;
(3).设定所述被测信号的载波周期为Tc,采样器的采样周期为Ts,通过(1)式计算局部求峰值的点数Np,其中是向上取整数运算符:
从射频已调波形数据序列SPF的第一个点开始,逐次取Np个点的序列,然后求这Np个点的极大值,第一段Np个点的序列的极大值记做Mp1,下一段的极大值记做Mp2,再下一段的极大值记做Mp3,第n段的极大值记做Mpn,一直取到整个波形SPF的结束,形成了一个新的序列SP,从而得到包络波形,如表4所示。显然每两个相邻的Mp值的时间间隔是Tp=NpTs。显然表4所述序列SP的持续时间LTP需要大于一个被测脉冲的周期。如图5所示,这种方法能够有效检测出射频已调脉冲的外围轮廓,即得到包络数据。
表4分段极大值序列Sp表
(4).直接检波得到的包络含有数字量化等原因造成的量化噪声,如图6和图7所示。这种噪声影响上升沿测试的准确度,甚至使得测量失效,为此要对所述包络波形作数字低通滤波,低通滤波器的带宽的初选值为Bp0,由(2)式计算得到:
其中TErise和TEfall分别是上升时间和下降时间的预期值,min()为求最小值函数,低通滤波器最终选择的带宽为Bp,由(3)式计算得到:
Bp=min(Bp0,fc) (3),
检波包络滤波后波形如图8所示,滤波后脉冲顶部局部放大图9所示,显然从时域波形看噪声已经得到有效抑制。
滤波器会造成上升时间和下降时间测量不确定度的增加,其贡献的不确定度分量是:
其中K是一个常数,较优地可以将K取为0.4。
(5).如果一个波形中含有多个脉冲,需要进行以下处理:(1)找出平顶脉冲顶部幅度的中位值。具体操作方法是先求出整个脉冲波形的幅度极大值Amax,然后找出波形中所有幅度大于kmaxAmax的点的幅度序列(kmax是一个接近1但是小于1的系数,较优地可以取0.95),然后求该幅度序列的中位值,记做Amax median。(2)逐个截取单个脉冲。在原波形序序列中找出所有幅度大于kmAmax median的点的位置(km是一个介于0到0.1之间的系数,较优地可以取0.05),显然这些位置形成了如表5所示的结构,每个群代表一个单独的脉冲,群内各个位置数相邻,后一个群的第一个位置数和前一个群的最后一个位置数是不连续的,存在差值,并且该差值等于脉冲宽度内的点数,利用该差值将各个群分离出来,从而得到一个单独的脉冲,如图10所示。(3)从单独脉冲中切分上升沿部分和下降沿部分。具体做法是:求出单独脉冲最大值的位置,以该位置为界限,左边的是上升沿对应的部分,右边是下降沿对应的部分,如图11所示。
表5位置序列
(6).对切分得到的上升沿部分和下降沿部分进行处理,求取上升时间、下降时间、脉冲宽度等参量。单个脉冲数据的脉冲参量,包括但不限于:下降时间、脉冲宽度等。一般上升时间指0.1Apeak到0.9Apeak之间的时间,其中Apeak指脉冲峰值。在一些特殊情况下,由于脉冲底部和平坦顶部噪声的作用,0.1Apeak和0.9Apeak这两个点不一定能准确定位,在这种情况下下可以先求kLApeak到kHApeak之间的时间trise0,kL是一个介于0.1到0.5之间的常数,如果有可能获得准确时间位置应该取得更接近0.1一些,kH是一个介于0.5到0.9之间的常数,如果有可能获得准确时间位置应该取得更接近0.9一些。然后可以利用式(4)估算0.1倍峰值到0.9倍峰值的上升时间trise:
脉冲的下降时间可以用类似方法求得,其他的脉冲宽度等参量也可以基于图11所示的截断波形求得。
实施案例
利用图3的实验装置进行实验,得到若干组实验结果,如表6~表8表所示。显然本发明方法测量值和设计值、基带脉冲测量结果有很好的吻合度,比之表1~表2的检波器测量结果,可信度更高。
表6 2GHz载波频率,上升下降时间设计值为100ns的已调脉冲测量结果
表7 2GHz载波频率,上升下降时间设计值为300ns的已调脉冲测量结果
表8 2GHz载波频率,上升下降时间设计值为500ns的已调脉冲测量结果
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种测量调制脉冲参数的方法,其特征在于:包括以下步骤:
(1).测量系统的硬件射频采样速率大于被测信号的载波频率fc;
(2).对所述被测信号的波形序列数据求绝对值;
(3).设定所述被测信号的载波周期为Tc,采样器的采样周期为Ts,通过(1)式计算局部求峰值的点数Np,其中是向上取整数运算符:
从射频已调波形数据序列SPF的第一个点开始,逐次取Np个点的序列,然后求这Np个点的极大值,第一段Np个点的序列的极大值记做Mp1,下一段的极大值记做Mp2,再下一段的极大值记做Mp3,第n段的极大值记做Mpn,一直取到整个波形SPF的结束,形成了一个新的序列SP,从而得到包络波形;
(4).对所述包络波形作数字低通滤波,低通滤波器的带宽的初选值为Bp0,由(2)式计算得到:
其中TErise和TEfall分别是上升时间和下降时间的预期值,min()为求最小值函数,低通滤波器最终选择的带宽为Bp,由(3)式计算得到:
Bp=min(Bp0,fc) (3);
(5).若经过低通滤波后的波形中含有多个脉冲,则进行以下处理:<1>找出平顶脉冲顶部幅度的中位值;<2>逐个获取单个脉冲;<3>从单个脉冲中切分上升沿部分和下降沿部分;
(6).对切分得到的所述上升沿部分和所述下降沿部分进行处理,求取上升时间、下降时间、脉冲宽度参量。
2.根据权利要求1所述的方法,其特征在于:其中步骤(1)中的采样速率为所述被测信号的载波频率fc的8倍以上。
3.根据权利要求1所述的方法,其特征在于:所述步骤(5)中的步骤<1>的具体操作为:先求出整个脉冲波形的幅度极大值Amax,然后找出波形中所有幅度大于kmaxAmax的点的幅度序列,然后求该幅度序列的中位值,记做Amaxmedian,其中,Kmax为小于1的系数。
4.根据权利要求3所述的方法,其特征在于:所述系数Kmax取0.95。
5.根据权利要求3所述的方法,其特征在于:所述步骤(5)中的步骤<2>的具体操作为:在所述被测信号的波形序列中找出所有幅值大于kmAmaxmedian的点,将这些点的位置形成各个群,群内各个位置数相邻,而后一个群的第一个位置数和前一个群的最后一个位置数是不连续的,存在差值,并且该差值等于脉冲宽度内的点数,利用该差值将各个群分离出来,从而得到单个脉冲,km为介于0到0.1之间的系数。
6.根据权利要求5所述的方法,其特征在于:所述km取0.05。
7.根据权利要求1所述的方法,其特征在于:所述步骤(6)中的所述上升时间取0.1Apeak到0.9Apeak,其中Apeak为脉冲峰值。
8.根据权利要求1所述的方法,其特征在于:所述采样器为数字示波器。
9.根据权利要求7所述的方法,其特征在于:当0.1Apeak和0.9Apeak这两个点不能准确定位时,先求kLApeak到kHApeak之间的时间trise0,其中kL是一个介于0.1到0.5之间的常数,kH是一个介于0.5到0.9之间的常数,然后利用式(4)估算0.1倍峰值到0.9倍峰值的上升时间trise:
10.根据权利要求5所述的方法,其特征在于:所述单个脉冲的脉冲参量包括下降时间、脉冲宽度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510487559.6A CN105277801B (zh) | 2015-08-10 | 2015-08-10 | 一种测量调制脉冲参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510487559.6A CN105277801B (zh) | 2015-08-10 | 2015-08-10 | 一种测量调制脉冲参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105277801A CN105277801A (zh) | 2016-01-27 |
CN105277801B true CN105277801B (zh) | 2018-08-07 |
Family
ID=55147162
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510487559.6A Active CN105277801B (zh) | 2015-08-10 | 2015-08-10 | 一种测量调制脉冲参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105277801B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106597122B (zh) * | 2017-01-10 | 2019-03-19 | 烟台大学 | 一种雷达及通信信号的脉冲宽度检测算法 |
CN108509377B (zh) * | 2018-03-19 | 2021-07-27 | 东南大学 | 一种基于沿特征提取的脉冲信号到达时间与脉宽估计方法 |
CN112666403A (zh) * | 2020-12-30 | 2021-04-16 | 上海微波技术研究所(中国电子科技集团公司第五十研究所) | 一种基于高频采样的脉冲调制信号参数测量系统及方法 |
CN114089290B (zh) * | 2021-11-12 | 2022-08-30 | 扬州宇安电子科技有限公司 | 一种快速高效测频引导噪声干扰方法 |
WO2024032432A1 (zh) * | 2022-08-08 | 2024-02-15 | 深圳市恒运昌真空技术有限公司 | 一种射频电源信号采集方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1538763A (zh) * | 2003-04-17 | 2004-10-20 | 华为技术有限公司 | 一种脉冲宽度调制方法及装置 |
CN201039103Y (zh) * | 2007-04-26 | 2008-03-19 | 华南理工大学 | 用于脉冲调幅波的检波电路 |
US7516032B2 (en) * | 2001-12-20 | 2009-04-07 | Stmicroelectronics Pvt. Ltd. | Resolution in measuring the pulse width of digital signals |
CN102664588A (zh) * | 2012-05-14 | 2012-09-12 | 中国航空工业集团公司北京长城计量测试技术研究所 | 用于频率调制和幅度调制信号的数字化解调方法及装置 |
CN104022761A (zh) * | 2014-05-23 | 2014-09-03 | 中国电子科技集团公司第四十一研究所 | 一种脉冲调制信号的测试方法 |
-
2015
- 2015-08-10 CN CN201510487559.6A patent/CN105277801B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7516032B2 (en) * | 2001-12-20 | 2009-04-07 | Stmicroelectronics Pvt. Ltd. | Resolution in measuring the pulse width of digital signals |
CN1538763A (zh) * | 2003-04-17 | 2004-10-20 | 华为技术有限公司 | 一种脉冲宽度调制方法及装置 |
CN201039103Y (zh) * | 2007-04-26 | 2008-03-19 | 华南理工大学 | 用于脉冲调幅波的检波电路 |
CN102664588A (zh) * | 2012-05-14 | 2012-09-12 | 中国航空工业集团公司北京长城计量测试技术研究所 | 用于频率调制和幅度调制信号的数字化解调方法及装置 |
CN104022761A (zh) * | 2014-05-23 | 2014-09-03 | 中国电子科技集团公司第四十一研究所 | 一种脉冲调制信号的测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105277801A (zh) | 2016-01-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105277801B (zh) | 一种测量调制脉冲参数的方法 | |
EP2605025B1 (en) | Frequency mask trigger with non uniform bandwidth segments | |
EP1672376B1 (en) | Real-time RF/IF signature trigger | |
CN106443203A (zh) | 一种脉冲信号检测系统及方法 | |
CN109617631B (zh) | 基于数字信道化瞬时参数测量的侦察系统自适应接收方法 | |
US20150066409A1 (en) | Real-time Oscilloscope For Generating a Fast Real-time Eye Diagram | |
US20200209281A1 (en) | Measurement apparatus and method for analyzing a waveform of a signal | |
US8762096B2 (en) | Trigger generation for digital modulation signal analysis | |
KR101765239B1 (ko) | 광대역 다채널 수신기에서 채널간 간섭에 의한 오탐지 쇼트 펄스 pdw 제거 방법 및 장치 | |
EP2985610A1 (en) | Method for determining a correlated waveform on a real time oscilloscope | |
JP5752481B2 (ja) | パッシブレーダ装置、誘導装置及び電波検出方法 | |
KR101207779B1 (ko) | 전자전 장비를 위한 신호 처리 장치 및 신호 처리 방법 | |
CN105049105B (zh) | 一种频率分集信号的频率提取方法 | |
Kubankova et al. | Algorithms of digital modulation classification and their verification | |
CN108957145B (zh) | 局域全闪定位网的地闪回击脉冲检出系统及方法 | |
US6853933B2 (en) | Method of identifying spectral impulses for Rj Dj separation | |
US11640523B2 (en) | Method for training a neural network, method for automatically characterizing a measurement signal, measurement apparatus and method for analyzing | |
US9395407B2 (en) | Enhanced dynamic range RF pulse measurement system | |
Recordon et al. | A comparator-based technique for identification of intentional electromagnetic interference attacks | |
US10712367B2 (en) | Method for analyzing a signal as well as measurement and analyzing device | |
Carobbi et al. | On the alleged inconsistency in cispr 16-1-1 performance tests for disturbance analyzers | |
US20240255572A1 (en) | Test and measurement instrument having tailored jitter compensation | |
Westenberger | Use of time domain methods for CISPR16 compliant EMI measurements | |
CN108007548A (zh) | 一种通过扫频诊断设备故障的方法 | |
Kubankova et al. | Methods of digital modulation recognition and their testing |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20211230 Address after: 100191 No. 40, Haidian District, Beijing, Xueyuan Road Patentee after: CHINA ACADEMY OF INFORMATION AND COMMUNICATIONS Address before: 100191 No. 52 Garden North Road, Beijing, Haidian District Patentee before: CHINA ACADEME OF TELECOMMUNICATION RESEARCH OF MIIT |
|
TR01 | Transfer of patent right |