CN105721078A - 一种通用群时延纹波自动化测试方法 - Google Patents
一种通用群时延纹波自动化测试方法 Download PDFInfo
- Publication number
- CN105721078A CN105721078A CN201610159382.1A CN201610159382A CN105721078A CN 105721078 A CN105721078 A CN 105721078A CN 201610159382 A CN201610159382 A CN 201610159382A CN 105721078 A CN105721078 A CN 105721078A
- Authority
- CN
- China
- Prior art keywords
- frequency
- test
- delay
- group delay
- group
- 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
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/364—Delay profiles
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18519—Operations control, administration or maintenance
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/14—Relay systems
- H04B7/15—Active relay systems
- H04B7/185—Space-based or airborne stations; Stations for satellite systems
- H04B7/1851—Systems using a satellite or space-based relay
- H04B7/18515—Transmission equipment in satellites or space-based relays
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Electromagnetism (AREA)
- Radio Relay Systems (AREA)
Abstract
一种通用群时延纹波自动化测试方法。测试信号源产生两路基带信号,一路送数字示波器直接采样作为参考,另一路对高频载波进行调频,此调频波通过卫星转发器后,由测试频谱仪变频到中频,再经过数字示波器采样,得出欲测试的频率范围内的群时延曲线数据,根据最小二乘法将得到的群时延曲线数据进行拟合,得出反映整体趋势的拟合曲线数据,比对拟合前后的两组曲线数据,找出指标频率范围内的最大差值,即群时延纹波。该方法使用“最小二乘法”拟合曲线,从总偏差最小的角度来取近似曲线,即从大量给定数据中找出规律,并构造一条曲线反映数据点的总趋势,消除其局部波动,从而得出反映全局特性的群时延纹波。
Description
技术领域
本发明适用于航天器地面测试领域,特别是对通信卫星转发器群时延纹波自动化测试,属于射频测试技术领域。
背景技术
群时延是描述传输系统相位频率特性的一个重要参量,群时延纹波是描述群时延的波动、起伏情况,是衡量连续信号通过传输系统后相位线性度的重要指标。
群时延纹波不仅决定了系统产生的信号传输时延的波动情况,而且与信号传输失真密切相关,对信号的传输质量有很大影响,是现在卫星通信系统的一项主要技术指标。一旦系统的群时延纹波特性不好,将会造成信号相位失真,降低信噪比,增加误码率。
传统通信卫星转发器通道群时延特性测试利用调制器对信号调制,经卫星转发器转发再由解调器解调后,送入鉴频器与参考信号进行比相,从而得到群时延特性数据,但对群时延曲线数据的波动、起伏等情况未作分析,因此无法深入了解卫星转发器通道的群时延纹波特性,从而无法全面了解卫星转发器通道对传输信号质量的影响。
随着通信卫星转发器技术的不断发展,对转发器通道的群时延性能要求越来越高,对其波动、起伏情况要求更为严格,就此提出一种群时延纹波特性,对此发明一种通用群时延纹波自动化测试方法,达到卫星转发器群时延纹波测试的目的。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种通用群时延纹波自动化测试方法,能够根据卫星转发器的特性进行重要参数的设置,减小由于固定参数引入的测试误差,同时对群时延曲线数据进行全局趋势的拟合,比对拟合前后两组曲线数据,从而得出卫星转发器群时延的波动、起伏特性,即群时延纹波,更好地描述转发器的通道性能。
本发明的技术方案:一种通用群时延纹波自动化测试方法,步骤如下:
(1)根据被测对象带宽,设置测试步进、调制频率,设置“最小二乘法”的拟合阶数;
(2)由信号源产生的一路基带信号直接送入数字示波器,另一路为上行调制信号送入卫星,卫星转发下来的下行调制信号进入频谱仪,频谱仪将下行调制信号变频到中频送入数字示波器;
(3)数字示波器对不同频率点的基带信号和下行调制信号分别采样,通过采样数据得出基带信号和下行调制信号的初始相位,比较两组初始相位,得出相位差,再除以调制信号角频率,计算得出该频率点的时延,从而获得被测对象带宽内所有测试频率点的群时延曲线数据;再经过频率中心值处理和群时延最小值处理,得到每个频点的频率对应值和每个频点的群时延最小值处理后的数值,即得出群时延相对数值;
(4)根据“最小二乘法”对群时延相对数值进行相应阶数拟合,得出拟合曲线数据;
(5)比对拟合前、后各测试频点的时延数据,找出指标带宽内的最大差值,即群时延纹波。
步骤(1)中测试步进B为被测对象测试带宽,测试点数根据测试要求设定;所述的拟合阶数设置为6。
步骤(3)中频率中心值处理的具体方法为:各测试频点的下行频率-下行中心频率,得到每个频点的频率对应值。
步骤(3)中群时延最小值处理的具体方法为:首先获得所有测试频率点中最小的群时延数据值,之后用每个测试频点群时延数据值-最小的群时延数据值,得到每个测试频点群时延最小化处理后的数值。
步骤(5)的具体方法为:将每个测试频点的群时延最小化处理后的数值与每个测试频点经“最小二乘法”拟合后的数据对应频点分别作差,并取其绝对值,得出每个测试频点的差值绝对值数据;在每个测试频点的差值绝对值数据中找出指标带宽内的最大值,即为指标带宽内的群时延纹波。
本发明与现有技术相比的优点在于:
(1)本发明可以对测试仪器的频率步进、调制频率和拟合阶数等重要参数进行设置,减小由于固定参数不合理造成的测试误差;
(2)本发明采用“最小二乘法”拟合曲线,不要求曲线完全通过所有已知数据点,而是从总偏差最小的角度来取近似曲线,即从大量给定数据中找出规律,并构造一条曲线反映数据点的总趋势,以消除其局部波动,从而得出反映全局特性的群时延纹波,更为准确地描述被测对象的群时延波动、起伏特性。
附图说明
图1是群时延纹波自动化测试设备连接图。
图2是群时延纹波自动化测试方法流程图。
具体实施方式
下面结合附图对本发明做进一步说明。
如图1所示,通用群时延纹波自动化测试方法是建立在包括测试信号源、测试频谱仪、数字示波器、测试软件在内的测试系统基础上。其中测试信号源提供两路信号,一路为基带信号(即调制频率信号)进入数字示波器作为参考信号,另一路用基带信号调制的射频信号进入卫星转发器。测试频谱仪将经过转发器后的射频调制信号变频到中频,送入数字示波器,数字示波器对两路信号进行数字采样,再经过测试软件特殊算法进行计算(数字解调)得到群时延曲线数据,之后按照“最小二乘法”进行拟合运算,比对拟合前后曲线数据,找出指标带宽内的最大差值,即群时延纹波。
如图2所示,通用群时延纹波自动化测试方法的详细实施步骤如下:
1、根据被测对象(卫星转发器)带宽B,通过软件设置测试频率步进b(测试频点数可根据测试要求设定)、调制频率,设置“最小二乘法”的拟合阶数N(N∈整数,常设为6)。
2、由信号源产生的一路基带信号直接送入数字示波器,另一路为上行调制信号送入卫星,卫星转发下来的下行调制信号进入频谱仪,频谱仪将下行调制信号变频到中频送入数字示波器。
3、数字示波器对不同频率点的基带信号和下行调制信号分别采样,通过采样数据得出基带信号和下行调制信号的初始相位,比较两组初始相位,得出相位差,除以调制信号角频率,计算出该频率点的时延,从而获得所有测试频率点的群时延曲线数据。(具体见航天器工程,第19卷、第1期,卫星群时延特性测量新方案的设计和实现)。再经过频率(中心值处理)和群时延(最小值处理),即得出群时延的相对数值。频率中心值处理的具体方法为:各测试频点的下行频率-下行中心频率,得到每个频点的频率对应值。群时延最小值处理的具体方法为:首先获得所有测试频率点中最小的群时延数据值,之后用每个测试频点群时延数据值-最小的群时延数据值,得到每个测试频点群时延最小化处理后的数值。
4、根据“最小二乘法”对群时延相对数值进行相应阶数拟合,得出拟合曲线数据。
5、将每个测试频点的群时延最小化处理后的数值与每个测试频点经“最小二乘法”拟合后的数据对应频点分别作差,并取其绝对值,得出每个测试频点的差值绝对值数据;在每个测试频点的差值绝对值数据中找出指标带宽的内的最大值,即为指标带宽内的群时延纹波。
实施例
以一路卫星转发器为例,上行中心频率为12700MHz,下行中心频率为10700MHz,带宽为48MHz,进行以下步骤:
1、按照文件要求,测试群时延特性需要左右各扩展4MHz,即测试带宽变为56MHz,通过测试软件设置频率步进为2MHz,即测试点数为29,调制频率为0.4MHz,设置“最小二乘法”拟合阶数为6阶。
2、信号源根据软件设置,产生一路基带信号(调制频率为0.4MHz)直接送入数字示波器,另一路为上行调制信号(载波为12700MHz,调制频率为0.4MHz)送入卫星,卫星转发下来的下行调制信号(载波为10700MHz,调制频率为0.4MHz)进入频谱仪,频谱仪将下行调制信号变频到中频(21.4MHz)送入数字示波器。
3、数字示波器对基带信号和下行调制信号分别采样,通过采样数据得出基带信号和下行调制信号的初始相位,比较两组初始相位,得出相位差,除以调制信号角频率,计算出该频点的时延,从而获得所有测试频点的群时延曲线数据。如表1所示,第1列为序号,第2列为下行频率,第3列为群时延数据,为方便数据处理,对下行频率进行中心值处理,对群时延进行最小值处理,生成第4列频率(中心值处理)数据和第5列群时延(最小值处理)数据。
4、按照“最小二乘法”对群时延曲线数据(即表1第4、5列数据)进行6阶拟合,拟合得出第6列群时延(二乘法拟合)数据。
5、将第5列群时延(最小值处理)数据和第6列群时延(二乘法拟合)数据对应频点分别作差,并取其绝对值,得出第7列差值绝对值数据,按照文件要求,找出第7列数据中指标36MHz带宽的内的最大值为1.6ns(第21行),该值就是指标带宽内的群时延纹波。
表1群时延纹波数据
本发明未详细描述内容为本领域技术人员公知技术。
Claims (5)
1.一种通用群时延纹波自动化测试方法,其特征在于步骤如下:
(1)根据被测对象带宽,设置测试步进、调制频率,设置“最小二乘法”的拟合阶数;
(2)由信号源产生的一路基带信号直接送入数字示波器,另一路为上行调制信号送入卫星,卫星转发下来的下行调制信号进入频谱仪,频谱仪将下行调制信号变频到中频送入数字示波器;
(3)数字示波器对不同频率点的基带信号和下行调制信号分别采样,通过采样数据得出基带信号和下行调制信号的初始相位,比较两组初始相位,得出相位差,再除以调制信号角频率,计算得出该频率点的时延,从而获得被测对象带宽内所有测试频率点的群时延曲线数据;再经过频率中心值处理和群时延最小值处理,得到每个频点的频率对应值和每个频点的群时延最小值处理后的数值,即得出群时延相对数值;
(4)根据“最小二乘法”对群时延相对数值进行相应阶数拟合,得出拟合曲线数据;
(5)比对拟合前、后各测试频点的时延数据,找出指标带宽内的最大差值,即群时延纹波。
2.根据权利要求1所述的一种通用群时延纹波自动化测试方法,其特征在于:步骤(1)中测试步进B为被测对象测试带宽,测试点数根据测试要求设定;所述的拟合阶数设置为6。
3.根据权利要求1所述的一种通用群时延纹波自动化测试方法,其特征在于:步骤(3)中频率中心值处理的具体方法为:各测试频点的下行频率-下行中心频率,得到每个频点的频率对应值。
4.根据权利要求1所述的一种通用群时延纹波自动化测试方法,其特征在于:步骤(3)中群时延最小值处理的具体方法为:首先获得所有测试频率点中最小的群时延数据值,之后用每个测试频点群时延数据值-最小的群时延数据值,得到每个测试频点群时延最小化处理后的数值。
5.根据权利要求1所述的一种通用群时延纹波自动化测试方法,其特征在于:步骤(5)的具体方法为:将每个测试频点的群时延最小化处理后的数值与每个测试频点经“最小二乘法”拟合后的数据对应频点分别作差,并取其绝对值,得出每个测试频点的差值绝对值数据;在每个测试频点的差值绝对值数据中找出指标带宽内的最大值,即为指标带宽内的群时延纹波。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610159382.1A CN105721078B (zh) | 2016-03-21 | 2016-03-21 | 一种通用群时延纹波自动化测试方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610159382.1A CN105721078B (zh) | 2016-03-21 | 2016-03-21 | 一种通用群时延纹波自动化测试方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105721078A true CN105721078A (zh) | 2016-06-29 |
CN105721078B CN105721078B (zh) | 2018-03-09 |
Family
ID=56158818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610159382.1A Active CN105721078B (zh) | 2016-03-21 | 2016-03-21 | 一种通用群时延纹波自动化测试方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105721078B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107294626A (zh) * | 2017-05-16 | 2017-10-24 | 上海卫星工程研究所 | 高效的高码率信号误码率测试方法 |
CN110351165A (zh) * | 2019-07-22 | 2019-10-18 | 中国电子科技集团公司第三十八研究所 | 一种数字阵列接收通道群时延起伏的测量方法及系统 |
CN110557190A (zh) * | 2019-08-25 | 2019-12-10 | 中国科学院国家授时中心 | 在轨卫星的转发器群时延测量方法 |
CN111812966A (zh) * | 2020-07-02 | 2020-10-23 | 北京航天飞行控制中心 | 基于多航天器的天地时差的确定方法、确定装置与处理器 |
CN116073891A (zh) * | 2023-03-13 | 2023-05-05 | 北京中天星控科技开发有限公司 | 一种卫星群时延测试方法、装置及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102325058A (zh) * | 2011-09-06 | 2012-01-18 | 北京空间飞行器总体设计部 | 一种变频系统群时延测试方法 |
-
2016
- 2016-03-21 CN CN201610159382.1A patent/CN105721078B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102325058A (zh) * | 2011-09-06 | 2012-01-18 | 北京空间飞行器总体设计部 | 一种变频系统群时延测试方法 |
Non-Patent Citations (2)
Title |
---|
沙海: "卫星导航系统传输信道的群时延测量方法研究与应用", 《中国优秀学位论文全文数据库 信息科技辑》 * |
黄坤超: "时延测试方法研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107294626A (zh) * | 2017-05-16 | 2017-10-24 | 上海卫星工程研究所 | 高效的高码率信号误码率测试方法 |
CN110351165A (zh) * | 2019-07-22 | 2019-10-18 | 中国电子科技集团公司第三十八研究所 | 一种数字阵列接收通道群时延起伏的测量方法及系统 |
CN110557190A (zh) * | 2019-08-25 | 2019-12-10 | 中国科学院国家授时中心 | 在轨卫星的转发器群时延测量方法 |
CN111812966A (zh) * | 2020-07-02 | 2020-10-23 | 北京航天飞行控制中心 | 基于多航天器的天地时差的确定方法、确定装置与处理器 |
CN111812966B (zh) * | 2020-07-02 | 2022-01-28 | 北京航天飞行控制中心 | 基于多航天器的天地时差的确定方法、确定装置与处理器 |
CN116073891A (zh) * | 2023-03-13 | 2023-05-05 | 北京中天星控科技开发有限公司 | 一种卫星群时延测试方法、装置及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105721078B (zh) | 2018-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105721078A (zh) | 一种通用群时延纹波自动化测试方法 | |
US11223431B2 (en) | Systems, methods, and devices for electronic spectrum management | |
CN102164013B (zh) | 电台共址干扰的emc综合预测分析系统和分析方法 | |
CN111381226B (zh) | 一种基于多频段接收信号强度的测距增强方法 | |
CN101969329B (zh) | 天线组阵中基于循环互相关的信号相位差估计装置与方法 | |
CN102638319B (zh) | 一种导航卫星二进制偏移载波信号的调制性能测试方法 | |
CN104601259A (zh) | 具有同相正交信号失衡估测及校正技术的无线通信接收器 | |
CN105656834A (zh) | 一种新型宽带接收机iq通道失配的数字校正方法 | |
CN109547128A (zh) | Massive MIMO端到端性能的快速测试方法 | |
CN108111208A (zh) | 一种灵活码本配置的毫米波系统角度域信道估计方法 | |
CN103869334B (zh) | Gnss空间信号畸变的自动识别与处理方法 | |
CN106842198A (zh) | 一种基于凸优化理论的子带融合成像算法 | |
CN106324598A (zh) | 一种多通道子带信号的均衡、同步方法及其系统 | |
CN105681232A (zh) | 一种基于共享信道和压缩感知的大规模mimo信道估计方法 | |
CN103188067B (zh) | 一种扩频系统的码片时钟频率偏差误差估计及校正的方法 | |
CN109633625B (zh) | 基于ofdm信号归一化一阶矩的单脉冲测角方法及跟踪应用 | |
EP3556043A1 (en) | Improved radio communications system and method based on time twisted waves | |
CN104967491B (zh) | 多通道幅相测试系统信号接收处理方法 | |
CN105025570A (zh) | 基于Alamouti编码的移动台定位的设计方法 | |
CN110149118A (zh) | 一种基于非均匀滤波器组的动态信道化方法 | |
CN106789796A (zh) | 一种基于正交数字包络检波的0/π调制测角方法 | |
CN114257280A (zh) | 一种基于空间扫描的自适应毫米波波束搜索方法 | |
CN103368891A (zh) | 信号处理装置及信号处理方法 | |
CN109302243B (zh) | 一种n端口网分的抗同频干扰扫描方法 | |
CN111953435A (zh) | 一种基于pcal信号的天线组阵设备链路标校方法 |
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 |