CN105699770B - 百万通道频谱仪的频点输出信号失真恢复方法 - Google Patents
百万通道频谱仪的频点输出信号失真恢复方法 Download PDFInfo
- Publication number
- CN105699770B CN105699770B CN201610068254.6A CN201610068254A CN105699770B CN 105699770 B CN105699770 B CN 105699770B CN 201610068254 A CN201610068254 A CN 201610068254A CN 105699770 B CN105699770 B CN 105699770B
- Authority
- CN
- China
- Prior art keywords
- frequency point
- signal
- channel
- distorted
- fft
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种百万通道频谱仪的频点输出信号失真恢复方法，具体为：1）搭建二维级联式FFT运算模块，分析获得频谱泄漏的频点上能量的分配比；2）通过信号相邻通道内泄漏信号叠加方式对输出信号通道进行信号还原。本发明弥补了级联FFT算法存在天生的不足，提供一种实用的方法来进行信号的恢复，使得该算法在输出频谱时是可行的。
Description
技术领域
本发明涉及百万通道频谱仪领域，尤其是一种应用于ＦＡＳＴ的百万通道频谱仪的二维级联FFT（Fast Fourier Transform）输出频点信号失真的恢复方法。
背景技术
接收机系统是射电望远镜的重要组成部分，是决定望远镜性能的关键因素。最近5年，超宽带技术有重要发展，世界各大射电望远镜都在建造或考虑建造超宽带接收机。我国500米口径球面射电望远镜（Five-hundred-meter Aperture Spherical radioTelescope,简称FAST ）将于2016年投入使用，届时是已知的70 MHz和3 GHz之间的低频无线电波段中灵敏度最高的单天线射电望远镜。为了使其能够充分发挥其口径大的性能，确保望远镜在世界上的领先，我国自主研发世界一流的超宽带通用型数字后端(ChinaReconfigurable Analog-Digital backend,简称CRANE)。宽带接收机研制出之后不仅可以直接应用于FAST，而且还用于中国的其他射电望远镜，如上海65 米，新疆110米等。
目前，在整个CRANE开发中，硬件部分――通过与中国科学院自动化研究所的合作，不仅开发了现场可编程门阵列（Field Programmable Gate Array，简称FPGA）运算板FDB（FAST Digital Backend），还开发了3Gsps，12bit的宽带高精度模数转换电路板（ADC）芯片（简称FAST3212ADC），并用此实现3GHz带宽的一次性覆盖。此超宽带通用型数字后端同时也包括开发模拟信号前端电路板（Analog Front-end Board，简称AFB）。固件部分，用基于FPGA运算板的多种信号处理算法来实现FDB对于数字信号处理的通用性。在射电望远镜频谱仪中，如果想要获得更高的频谱分辨率，就尽可能的提高 FFT 通道数。现有通用的CASPER组 FFT通道数较少，处理速度不够快，若应用在FDB上，无法充分利用 FAST 所提供的超高性能。
发明内容
针对现有技术存在的问题，本发明的目的在于提供一种百万通道频谱仪的频点输出信号失真恢复方法，具体为：
1）搭建二维级联式FFT运算模块，分析获得频谱泄漏的频点上能量的分配比；
2）通过信号相邻通道内泄漏信号叠加方式对输出信号通道进行信号还原。
进一步，搭建所述二维级联式FFT运算模块，具体为：通过级联两个小点数的FFT来进行实现大点数的FFT，通过直接搭载两片Virtex-6芯片来满足运算要求，每一块芯片分别加载一个子FFT运算单元。
进一步，所述频谱泄漏的频点上能量分配比的分析方法为：向二维级联FFT运算模块中输入单频率原始信号，经过处理，得到失真的信号频率，绘制该失真信号频谱图，得到频谱泄漏后各频点通道内的能量分配情况。
进一步，所述步骤2）具体为：
A）确定能量泄漏在所需频点周围的分配比例，从而确定叠加通道数；
B）将泄漏到临近各个通道频率处的线谱相加，左右各5个通道进行第一次相加；
C）若得到的频率能量不能够增加3dB以上，则再分别在左右增加5个通道进行叠加；
D）直至得到的频率能量增加值为大于等于3dB，从而恢复该频率点处真实的信号能量值。
本发明弥补了级联FFT算法存在天生的不足，提供一种实用的方法来进行信号的恢复，使得该算法在输出频谱时是可行的。
附图说明
图1为本发明流程框架图；
图2为二维FFT考虑通道偏离后频谱泄漏通道情况；
图3为图2中局部放大图；
图4为通道偏离20造成能量泄漏情况下二维FFT 频点处的信号能量恢复值；
图5为图4中局部放大图；
图6为通道偏离32造成能量泄漏情况下二维FFT 频点处的信号能量恢复值；
图7为图6中局部放大图。
具体实施方式
下面，参考附图，对本发明进行更全面的说明，附图中示出了本发明的示例性实施例。然而，本发明可以体现为多种不同形式，并不应理解为局限于这里叙述的示例性实施例。而是，提供这些实施例，从而使本发明全面和完整，并将本发明的范围完全地传达给本领域的普通技术人员。
如图1所示，本发明一种百万通道频谱仪的频点输出信号失真恢复方法，具体为：
1）搭建二维级联式FFT运算模块，分析获得频谱泄漏的频点上能量的分配比；
2）通过信号相邻通道内泄漏信号叠加方式对输出信号通道进行信号还原。
本发明优化了FFT，研究级联（cascade） FFT（两个小点数FFT），步骤1）中搭建所述二维级联式FFT运算模块，具体为：通过级联两个小点数的FFT 来进行实现大点数的FFT。FDB 的硬件设计中通过直接搭载两片Virtex-6芯片来满足运算要求，每一块芯片分别加载一个子FFT运算单元。一级FFT运算单元的分辨率可能会造成最终频点输出时信号失真，
本发明通过在Matlab 上进行级联FFT 算法模拟，研究信号失真的情况，解决对于信号失真情况下的信号恢复的问题。
搭建二维级联式FFT运算模块，分析获得频谱泄漏的频点上能量的分配比，通过信号相邻通道内泄漏信号叠加方式对输出信号通道进行信号还原。
1）分析频谱泄漏的频点上能量的分配比：
利用Matlab模拟分析，得出在预期大小的二维级联FFT运算单元处理后，频谱泄漏的频点周围能量的分配情况，决定下一步选择的叠加通道数。
具体为：向二维级联FFT运算模块中输入单频率原始信号，经过处理，得到失真的信号频率，绘制该失真信号频谱图，得到频谱泄漏后各频点通道内的能量分配情况。
图2给出了有信号通道偏离时，频谱泄漏的频点周围能量的分配情况，由黑色加号点给出，它反映出有频谱泄漏的频点上能量的分配比（未归一化），真实的通道能量由对应的蓝色点给出。图3为其局部放大图。
2）选择将泄漏到临近通道频率处的线谱相加，进而恢复该频率点处真实的信号能量值。
该步骤具体为：
A）确定能量泄漏在所需频点周围的各个频点，从而确定叠加通道数；
B）将泄漏到临近各个通道频率处的线谱相加，左右各5个通道进行第一次相加；
C）若得到的频率能量不能够增加3dB以上，则再分别在左右增加5个通道进行叠加；
D）直至得到的频率能量增加值为大于等于3dB，从而恢复该频率点处真实的信号能量值。
例如，我们分别选取了二级级联FFT运算单元处理后，以最近的整数频点起，第20，32个通道处的频率（即通道偏离为20，32）。如图4所示，通道偏离20造成能量泄漏情况下二维FFT 频点处的信号能量恢复值，由最上方粉色点给出。图5为图4的局部放大图，清晰反映出能量恢复后点的位置；如图6所示，通道偏离32造成能量泄漏情况下二维FFT 频点处的信号能量恢复值，由最上方粉色点给出。图7此为图6的局部放大图，清晰反映出能量恢复后点的位置。上述附图分别给出了通道偏离为20，32时，频谱泄漏后的信号恢复分析：图6、图7可以看出，相邻通道内的泄漏能量相加后，信号能量增加了约2dB。理论上应该是增加3dB，我们只取了周围40条信号通道的值进行相加。图中能量值最大的粉色点是恢复后的真实信号强度。通道偏离大于32时，靠近下一个二级FFT运算单元整数频点，所以具有相对称的性质，真实信号通道的能量会向之后的频点泄漏越来越多的能量。
本方法解决了以下关键问题：
1）分析算法在实际应用过程中带来的频谱失真，掌握影响频谱正确输出因素
要实现二维级联FFT算法在频谱仪中真正能够利用，必须先清楚地得到算法对频谱输出正确性的影响。在进行二维级联FFT运算时，信号会因为第一级FFT运算单元输出长度而造成频谱泄漏。此二维级联FFT达到的频谱分辨率为1KHz，大小为N1×N2。在经过第一级FFT运算后，得到的“分辨率”为1K×N1Hz,此时可能会导致一级运算单元输出失真，即一级输出信号位置落在N1的非整数倍，最终会在原有精度上造成一定程度的频谱泄漏。通过Matlab模拟，量化性的得到频谱失真特性，并绘制成图，确定后续进行信号恢复的基础以及信号恢复算法的选定。
2）发展实用的信号恢复方法，以实现二维级联FFT的可行性。
我们最终是需要将二维级联FFT应用到硬件上搭建百万通道频谱仪的，级联FFT算法存在天生的不足，故必须发展出一种实用的方法来进行信号的恢复，使得该算法在输出频谱时是可行的。
Claims (4)
1.一种百万通道频谱仪的频点输出信号失真恢复方法，其特征在于，该方法具体为：
1)搭建二维级联式FFT运算模块，分析获得频谱泄漏的频点上能量的分配比；
2)通过信号相邻通道内泄漏信号叠加方式对输出信号通道进行信号还原。
2.如权利要求1所述的百万通道频谱仪的频点输出信号失真恢复方法，其特征在于，搭建所述二维级联式FFT运算模块，具体为：通过级联两个小点数的FFT来进行实现大点数的FFT，通过直接搭载两片Virtex-6芯片来满足运算要求，每一块芯片分别加载一个子FFT运算单元。
3.如权利要求1所述的百万通道频谱仪的频点输出信号失真恢复方法，其特征在于，所述频谱泄漏的频点上能量分配比的分析方法为：向二维级联FFT运算模块中输入单频率原始信号，经过处理，得到失真信号的频率，绘制该失真信号的频谱图，得到频谱泄漏后各频点通道内的能量分配情况。
4.如权利要求1所述的百万通道频谱仪的频点输出信号失真恢复方法，其特征在于，所述步骤2)具体为：
A)确定能量泄漏在所需频点周围的分配比例，从而确定叠加通道数；
B)将泄漏到临近各个通道频率处的线谱相加，左右各5个通道进行第一次相加；
C)若得到的频率能量不能够增加3dB以上，则再分别在左右增加5个通道进行叠加；
D)直至得到的频率能量增加值为大于等于3dB，从而恢复该频点处真实的信号能量值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610068254.6A CN105699770B (zh) | 2016-02-01 | 2016-02-01 | 百万通道频谱仪的频点输出信号失真恢复方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610068254.6A CN105699770B (zh) | 2016-02-01 | 2016-02-01 | 百万通道频谱仪的频点输出信号失真恢复方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105699770A CN105699770A (zh) | 2016-06-22 |
| CN105699770B true CN105699770B (zh) | 2019-01-08 |
Family
ID=56228925
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610068254.6A Active CN105699770B (zh) | 2016-02-01 | 2016-02-01 | 百万通道频谱仪的频点输出信号失真恢复方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105699770B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4918381A (en) * | 1989-04-06 | 1990-04-17 | Hewlett-Packard Company | Automated method for determining total harmonic distortion |
| WO2006082681A1 (ja) * | 2005-02-01 | 2006-08-10 | Hitachi Kokusai Electric Inc. | スペクトル解析方法、歪検出装置、歪補償増幅装置 |
| US20080019464A1 (en) * | 2002-04-22 | 2008-01-24 | Cognio, Inc. | Signal Pulse Detection Scheme for Use in Real-Time Spectrum Analysis |
| CN101639499A (zh) * | 2009-09-01 | 2010-02-03 | 中国电子科技集团公司第四十一研究所 | 基波跟踪的小信号谐波失真的测量装置及方法 |
| CN102662108A (zh) * | 2012-05-14 | 2012-09-12 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种正弦波局域失真的测量方法及装置 |
| CN102998528A (zh) * | 2012-12-12 | 2013-03-27 | 湖南大学 | 动态正弦畸变信号中频率分量的提取方法及其实现装置 |
| CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
| CN104391282A (zh) * | 2014-11-24 | 2015-03-04 | 天津大学 | 通过频谱校正来改善成像质量的方法 |
-
2016
- 2016-02-01 CN CN201610068254.6A patent/CN105699770B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4918381A (en) * | 1989-04-06 | 1990-04-17 | Hewlett-Packard Company | Automated method for determining total harmonic distortion |
| US20080019464A1 (en) * | 2002-04-22 | 2008-01-24 | Cognio, Inc. | Signal Pulse Detection Scheme for Use in Real-Time Spectrum Analysis |
| WO2006082681A1 (ja) * | 2005-02-01 | 2006-08-10 | Hitachi Kokusai Electric Inc. | スペクトル解析方法、歪検出装置、歪補償増幅装置 |
| CN101639499A (zh) * | 2009-09-01 | 2010-02-03 | 中国电子科技集团公司第四十一研究所 | 基波跟踪的小信号谐波失真的测量装置及方法 |
| CN102662108A (zh) * | 2012-05-14 | 2012-09-12 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种正弦波局域失真的测量方法及装置 |
| CN102998528A (zh) * | 2012-12-12 | 2013-03-27 | 湖南大学 | 动态正弦畸变信号中频率分量的提取方法及其实现装置 |
| CN103308766A (zh) * | 2013-05-15 | 2013-09-18 | 湖南大学 | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 |
| CN104391282A (zh) * | 2014-11-24 | 2015-03-04 | 天津大学 | 通过频谱校正来改善成像质量的方法 |
Non-Patent Citations (4)
| Title |
|---|
| A 6 × 320-MHz 1024-channel FFT cross-spectrum analyzer for radio astronomy;Y. Chikada et.al;《Proceedings of the IEEE》;19870930;第75卷(第9期);第1203-1208页 * |
| Data processing and synchronization of the digital backend for FAST;Yu Xinying et.al;《2015 IEEE International Symposium on Precision Clock Synchronization for Measurement, Control, and Communication (ISPCS)》;20151016;第48页 * |
| FFT 分析中一种减少频谱泄漏的方法;林俊武 等;《三明学院学报》;20070630;第24卷(第2期);第166-169页 * |
| 离散频谱泄漏抑制方法研究进展;杨辉跃 等;《计算机应用研究》;20150630;第32卷(第6期);第1608-1013页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105699770A (zh) | 2016-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105158735B (zh) | 基于压缩采样阵列的空频二维谱估计方法 | |
| CN112054976A (zh) | 一种超宽带多信道信号并行处理分析方法及系统 | |
| Correia et al. | OFDM-like high order backscatter modulation | |
| CN105699770B (zh) | 百万通道频谱仪的频点输出信号失真恢复方法 | |
| CN117434426B (zh) | 一种开关电容滤波器的测试方法、系统及装置 | |
| CN105372494A (zh) | 2GHz带宽实时FFT频谱仪系统 | |
| CN106031046A (zh) | 基于压缩感知的信号处理方法及装置 | |
| CN103001712A (zh) | 一种多频段、多信号二阶互调测试系统 | |
| CN109001554A (zh) | 抗扰度试验多频率信号发生器 | |
| CN111030765A (zh) | 一种能识别镜频信号的外差扫描式频谱分析系统 | |
| CN111082834A (zh) | 一种基于啁啾信号正交解调的射频时延快速测量装置 | |
| CN111999587B (zh) | 一种调制域网络参数测试的迭代矢量误差修正方法 | |
| CN203479903U (zh) | 频谱分析仪通用跟踪源系统 | |
| Ma et al. | A novel real-time Fourier and inverse Fourier transforming system based on non-uniform coupled-line phaser | |
| CN109525265A (zh) | 基于欠采样补偿的宽带测量装置及其方法 | |
| CN104967491A (zh) | 多通道幅相测试系统信号接收处理方法 | |
| CN104270209B (zh) | 基于不同校准平面的rru驻波比的检测方法和装置 | |
| CN105021900B (zh) | 一种用于多通道测量的通道捷变差分测量方法 | |
| CN106301469B (zh) | 一种同轴电缆的幅相稳定性测试方法 | |
| CN104683041A (zh) | 射频元器件测试方法 | |
| CN109302243B (zh) | 一种n端口网分的抗同频干扰扫描方法 | |
| CN105186089A (zh) | 一种小型化大频比微波双频功分器 | |
| CN106603069A (zh) | 一种甚低频无源相位合成网络 | |
| CN115575957B (zh) | 一种毫米波大气臭氧成分探测仪 | |
| Kang et al. | Electromagnetic signal intelligent identification based on radio frequency fingerprints |
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 |