CN104133404A - 一种信号处理方法及装置 - Google Patents
一种信号处理方法及装置 Download PDFInfo
- Publication number
- CN104133404A CN104133404A CN201410351504.8A CN201410351504A CN104133404A CN 104133404 A CN104133404 A CN 104133404A CN 201410351504 A CN201410351504 A CN 201410351504A CN 104133404 A CN104133404 A CN 104133404A
- Authority
- CN
- China
- Prior art keywords
- signal
- amplitude
- default
- spectral line
- algorithm
- 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
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本申请提供了一种信号处理方法,包括:依据预设的调理条件对原始信号调理,得到模拟信号;将模拟信号进行A/D转换,得到数字信号;依据FIR算法和FFT算法对数字信号进行频谱分析,得到第一信号;依据预设修正算法对第一信号的频率和幅值修正处理,得到第二信号;依据预设补偿规则,对第二信号的幅值补偿,得到第三信号并输出至计算设备。采用该方法,在根据FIR算法和FFT算法对数字信号进行频谱分析后,进行频率和幅值的修正处理、对幅值进行补偿。对信号进行修正处理,消减了FFT算法的频谱泄漏现象,同时缓解了FFT算法的栅栏效应;而对信号进行幅值补偿,弥补了FIR带通波动和过渡带衰减的问题,提高了信号处理结果的精度。
Description
技术领域
本发明涉及信号处理领域,更具体的说,是涉及一种信号处理方法及装置。
背景技术
伴随着大规模集成电路的发展和控制算法的深入研究,数字信号处理以其精度高、可靠性高、灵活性大、易于大规模集成等优点,在各个领域均得到广泛应用。
数字信号处理是采用数值计算的方法对数字序列进行处理,将信号变换成符合需要的某种形式。如,对数字信号进行滤波以限制其他频带或滤除噪声和干扰,或将他们与其他信号进行分离,对信号进行频谱分析以了解信号的频谱组成,进而对信号进行识别。
现有技术中,FIR(Finite impact Response,有限长冲激响应)算法和FFT(Fast Fourier Transform,快速傅立叶变换)算法为频谱分析的主要方法。如图1所示的传统的信号处理流程,原始信号经由信号调理、A/D(Analog-to-digital,模拟/数字)转换后,送至数字信号处理单元进行频谱分析,然后计算。其中,该频谱分析环节由FIR和FFT组成。
由于低阶FIR存在带通波动大和过渡带衰减的问题,造成不同频率点的信号衰减系数不同,会造成临界频点的信号误判。而FFT算法在应用时存在影响测量准确性的频谱泄漏问题。
因此,采用现有技术中的方法在进行信号处理时,容易出现信号处理结果精度较低的问题。
发明内容
有鉴于此,本发明提供了一种信号处理方法及装置,解决了现有技术中信号处理结果精度较低的问题。
为实现上述目的,本发明提供如下技术方案:
一种信号处理方法,包括:
依据预设的调理条件对原始信号进行调理,得到模拟信号;
将所述模拟信号进行模拟/数字A/D转换,得到数字信号;
依据有限长冲激响应FIR算法和快速傅立叶变换FFT算法对所述数字信号进行频谱分析,得到第一信号;
依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
上述的方法,优选的,所述修正算法包括FFT差值修正算法,所述依据预设修正算法对所述第一信号的频率和幅值进行修正处理具体包括:
依据预设的FFT插值算法计算出所述第二信号的离散频谱;
选择所述离散频谱中预设个数的谱线;
计算各个谱线对应的谱线幅值;
依据预设的窗函数和所述谱线幅值、以及预设的拟合条件计算得到修正信息,并依据所述修正信息对所述第一信号的频率和幅值进行修正,得到第二信号。
上述的方法,优选的,所述依据预设补偿规则,对所述第二信号的幅值补偿具体包括:
从预设的FIR系数补偿表中查找与所述第二信号的频率对应的补偿信息;
依据所述补偿信息对所述第二信号的幅值进行查表补偿,得到第三信号。
上述的方法,优选的,所述选择所述离散频谱中预设个数的谱线包括:
确定所述离散谱线中幅值最大的谱线;
选择与所述最大谱线相邻的所述预设个数中剩余的谱线。
上述的方法,优选的,所述依据预设的窗函数和所述谱线幅值,以及预设的拟合条件,对所述第一信号的频率和幅值进行修正包括:
依据所述幅值分别构造第一函数和第二函数;
依据所述预设的窗函数、MATLAB的拟合函数,分别对所述第一函数和第二函数拟合,得到所述第一信号的实际频率值函数和实际幅值函数的多项式拟合公式,并利用所述多项式拟合公式对所述第一信号的频率和幅值进行修正处理。
上述的方法,优选的,所述窗函数采用布莱克曼窗函数。
上述的方法,优选的,所述确定所述离散谱线中幅值最大的谱线包括:
以所述预设位置为准,对所述离散谱线采用冒泡排序法排序,得到幅值最大的谱线。
上述的方法,优选的,所述预设个数包括:2个、3个或5个。
一种信号处理装置,包括:
调理模块,用于依据预设的调理条件对原始信号进行调理,得到模拟信号;
转换模块,用于将所述模拟信号进行A/D转换,得到数字信号;
频谱分析模块,用于依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;
修正模块,用于依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
补偿模块,用于依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
上述的装置,优选的,所述修正算法包括FFT差值修正算法,所述修正模块具体用于:
依据预设的FFT插值算法计算出所述第二信号的离散频谱;
选择所述离散频谱中预设个数的谱线;
计算各个谱线对应的谱线幅值;
依据预设的窗函数和所述谱线幅值、以及预设的拟合条件计算得到修正信息,并依据所述修正信息对所述第一信号的频率和幅值进行修正,得到第二信号。
经由上述的技术方案可知,与现有技术相比,本发明提供了一种信号处理方法,包括:依据预设的调理条件对原始信号进行调理,得到模拟信号;将所述模拟信号进行A/D转换,得到数字信号;依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。采用该方法,在根据FIR算法和FFT算法对数字信号进行频谱分析完成后,继续进行频率和幅值的修正处理,以及对幅值进行补偿。对信号进行修正处理,消减了FFT算法的频谱泄漏现象,同时缓解了FFT算法的栅栏效应;而对信号进行幅值补偿,弥补了FIR带通波动和过渡带衰减的问题,提高了信号处理结果的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为传统的信号处理流程图;
图2为本申请提供的一种信号处理方法实施例1的流程图;
图3为本申请提供的一种信号处理方法实施例2的流程图;
图4为本申请提供的一种信号处理方法实施例2中步骤S307的具体流程图;
图5为本申请提供的一种信号处理方法实施例3的流程图;
图6为本申请提供的一种信号处理方法实施例3中示例的响应图;
图7为本申请提供的一种信号处理方法实施例3中示例的频率与幅值衰减系数之间的关系图;
图8为本申请提供的一种信号处理装置实施例1的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图2所示的本申请实施例1提供的一种信号处理方法的流程图,包括:
步骤S201:依据预设的调理条件对原始信号进行调理,得到模拟信号;
其中,该调理条件可以为后续的A/D转换结构输入信号的条件。
具体的,该调理过程包括:功率放大、滤波、线性化补偿、隔离和/或保护等,调理的过程可根据实际情况进行设定,本实施例不做赘述。
步骤S202:将所述模拟信号进行A/D转换,得到数字信号;
其中,由于后续的FIR算法和FFT算法需要依据数字格式信号,而经过调理后得到的为模拟格式信号,因此,对该模拟信号进行A/D转换,使该信号转换为模拟信号。
步骤S203:依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;
其中,依据FIR算法和FFT算法依次对数字信号进行频谱分析,得到第一信号。
具体的,该频谱分析过程为:采用FFT算法实现对数字信号的变换,采用FIR算法实现对数字信号的滤波。
但是,采用FFT算法,该变换过程中容易出现频谱泄漏现象;而采用FIR算法的滤波存在带通波动大和过渡带衰减的问题。
因此,在后续的步骤S204-205,以对上述问题进行弥补。
步骤S204:依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
其中,在进行信号处理之前,预设修正算法。
具体的,根据该修正算法,对第一信号的频率和幅值分别进行修正,以克服FFT算法进行数字变换时出现的频谱泄漏的问题。
具体的修正过程在后续的实施例中进行说明,本实施例不做赘述。
步骤S205:依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
其中,在进行信号处理之前,预设补偿规则。
具体的,根据该补偿规则,对经过修正的第二信号进行幅值补偿,以克服采用FIR算法滤波导致的带通波动大和过渡带衰减的问题。
其中,该第三信号输出给后续的计算设备,以进行后续的计算过程。
需要说明的是,本申请中对数字信号进行频谱分析时,采用的为FIR算法和FFT算法,实际实施中,在进行频谱分析时,也可只采用FIR算法和FFT算法中的一种,相应的,当采用的为上述的两种算法时,执行步骤S204和205,对频谱分析完成的信号进行频率、幅值的修正,以及幅值补偿;当采用的只是FIR算法时,只执行步骤S205,对频谱分析完成的信号进行幅值补偿;当采用的只是FFT算法时,只执行步骤S204,对频谱分析完成的信号进行频率、幅值的修正。
当然,在实际实施中,在进行频谱分析时,还可只采用FIR算法和FFT算法外的其他算法,当频谱分析的结果也有带通波动大和过渡带衰减的问题时,则执行步骤S205;当频谱分析的结果也有频谱泄漏现象时,则执行步骤S204。
需要说明的是,实际实施中,本申请中的方法可移植到嵌入式软件或DSP(digital signal processor,数字信号处理器),将原始信号输入该嵌入式软件或DSP中实现该方法。
综上,本实施例提供了一种信号处理方法,包括:依据预设的调理条件对原始信号进行调理,得到模拟信号;将所述模拟信号进行A/D转换,得到数字信号;依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。采用该方法,在根据FIR算法和FFT算法对数字信号进行频谱分析完成后,继续进行频率和幅值的修正处理,以及对幅值进行补偿。对信号进行修正处理,消减了FFT算法的频谱泄漏现象,同时缓解了FFT算法的栅栏效应;而对信号进行幅值补偿,弥补了FIR带通波动和过渡带衰减的问题,提高了信号处理结果的精度。
在具体实现中,该修正算法可采用FFT差值修正算法。
如图3,示出的本申请提供的一种信号处理方法实施例2的流程图,包括:
步骤S301:依据预设的调理条件对原始信号进行调理,得到模拟信号;
步骤S302:将所述模拟信号进行A/D转换,得到数字信号;
步骤S303:依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;
其中,步骤S301-303与实施例1中的步骤S201-203一致,本实施例不做赘述。
步骤S304:依据预设的FFT插值算法计算出所述第一信号的离散频谱;
具体的,依据预设的第一采样频率对第一信号进行采样,得到离散信号;利用傅里叶变换和预设的窗函数,将该离散信号的时域信号转换为频域信号,得到该第一信号的连续频谱;对该连续频谱继续采样得到离散的若干条谱线,该采样频率与第一采样频率相关。
例如,该第一信号为频率f0、幅值为A、相角为的无限长正弦波信号x(t),该第一信号的表达式为:
以采样频率fs对x(t)采样后,得到如下形式的离散信号:
采用的窗函数的时域形式为w(n),其连续频谱为W(2πf),即
利用傅里叶变换,将是与信号红钻和为频域,得到x(t)信号的连续频谱为
忽略负频点处谱峰旁瓣影响,得到
对该链路频谱以采样频率进行采样,其中,N为数据长度,得到
该采样得到的算式(6)即为离散的若干条谱线。
步骤S305:选择所述离散频谱中预设个数的谱线;
具体的,包括:确定所述离散谱线中幅值最大的谱线;选择与所述最大谱线相邻的所述预设个数中剩余的谱线。
其中确定所述离散谱线中幅值最大的谱线的过程具体为:以所述预设位置为准,对所述离散谱线采用冒泡排序法排序,得到幅值最大的谱线。
其中,所述预设个数包括:2个、3个或5个。
其中,当为2个时,选取的谱线包括:幅值最大的谱线和该谱线附近的次大谱线;当为3个时,选取的谱线包括:幅值最大的谱线和该谱线前后相邻的2个谱线;当为5个时,选取的谱线包括:幅值最大的谱线和该谱线前后相邻的4个谱线,前后各两个谱线。
本实施例中以选取3条谱线为例进行说明。
设k0可能部位整数,其附近最大的谱线为k1,k1左侧谱线为k2、右侧谱线为k3,选择谱线为k1、k2和k3。
步骤S306:计算各个谱线对应的谱线幅值;
其中,根据步骤S305中选择的谱线计算该谱线与预设位置的偏移距离值,并且依据采样得到的离散公式计算得到选择的谱线对应的幅值。
例如,选择谱线为k1、k2和k3,各条谱线的幅值分别为 引入参数α=k0-k1。
由于k1为k0附近最大的谱线,所以α的取值范围为[-0.5,0.5],反应了k0偏离k1的距离,得到k1=k0-α、k2=k0-α-1、k2=k0-α+1,代入式(6)得
步骤S307:依据预设的窗函数和所述谱线幅值、以及预设的拟合条件计算得到修正信息,并依据所述修正信息对所述第一信号的频率和幅值进行修正,得到第二信号;
其中,依据上述的谱线以及谱线的幅值、结合预设的窗函数,计算得到一组与待处理的第以信号对应的参数组,将该参数组与预设的拟合公式进行拟合,得到幅值、频率的修正信息,根据该修正信息对第一信号的幅值和频率进行修正,得到第二信号。
如图4所示,该步骤S307的具体流程图,包括:
步骤S401:依据所述幅值分别构造第一函数和第二函数;
继续上例,依据谱线k1、k2和k3的幅值构造函数1和函数2:
函数1:
函数2:
步骤S402:依据所述预设的窗函数、MATLAB的拟合函数,分别对所述第一函数和第二函数拟合,得到所述第一信号的实际频率值函数和实际幅值函数的多项式拟合公式,并利用所述多项式拟合公式对所述第一信号的频率和幅值进行修正处理。
其中,本实施例中该窗函数采用布莱克曼窗函数。
令拟合出g(x)的反函数,即
α=g-1(β)=a0+a1β+a2β2+...+anβn (12)
其中,a0、a1、a2...an可取值任意实数。
具体可利用MATLAB(MATrix LABoratory,矩阵实验室)相关多项式插值工具拟合g(x)的反函数。
由此可知,第一信号的实际频率值为
f0=k0Δf=(k1+α)Δf (13)
同理,利用多项式插值拟合出式(11)中ρ(x)的反函数,即
A=N-1(y1+y2+y3)ρ-1(α)=N-1(y1+y2+y3)(b0+b1α+b2α2+...+bnαn) (14)
其中,b0、b1、b2...bn可取值任意实数。
其中,布莱克曼窗函数的时域表达式为
由于被处理的第一信号的频率一致,则α可计算得到,β也可计算得到,同理,也可计算得到。
变化采样频率,如f=f+0.001,重复执行上述步骤中对x(t)进行采样,以及后需的计算过程,最终得出一组α、β、γ值。
利用MATLAB的拟合函数采用polyfit,xx=polyfitg-1(β,α,3)、yy=polyfitg-1(α,γ,3)。
多项式拟合函数p=polyfitg-1(x,y,n)中,n表示多项式的最高阶数,x、y为将要拟合的数据,p为拟合得到的多项式系数p=[p1,p2,...,pn,pn+1],并且满足y=p1xn+p2xn-1+...+pnx+pn+1。
结合上述的拟合函数,得到式(12)的多项式拟合公式(频率值的拟合量)
α=1.9458-2.6332β+0.7089β2-0.0221β3 (16)
结合上述的拟合函数,得到式(14)的多项式拟合公式(幅值的拟合量)
A=N-1(y1+y2+y3)(2.1770+0.0006α+0.5631α2-0.0105α3) (17)
依据式(16)和(17)中的拟合量,对第一信号中的幅值A和频率f0进行补偿,得到第二信号。
需要说明的是,本实施例中采用了多项式逼近的方式,FFT修正算法数据计算量小,简化了计算过程。
并且,在进行补偿过程中,采用了采样的方法,降低了***的开销。
步骤S308:依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
其中,步骤S308与实施例1中的步骤S205一致,本实施例不做赘述。
需要说明的是,本实施例中的窗函数采用的为布莱克曼窗函数,当然,并不限定于此,实际实施中,还可采用矩形窗函数、汉宁窗函数或者海明窗函数,本实施例中不再详细解释。
综上,本实施例提供的一种信息处理方法,采用FFT差值修正算法对第一信号的幅值和频率进行修正,得到第二信号。由于该FFT差值修正算法,与前面步骤的对所述数字信号进行频谱分析采用的FFT算法相应,能够具有针对性的对采用FFT算法进行频谱分析时出现的频谱泄漏的问题。
在具体实现中,该可采用FIR系数补偿表对信号的幅值进行补偿。
如图5,示出的本申请提供的一种信号处理方法实施例3的流程图,包括:
步骤S501:依据预设的调理条件对原始信号进行调理,得到模拟信号;
步骤S502:将所述模拟信号进行A/D转换,得到数字信号;
步骤S503:依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;
步骤S504:依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
其中,步骤S501-504与实施例1中的步骤S201-204一致,本实施例不做赘述。
步骤S505:从预设的FIR系数补偿表中查找与所述第二信号的频率对应的补偿信息;
其中,该FIR系数补偿表中包含频率和衰减系数的对应关系。
具体的,根据该第二信号的频率,从该FIR系数补偿表中查找对应的补偿信息,即衰减系数。
需要说明的是,该系数补偿表可以为统一设置在计算设备中的,也可为预设在其他存储设备中,本实施例中不做限制。
步骤S506:依据所述补偿信息对所述第二信号的幅值进行查表补偿,得到第三信号并输出至计算设备。
其中,根据查找到的补偿信息及衰减系数,对第二信号的幅值进行补偿,以弥补采用FIR算法进行带通滤波带来的带通波动大和过渡带衰减。
下面介绍补偿表的建立过程。
以采样频率为500Hz、阶数为128、2~12Hz带通滤波为例进行详细说明。
利用MATLAB软件按照上述要求生成一组滤波系数;利用MATLAB画出该类别系数对应的频率响应图,该响应图如图6;将图6中频率响应横轴数据由归一化频率转换成实际频率,并且只关注过渡带及带通内的频率分量,纵轴数据转换成衰减系数,得到图7所示的频率与幅值衰减系数之间的关系图;观察图7可知,0Hz分量衰减为原始分量的0.5倍,7Hz分量衰减为原始分量的1倍,14Hz分量衰减至原始分量的0.25倍。显见,该滤波器通带阻带波动大,若不做处理,难以满足工程需要。将上述频率响应变化成表格的形式,以使完成FFT计算结束后对照频率值进行幅值的修正。
综上,本实施例提供的一种信息处理方法,采用FIR系数补偿表对信号的幅值进行补偿,得到第三信号。该FIR系数补偿表为预设的表格,当进行补偿时,只需根据待补偿信号的频率,从该补偿表中进行查找即可,数据计算量小,并且该FIR系数补偿表,与前面步骤的对所述数字信号进行频谱分析采用的FIR算法相应,能够具有针对性的对采用FIR算法进行频谱分析时出现的带通波动大和过渡带衰减的问题。
与上述本申请提供的一种信号处理方法实施例相对应的,本申请还提供了一种信号处理装置实施例。
如图8,示出了本申请提供的一种信号处理装置实施例1的结构示意图,包括:调理模块801、转换模块802、频谱分析模块803、修正模块804和补偿模块805;
调理模块801,用于依据预设的调理条件对原始信号进行调理,得到模拟信号;
其中,该调理模块801的调理条件可以为后续的A/D转换结构输入信号的条件。
具体的,该调理过程包括:功率放大、滤波、线性化补偿、隔离和/或保护等,调理的过程可根据实际情况进行设定,本实施例不做赘述。
转换模块802,用于将所述模拟信号进行A/D转换,得到数字信号;
其中,由于后续的FIR算法和FFT算法需要依据数字格式信号,而经过调理后得到的为模拟格式信号,因此,转换模块802对该模拟信号进行A/D转换,使该信号转换为模拟信号。
频谱分析模块803,用于依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;
其中,频谱分析模块803依据FIR算法和FFT算法依次对数字信号进行频谱分析,得到第一信号。
具体的,该频谱分析模块803频谱分析过程为:采用FFT算法实现对数字信号的变换,采用FIR算法实现对数字信号的滤波。
但是,采用FFT算法,该变换过程中容易出现频谱泄漏现象;而采用FIR算法的滤波存在带通波动大和过渡带衰减的问题。
因此,在后续的修正模块804和补偿模块805,以对上述问题进行弥补。
修正模块804,用于依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
其中,在进行信号处理之前,对修正模块804预设修正算法。
具体的,修正模块804根据该修正算法,对第一信号的频率和幅值分别进行修正,以克服FFT算法进行数字变换时出现的频谱泄漏的问题。
其中,该第三信号输出给后续的计算设备,以进行后续的计算过程。
其中,所述修正算法包括FFT差值修正算法,该修正模块804具体用于:依据预设的FFT插值算法计算出所述第二信号的离散频谱;选择所述离散频谱中预设个数的谱线;计算各个谱线与预设位置的偏移距离值,以及所述偏移距离值对应的谱线幅值;依据预设的窗函数和所述谱线幅值、以及预设的拟合条件计算得到修正信息,并依据所述修正信息对所述第一信号的频率和幅值进行修正,得到第二信号。
具体的修正过程在上述的方法实施例中已进行说明,本实施例不做赘述。
补偿模块805,用于依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
其中,在进行信号处理之前,对补偿模块805预设补偿规则。
具体的,补偿模块805根据该补偿规则,对经过修正的第二信号进行幅值补偿,以克服采用FIR算法滤波导致的带通波动大和过渡带衰减的问题。
需要说明的是,本申请中对数字信号进行频谱分析时,采用的为FIR算法和FFT算法,实际实施中,在进行频谱分析时,也可只采用FIR算法和FFT算法中的一种,相应的,当采用的为上述的两种算法时,触发修正模块804和补偿模块805,对频谱分析完成的信号进行频率、幅值的修正,以及幅值补偿;当采用的只是FIR算法时,只触发补偿模块805,对频谱分析完成的信号进行幅值补偿;当采用的只是FFT算法时,只触发修正模块804,对频谱分析完成的信号进行频率、幅值的修正。
当然,在实际实施中,在进行频谱分析时,还可只采用FIR算法和FFT算法外的其他算法,当频谱分析的结果也有带通波动大和过渡带衰减的问题时,则触发补偿模块805;当频谱分析的结果也有频谱泄漏现象时,则触发修正模块804。
需要说明的是,实际实施中,本申请中的方法可移植到嵌入式软件或DSP,将原始信号输入该嵌入式软件或DSP中实现该方法。
综上,本实施例提供了一种信号处理装置,包括:调理模块,用于依据预设的调理条件对原始信号进行调理,得到模拟信号;转换模块,用于将所述模拟信号进行A/D转换,得到数字信号;频谱分析模块,用于依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;修正模块,用于依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;补偿模块,用于依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。。采用该装置,在根据FIR算法和FFT算法对数字信号进行频谱分析完成后,继续进行频率和幅值的修正处理,以及对幅值进行补偿。对信号进行修正处理,消减了FFT算法的频谱泄漏现象,同时缓解了FFT算法的栅栏效应;而对信号进行幅值补偿,弥补了FIR带通波动和过渡带衰减的问题,提高了信号处理结果的精度。
Claims (10)
1.一种信号处理方法,其特征在于,包括:
依据预设的调理条件对原始信号进行调理,得到模拟信号;
将所述模拟信号进行模拟/数字A/D转换,得到数字信号;
依据有限长冲激响应FIR算法和快速傅立叶变换FFT算法对所述数字信号进行频谱分析,得到第一信号;
依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
2.根据权利要求1所述的方法,其特征在于,所述修正算法包括FFT差值修正算法,所述依据预设修正算法对所述第一信号的频率和幅值进行修正处理具体包括:
依据预设的FFT插值算法计算出所述第二信号的离散频谱;
选择所述离散频谱中预设个数的谱线;
计算各个谱线对应的谱线幅值;
依据预设的窗函数和所述谱线幅值、以及预设的拟合条件计算得到修正信息,并依据所述修正信息对所述第一信号的频率和幅值进行修正,得到第二信号。
3.根据权利要求1所述的方法,其特征在于,所述依据预设补偿规则,对所述第二信号的幅值补偿具体包括:
从预设的FIR系数补偿表中查找与所述第二信号的频率对应的补偿信息;
依据所述补偿信息对所述第二信号的幅值进行查表补偿,得到第三信号。
4.根据权利要求2所述的方法,其特征在于,所述选择所述离散频谱中预设个数的谱线包括:
确定所述离散谱线中幅值最大的谱线;
选择与所述最大谱线相邻的所述预设个数中剩余的谱线。
5.根据权利要求2所述的方法,其特征在于,所述依据预设的窗函数和所述谱线幅值,以及预设的拟合条件,对所述第一信号的频率和幅值进行修正包括:
依据所述幅值分别构造第一函数和第二函数;
依据所述预设的窗函数、矩阵实验室MATLAB的拟合函数,分别对所述第一函数和第二函数拟合,得到所述第一信号的实际频率值函数和实际幅值函数的多项式拟合公式,并利用所述多项式拟合公式对所述第一信号的频率和幅值进行修正处理。
6.根据权利要求2所述的方法,其特征在于,所述窗函数采用布莱克曼窗函数。
7.根据权利要求4所述的方法,其特征在于,所述确定所述离散谱线中幅值最大的谱线包括:
以所述预设位置为准,对所述离散谱线采用冒泡排序法排序,得到幅值最大的谱线。
8.根据权利要求2所述的方法,其特征在于,所述预设个数包括:2个、3个或5个。
9.一种信号处理装置,其特征在于,包括:
调理模块,用于依据预设的调理条件对原始信号进行调理,得到模拟信号;
转换模块,用于将所述模拟信号进行A/D转换,得到数字信号;
频谱分析模块,用于依据FIR算法和FFT算法对所述数字信号进行频谱分析,得到第一信号;
修正模块,用于依据预设修正算法对所述第一信号的频率和幅值进行修正处理,得到第二信号;
补偿模块,用于依据预设补偿规则,对所述第二信号的幅值补偿,得到第三信号并输出至计算设备。
10.根据权利要求9所述的装置,其特征在于,所述修正算法包括FFT差值修正算法,所述修正模块具体用于:
依据预设的FFT插值算法计算出所述第二信号的离散频谱;
选择所述离散频谱中预设个数的谱线;
计算各个谱线对应的谱线幅值;
依据预设的窗函数和所述谱线幅值、以及预设的拟合条件计算得到修正信息,并依据所述修正信息对所述第一信号的频率和幅值进行修正,得到第二信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410351504.8A CN104133404B (zh) | 2014-07-23 | 2014-07-23 | 一种信号处理方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410351504.8A CN104133404B (zh) | 2014-07-23 | 2014-07-23 | 一种信号处理方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104133404A true CN104133404A (zh) | 2014-11-05 |
CN104133404B CN104133404B (zh) | 2016-09-07 |
Family
ID=51806124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410351504.8A Active CN104133404B (zh) | 2014-07-23 | 2014-07-23 | 一种信号处理方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104133404B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104914306A (zh) * | 2014-12-03 | 2015-09-16 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号幅值测量方法 |
CN104914308A (zh) * | 2014-12-03 | 2015-09-16 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号相位测量方法 |
CN104931777A (zh) * | 2014-12-03 | 2015-09-23 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号频率测量方法 |
CN105137175A (zh) * | 2015-08-06 | 2015-12-09 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
CN105938508A (zh) * | 2016-05-13 | 2016-09-14 | 中国水利水电科学研究院 | 一种精确计算振动或压力脉动信号频率及幅值的方法 |
CN106597098A (zh) * | 2016-11-09 | 2017-04-26 | 深圳市鼎阳科技有限公司 | 一种频谱分析仪的数据处理方法和装置 |
CN107656139A (zh) * | 2017-09-19 | 2018-02-02 | 南京南瑞继保电气有限公司 | 一种信号确定方法、装置、设备和计算机可读存储介质 |
CN109426177A (zh) * | 2017-08-24 | 2019-03-05 | 中车株洲电力机车研究所有限公司 | 一种列车主动降噪方法及装置 |
CN113673091A (zh) * | 2021-07-26 | 2021-11-19 | 哈尔滨工业大学 | 一种机载嵌入式实时大点数时序信号处理方法 |
CN113962122A (zh) * | 2021-09-30 | 2022-01-21 | 北京智芯仿真科技有限公司 | 集成电路全波电磁仿真低频基准频点确定方法及*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1033614C (zh) * | 1992-03-16 | 1996-12-18 | Rca汤姆森许可公司 | 视频信号均衡器 |
CN1309295A (zh) * | 2000-12-04 | 2001-08-22 | 浙江大学 | 用于药物筛选的多功能细胞微生理计 |
CN1573930A (zh) * | 2003-06-12 | 2005-02-02 | 阿尔派株式会社 | 语音修正装置 |
JP2005055412A (ja) * | 2003-08-01 | 2005-03-03 | High Speed Signal Processing Laboratory Inc | Fft型周波数解析装置。 |
CN1742322A (zh) * | 2003-01-24 | 2006-03-01 | 索尼爱立信移动通讯股份有限公司 | 噪声减小和视听语音活动检测 |
CN103926456A (zh) * | 2014-04-10 | 2014-07-16 | 宁波恒力达科技有限公司 | 基于改进fft的闪变值计算方法及逆变器 |
-
2014
- 2014-07-23 CN CN201410351504.8A patent/CN104133404B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1033614C (zh) * | 1992-03-16 | 1996-12-18 | Rca汤姆森许可公司 | 视频信号均衡器 |
CN1309295A (zh) * | 2000-12-04 | 2001-08-22 | 浙江大学 | 用于药物筛选的多功能细胞微生理计 |
CN1742322A (zh) * | 2003-01-24 | 2006-03-01 | 索尼爱立信移动通讯股份有限公司 | 噪声减小和视听语音活动检测 |
CN1573930A (zh) * | 2003-06-12 | 2005-02-02 | 阿尔派株式会社 | 语音修正装置 |
JP2005055412A (ja) * | 2003-08-01 | 2005-03-03 | High Speed Signal Processing Laboratory Inc | Fft型周波数解析装置。 |
CN103926456A (zh) * | 2014-04-10 | 2014-07-16 | 宁波恒力达科技有限公司 | 基于改进fft的闪变值计算方法及逆变器 |
Non-Patent Citations (1)
Title |
---|
谢磊,等: "有源滤波器补偿方式的比较与分析", 《电工技术杂志》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104931777B (zh) * | 2014-12-03 | 2018-03-02 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号频率测量方法 |
CN104914308A (zh) * | 2014-12-03 | 2015-09-16 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号相位测量方法 |
CN104931777A (zh) * | 2014-12-03 | 2015-09-23 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号频率测量方法 |
CN104914306A (zh) * | 2014-12-03 | 2015-09-16 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号幅值测量方法 |
CN104914306B (zh) * | 2014-12-03 | 2017-11-03 | 许继集团有限公司 | 一种基于两条dft复数谱线的信号幅值测量方法 |
CN105137175A (zh) * | 2015-08-06 | 2015-12-09 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
CN105137175B (zh) * | 2015-08-06 | 2017-11-07 | 湖南大学 | 一种加窗短时傅里叶变换三点插值动态频率测量方法 |
CN105938508A (zh) * | 2016-05-13 | 2016-09-14 | 中国水利水电科学研究院 | 一种精确计算振动或压力脉动信号频率及幅值的方法 |
CN105938508B (zh) * | 2016-05-13 | 2019-04-02 | 中国水利水电科学研究院 | 一种精确计算振动或压力脉动信号频率及幅值的方法 |
CN106597098B (zh) * | 2016-11-09 | 2019-11-12 | 深圳市鼎阳科技有限公司 | 一种频谱分析仪的数据处理方法和装置 |
CN106597098A (zh) * | 2016-11-09 | 2017-04-26 | 深圳市鼎阳科技有限公司 | 一种频谱分析仪的数据处理方法和装置 |
CN109426177A (zh) * | 2017-08-24 | 2019-03-05 | 中车株洲电力机车研究所有限公司 | 一种列车主动降噪方法及装置 |
CN107656139A (zh) * | 2017-09-19 | 2018-02-02 | 南京南瑞继保电气有限公司 | 一种信号确定方法、装置、设备和计算机可读存储介质 |
CN113673091A (zh) * | 2021-07-26 | 2021-11-19 | 哈尔滨工业大学 | 一种机载嵌入式实时大点数时序信号处理方法 |
CN113673091B (zh) * | 2021-07-26 | 2024-07-16 | 哈尔滨工业大学 | 一种机载嵌入式实时大点数时序信号处理方法 |
CN113962122A (zh) * | 2021-09-30 | 2022-01-21 | 北京智芯仿真科技有限公司 | 集成电路全波电磁仿真低频基准频点确定方法及*** |
CN113962122B (zh) * | 2021-09-30 | 2022-05-31 | 北京智芯仿真科技有限公司 | 集成电路全波电磁仿真低频基准频点确定方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN104133404B (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104133404A (zh) | 一种信号处理方法及装置 | |
CA2819661C (en) | Correction of pressure signals measured during supersonic wind tunnel testing | |
CN106679741A (zh) | 基于涡街流量计抗干扰信号的处理方法及*** | |
CN102043091B (zh) | 数字化高精度相位检测器 | |
CN106199185B (zh) | 一种基于连续对数扫频的线性脉冲响应测量方法及*** | |
CN105445624A (zh) | 结合小波变换及曲线拟合的电缆故障定位方法 | |
CN103399589B (zh) | 一种电液加速度伺服***随机振动控制方法 | |
CN110807524A (zh) | 单通道信号盲源分离幅度校正方法 | |
CN108459087A (zh) | 一种应用于板结构损伤检测的多模态Lamb波模态分离方法 | |
CN105307095A (zh) | 一种基于fft的高分辨率音频频率测量方法 | |
Amundsen et al. | Elimination of temporal dispersion from the finite-difference solutions of wave equations in elastic and anelastic models | |
Ghoreishi et al. | Adomian decomposition method (ADM) for nonlinear wave-like equations with variable coefficient | |
Hao et al. | Urban hazardous chemicals pipeline leakage positioning method based on CELMD-MCKD | |
CN105425298A (zh) | 一种消除有限差分正演过程中数值频散的方法和装置 | |
Mboup | Parameter estimation via differential algebra and operational culculus | |
CN105115552B (zh) | 基于侧音测相的超声波流量测量方法及装置 | |
CN104484545B (zh) | 一种基于方差贡献率的泄流结构振动响应动态融合方法 | |
CN109506762A (zh) | 基于滤波器的水听器接收信号修正方法 | |
CN104931777A (zh) | 一种基于两条dft复数谱线的信号频率测量方法 | |
CN105865318A (zh) | 应用于lvdt位移传感器的非相敏数字解调***及方法 | |
CN103323063B (zh) | 超声波流量计及其时间差测量方法 | |
CN104501832A (zh) | 一种改进型实用惯性传感器降噪装置 | |
CN112505413B (zh) | 一种时频分析方法和*** | |
CN103399590B (zh) | 一种电液加速度伺服***的随机振动控制方法 | |
RU2489687C1 (ru) | Устройство для измерения акустического сопротивления твердых материалов |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP01 | Change in the name or title of a patent holder | ||
CP01 | Change in the name or title of a patent holder |
Address after: The age of 412001 in Hunan Province, Zhuzhou Shifeng District Road No. 169 Patentee after: ZHUZHOU CRRC TIMES ELECTRIC Co.,Ltd. Address before: The age of 412001 in Hunan Province, Zhuzhou Shifeng District Road No. 169 Patentee before: ZHUZH CSR TIMES ELECTRIC Co.,Ltd. |