CN107493118B - 信号获取方法及装置 - Google Patents

信号获取方法及装置 Download PDF

Info

Publication number
CN107493118B
CN107493118B CN201710783735.XA CN201710783735A CN107493118B CN 107493118 B CN107493118 B CN 107493118B CN 201710783735 A CN201710783735 A CN 201710783735A CN 107493118 B CN107493118 B CN 107493118B
Authority
CN
China
Prior art keywords
signal
fft
narrow
band
noise signal
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
Application number
CN201710783735.XA
Other languages
English (en)
Other versions
CN107493118A (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.)
Changsha Haige Beidou Information Technology Co Ltd
Original Assignee
Changsha Haige Beidou Information Technology Co Ltd
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 Changsha Haige Beidou Information Technology Co Ltd filed Critical Changsha Haige Beidou Information Technology Co Ltd
Priority to CN201710783735.XA priority Critical patent/CN107493118B/zh
Publication of CN107493118A publication Critical patent/CN107493118A/zh
Application granted granted Critical
Publication of CN107493118B publication Critical patent/CN107493118B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/71Interference-related aspects the interference being narrowband interference
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/0082Monitoring; Testing using service channels; using auxiliary channels
    • H04B17/0087Monitoring; Testing using service channels; using auxiliary channels using auxiliary channels or channel simulators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/20Monitoring; Testing of receivers
    • H04B17/29Performance testing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Noise Elimination (AREA)

Abstract

本发明实施例提供一种信号获取方法及装置,属于数据处理领域。所述方法包括:获取高斯伪随机噪声信号X;对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft;在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft;对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y,从而本方法通过在频域对信号X_fft进行幅值加权可以有效提高获取的窄带高斯噪声信号Y的带外抑制度,减少窄带高斯噪声信号Y的过渡带带宽,并且,通过该方法获取窄带高斯噪声信号Y不需要复杂大量的硬件资源,简单、容易实现。

Description

