CN112305511A - 一种雷达接收机数字基带信号录取分析系统和方法 - Google Patents

一种雷达接收机数字基带信号录取分析系统和方法 Download PDF

Info

Publication number
CN112305511A
CN112305511A CN202011031659.5A CN202011031659A CN112305511A CN 112305511 A CN112305511 A CN 112305511A CN 202011031659 A CN202011031659 A CN 202011031659A CN 112305511 A CN112305511 A CN 112305511A
Authority
CN
China
Prior art keywords
power
signal
noise
formula
harmonic
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
Application number
CN202011031659.5A
Other languages
English (en)
Other versions
CN112305511B (zh
Inventor
陈利彬
吴来萍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Radio Measurement
Original Assignee
Beijing Institute of Radio Measurement
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Radio Measurement filed Critical Beijing Institute of Radio Measurement
Priority to CN202011031659.5A priority Critical patent/CN112305511B/zh
Publication of CN112305511A publication Critical patent/CN112305511A/zh
Application granted granted Critical
Publication of CN112305511B publication Critical patent/CN112305511B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • G01S7/4021Means for monitoring or calibrating of parts of a radar system of receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

本发明的一个实施例公开了一种雷达接收机数字基带信号录取分析系统和方法,包括:计算机设备和数据录取装置;其中,所述计算机设备响应于用户的输入而生成控制指令;所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;所述计算机设备对所述数据进行分析。

Description

一种雷达接收机数字基带信号录取分析系统和方法
技术领域
本发明涉及通信领域。更具体地,涉及一种雷达接收机数字基带信号录取分析系统和方法。
背景技术
现代雷达系统的接收机通常采用光纤作为数字基带信号的实时传输方式。为了对于接收机性能进行有效评估,需要采用数字基带信号分析方法完成接收机性能指标的验证。而传统的台式仪器设备如示波器、频谱仪等,由于不具备光纤接口也不具备满足雷达接收机需求的数字基带信号处理分析能力,无法在接收机性能测试中进行应用,导致测试验证手段存在缺失。
发明内容
本发明的目的在于提供一种雷达接收机数字基带信号录取分析系统和方法,以解决现有技术存在的问题中的至少一个。
为达到上述目的,本发明采用下述技术方案:
第一方面,本发明提供一种雷达接收机数字基带信号录取分析系统,包括:
计算机设备和数据录取装置;
其中,所述计算机设备响应于用户的输入而生成控制指令;
所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;
所述计算机设备对所述数据进行分析。
在一个具体实施方式中,所述数据录取装置包括:
FPGA核心控制器、同步触发连接器、光电收发器模块、第一四通道指示灯、第二四通道指示灯、隔离器和存储器;
其中,
所述光电收发器模块包括:
第一光电收发器、第二光电收发器、第三光电收发器和第四光电收发器;
所述光电收发器模块用于接收外部被测接收机输出的光信号并将所述光信号转换为串行电信号,或将FPGA核心控制器输出的串行电信号转换为光信号;
所述同步触发连接器用于接收外部接收机发出的同步信号;
所述FPGA核心控制器用于接收所述同步信号和所述串行电信号,并对所述同步信号和所述串行电信号进行处理,将处理后得到的数据发送给计算机设备;
所述第一四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于发送状态;
所述第二四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于接收链路成功连接状态;
所述隔离器用于对所述同步触发连接器输出的同步信号进行光耦隔离,并将所述进行隔离后的同步信号输出至FPGA核心控制器;
所述存储器在所述FPGA核心控制器的控制下进行数据的读写,完成并行接收数据缓存。
在一个具体实施方式中,所述计算机设备包括指标分析模块;
所述指标分析模块被配置为:
对接收机打包的信号进行拆分,获取用户所选的一路数字基带信号的数据;
按照用户所设置的一路数字基带信号所需的FFT点数N对于所述数据进行加窗和FFT运算,以获得功率谱线,其中,N为用户设置的一路数字基带信号所需的FFT点数,记|Y(k)|2为运算后的第k点谱线功率,P为所加窗函数的阶数,L为所加窗函数带来的功率损失百分比,S为杂散数和H为谐波总数;
在0至
Figure BDA0002703855010000031
数字基带信号范围内对于所述谱线功率|Y(k)|2求最大值,得到最大值对应的位置为基波信号的位置I1,基波信号取值范围B1为I1-P至I1+P;二次谐波取值范围B2为2×I1-P至2×I1+P,同理可得到三次至六次谐波的取值范围B3~B6;在P至
Figure BDA0002703855010000032
且不包含B1~B6的数字基带信号范围内,对于|Y(k)|2求最大值,得到最大值对应的位置为最大杂散信号的位置Is,最大杂散信号B7的取值范围为Is-P至Is+P;
当最大杂散的谱线功率小于二次至六次谐波谱线功率最小值时,所述杂散数S等于0和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(1),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000033
当最大杂散的谱线功率大于或者等于二次至六次谐波谱线功率最小值时,所述杂散数S等于1和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(2),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000041
公式(1)和(2)中,N为用户设置的一路数字基带信号所需的FFT点数,|Y(k)|2为运算后的第k点谱线的功率,P为所加窗函数的阶数,S为杂散数和H为谐波总数;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、噪底指标PNoise_floor、在0至
Figure BDA0002703855010000042
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置I1和所加窗函数的阶数P代入公式(3),得到实测信号功率PMain
Figure BDA0002703855010000043
式中,L为所加窗函数带来的功率损失百分比、PNoise_floor为噪底指标,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率,I1为在0至
Figure BDA0002703855010000044
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、所加窗函数的阶数P和在P至
Figure BDA0002703855010000045
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置Is代入到公式(4)得到最差杂散WoSpur,即杂散信号功率PSpur
Figure BDA0002703855010000046
式中,L为所加窗函数带来的功率损失百分比、P为所加窗函数的阶数、|Y(k)|2为运算后的第k点谱线的功率,Is为在P至
Figure BDA0002703855010000051
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置;
将所述用户设置的一路数字基带信号所需的FFT点数N和噪底指标PNoise_floor代入公式(5),得到噪声总功率PNoise,将所述实测信号功率PMain和噪声总功率PNoise代入公式(6),得到信噪比SNR;
Figure BDA0002703855010000052
式中,N为用户设置的一路数字基带信号所需的FFT点数,PNoise_floor为噪底指标;
Figure BDA0002703855010000053
式中,PMain为实测信号功率,PNoise为噪声总功率;
将所述所加窗函数带来的功率损失百分比L、在0至
Figure BDA0002703855010000054
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置I1、运算后的第k点谱线的功率|Y(k)|2和所加窗函数的阶数P代入公式(7),得到所述二次至六次谐波总功率PSumHar;将所述实测信号功率PMain,噪声总功率PNoise和二次至六次谐波总功率PSumHar代入公式(8),得到信纳比SINAD;
Figure BDA0002703855010000055
式中,L为所加窗函数带来的功率损失百分比、I1为在0至
Figure BDA0002703855010000061
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率;
Figure BDA0002703855010000062
式中,PMain为实测信号功率,PNoise为噪声总功率,PSumHar为二次至六次谐波总功率;
将所述信纳比SINAD代入公式(9)得到有效位ENOB,
Figure BDA0002703855010000063
式中,SINAD为信纳比;
将所述二次至六次谐波总功率PSumHar代入公式(10)得到二次至六次谐波总功率的平均值
Figure BDA0002703855010000064
将所述实测信号功率PMain、二次至六次谐波总功率的平均值
Figure BDA0002703855010000065
和杂散信号功率PSpur代入公式(11),得到无杂散动态范围SFDR;
Figure BDA0002703855010000066
式中,PSumHar为二次至六次谐波总功率;
Figure BDA0002703855010000067
式中,PMain为实测信号功率、
Figure BDA0002703855010000068
为二次至六次谐波总功率的平均值,PSpur为杂散信号功率;
将所述二次至六次谐波总功率PSumHar和实测信号功率PMain代入公式(12),得到总谐波失真THD;
Figure BDA0002703855010000071
式中,PSumHar为二次至六次谐波总功率,PMain为实测信号功率。
在一个具体实施方式中,所述FPGA核心控制器包括:高速数据收发模块、DDR控制模块、收发控制模块、总线接口模块、发送缓存模块、触发管理模块和指示灯管理模块;
所述高速数据收发模块的接收部分用于接收所述第一至第四光电收发器输出的串行电信号,完成串-并转换及解码,再将并行接收数据发送至所述DDR控制模块进行缓存;所述高速数据收发模块的发送部分用于将所述发送缓存模块缓存的并行发送数据进行编码和并-串转换,再将串行发送数据通过第一至第四光电收发器发送出去;
所述DDR控制模块用于对所述存储器读写控制,完成并行接收数据缓存,将所述并行接收数据发送给所述收发控制模块;
所述收发控制模块用于完成所述总线接口模块和DDR控制模块之间的数据缓冲,接收计算机设备通过所述总线接口模块下发的数据,通过所述DDR控制模块,从所述存储器中读取数据,并将数据通过所述总线接口模块发送到计算机设备;
发送缓存模块在所述所述收发控制模块的控制下按照外部同步信号的节拍发送数据,用于缓存等待发送的并行数据;
触发管理模块用于接收所述隔离后的同步信号,并将所述同步触发信号发送给所述收发控制模块;
指示灯管理模块在所述收发控制模块的指示下对所述第一四通道指示灯和第二四通道指示灯的工作状态进行控制。
在一个具体实施方式中,所述第k点谱线功率|Y(k)|2,k=0,…,N-1。
在一个具体实施方式中,所述计算机设备还包括:
人机交互模块、总线驱动模块、逻辑控制及数据存储模块和显示器;
所述计算机设备的人机交互模块响应于用户的输入而生成控制指令;
所述总线驱动模块用于接收所述总线接口模块发送的数据,并将数据发送到所述逻辑控制及数据存储模块;
所述逻辑控制及数据存储模块用于对接收到的数据进行存储;
所述显示器用于将所述计算机设备分析后的数据以图象的方式显示出来。
在一个具体实施方式中,所述总线接口模块为:PCIe总线接口模块、PCI总线接口模块、USB总线接口模块或网络总线接口模块。
第二方面,本发明提供一种进行数字基带信号录取分析的方法,包括:
所述计算机设备响应于用户的输入而生成控制指令;
所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;
所述计算机设备对所述数据进行分析。
在一个具体实施方式中,所述数据录取装置包括:FPGA核心控制器、同步触发连接器、光电收发器模块、第一四通道指示灯、第二四通道指示灯、隔离器和存储器;
其中,
所述光电收发器模块包括:
第一光电收发器、第二光电收发器、第三光电收发器和第四光电收发器;
所述方法还包括:
所述光电收发器模块用于接收外部被测接收机输出的光信号并将所述光信号转换为串行电信号,或将FPGA核心控制器输出的串行电信号转换为光信号;
所述同步触发连接器用于接收外部接收机发出的同步信号;
所述FPGA核心控制器用于接收所述同步信号和所述串行电信号,并对所述同步信号和所述串行电信号进行处理,将处理后得到的数据发送给计算机设备;
所述第一四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于发送状态;
所述第二四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于接收链路成功连接状态;
所述隔离器用于对所述同步触发连接器输出的同步信号进行光耦隔离,并将所述进行隔离后的同步信号输出至FPGA核心控制器;
所述存储器在所述FPGA核心控制器的控制下进行数据的读写,完成并行接收数据缓存。
在一个具体实施方式中,所述计算机设备包括指标分析模块;
所述方法还包括:
所述指标分析模块被配置为:
对接收机打包的信号进行拆分,获取用户所选的一路数字基带信号的数据;
按照用户所设置的一路数字基带信号所需的FFT点数N对于所述数据进行加窗和FFT运算,以获得功率谱线,其中,N为用户设置的一路数字基带信号所需的FFT点数,记|Y(k)|2为运算后的第k点谱线功率,P为所加窗函数的阶数,L为所加窗函数带来的功率损失百分比,S为杂散数和H为谐波总数;
在0至
Figure BDA0002703855010000101
数字基带信号范围内对于所述谱线功率|Y(k)|2求最大值,得到最大值对应的位置为基波信号的位置I1,基波信号取值范围B1为I1-P至I1+P;二次谐波取值范围B2为2×I1-P至2×I1+P,同理可得到三次至六次谐波的取值范围B3~B6;在P至
Figure BDA0002703855010000102
且不包含B1~B6的数字基带信号范围内,对于|Y(k)|2求最大值,得到最大值对应的位置为最大杂散信号的位置Is,最大杂散信号B7的取值范围为Is-P至Is+P;
当最大杂散的谱线功率小于二次至六次谐波谱线功率最小值时,所述杂散数S等于0和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(1),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000103
当最大杂散的谱线功率大于或者等于二次至六次谐波谱线功率最小值时,所述杂散数S等于1和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(2),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000111
公式(1)和(2)中,N为用户设置的一路数字基带信号所需的FFT点数,|Y(k)|2为运算后的第k点谱线的功率,P为所加窗函数的阶数,S为杂散数和H为谐波总数;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、噪底指标PNoise_floor、在0至
Figure BDA0002703855010000112
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置I1和所加窗函数的阶数P代入公式(3),得到实测信号功率PMain
Figure BDA0002703855010000113
式中,L为所加窗函数带来的功率损失百分比、PNoise_floor为噪底指标,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率,I1为在0至
Figure BDA0002703855010000114
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、所加窗函数的阶数P和在P至
Figure BDA0002703855010000115
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置Is代入到公式(4)得到最差杂散WoSpur,即杂散信号功率PSpur
Figure BDA0002703855010000116
式中,L为所加窗函数带来的功率损失百分比、P为所加窗函数的阶数、|Y(k)|2为运算后的第k点谱线的功率,Is为在P至
Figure BDA0002703855010000121
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置;
将所述用户设置的一路数字基带信号所需的FFT点数N和噪底指标PNoise_floor代入公式(5),得到噪声总功率PNoise,将所述实测信号功率PMain和噪声总功率PNoise代入公式(6),得到信噪比SNR;
Figure BDA0002703855010000122
式中,N为用户设置的一路数字基带信号所需的FFT点数,PNoise_floor为噪底指标;
Figure BDA0002703855010000123
式中,PMain为实测信号功率,PNoise为噪声总功率;
将所述所加窗函数带来的功率损失百分比L、在0至
Figure BDA0002703855010000124
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置I1、运算后的第k点谱线的功率|Y(k)|2和所加窗函数的阶数P代入公式(7),得到所述二次至六次谐波总功率PSumHar;将所述实测信号功率PMain,噪声总功率PNoise和二次至六次谐波总功率PSumHar代入公式(8),得到信纳比SINAD;
Figure BDA0002703855010000125
式中,L为所加窗函数带来的功率损失百分比、I1为在0至
Figure BDA0002703855010000131
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率;
Figure BDA0002703855010000132
式中,PMain为实测信号功率,PNoise为噪声总功率,PSumHar为二次至六次谐波总功率;
将所述信纳比SINAD代入公式(9)得到有效位ENOB,
Figure BDA0002703855010000133
式中,SINAD为信纳比;
将所述二次至六次谐波总功率PSumHar代入公式(10)得到二次至六次谐波总功率的平均值
Figure BDA0002703855010000134
将所述实测信号功率PMain、二次至六次谐波总功率的平均值
Figure BDA0002703855010000135
和杂散信号功率PSpur代入公式(11),得到无杂散动态范围SFDR;
Figure BDA0002703855010000136
式中,PSumHar为二次至六次谐波总功率;
Figure BDA0002703855010000137
式中,PMain为实测信号功率、
Figure BDA0002703855010000138
为二次至六次谐波总功率的平均值,PSpur为杂散信号功率;
将所述二次至六次谐波总功率PSumHar和实测信号功率PMain代入公式(12),得到总谐波失真THD;
Figure BDA0002703855010000141
式中,PSumHar为二次至六次谐波总功率,PMain为实测信号功率。
本发明的有益效果如下:
本发明提供了一种雷达接收机数字基带信号录取分析系统和方法,可在测试过程中完成雷达接收机数字基带信号的接收、录取和处理分析,为雷达接收机指标的评估提供有效依据。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出根据本发明一个实施例的雷达接收机数字基带信号录取分析系统的示意图。
图2示出根据本发明一个实施例的数字基带信号录取分析系统主界面截图。
图3示出根据本发明一个实施例的数字基带信号录取分析系统实时分析界面截图。
图4示出根据本发明一个实施例的数字基带信号录取分析系统数据回放界面截图。
图5示出利用本申请的一个实施例的数字基带信号录取分析系统进行数字基带信号录取分析的方法流程图。
具体实施方式
为使本发明的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
第一实施例
本发明公开了一种雷达接收机数字基带信号录取分析系统,该系统的示意图如图1所示,该系统包括:计算机设备和数据录取装置;
所述计算机设备响应于用户的输入而生成控制指令;
所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;
所述计算机设备对所述数据进行分析。
在一个具体实施例中,所述数据录取装置包括:
FPGA核心控制器、同步触发连接器、光电收发器模块、第一四通道指示灯、第二四通道指示灯、隔离器和存储器;
其中,
所述光电收发器模块包括:
第一光电收发器、第二光电收发器、第三光电收发器和第四光电收发器;
所述光电收发器模块用于接收外部被测接收机输出的光信号并将所述光信号转换为串行电信号,或将FPGA核心控制器输出的串行电信号转换为光信号;
所述同步触发连接器用于接收外部接收机发出的同步信号;
所述FPGA核心控制器用于接收所述同步信号和所述串行电信号,并对所述同步信号和所述串行电信号进行处理,将处理后得到的数据发送给计算机设备;
所述第一四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于发送状态;
所述第二四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于接收链路成功连接状态;
所述隔离器用于对所述同步触发连接器输出的同步信号进行光耦隔离,对被测设备进行有效的防护,防止外部接收机电源和本发明所述系统电源之间的相互干扰,并将所述进行隔离后的同步信号输出至FPGA核心控制器;
所述存储器在所述FPGA核心控制器的控制下进行数据的读写,完成并行接收数据缓存。
其中,在一个具体实施例中,所述FPGA核心控制器包括:
高速数据收发模块、DDR控制模块、收发控制模块、总线接口模块、发送缓存模块、触发管理模块和指示灯管理模块;
所述高速数据收发模块的接收部分用于接收所述第一至第四光电收发器输出的串行电信号,完成串-并转换及解码,再将并行接收数据发送至所述DDR控制模块进行缓存;所述高速数据收发模块的发送部分用于将所述发送缓存模块缓存的并行发送数据进行编码和并-串转换,再将串行发送数据通过第一至第四光电收发器发送出去;
所述DDR控制模块用于对所述存储器读写控制,完成并行接收数据缓存,将所述并行接收数据发送给所述收发控制模块;
所述收发控制模块用于完成所述总线接口模块和DDR控制模块之间的数据缓冲,接收计算机设备通过所述总线接口模块下发的数据,通过所述DDR控制模块,从所述存储器中读取数据,并将数据通过所述总线接口模块发送到计算机设备;
发送缓存模块在所述所述收发控制模块的控制下按照外部同步信号的节拍发送数据,用于缓存等待发送的并行数据;
触发管理模块用于接收所述隔离后的同步信号,并将所述同步触发信号发送给所述收发控制模块;
指示灯管理模块在所述收发控制模块的指示下对所述第一四通道指示灯和第二四通道指示灯的工作状态进行控制;
总线接口模块用于实现FPGA核心控制器和计算机设备之间的复杂通信。
在一个具体实施例中,总线接口模块主要为:PCIe总线接口模块、PCI总线接口模块、USB总线接口模块或网络总线接口模块。
在一个具体实施例中,计算机设备包括:总线驱动模块、逻辑控制及数据存储模块、指标分析模块、人机交互模块和显示器。
在一个具体实施例中,同步触发连接器输出端与隔离器输入端相连,隔离器输出端与FPGA核心控制器中的触发管理模块输入端相连,FPGA核心控制器中的指示灯管理模块分别于第一四通道指示灯、第二四通道指示灯输入端相连,FPGA核心控制器中的高速数据收发模块与第一至第四光电收发器双向连接,FPGA核心控制器中DDR控制模块与存储器双向连接,FPGA核心控制器中的总线接口模块与计算机设备双向连接。
在一个具体实施例中,FPGA核心控制器中,触发管理模块的同步触发信号输出端与收发控制模块的同步触发信号输入端连接。收发控制模块的状态指示输出端与指示灯管理模块的状态指示输入端连接。收发控制模块的并行发送数据输出端与发送缓存模块的并行发送数据输入端连接,发送缓存模块的并行发送数据输出端与高速数据收发模块的并行发送数据输入端连接。高速数据收发模块的并行接收数据输出端与DDR控制模块的并行接收数据输入端连接,DDR控制模块的并行接收数据输出端与收发控制模块的并行接收数据输入端连接,收发控制模块的接收缓存控制信号与DDR控制模块的接收缓存控制信号双向连接。收发控制模块的控制数据接口与总线接口模块的控制数据接口双向连接,总线接口模块的总线数据接口与计算机设备中的总线驱动模块双向连接。
在一个具体实施例中,所述计算机设备的人机交互模块响应于用户的输入而生成控制指令;所述总线驱动模块用于接收所述总线接口模块发送的数据,并将数据发送到所述逻辑控制及数据存储模块;所述逻辑控制及数据存储模块用于对接收到的数据进行存储;指标分析模块用于对采集的数据进行指标分析;所述显示器用于将所述计算机设备分析后的数据以图象的方式显示出来。
在一个具体实施例中,所述指标分析模块被配置为:
对接收机打包的信号进行拆分,获取用户所选的一路数字基带信号的数据;
按照用户所设置的一路数字基带信号所需的FFT点数N对于所述数据进行加窗和FFT运算,以获得功率谱线,其中,N为用户设置的一路数字基带信号所需的FFT点数,记|Y(k)|2为运算后的第k(k=0,…,N-1)点谱线功率,P为所加窗函数的阶数,L为所加窗函数带来的功率损失百分比,S为杂散数和H为谐波总数;
在0至
Figure BDA0002703855010000181
数字基带信号范围内对于所述谱线功率|Y(k)|2求最大值,得到最大值对应的位置为基波信号的位置I1,基波信号取值范围B1为I1-P至I1+P;二次谐波取值范围B2为2×I1-P至2×I1+P,同理可得到三次至六次谐波的取值范围B3~B6;在P至
Figure BDA0002703855010000182
且不包含B1~B6的数字基带信号范围内,对于|Y(k)|2求最大值,得到最大值对应的位置为最大杂散信号的位置Is,最大杂散信号B7的取值范围为Is-P至Is+P;
当最大杂散的谱线功率小于二次至六次谐波谱线功率最小值时,所述杂散数S等于0和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(1),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000191
当最大杂散的谱线功率大于或者等于二次至六次谐波谱线功率最小值时,所述杂散数S等于1和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(2),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000192
公式(1)和(2)中,N为用户设置的一路数字基带信号所需的FFT点数,|Y(k)|2为运算后的第k点谱线的功率,P为所加窗函数的阶数,S为杂散数和H为谐波总数;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、噪底指标PNoise_floor、在0至
Figure BDA0002703855010000193
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置I1和所加窗函数的阶数P代入公式(3),得到实测信号功率PMain
Figure BDA0002703855010000201
式中,L为所加窗函数带来的功率损失百分比、PNoise_floor为噪底指标,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率,I1为在0至
Figure BDA0002703855010000202
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、所加窗函数的阶数P和在P至
Figure BDA0002703855010000203
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置Is代入到公式(4)得到最差杂散WoSpur,即杂散信号功率PSpur
Figure BDA0002703855010000204
式中,L为所加窗函数带来的功率损失百分比、P为所加窗函数的阶数、|Y(k)|2为运算后的第k点谱线的功率,Is为在P至
Figure BDA0002703855010000205
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置;
将所述用户设置的一路数字基带信号所需的FFT点数N和噪底指标PNoise_floor代入公式(5),得到噪声总功率PNoise,将所述实测信号功率PMain和噪声总功率PNoise代入公式(6),得到信噪比SNR;
Figure BDA0002703855010000206
式中,N为用户设置的一路数字基带信号所需的FFT点数,PNoise_floor为噪底指标;
Figure BDA0002703855010000211
式中,PMain为实测信号功率,PNoise为噪声总功率;
将所述所加窗函数带来的功率损失百分比L、在0至
Figure BDA0002703855010000212
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置I1、运算后的第k点谱线的功率|Y(k)|2和所加窗函数的阶数P代入公式(7),得到所述二次至六次谐波总功率PSumHar;将所述实测信号功率PMain,噪声总功率PNoise和二次至六次谐波总功率PSumHar代入公式(8),得到信纳比SINAD;
Figure BDA0002703855010000213
式中,L为所加窗函数带来的功率损失百分比、I1为在0至
Figure BDA0002703855010000214
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率;
Figure BDA0002703855010000215
式中,PMain为实测信号功率,PNoise为噪声总功率,PSumHar为二次至六次谐波总功率;
将所述信纳比SINAD代入公式(9)得到有效位ENOB,
Figure BDA0002703855010000221
式中,SINAD为信纳比;
将所述二次至六次谐波总功率PSumHar代入公式(10)得到二次至六次谐波总功率的平均值
Figure BDA0002703855010000222
将所述实测信号功率PMain、二次至六次谐波总功率的平均值
Figure BDA0002703855010000223
和杂散信号功率PSpur代入公式(11),得到无杂散动态范围SFDR;
Figure BDA0002703855010000224
式中,PSumHar为二次至六次谐波总功率;
Figure BDA0002703855010000225
式中,PMain为实测信号功率、
Figure BDA0002703855010000226
为二次至六次谐波总功率的平均值,PSpur为杂散信号功率;
将所述二次至六次谐波总功率PSumHar和实测信号功率PMain代入公式(12),得到总谐波失真THD;
Figure BDA0002703855010000227
式中,PSumHar为二次至六次谐波总功率,PMain为实测信号功率。
在一个具体实施例中,计算机设备中,总线驱动模块、逻辑控制及数据存储模块、指标分析模块和人机交互模块依次双向连接,人机交互模块与显示器双向连接。
本发明公开了一种雷达接收机数字基带信号录取分析系统,解决了传统的台式仪器设备不具备光纤接口、且不具备基带数据处理能力,无法通过数字基带分析方法完成接收机性能指标测试的问题。
第二实施例
如图5所示,本发明公开了一种进行数字基带信号录取分析的方法,包括:
所述计算机设备响应于用户的输入而生成控制指令;所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;所述计算机设备对所述数据进行分析。
在一个具体实施例中,
所述数据录取装置包括:FPGA核心控制器、同步触发连接器、光电收发器模块、第一四通道指示灯、第二四通道指示灯、隔离器和存储器;
其中,
所述光电收发器模块包括:
第一光电收发器、第二光电收发器、第三光电收发器和第四光电收发器;
所述方法还包括:
所述光电收发器模块用于接收外部被测接收机输出的光信号并将所述光信号转换为串行电信号,或将FPGA核心控制器输出的串行电信号转换为光信号;
所述同步触发连接器用于接收外部接收机发出的同步信号;
所述FPGA核心控制器用于接收所述同步信号和所述串行电信号,并对所述同步信号和所述串行电信号进行处理,将处理后得到的数据发送给计算机设备;
所述第一四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于发送状态;
所述第二四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于接收链路成功连接状态;
所述隔离器用于对所述同步触发连接器输出的同步信号进行光耦隔离,并将所述进行隔离后的同步信号输出至FPGA核心控制器;
所述存储器在所述FPGA核心控制器的控制下进行数据的读写,完成并行接收数据缓存。
在一个具体实施例中,进行雷达接收机数字基带信号录取分析方法的流程如下:
被测接收机与本发明所述系统通过完成光纤、同步触发信号连接并正常通电工作后,链路连接成功的光电收发器所对应的第二四通道指示灯点亮,打开图2所示的数字基带信号录取分析系统主界面软件,选择发送或接收光纤通道、设置同步和定时方式后,可完本发明所述系统的功能控制。本发明所述系统包含两种指标分析功能,分别为实时分析和数据回放,分别如图3和图4所示。实时分析功能为将接收到的数据写入计算机内存后直接进行数据指标分析;而数据回放则需要在图4中选择需要回访的文件名称和路径,再将文件读入计算机内存后进行数据指标分析。
所述计算机设备响应于用户的输入而生成控制指令,在一个具体实施例中,实时分析功能的实现,用户需要完成如图2所示的界面截屏中采样频率、有效位数、满量程幅度、码制、FFT点数、光纤通道(对应第一至第四光电收发器接收通道)、信号类型、信号数(接收机通常会将几个接收通道的数据打包在一个光纤通道上,本选项用于指示每个光电收发器所接收到的基带数据具体打包了几个接收通道的数字基带信号)、第几路信号(用于信号数中的一个通道完成信号的时域、频域波形现实,并进行指标分析)、窗选择的设置,然后单击“显示”按钮,即可开始对接收到的数据进行分析;如需要停止数据分析,则单击“停止”按钮;如需要退出实施分析功能,则单击“退出”按钮。
在一个具体实施例中,在进行指标分析前,所述计算机设备包括指标分析模块;
所述方法还包括:所述指标分析模块被配置为:
按照光纤通道、信号数和第几路信号对接收机打包的信号进行拆分,获取用户所选的一路数字基带信号的数据;
按照用户所设置的一路数字基带信号所需的FFT点数N对于所述数据进行加窗和FFT运算,以获得功率谱线,其中,N为用户设置的一路数字基带信号所需的FFT点数,记|Y(k)|2为运算后的第k(k=0,…,N-1)点谱线功率,P为所加窗函数的阶数,L为所加窗函数带来的功率损失百分比,S为杂散数和H为谐波总数;
在一个具体实施例中,噪底:在0至
Figure BDA0002703855010000251
数字基带信号范围内对于|Y(k)|2求最大值,最大值对应的位置为基波信号的位置I1,则基波信号取值范围波B1为I1-P至I1+P;二次谐波取值范围B2为2×I1-P至2×I1+P,同理可得到三次至六次谐波的取值范围B3~B6;在P至
Figure BDA0002703855010000252
且不包含B1~B6的数字基带信号范围内,对于|Y(k)|2求最大值,最大值对应的位置为最大杂散信号的位置Is,则最大杂散信号B7的取值范围为Is-P至Is+P,可以按照后续的计算公式完成噪底指标PNoise_floor的分析计算。
其中,在一个具体实施例中,当最大杂散的谱线功率小于二次至六次谐波谱线功率最小值时,所述杂散数S等于0和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(1),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000261
反之,当最大杂散的谱线功率大于或者等于二次至六次谐波谱线功率最小值时,所述杂散数S等于1和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(2),对噪底指标PNoise_floor进行分析计算;
Figure BDA0002703855010000262
公式(1)和(2)中,N为用户设置的一路数字基带信号所需的FFT点数,|Y(k)|2为运算后的第k点谱线的功率,P为所加窗函数的阶数,S为杂散数和H为谐波总数;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、噪底指标PNoise_floor、在0至
Figure BDA0002703855010000263
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置I1和所加窗函数的阶数P代入公式(3),得到实测信号功率PMain
Figure BDA0002703855010000264
式中,L为所加窗函数带来的功率损失百分比、PNoise_floor为噪底指标,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率,I1为在0至
Figure BDA0002703855010000271
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、所加窗函数的阶数P和在P至
Figure BDA0002703855010000272
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置Is代入到公式(4)得到最差杂散WoSpur,即杂散信号功率PSpur
Figure BDA0002703855010000273
式中,L为所加窗函数带来的功率损失百分比、P为所加窗函数的阶数、|Y(k)|2为运算后的第k点谱线的功率,Is为在P至
Figure BDA0002703855010000274
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置;
将所述用户设置的一路数字基带信号所需的FFT点数N和噪底指标PNoise_floor代入公式(5),得到噪声总功率PNoise,将所述实测信号功率PMain和噪声总功率PNoise代入公式(6),得到信噪比SNR;
Figure BDA0002703855010000275
式中,N为用户设置的一路数字基带信号所需的FFT点数,PNoise_floor为噪底指标;
Figure BDA0002703855010000281
式中,PMain为实测信号功率,PNoise为噪声总功率;
将所述所加窗函数带来的功率损失百分比L、在0至
Figure BDA0002703855010000282
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置I1、运算后的第k点谱线的功率|Y(k)|2和所加窗函数的阶数P代入公式(7),得到所述二次至六次谐波总功率PSumHar;将所述实测信号功率PMain,噪声总功率PNoise和二次至六次谐波总功率PSumHar代入公式(8),得到信纳比SINAD;
Figure BDA0002703855010000283
式中,L为所加窗函数带来的功率损失百分比、I1为在0至
Figure BDA0002703855010000284
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率;
Figure BDA0002703855010000285
式中,PMain为实测信号功率,PNoise为噪声总功率,PSumHar为二次至六次谐波总功率;
将所述信纳比SINAD代入公式(9)得到有效位ENOB,
Figure BDA0002703855010000286
式中,SINAD为信纳比;
将所述二次至六次谐波总功率PSumHar代入公式(10)得到二次至六次谐波总功率的平均值
Figure BDA0002703855010000291
将所述实测信号功率PMain、二次至六次谐波总功率的平均值
Figure BDA0002703855010000292
和杂散信号功率PSpur代入公式(11),得到无杂散动态范围SFDR;
Figure BDA0002703855010000293
式中,PSumHar为二次至六次谐波总功率;
Figure BDA0002703855010000294
式中,PMain为实测信号功率、
Figure BDA0002703855010000295
为二次至六次谐波总功率的平均值,PSpur为杂散信号功率;
将所述二次至六次谐波总功率PSumHar和实测信号功率PMain代入公式(12),得到总谐波失真THD;
Figure BDA0002703855010000296
式中,PSumHar为二次至六次谐波总功率,PMain为实测信号功率。
本发明公开了一种雷达接收机数字基带信号录取分析方法,解决了传统的台式仪器设备不具备光纤接口、且不具备基带数据处理能力,无法通过数字基带分析方法完成接收机性能指标测试的问题。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。

Claims (10)

1.一种雷达接收机数字基带信号录取分析系统,其特征在于,包括:
计算机设备和数据录取装置;
其中,所述计算机设备响应于用户的输入而生成控制指令;
所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;
所述计算机设备对所述数据进行分析。
2.根据权利要求1所述的系统,其特征在于,
所述数据录取装置包括:
FPGA核心控制器、同步触发连接器、光电收发器模块、第一四通道指示灯、第二四通道指示灯、隔离器和存储器;
其中,
所述光电收发器模块包括:
第一光电收发器、第二光电收发器、第三光电收发器和第四光电收发器;
所述光电收发器模块用于接收外部被测接收机输出的光信号并将所述光信号转换为串行电信号,或将FPGA核心控制器输出的串行电信号转换为光信号;
所述同步触发连接器用于接收外部接收机发出的同步信号;
所述FPGA核心控制器用于接收所述同步信号和所述串行电信号,并对所述同步信号和所述串行电信号进行处理,将处理后得到的数据发送给计算机设备;
所述第一四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于发送状态;
所述第二四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于接收链路成功连接状态;
所述隔离器用于对所述同步触发连接器输出的同步信号进行光耦隔离,并将所述进行隔离后的同步信号输出至FPGA核心控制器;
所述存储器在所述FPGA核心控制器的控制下进行数据的读写,完成并行接收数据缓存。
3.根据权利要求1所述的系统,其特征在于,
所述计算机设备包括指标分析模块;
所述指标分析模块被配置为:
对接收机打包的信号进行拆分,获取用户所选的一路数字基带信号的数据;
按照用户所设置的一路数字基带信号所需的FFT点数N对于所述数据进行加窗和FFT运算,以获得功率谱线,其中,N为用户设置的一路数字基带信号所需的FFT点数,记|Y(k)|2为运算后的第k点谱线功率,P为所加窗函数的阶数,L为所加窗函数带来的功率损失百分比,S为杂散数和H为谐波总数;
在0至
Figure FDA0002703853000000021
数字基带信号范围内对于所述谱线功率|Y(k)|2求最大值,得到最大值对应的位置为基波信号的位置I1,基波信号取值范围B1为I1-P至I1+P;二次谐波取值范围B2为2×I1-P至2×I1+P,同理可得到三次至六次谐波的取值范围B3~B6;在P至
Figure FDA0002703853000000022
且不包含B1~B6的数字基带信号范围内,对于|Y(k)|2求最大值,得到最大值对应的位置为最大杂散信号的位置Is,最大杂散信号B7的取值范围为Is-P至Is+P;
当最大杂散的谱线功率小于二次至六次谐波谱线功率最小值时,所述杂散数S等于0和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(1),对噪底指标PNoise_floor进行分析计算;
Figure FDA0002703853000000031
当最大杂散的谱线功率大于或者等于二次至六次谐波谱线功率最小值时,所述杂散数S等于1和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(2),对噪底指标PNoise_floor进行分析计算;
Figure FDA0002703853000000032
公式(1)和(2)中,N为用户设置的一路数字基带信号所需的FFT点数,|Y(k)|2为运算后的第k点谱线的功率,P为所加窗函数的阶数,S为杂散数和H为谐波总数;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、噪底指标PNoise_floor、在0至
Figure FDA0002703853000000033
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置I1和所加窗函数的阶数P代入公式(3),得到实测信号功率PMain
Figure FDA0002703853000000034
式中,L为所加窗函数带来的功率损失百分比、PNoise_floor为噪底指标,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率,I1为在0至
Figure FDA0002703853000000041
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、所加窗函数的阶数P和在P至
Figure FDA0002703853000000042
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置Is代入到公式(4)得到最差杂散WoSpur,即杂散信号功率PSpur
Figure FDA0002703853000000043
式中,L为所加窗函数带来的功率损失百分比、P为所加窗函数的阶数、|Y(k)|2为运算后的第k点谱线的功率,Is为在P至
Figure FDA0002703853000000044
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置;
将所述用户设置的一路数字基带信号所需的FFT点数N和噪底指标PNoise_floor代入公式(5),得到噪声总功率PNoise,将所述实测信号功率PMain和噪声总功率PNoise代入公式(6),得到信噪比SNR;
Figure FDA0002703853000000045
式中,N为用户设置的一路数字基带信号所需的FFT点数,PNoise_floor为噪底指标;
Figure FDA0002703853000000051
式中,PMain为实测信号功率,PNoise为噪声总功率;
将所述所加窗函数带来的功率损失百分比L、在0至
Figure FDA0002703853000000052
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置I1、运算后的第k点谱线的功率|Y(k)|2和所加窗函数的阶数P代入公式(7),得到所述二次至六次谐波总功率PSumHar;将所述实测信号功率PMain,噪声总功率PNoise和二次至六次谐波总功率PSumHar代入公式(8),得到信纳比SINAD;
Figure FDA0002703853000000053
式中,L为所加窗函数带来的功率损失百分比、I1为在0至
Figure FDA0002703853000000054
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率;
Figure FDA0002703853000000055
式中,PMain为实测信号功率,PNoise为噪声总功率,PSumHar为二次至六次谐波总功率;
将所述信纳比SINAD代入公式(9)得到有效位ENOB,
Figure FDA0002703853000000056
式中,SINAD为信纳比;
将所述二次至六次谐波总功率PSumHar代入公式(10)得到二次至六次谐波总功率的平均值
Figure FDA0002703853000000061
将所述实测信号功率PMain、二次至六次谐波总功率的平均值
Figure FDA0002703853000000062
和杂散信号功率PSpur代入公式(11),得到无杂散动态范围SFDR;
Figure FDA0002703853000000063
式中,PSumHar为二次至六次谐波总功率;
Figure FDA0002703853000000064
式中,PMain为实测信号功率、
Figure FDA0002703853000000065
为二次至六次谐波总功率的平均值,PSpur为杂散信号功率;
将所述二次至六次谐波总功率PSumHar和实测信号功率PMain代入公式(12),得到总谐波失真THD;
Figure FDA0002703853000000066
式中,PSumHar为二次至六次谐波总功率,PMain为实测信号功率。
4.根据权利要求2所述的系统,其特征在于,所述FPGA核心控制器包括:
高速数据收发模块、DDR控制模块、收发控制模块、总线接口模块、发送缓存模块、触发管理模块和指示灯管理模块;
所述高速数据收发模块的接收部分用于接收所述第一至第四光电收发器输出的串行电信号,完成串-并转换及解码,再将并行接收数据发送至所述DDR控制模块进行缓存;所述高速数据收发模块的发送部分用于将所述发送缓存模块缓存的并行发送数据进行编码和并-串转换,再将串行发送数据通过第一至第四光电收发器发送出去;
所述DDR控制模块用于对所述存储器读写控制,完成并行接收数据缓存,将所述并行接收数据发送给所述收发控制模块;
所述收发控制模块用于完成所述总线接口模块和DDR控制模块之间的数据缓冲,接收计算机设备通过所述总线接口模块下发的数据,通过所述DDR控制模块,从所述存储器中读取数据,并将数据通过所述总线接口模块发送到计算机设备;
发送缓存模块在所述所述收发控制模块的控制下按照外部同步信号的节拍发送数据,用于缓存等待发送的并行数据;
触发管理模块用于接收所述隔离后的同步信号,并将所述同步触发信号发送给所述收发控制模块;
指示灯管理模块在所述收发控制模块的指示下对所述第一四通道指示灯和第二四通道指示灯的工作状态进行控制。
5.根据权利要求3所述的系统,其特征在于,
所述第k点谱线功率|Y(k)|2,k=0,…,N-1。
6.根据权利要求3所述的系统,其特征在于,
所述计算机设备还包括:
人机交互模块、总线驱动模块、逻辑控制及数据存储模块和显示器;
所述计算机设备的人机交互模块响应于用户的输入而生成控制指令;
所述总线驱动模块用于接收所述总线接口模块发送的数据,并将数据发送到所述逻辑控制及数据存储模块;
所述逻辑控制及数据存储模块用于对接收到的数据进行存储;
所述显示器用于将所述计算机设备分析后的数据以图象的方式显示出来。
7.根据权利要求4所述的系统,其特征在于,
所述总线接口模块为:PCIe总线接口模块、PCI总线接口模块、USB总线接口模块或网络总线接口模块。
8.一种利用权利要求1-7中任一项所述的系统进行数字基带信号录取分析的方法,其特征在于,包括:
所述计算机设备响应于用户的输入而生成控制指令;
所述数据录取装置根据所述控制指令采集外部被测接收机输出的光信号以及同步信号,并生成数据;
所述计算机设备对所述数据进行分析。
9.根据权利要求8所述的方法,其特征在于,
所述数据录取装置包括:FPGA核心控制器、同步触发连接器、光电收发器模块、第一四通道指示灯、第二四通道指示灯、隔离器和存储器;
其中,
所述光电收发器模块包括:
第一光电收发器、第二光电收发器、第三光电收发器和第四光电收发器;
所述方法还包括:
所述光电收发器模块用于接收外部被测接收机输出的光信号并将所述光信号转换为串行电信号,或将FPGA核心控制器输出的串行电信号转换为光信号;
所述同步触发连接器用于接收外部接收机发出的同步信号;
所述FPGA核心控制器用于接收所述同步信号和所述串行电信号,并对所述同步信号和所述串行电信号进行处理,将处理后得到的数据发送给计算机设备;
所述第一四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于发送状态;
所述第二四通道指示灯用于指示所述FPGA核心控制器控制下的第一至第四光电收发器是否处于接收链路成功连接状态;
所述隔离器用于对所述同步触发连接器输出的同步信号进行光耦隔离,并将所述进行隔离后的同步信号输出至FPGA核心控制器;
所述存储器在所述FPGA核心控制器的控制下进行数据的读写,完成并行接收数据缓存。
10.根据权利要求9所述的方法,其特征在于,
所述计算机设备包括指标分析模块;
所述方法还包括:
所述指标分析模块被配置为:
对接收机打包的信号进行拆分,获取用户所选的一路数字基带信号的数据;
按照用户所设置的一路数字基带信号所需的FFT点数N对于所述数据进行加窗和FFT运算,以获得功率谱线,其中,N为用户设置的一路数字基带信号所需的FFT点数,记|Y(k)|2为运算后的第k点谱线功率,P为所加窗函数的阶数,L为所加窗函数带来的功率损失百分比,S为杂散数和H为谐波总数;
在0至
Figure FDA0002703853000000091
数字基带信号范围内对于所述谱线功率|Y(k)|2求最大值,得到最大值对应的位置为基波信号的位置I1,基波信号取值范围B1为I1-P至I1+P;二次谐波取值范围B2为2×I1-P至2×I1+P,同理可得到三次至六次谐波的取值范围B3~B6;在P至
Figure FDA0002703853000000092
且不包含B1~B6的数字基带信号范围内,对于|Y(k)|2求最大值,得到最大值对应的位置为最大杂散信号的位置Is,最大杂散信号B7的取值范围为Is-P至Is+P;
当最大杂散的谱线功率小于二次至六次谐波谱线功率最小值时,所述杂散数S等于0和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(1),对噪底指标PNoise_floor进行分析计算;
Figure FDA0002703853000000101
当最大杂散的谱线功率大于或者等于二次至六次谐波谱线功率最小值时,所述杂散数S等于1和谐波总数H等于5:
将所述用户设置的一路数字基带信号所需的FFT点数N,运算后的第k点谱线功率|Y(k)|2,所加窗函数的阶数P,杂散数S,谐波总数H代入公式(2),对噪底指标PNoise_floor进行分析计算;
Figure FDA0002703853000000102
公式(1)和(2)中,N为用户设置的一路数字基带信号所需的FFT点数,|Y(k)|2为运算后的第k点谱线的功率,P为所加窗函数的阶数,S为杂散数和H为谐波总数;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、噪底指标PNoise_floor、在0至
Figure FDA0002703853000000111
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置I1和所加窗函数的阶数P代入公式(3),得到实测信号功率PMain
Figure FDA0002703853000000112
式中,L为所加窗函数带来的功率损失百分比、PNoise_floor为噪底指标,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率,I1为在0至
Figure FDA0002703853000000113
数字基带信号范围内|Y(k)|2最大值对应的基波信号的位置;
将所述所加窗函数带来的功率损失百分比L、运算后的第k点谱线的功率|Y(k)|2、所加窗函数的阶数P和在P至
Figure FDA0002703853000000114
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置Is代入到公式(4)得到最差杂散WoSpur,即杂散信号功率PSpur
Figure FDA0002703853000000115
式中,L为所加窗函数带来的功率损失百分比、P为所加窗函数的阶数、|Y(k)|2为运算后的第k点谱线的功率,Is为在P至
Figure FDA0002703853000000116
且不包含B1~B6的数字基带信号范围内,|Y(k)|2最大值对应的位置为最大杂散信号的位置;
将所述用户设置的一路数字基带信号所需的FFT点数N和噪底指标PNoise_floor代入公式(5),得到噪声总功率PNoise,将所述实测信号功率PMain和噪声总功率PNoise代入公式(6),得到信噪比SNR;
Figure FDA0002703853000000121
式中,N为用户设置的一路数字基带信号所需的FFT点数,PNoise_floor为噪底指标;
Figure FDA0002703853000000122
式中,PMain为实测信号功率,PNoise为噪声总功率;
将所述所加窗函数带来的功率损失百分比L、在0至
Figure FDA0002703853000000123
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置I1、运算后的第k点谱线的功率|Y(k)|2和所加窗函数的阶数P代入公式(7),得到所述二次至六次谐波总功率PSumHar;将所述实测信号功率PMain,噪声总功率PNoise和二次至六次谐波总功率PSumHar代入公式(8),得到信纳比SINAD;
Figure FDA0002703853000000124
式中,L为所加窗函数带来的功率损失百分比、I1为在0至
Figure FDA0002703853000000125
数字基带信号范围内|Y(k)|2最大值对应基波信号的位置,P为所加窗函数的阶数,|Y(k)|2为运算后的第k点谱线的功率;
Figure FDA0002703853000000126
式中,PMain为实测信号功率,PNoise为噪声总功率,PSumHar为二次至六次谐波总功率;
将所述信纳比SINAD代入公式(9)得到有效位ENOB,
Figure FDA0002703853000000131
式中,SINAD为信纳比;
将所述二次至六次谐波总功率PSumHar代入公式(10)得到二次至六次谐波总功率的平均值
Figure FDA0002703853000000132
将所述实测信号功率PMain、二次至六次谐波总功率的平均值
Figure FDA0002703853000000133
和杂散信号功率PSpur代入公式(11),得到无杂散动态范围SFDR;
Figure FDA0002703853000000134
式中,PSumHar为二次至六次谐波总功率;
Figure FDA0002703853000000135
式中,PMain为实测信号功率、
Figure FDA0002703853000000136
为二次至六次谐波总功率的平均值,PSpur为杂散信号功率;
将所述二次至六次谐波总功率PSumHar和实测信号功率PMain代入公式(12),得到总谐波失真THD;
Figure FDA0002703853000000137
式中,PSumHar为二次至六次谐波总功率,PMain为实测信号功率。
CN202011031659.5A 2020-09-27 2020-09-27 一种雷达接收机数字基带信号录取分析系统和方法 Active CN112305511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011031659.5A CN112305511B (zh) 2020-09-27 2020-09-27 一种雷达接收机数字基带信号录取分析系统和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011031659.5A CN112305511B (zh) 2020-09-27 2020-09-27 一种雷达接收机数字基带信号录取分析系统和方法

Publications (2)

Publication Number Publication Date
CN112305511A true CN112305511A (zh) 2021-02-02
CN112305511B CN112305511B (zh) 2024-05-07

Family

ID=74488498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011031659.5A Active CN112305511B (zh) 2020-09-27 2020-09-27 一种雷达接收机数字基带信号录取分析系统和方法

Country Status (1)

Country Link
CN (1) CN112305511B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113805148A (zh) * 2021-08-18 2021-12-17 北京无线电测量研究所 一种多功能数字tr组件光纤测试装置、系统和方法
CN116955258A (zh) * 2023-09-19 2023-10-27 成都立思方信息技术有限公司 一种可灵活连接的触发连接器、信号采集控制设备及系统

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201732160U (zh) * 2010-05-05 2011-02-02 北京航空航天大学 一种双接口雷达数据记录仪
CN102170318A (zh) * 2011-01-11 2011-08-31 南京国睿安泰信科技股份有限公司 用于接收机性能测试的频谱分析算法
US20130016761A1 (en) * 2011-07-15 2013-01-17 Renesas Mobile Corporation Transmitter with a Variable Sampling Rate
CN105676197A (zh) * 2016-03-30 2016-06-15 西安电子科技大学 一种阵列雷达回波模拟器及其模拟方法
CN105785335A (zh) * 2016-03-28 2016-07-20 电子科技大学 一种基于cPCI的数字阵接收通道性能自动测试系统
CN106406751A (zh) * 2016-06-03 2017-02-15 中北大学 基于多通道lvds接口的高速大容量i/q数据记录仪
CN108051786A (zh) * 2017-10-30 2018-05-18 北京航天福道高技术股份有限公司 一种宽带目标模拟器验证平台及验证方法
CN108333566A (zh) * 2018-02-07 2018-07-27 哈尔滨工业大学 一种基于zynq的便携式雷达测试系统及测试方法
CN110632568A (zh) * 2019-11-05 2019-12-31 中国科学院电子学研究所 合成孔径雷达实时成像处理器的测试信号源
CN111505593A (zh) * 2020-04-30 2020-08-07 北京无线电测量研究所 一种频综综合测试系统及测试方法
CN111562552A (zh) * 2020-04-30 2020-08-21 北京无线电测量研究所 一种地面雷达同步网络主控装置及控制方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201732160U (zh) * 2010-05-05 2011-02-02 北京航空航天大学 一种双接口雷达数据记录仪
CN102170318A (zh) * 2011-01-11 2011-08-31 南京国睿安泰信科技股份有限公司 用于接收机性能测试的频谱分析算法
US20130016761A1 (en) * 2011-07-15 2013-01-17 Renesas Mobile Corporation Transmitter with a Variable Sampling Rate
CN105785335A (zh) * 2016-03-28 2016-07-20 电子科技大学 一种基于cPCI的数字阵接收通道性能自动测试系统
CN105676197A (zh) * 2016-03-30 2016-06-15 西安电子科技大学 一种阵列雷达回波模拟器及其模拟方法
CN106406751A (zh) * 2016-06-03 2017-02-15 中北大学 基于多通道lvds接口的高速大容量i/q数据记录仪
CN108051786A (zh) * 2017-10-30 2018-05-18 北京航天福道高技术股份有限公司 一种宽带目标模拟器验证平台及验证方法
CN108333566A (zh) * 2018-02-07 2018-07-27 哈尔滨工业大学 一种基于zynq的便携式雷达测试系统及测试方法
CN110632568A (zh) * 2019-11-05 2019-12-31 中国科学院电子学研究所 合成孔径雷达实时成像处理器的测试信号源
CN111505593A (zh) * 2020-04-30 2020-08-07 北京无线电测量研究所 一种频综综合测试系统及测试方法
CN111562552A (zh) * 2020-04-30 2020-08-21 北京无线电测量研究所 一种地面雷达同步网络主控装置及控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
WALT KESTER: "Understand SINAD, ENOB, SNR, THD, THD + N, and SFDR so You Don\'t Get Lost in the Noise Floor", 《MT-003 TUTORIAL》, pages 1 - 8 *
黄庆等: "一种基于时-频域计算的ADC动态参数测试方法", 《微电子学》, pages 782 - 786 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113805148A (zh) * 2021-08-18 2021-12-17 北京无线电测量研究所 一种多功能数字tr组件光纤测试装置、系统和方法
CN113805148B (zh) * 2021-08-18 2023-11-28 北京无线电测量研究所 一种多功能数字tr组件光纤测试装置、系统和方法
CN116955258A (zh) * 2023-09-19 2023-10-27 成都立思方信息技术有限公司 一种可灵活连接的触发连接器、信号采集控制设备及系统
CN116955258B (zh) * 2023-09-19 2023-11-28 成都立思方信息技术有限公司 一种可灵活连接的触发连接器、信号采集控制设备及系统

Also Published As

Publication number Publication date
CN112305511B (zh) 2024-05-07

Similar Documents

Publication Publication Date Title
CN105785335B (zh) 一种基于cPCI的数字阵接收通道性能自动测试系统
CN112305511A (zh) 一种雷达接收机数字基带信号录取分析系统和方法
CN108563144B (zh) 一种弹载雷达信号处理半实物仿真测试系统
US7263286B2 (en) Fast testing system for optical transceiver and testing method thereof
CN203038273U (zh) 多功能数据采集卡
US20140362901A1 (en) Methods and systems for providing optimum decision feedback equalization of high-speed serial data links
US8913517B2 (en) Method for measuring sensitivity of data packet signal receiver
CN101193324A (zh) 机顶盒集成测试系统
CN112035302B (zh) 一种总线数据的实时监控分析方法、装置及系统
CN109783418A (zh) 一种宽频带信号高速采集与波形存储分析系统
CN110221261B (zh) 一种雷达波形产生模块测试分析方法及装置
CN105067894A (zh) 混频器变频损耗的测试方法和系统
CN114184941A (zh) 一种针对含有adc芯片的模块的测试系统及方法
CN113938191A (zh) 一种对光模块进行参数测试的方法及装置
US6356220B1 (en) A/D converting device, A/D converting method and optical pulse testing apparatus
CN208971519U (zh) 多通道光模块误码率测试系统
CN101998135A (zh) 移动电视信号采集及播放系统、控制方法
US20230070298A1 (en) Real-equivalent-time clock recovery for a nearly-real-time real-equivalent-time oscilloscope
CN113806277B (zh) 基于srio协议的fpga与dsp的数据传输系统
CN101738308A (zh) 光缆测试数据采集方法、采集装置及系统
CN114371392A (zh) 触控数据检测装置、触控芯片及电子设备
WO2021115445A1 (zh) 一种光模块劣化测试方法、系统、设备和存储介质
CN201345098Y (zh) 一种同时读写存储器的数据采集装置
CN107144751A (zh) 一种多信道矢量网络参数分析系统及方法
CN219420769U (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