信号获取方法及装置
技术领域
本发明涉及数据处理领域,具体而言,涉及一种信号获取方法及装置。
背景技术
窄带高斯噪声信号,在测试通信系统的性能时,发挥着重要作用。窄带信号可以用于干扰机产生针对特定频点的压制式干扰信号,以测试无线接收机的抗窄带干扰能力,或用于信道模拟,在中频与信号合成具有特定载噪比的测试信号等。上述应用中,往往要求高斯窄带信号本身具有良好的频谱特性,如过渡带陡峭、较高的带外抑制水平等,以确保测试的准确性。
通常为了测量通信系统在一定信噪比条件下的通信指标,加入噪声的强度以及一定带宽下噪声信号本身的带外底噪水平,都需要可控,所以信号的产生大多是在数字域进行,然后再通过DAC变换到模拟域。目前数字域产生高斯窄带噪声信号的方法,主要是通过FIR、IIR等带通滤波器,得到所需的高斯窄带噪声信号。这种基于滤波器的高斯窄带噪声信号产生方法,随着所需高斯窄带噪声信号的过渡带陡峭程度的提高、同时兼顾很高的带外抑制水平,其滤波器阶数将难以承受,需要大量的硬件资源,工程实现难度大。比如,采样率20MHz下,一个带宽2MHz,过渡带50KHz,带外抑制65dB的窄带高斯噪声信号,所需FIR滤波器阶数将达到800阶以上。使用IIR椭圆带通滤波器,仍需要10组20阶的高阶滤波器,数字实现如此高阶的滤波器存在有效位和数据截断等问题。而使用Matlab等软件方法事先由高阶滤波器产生好窄带噪声信号,然后截取一定个数的信号点存入硬件ROM,硬件通过查表输出该信号。由于截取的信号点数有限,其输出信号在短时间内即会整体呈现出周期性,比如以20MHz采样率为例,1s即对应20000000个点,需要大量ROM空间,且其信号整体的伪随机特性只能持续1s,在对信号随机性要求高的场合不适用。
发明内容
有鉴于此,本发明实施例的目的在于提供一种信号获取方法及装置,以改善上述问题。
第一方面,本发明实施例提供了一种信号获取方法,应用于一电子设备中,所述方法包括:获取高斯伪随机噪声信号X;对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft;在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft;对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y。
进一步地,对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft,包括:对所述高斯伪随机噪声信号X进行N点傅里叶变换,获得变换后的信号X_fft,其中,N>=fs/WS,fs为所述窄带高斯噪声信号Y的采样频率,W为所述窄带高斯噪声信号Y的带宽,WS为所述窄带高斯噪声信号Y的过渡带宽。
进一步地,在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft,包括:取索引n为1到N,在频域,对所述信号X_fft进行N点遍历,若检测到当前遍历点的位置索引n<[(f0-W)/(fs/N)]或者n>[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft取值为0,若检测到当前遍历点的位置索引n>=[(f0-W)/(fs/N)]且n<=[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft保持不变,其中,f0为中心频率。
进一步地,对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y,包括:对所述信号Y_fft进行傅里叶反变换,获得经反变换信号;保留所述经反变换信号的实部,以获取窄带高斯噪声信号Y。
进一步地,对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y之后,所述方法还包括:将所述窄带高斯噪声信号Y中的每个点除以一预设常数,以获取降低了所述窄带高斯噪声信号Y的带外抑制度后的信号Z1;或将所述窄带高斯噪声信号Y中的每个点乘以一预设常数,以获取提高了所述窄带高斯噪声信号Y的带外抑制度后的信号Z2。
第二方面,本发明实施例提供了一种信号获取装置,运行一电子设备中,所述装置包括:伪随机噪声获取模块,用于获取高斯伪随机噪声信号X;傅里叶变换模块,用于对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft;幅值加权模块,用于在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft;傅里叶反变换模块,用于对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y。
进一步地,所述傅里叶变换模块,具体用于对所述高斯伪随机噪声信号X进行N点傅里叶变换,获得变换后的信号X_fft,其中,N>=fs/WS,fs为所述窄带高斯噪声信号Y的采样频率,W为所述窄带高斯噪声信号Y的带宽,WS为所述窄带高斯噪声信号Y的过渡带宽。
进一步地,所述幅值加权模块,具体用于取索引n为1到N,在频域,对所述信号X_fft进行N点遍历,若检测到当前遍历点的位置索引n<[(f0-W)/(fs/N)]或者n>[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft取值为0,若检测到当前遍历点的位置索引n>=[(f0-W)/(fs/N)]且n<=[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft保持不变,其中,f0为中心频率。
进一步地,所述傅里叶反变换模块包括:反变换信号获取单元,用于对所述信号Y_fft进行傅里叶反变换,获得经反变换信号;窄带高斯噪声信号获取单元,用于保留所述经反变换信号的实部,以获取窄带高斯噪声信号Y。
进一步地,所述装置还包括:第一带外抑制模块,用于将所述窄带高斯噪声信号Y中的每个点除以一预设常数,以获取降低了所述窄带高斯噪声信号Y的带外抑制度后的信号Z1;第二带外抑制模块,用于将所述窄带高斯噪声信号Y中的每个点乘以一预设常数,以获取提高了所述窄带高斯噪声信号Y的带外抑制度后的信号Z2。
本发明实施例的有益效果是:
本发明实施例提供一种信号获取方法及装置,首先获取高斯伪随机噪声信号X,对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft,然后在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft,再对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y,从而本方法通过在频域对信号X_fft进行幅值加权可以有效提高获取的窄带高斯噪声信号Y的带外抑制度,减少窄带高斯噪声信号Y的过渡带带宽,并且,通过该方法获取窄带高斯噪声信号Y不需要复杂大量的硬件资源,简单、容易实现。
本发明的其他特征和优点将在随后的说明书阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明实施例了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1示出了一种可应用于本申请实施例中的电子设备的结构框图;
图2为本发明实施例提供的一种信号获取方法的流程图;
图3为本发明实施例提供的另一种信号获取方法的流程图;
图4为本发明实施例提供的一种信号获取装置的结构框图。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本发明的描述中,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
请参照图1,图1示出了一种可应用于本申请实施例中的电子设备100的结构框图。电子设备100可以包括信号获取装置、存储器101、存储控制器102、处理器103、外设接口104、输入输出单元105、音频单元106、显示单元107。
所述存储器101、存储控制器102、处理器103、外设接口104、输入输出单元105、音频单元106、显示单元107各元件相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述信号获取装置包括至少一个可以软件或固件(firmware)的形式存储于所述存储器101中或固化在所述信号获取装置的操作系统(operating system,OS)中的软件功能模块。所述处理器103用于执行存储器101中存储的可执行模块,例如所述信号获取装置包括的软件功能模块或计算机程序。
其中,存储器101可以是,但不限于,随机存取存储器(Random Access Memory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(Programmable Read-OnlyMemory,PROM),可擦除只读存储器(Erasable Programmable Read-Only Memory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory,EEPROM)等。其中,存储器101用于存储程序,所述处理器103在接收到执行指令后,执行所述程序,前述本发明实施例任一实施例揭示的流过程定义的服务器所执行的方法可以应用于处理器103中,或者由处理器103实现。
处理器103可以是一种集成电路芯片,具有信号的处理能力。上述的处理器103可以是通用处理器,包括中央处理器(Central Processing Unit,简称CPU)、网络处理器(Network Processor,简称NP)等;还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器103也可以是任何常规的处理器等。
所述外设接口104将各种输入/输出装置耦合至处理器103以及存储器101。在一些实施例中,外设接口104,处理器103以及存储控制器102可以在单个芯片中实现。在其他一些实例中,他们可以分别由独立的芯片实现。
输入输出单元105用于提供给用户输入数据实现用户与所述服务器(或本地终端)的交互。所述输入输出单元105可以是,但不限于,鼠标和键盘等。
音频单元106向用户提供音频接口,其可包括一个或多个麦克风、一个或者多个扬声器以及音频电路。
显示单元107在所述电子设备100与用户之间提供一个交互界面(例如用户操作界面)或用于显示图像数据给用户参考。在本实施例中,所述显示单元107可以是液晶显示器或触控显示器。若为触控显示器,其可为支持单点和多点触控操作的电容式触控屏或电阻式触控屏等。支持单点和多点触控操作是指触控显示器能感应到来自该触控显示器上一个或多个位置处同时产生的触控操作,并将该感应到的触控操作交由处理器103进行计算和处理。
所述外设接口104将各种输入/输入装置耦合至处理器103以及存储器101。在一些实施例中,外设接口104,处理器103以及存储控制器102可以在单个芯片中实现。在其他一些实例中,他们可以分别由独立的芯片实现。
输入输出单元105用于提供给用户输入数据实现用户与处理终端的交互。所述输入输出单元105可以是,但不限于,鼠标和键盘等。
可以理解,图1所示的结构仅为示意,所述电子设备100还可包括比图1中所示更多或者更少的组件,或者具有与图1所示不同的配置。图1中所示的各组件可以采用硬件、软件或其组合实现。
请参照图2,图2为本发明实施例提供的一种信号获取方法的流程图,所述方法应用于上述电子设备中,所述方法具体包括如下步骤:
步骤S110:获取高斯伪随机噪声信号X。
本实施例中可以采用m序列和查表法来获取高斯伪随机噪声信号X,其具体过程为:通过两个相互独立的均匀分布(如m序列)的伪随机数发生器1和发生器2分别产生随机数流r1和r2;再以r1为地址查询1n自然对数表,得打x1,以r2为地址查询正弦表,得到x2;然后将x1和x2相乘即可得到高斯伪随机噪声信号X。
步骤S120:对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft。
对所述高斯伪随机噪声信号X进行N点傅里叶变换,其中,N>=fs/WS,fs为所述窄带高斯噪声信号Y的采样频率,W为所述窄带高斯噪声信号Y的带宽,WS为所述窄带高斯噪声信号Y的过渡带宽。
步骤S130:在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft。
高斯伪随机噪声信号X进行傅里叶变换的点数为N,其高斯伪随机噪声信号X的中心频率为f0,其满足f0+w/2<fs/2。
在频域,对所述信号X_fft进行幅值加权,也就是保留窄带高斯噪声信号通带频点内的幅值,其他带外频点对应的幅值置0,即加权系数为1或0,从而可以不使用乘法器即可实现。当然,也可以通过乘法器的使用得到任意的加权系数,不限于1或0。
对所述信号X_fft进行幅值加权的过程为:取索引n为1到N,在频域,对所述信号X_fft进行N点遍历,若检测到当前遍历点的位置索引n<[(f0-W)/(fs/N)]或者n>[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft取值为0,即X_fft(n)=0,若检测到当前遍历点的位置索引n>=[(f0-W)/(fs/N)]且n<=[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft保持不变,即X_fft(n)=X_fft(n),其中,符号[]表示四舍五入取整。
当1到N个点遍历完成,则对所述信号X_fft的幅值加权过程也完成。
需要说明的是,也可以在时域对所述信号X_fft进行幅值加权。
步骤S140:对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y。
对所述信号Y_fft进行傅里叶反变换,获得经反变换信号,保留所述经反变换信号的实部,从而获取窄带高斯噪声信号Y。
请参照图3,还包括,步骤S150:将所述窄带高斯噪声信号Y中的每个点除以一预设常数,以获取降低了所述窄带高斯噪声信号Y的带外抑制度后的信号Z1。
为了降低窄带高斯信号Y的带外抑制度,将信号Y中的每个点均除以一预设常数,得到信号Z1,该预设常数可以为2的幂次,当然也可以是其他数,然后通过将数据右移实现,以节省除法器,每右移一位,带外抑制度降低约6dB。
步骤S160:将所述窄带高斯噪声信号Y中的每个点乘以一预设常数,以获取提高了所述窄带高斯噪声信号Y的带外抑制度后的信号Z2。
为了提高窄带高斯信号Y的带外抑制度,将信号Y中的每个点分别乘以一个预设常数,得到信号Z2,该预设常数可以为2的幂次,当然也可以是其他数,然后通过将数据左移实现,以节省乘法器,每左移一位,带外抑制度提高约6dB。
信号Z1和信号Z2即是信号Y调整带外抑制度后的结果,Z1和Z2可以作为数字形式的窄带高斯噪声信号输出。
下面通过以具体实施例对上述方法的实现过程进行具体说明。
以将该信号获取方法在干扰机信号源上的应用为例。
干扰机信号源窄带高斯信号发生器的数字部分,可以采用现场可编程门阵列(FPGA)实现,数模转换部分可以使用高性能的14为的数模转换器。
其中,干扰机信号源产生中心频率为f0=15.5MHz、带宽W=2MHz、过渡带宽小于WS<50KHz且带外抑制度优于65dB的窄带高斯噪声信号。
在FPGA设计了如下模块,包括高斯伪随机噪声产生模块、FFT双缓冲模块、N点FFT模块、幅值加权模块、N点IFFT模块、IFFT输出双缓冲模块、底噪抑制调制模块及主控状态机模块。
首先选择系统工作频率fs=40MHz,可以计算得到将高斯伪噪声进行傅里叶变换的点数N应大于fs/50000=800,实际的窄带高斯噪声信号的过渡带宽为WS=40KHz,N>=fs/WS=1024,所以可取傅里叶变换点数N=1024。
然后将高斯伪随机噪声产生模块产生的14bit量化的高斯伪随机噪声信号X并写入FFT双缓冲模块,当FFT双缓冲模块中的数据满1024点时,取出作为一帧送往N点FFT模块,进行1024点傅里叶变换,得到变换后的信号X_fft。
再将信号X_fft通过幅值加权模块进行幅值加权获得信号Y_fft,也就是,取索引n=1到1024,逐个遍历X_fft中的1024个点,根据计算得到n<[(f0-W)/(fs/N)]=345,或者n>[(f0+W)/(fs/N)]=450,所以,令n<345或n>450索引位置对应的数据点X_fft(n)=0,其他索引位置对应的数据点不变,并将获得的信号Y_fft再进行傅里叶反变换。
信号Y_fft经过N点IFFT模块进行傅里叶反变换,获得经反变换信号,然后保留所述经反变换信号的实部,以获取窄带高斯噪声信号Y。还可以将经反变换信号有效量化位宽为12位,并写入IFFT输出双缓冲模块,在通过底噪抑制调制模块进行带宽抑制时,从IFFT输出双缓冲模块中读取数据,并右移2位,也就是最高2位补零,得到14位量化的数字窄带高斯噪声信号Z1;或者,从IFFT输出双缓冲模块中读取数据,并左移2位,也就是最低2位补零,得到14位量化的数字窄带高斯噪声信号Z2。
将信号Z1或Z2输入数模转换器得到中心频率f0=15.5MHz、带宽W=2MHz、过渡带宽WS<50KHz,且带外抑制优于80dB的窄带高斯噪声信号。
另外,系统中的主控机对FFT双缓冲的输入和底噪抑制模块从IFFT输出缓冲中读取数据的过程进行控制,确保进入数模转换器的数据是连续的。同时,对高斯伪随机噪声产生模块中的m序列多项式抽头进行重配置,使得噪声信号在长的观测时间上仍然满足随机性,提高测试的可靠性。
所以本实施例中的信号获取方法,提高FFT的点数即可减少窄带高斯噪声信号的过渡带带宽,而增加数模转换器的有效位,即可提高窄带高斯噪声信号的带外抑制度,二者可兼得。
请参照图4,图4为本发明实施例提供的一种信号获取装置的结构框图,所述装置运行于一电子设备中,所述装置包括:伪随机噪声获取模块、傅里叶变换模块、幅值加权模块及傅里叶反变换模块。
伪随机噪声获取模块,用于获取高斯伪随机噪声信号X。
傅里叶变换模块,用于对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft。
幅值加权模块,用于在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft。
傅里叶反变换模块,用于对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y。
作为一种方式,所述傅里叶变换模块,具体用于对所述高斯伪随机噪声信号X进行N点傅里叶变换,获得变换后的信号X_fft,其中,N>=fs/WS,fs为所述窄带高斯噪声信号Y的采样频率,W为所述窄带高斯噪声信号Y的带宽,WS为所述窄带高斯噪声信号Y的过渡带宽。
作为一种方式,所述幅值加权模块,具体用于取索引n为1到N,在频域,对所述信号X_fft进行N点遍历,若检测到当前遍历点的位置索引n<[(f0-W)/(fs/N)]或者n>[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft取值为0,若检测到当前遍历点的位置索引n>=[(f0-W)/(fs/N)]且n<=[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft保持不变,其中,f0为中心频率。
作为一种方式,所述傅里叶反变换模块包括:
反变换信号获取单元,用于对所述信号Y_fft进行傅里叶反变换,获得经反变换信号。
窄带高斯噪声信号获取单元,用于保留所述经反变换信号的实部,以获取窄带高斯噪声信号Y。
作为一种方式,所述装置还包括:
第一带外抑制模块,用于将所述窄带高斯噪声信号Y中的每个点除以一预设常数,以获取降低了所述窄带高斯噪声信号Y的带外抑制度后的信号Z1。
第二带外抑制模块,用于将所述窄带高斯噪声信号Y中的每个点乘以一预设常数,以获取提高了所述窄带高斯噪声信号Y的带外抑制度后的信号Z2。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的装置的具体工作过程,可以参考前述方法中的对应过程,在此不再过多赘述。
综上所述,本发明实施例提供一种信号获取方法及装置,首先获取高斯伪随机噪声信号X,对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft,然后在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft,再对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y,从而本方法通过在频域对信号X_fft进行幅值加权可以有效提高获取的窄带高斯噪声信号Y的带外抑制度,减少窄带高斯噪声信号Y的过渡带带宽,并且,通过该方法获取窄带高斯噪声信号Y不需要复杂大量的硬件资源,简单、容易实现。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

Claims (6)

1.一种信号获取方法,应用于一电子设备中,其特征在于,所述方法包括:
获取高斯伪随机噪声信号X;
对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft;
在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft;
对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y;
将所述窄带高斯噪声信号Y中的每个点除以一预设常数,以获取降低了所述窄带高斯噪声信号Y的带外抑制度后的信号Z1;或
将所述窄带高斯噪声信号Y中的每个点乘以一预设常数,以获取提高了所述窄带高斯噪声信号Y的带外抑制度后的信号Z2;
其中,所述预设常数大于1;
在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft,包括:
取索引n为1到N,在频域,对所述信号X_fft进行N点遍历,若检测到当前遍历点的位置索引n<[(f0-W)/(fs/N)]或者n>[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft取值为0,若检测到当前遍历点的位置索引n>=[(f0-W)/(fs/N)]且n<=[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft保持不变,其中,f0为中心频率。
2.根据权利要求1所述的方法,其特征在于,对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft,包括:
对所述高斯伪随机噪声信号X进行N点傅里叶变换,获得变换后的信号X_fft,其中,N>=fs/WS,fs为所述窄带高斯噪声信号Y的采样频率,W为所述窄带高斯噪声信号Y的带宽,WS为所述窄带高斯噪声信号Y的过渡带宽。
3.根据权利要求1所述的方法,其特征在于,对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y,包括:
对所述信号Y_fft进行傅里叶反变换,获得经反变换信号;
保留所述经反变换信号的实部,以获取窄带高斯噪声信号Y。
4.一种信号获取装置,运行一电子设备中,其特征在于,所述装置包括:
伪随机噪声获取模块,用于获取高斯伪随机噪声信号X;
傅里叶变换模块,用于对所述高斯伪随机噪声信号X进行傅里叶变换,获得变换后的信号X_fft;
幅值加权模块,用于在频域,对所述信号X_fft进行幅值加权,获得信号Y_fft;
傅里叶反变换模块,用于对所述信号Y_fft进行傅里叶反变换,获取窄带高斯噪声信号Y;
第一带外抑制模块,用于将所述窄带高斯噪声信号Y中的每个点除以一预设常数,以获取降低了所述窄带高斯噪声信号Y的带外抑制度后的信号Z1;
第二带外抑制模块,用于将所述窄带高斯噪声信号Y中的每个点乘以一预设常数,以获取提高了所述窄带高斯噪声信号Y的带外抑制度后的信号Z2;
其中,所述预设常数大于1;
其中,所述幅值加权模块,具体用于取索引n为1到N,在频域,对所述信号X_fft进行N点遍历,若检测到当前遍历点的位置索引n<[(f0-W)/(fs/N)]或者n>[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft取值为0,若检测到当前遍历点的位置索引n>=[(f0-W)/(fs/N)]且n<=[(f0+W)/(fs/N)],则将该索引位置对应的所述信号X_fft保持不变,其中,f0为中心频率。
5.根据权利要求4所述的装置,其特征在于,所述傅里叶变换模块,具体用于对所述高斯伪随机噪声信号X进行N点傅里叶变换,获得变换后的信号X_fft,其中,N>=fs/WS,fs为所述窄带高斯噪声信号Y的采样频率,W为所述窄带高斯噪声信号Y的带宽,WS为所述窄带高斯噪声信号Y的过渡带宽。
6.根据权利要求4所述的装置,其特征在于,所述傅里叶反变换模块包括:
反变换信号获取单元,用于对所述信号Y_fft进行傅里叶反变换,获得经反变换信号;
窄带高斯噪声信号获取单元,用于保留所述经反变换信号的实部,以获取窄带高斯噪声信号Y。
CN201710783735.XA 2017-09-01 2017-09-01 信号获取方法及装置 Active CN107493118B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710783735.XA CN107493118B (zh) 2017-09-01 2017-09-01 信号获取方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710783735.XA CN107493118B (zh) 2017-09-01 2017-09-01 信号获取方法及装置

Publications (2)

Publication Number Publication Date
CN107493118A CN107493118A (zh) 2017-12-19
CN107493118B true CN107493118B (zh) 2020-01-24

Family

ID=60651398

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710783735.XA Active CN107493118B (zh) 2017-09-01 2017-09-01 信号获取方法及装置

Country Status (1)

Country Link
CN (1) CN107493118B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116388696B (zh) * 2023-06-06 2023-08-15 苏州华兴源创科技股份有限公司 噪声信号的生成方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1337809A (zh) * 2000-08-15 2002-02-27 华为技术有限公司 一种产生窄带数字高频噪声信号的方法及装置
CN1667432A (zh) * 2004-03-11 2005-09-14 株式会社东芝 加权函数生成方法、基准信号生成方法、传输信号生成装置、信号处理装置及天线
KR20060020517A (ko) * 2004-08-31 2006-03-06 학교법인 동의학원 오버샘플링식 오에프디엠 변복조 방식 시스템의송수신장치 및 방법
CN101552758A (zh) * 2008-04-02 2009-10-07 扬智科技股份有限公司 正交频分复用系统中精符号时序同步方法
CN101651463A (zh) * 2008-08-13 2010-02-17 大唐移动通信设备有限公司 一种基带信号窄带干扰的抑制方法及装置
CN106249208A (zh) * 2016-07-11 2016-12-21 西安电子科技大学 基于分数阶傅里叶变换的噪声调幅干扰下信号检测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101771643B (zh) * 2008-12-29 2012-10-03 中国移动通信集团公司 一种信号处理的方法和设备
CN101977066B (zh) * 2010-10-21 2013-08-28 上海华为技术有限公司 窄带干扰的抑制方法、装置及网络系统
CN106340292B (zh) * 2016-09-08 2019-08-20 河海大学 一种基于连续噪声估计的语音增强方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1337809A (zh) * 2000-08-15 2002-02-27 华为技术有限公司 一种产生窄带数字高频噪声信号的方法及装置
CN1667432A (zh) * 2004-03-11 2005-09-14 株式会社东芝 加权函数生成方法、基准信号生成方法、传输信号生成装置、信号处理装置及天线
KR20060020517A (ko) * 2004-08-31 2006-03-06 학교법인 동의학원 오버샘플링식 오에프디엠 변복조 방식 시스템의송수신장치 및 방법
CN101552758A (zh) * 2008-04-02 2009-10-07 扬智科技股份有限公司 正交频分复用系统中精符号时序同步方法
CN101651463A (zh) * 2008-08-13 2010-02-17 大唐移动通信设备有限公司 一种基带信号窄带干扰的抑制方法及装置
CN106249208A (zh) * 2016-07-11 2016-12-21 西安电子科技大学 基于分数阶傅里叶变换的噪声调幅干扰下信号检测方法

Also Published As

Publication number Publication date
CN107493118A (zh) 2017-12-19

Similar Documents

Publication Publication Date Title
KR101883066B1 (ko) 레이더 및 소나 애플리케이션들에서 고정 소수점 고속 푸리에 변환들을 스케일링하는 방법
Mert et al. Detrended fluctuation thresholding for empirical mode decomposition based denoising
Gedalyahu et al. Multichannel sampling of pulse streams at the rate of innovation
JP4081526B2 (ja) 信号をサンプリングおよび/または再構成するためのサンプリング方法、再構成方法および機器
US9946687B2 (en) Fixed-point high dynamic range fast fourier transforms
CN101867387A (zh) 低于奈奎斯特速率采样下信号重构技术方案
JP2014077791A (ja) 複数信号の相互関係を示す方法及び試験測定装置
CN107493118B (zh) 信号获取方法及装置
US20080126459A1 (en) Method and device for generating a filter coefficient in real time
US20150123721A1 (en) System and Method for Gaussian Random Noise Generation
Looney et al. Measuring phase synchrony using complex extensions of EMD
CN107710624A (zh) 在非均匀采样模数转换器中增加采样
US10622971B2 (en) Method and device for adjusting passband width of filter
CN110755055A (zh) 一种用于确定脉搏波形的波形评价信息的方法与设备
KR20210035272A (ko) 위상 측정 정확도 향상을 위한 방법 및 장치
JP2009030983A (ja) 波形発生装置、波形生成装置、試験装置およびプログラム
Fleischmann et al. Implementation of a cross‐spectrum FFT analyzer for a phase‐noise test system in a low‐cost FPGA
JP4512835B2 (ja) コンスタレーション表示方法及び装置
Lainscsek et al. Analytical derivation of nonlinear spectral effects and 1/f scaling artifact in signal processing of real-world data
US9537690B1 (en) Method and apparatus for extraction of baseband waveform from amplitude modulated signal via time domain sampling
CN116994603A (zh) 一种音频信号的检测方法、装置及计算设备
JP7371868B2 (ja) 較正装置、変換装置、較正方法、および較正プログラム
US20230185872A1 (en) Fast Fourier Transforms With Incomplete Input Data Replacement
RU78615U1 (ru) Устройство для определения несущей частоты искаженного электрического сигнала
Lowdermilk et al. Finite arithmetic considerations for the FFT implemented in FPGA-based embedded processors in synthetic instruments

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