[go: up one dir, main page]

CN109031438A - Anti-interference method and system for multi-channel receiver - Google Patents

Anti-interference method and system for multi-channel receiver Download PDF

Info

Publication number
CN109031438A
CN109031438A CN201810623258.5A CN201810623258A CN109031438A CN 109031438 A CN109031438 A CN 109031438A CN 201810623258 A CN201810623258 A CN 201810623258A CN 109031438 A CN109031438 A CN 109031438A
Authority
CN
China
Prior art keywords
signal
processing
channel receiver
filtering
filter
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.)
Pending
Application number
CN201810623258.5A
Other languages
Chinese (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.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
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 Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN201810623258.5A priority Critical patent/CN109031438A/en
Publication of CN109031438A publication Critical patent/CN109031438A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Noise Elimination (AREA)

Abstract

The invention discloses a kind of anti-interference methods and system for multi-channel receiver, which comprises synchronizes superposition processing to multi-channel receiver received signal;Comb filtering processing is carried out to the signal after the synchronized superposition processing;α-Trimed filtering processing is carried out to the signal handled through comb filtering.The anti-interference method and system of multi-channel receiver described in the embodiment of the present invention, by introducing three kinds of filtering techniques, a variety of interference in signal can be filtered out, realizes and broadband weak signal is acquired and handles art, improve the accuracy of electrical prospecting method and electrical measuring instrument,.

Description

用于多通道接收器的抗干扰方法及系统Anti-interference method and system for multi-channel receiver

技术领域technical field

本发明涉及检测领域,特别涉及一种用于多通道接收器的抗干扰方法及系统。The invention relates to the detection field, in particular to an anti-jamming method and system for a multi-channel receiver.

背景技术Background technique

电法勘探是矿产资源勘查的有效手段之一,其种类繁多、适应性强,被广泛应用于深部构造探测、矿产资源勘探以及水文及工程勘察等领域。Electrical prospecting is one of the effective means of mineral resource exploration. It has various types and strong adaptability, and is widely used in deep structure detection, mineral resource exploration, hydrological and engineering exploration and other fields.

根据场源的性质,电法勘探可分为天然源方法和人工源方法。天然源方法中的大地电磁法(Magnetotelluric,MT)是目前探测深度最大的电法勘探方法。可控源音频大地电磁法(Controlled Source Audio-frequency Magnetotelluric, CSAMT)以人工场源代替天然场源,沿用MT的观测方式,但是克服了MT 场源随机性的缺点,信号强度也大为提高,除探测深度较MT小以外,其工作效率、精度以及横纵向分辨率都明显提高。According to the nature of field sources, electrical prospecting can be divided into natural source methods and artificial source methods. Among the natural source methods, the magnetotelluric method (Magnetotelluric, MT) is currently the electrical prospecting method with the largest detection depth. Controlled Source Audio-frequency Magnetotelluric (CSAMT) replaces natural field sources with artificial field sources, and follows the observation method of MT, but overcomes the shortcomings of randomness of MT field sources, and greatly improves the signal strength. Except that the detection depth is smaller than that of MT, its work efficiency, precision and horizontal and vertical resolution are all significantly improved.

现有的电法勘探方法中及电法勘探仪器通常采用多通道接收器,但是现有的电法勘探方法的精确度仍然不高,原因是外部环境产生的干扰较大,难以提取有用信号。外部干扰主要有由环境产生的随机干扰、电力线产生的工频干扰以及电火花等产生的尖峰干扰等。现有技术中的抗干扰方法还有很大的不足,干扰过大,导致电法勘探精确度不高。Existing electrical prospecting methods and electrical prospecting instruments usually use multi-channel receivers, but the accuracy of the existing electrical prospecting methods is still not high, because the interference generated by the external environment is large, and it is difficult to extract useful signals. External interference mainly includes random interference generated by the environment, power frequency interference generated by power lines, and spike interference generated by electric sparks. The anti-jamming method in the prior art still has great deficiencies, and the interference is too large, resulting in low accuracy of electrical prospecting.

发明内容Contents of the invention

有鉴于此,本申请实施例提供了一种用于多通道接收器的抗干扰方法及系统,能滤除信号中的多种干扰,提高电法勘探方法及电法勘探仪器的精度,本发明实施例一方面提供了一种用于多通道接收器的抗干扰方法,所述方法对多通道接收器接收的信号进行同步叠加处理;对所述经同步叠加处理过后的信号进行梳状滤波处理;对所述经梳状滤波处理的信号进行α-Trimed滤波处理。In view of this, the embodiment of the present application provides an anti-interference method and system for a multi-channel receiver, which can filter out various interferences in the signal, and improve the accuracy of the electrical prospecting method and the electrical prospecting instrument. On the one hand, the embodiment provides an anti-jamming method for a multi-channel receiver, the method performs synchronous superposition processing on signals received by the multi-channel receiver; comb filter processing is performed on the signal after synchronous superposition processing ; performing α-Trimed filter processing on the comb filter processed signal.

可选的,所述对多通道接收器接收的信号进行同步叠加处理,包括:将多通道接收器接收的信号分为多个信号段;将每个信号段中的信号分为多个信号节;将每个信号段中的多个信号节进行叠加,获取所述信号段的幅度与相位信息;将多个叠加后的信号段进行叠加,获取所述多通道接收的信号的幅度与相位信息。Optionally, performing synchronous superposition processing on the signal received by the multi-channel receiver includes: dividing the signal received by the multi-channel receiver into multiple signal segments; dividing the signal in each signal segment into multiple signal segments ; Superimpose multiple signal sections in each signal segment to obtain amplitude and phase information of the signal segment; superimpose multiple superimposed signal segments to obtain amplitude and phase information of the multi-channel received signal .

可选的,对所述经同步叠加处理过后的信号进行梳状滤波处理,包括:对多通道接收器接收的信号进行延时,成为延时信号;将所述多通道接收器接收的信号与所述延时信号叠加,以产生相位抵消。Optionally, performing comb filter processing on the signal after synchronous superposition processing includes: delaying the signal received by the multi-channel receiver to become a delayed signal; combining the signal received by the multi-channel receiver with The delayed signals are added to produce phase cancellation.

可选的,对所述经同步叠加处理过后的信号进行梳状滤波处理,包括:使用带有乘加器的FPGA以实现对多通道接收器接收的信号进行梳状滤波处理。Optionally, performing comb filter processing on the synchronously superimposed signal includes: using an FPGA with a multiplier-adder to implement comb filter processing on signals received by the multi-channel receiver.

可选的,对所述经同步叠加处理过后的信号进行梳状滤波处理,包括:使用多个带有乘加器的FPGA级联,以实现对多通道接收器接收的信号进行梳状滤波处理。Optionally, performing comb filter processing on the signal after synchronous superposition processing includes: using multiple FPGA cascades with multiplier-adders to implement comb filter processing on signals received by the multi-channel receiver .

可选的,对所述经梳状滤波处理的信号进行α-Trimed滤波处理,包括:设置窗口长度N和参数α,其中α表示被修整数据样本的百分比;根据窗口长度N和参数α,对多通道接收器接收的信号进行处理;将移动窗口一个元素,对多通道接收器接收的信号进行处理。Optionally, performing α-Trimed filter processing on the signal processed by comb filtering, including: setting window length N and parameter α, wherein α represents the percentage of trimmed data samples; according to window length N and parameter α, for The signal received by the multi-channel receiver is processed; the window is moved by one element to process the signal received by the multi-channel receiver.

可选的,所述根据窗口长度N和参数α,对多通道接收器接收的信号进行处理,包括:根据窗口长度N和参数α,滤除多余数据;对剩余数据做平均处理。Optionally, the processing the signal received by the multi-channel receiver according to the window length N and the parameter α includes: filtering out redundant data according to the window length N and the parameter α; performing averaging processing on the remaining data.

可选的,对所述经梳状滤波处理的信号进行α-Trimed滤波处理之前,所述方法还包括:获取多通道接收器接收的信号对应的数值;在信号头和信号尾分别添加附加信号,且添加的附加信号为对称添加;对所述多通道接收器接收的信号对应的数值按照从小到大排序。Optionally, before performing α-Trimed filtering on the comb-filtered signal, the method further includes: acquiring a value corresponding to the signal received by the multi-channel receiver; adding additional signals at the beginning and end of the signal respectively , and the added additional signal is a symmetrical addition; the values corresponding to the signals received by the multi-channel receiver are sorted from small to large.

可选的,在所述信号头和信号尾分别添加附加信号,包括:若所述窗口长度N为偶数,则在信号头和信号尾各添加N/2个数据;若所述窗口长度N为奇数,则在信号头和信号尾各添加(N-1)/2个数据。Optionally, adding additional signals to the signal header and the signal tail respectively, including: if the window length N is an even number, adding N/2 data to the signal header and the signal tail respectively; if the window length N is If the number is odd, add (N-1)/2 data to the signal head and signal tail respectively.

本发明实施例另一方面,提供了一种用于多通道接收器的抗干扰系统,所述系统包括:同步叠加滤波器,用于对多通道接收器接收的信号进行同步叠加处理;梳状滤波器,用于对所述经同步叠加处理过后的信号进行梳状滤波处理;α-Trimed滤波器,用于对所述经梳状滤波处理的信号进行α-Trimed滤波处理。Another aspect of the embodiment of the present invention provides an anti-jamming system for a multi-channel receiver, the system includes: a synchronous superposition filter for synchronous superposition processing on signals received by the multi-channel receiver; A filter is used to perform comb filter processing on the signal processed by synchronous superposition; an α-Trimed filter is used to perform α-Trimed filter processing on the signal processed by comb filter.

本发明实施例所述的多通道接收器的抗干扰方法和系统,通过引入三种滤波技术,能滤除信号中的多种干扰,实现了对宽频带弱信号采集与处理术,提高电法勘探方法及电法勘探仪器的精确度。The anti-jamming method and system of the multi-channel receiver described in the embodiment of the present invention can filter out various kinds of interference in the signal by introducing three kinds of filtering techniques, realize the acquisition and processing of wide-band weak signals, and improve the electrical method. Accuracy of prospecting methods and electrical prospecting instruments.

附图说明Description of drawings

为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings that need to be used in the description of the embodiments. Obviously, the drawings in the following description are some embodiments of the present invention. For those skilled in the art, other drawings can also be obtained based on these drawings without creative effort.

图1为本发明实施例提供的用于多通道接收器的抗干扰方法的流程图;FIG. 1 is a flowchart of an anti-jamming method for a multi-channel receiver provided by an embodiment of the present invention;

图2所示为本发明实施例的对多通道接收器接收的信号进行同步叠加处理的示意图;FIG. 2 is a schematic diagram of performing synchronous superposition processing on signals received by a multi-channel receiver according to an embodiment of the present invention;

图3所示为本发明实施例的FPGA乘加器的示意图;Fig. 3 shows the schematic diagram of the FPGA multiply-accumulator of the embodiment of the present invention;

图4所示为本发明实施例的梳状滤波器的状态机示意图;FIG. 4 is a schematic diagram of a state machine of a comb filter according to an embodiment of the present invention;

图5所示为本发明实施例的一维信号延展示意图;FIG. 5 is a schematic diagram of one-dimensional signal extension according to an embodiment of the present invention;

图6所示为本发明实施例的二维信号延展示意图;FIG. 6 is a schematic diagram of two-dimensional signal extension according to an embodiment of the present invention;

图7为本发明实施例的α-trimmed滤波器滤波的流程示意图。FIG. 7 is a schematic flowchart of α-trimmed filter filtering according to an embodiment of the present invention.

具体实施方式Detailed ways

为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In order to make the purpose, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments It is only some embodiments of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.

在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“设置”、“耦合”、“连接”应做广义理解,例如,“连接”可以是直接连接,也可以通过中间媒介间接连接,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。In the description of the present invention, it should also be noted that, unless otherwise specified and limited, the terms "set", "coupled" and "connected" should be interpreted in a broad sense, for example, "connected" can be a direct connection, or It can be connected indirectly through an intermediary, or it can be an internal communication between two components. Those of ordinary skill in the art can understand the specific meanings of the above terms in the present invention according to specific situations.

在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。In this document, relational terms such as first and second etc. are used only to distinguish one entity or operation from another without necessarily requiring or implying any such relationship between these entities or operations. Actual relationship or sequence. Furthermore, the term "comprises", "comprises" or any other variation thereof is intended to cover a non-exclusive inclusion such that a process, method, article or apparatus comprising a set of elements includes not only those elements, but also includes elements not expressly listed. other elements of or also include elements inherent in such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article or apparatus comprising said element.

本发明实施例公开了一种用于多通道接收器的抗干扰方法,所述方法包括:The embodiment of the present invention discloses an anti-jamming method for a multi-channel receiver, the method comprising:

S110,对多通道接收器接收的信号进行同步叠加处理;S110, performing synchronous superposition processing on signals received by the multi-channel receiver;

S120,对所述经同步叠加处理过后的信号进行梳状滤波处理;S120, performing comb filter processing on the signal after synchronous superposition processing;

S130,对所述经梳状滤波处理的信号进行α-Trimed滤波处理。S130. Perform α-Trimed filtering on the comb-filtered signal.

叠加平均处理技术是一种有效的消噪信号处理方法,本发明实施例利用系统高精度同步时钟,实现硬件电路的叠加滤波技术。具体的,如图2所示,所述对多通道接收器接收的原始信号进行同步叠加处理,包括:The superimposed average processing technology is an effective noise-removing signal processing method. The embodiment of the present invention utilizes the high-precision synchronous clock of the system to realize the superimposed filtering technology of the hardware circuit. Specifically, as shown in Figure 2, the synchronous superposition processing of the original signal received by the multi-channel receiver includes:

S111,将多通道接收器接收的原始信号分为多个信号段;S111, divide the original signal received by the multi-channel receiver into multiple signal segments;

S112,将每个信号段中的信号分为多个信号节;S112. Divide the signal in each signal segment into multiple signal segments;

S113,将每个信号段中的多个信号节进行叠加,获取所述信号段的幅度与相位信息;S113. Superimpose multiple signal segments in each signal segment to obtain amplitude and phase information of the signal segment;

S114,将多个叠加后的信号段进行叠加,获取所述多通道接收的信号的幅度与相位信息。S114. Superimpose multiple superimposed signal segments, and acquire amplitude and phase information of signals received by the multi-channels.

其中对信号进行叠加的过程如下,在原始输入信号中,既包含有有用信号也包含有随机干扰,:The process of superimposing signals is as follows. In the original input signal, both useful signals and random interference are included:

设原始输入信号是f(t),可表示为f(t)=s(t)+n(t)。Suppose the original input signal is f(t), which can be expressed as f(t)=s(t)+n(t).

其中,s(t)为振幅恒定的周期信号,n(t)为随机干扰。Among them, s(t) is a periodic signal with constant amplitude, and n(t) is random interference.

随机干扰n(t)服从高斯分布,其均值为零,方差为δ2,也称作高斯白噪声。根据高斯白噪声的性质可知此噪声的功率为δ2,即输入信号的信噪比为:The random disturbance n(t) obeys Gaussian distribution, its mean is zero, and its variance is δ 2 , which is also called Gaussian white noise. According to the properties of Gaussian white noise, it can be known that the power of this noise is δ 2 , that is, the signal-to-noise ratio of the input signal is:

SNRin=S/δ2 SNR in = S/δ 2

S为总功率,进行m次叠加平均后,输出信号为:S is the total power. After m times of stacking and averaging, the output signal is:

式中,T为时间间隔,f(tk+iT)为第i点的第K次样值,n(t)为噪声信号有效值。对输入信号进行上述叠加平均处理,其中有效信号经m次叠加平均后其值不变,而噪声叠加平均后变为:In the formula, T is the time interval, f(tk+iT) is the Kth sample value of the i-th point, and n(t) is the effective value of the noise signal. The above-mentioned superposition and average processing is performed on the input signal, and the value of the effective signal remains unchanged after m times of superposition and averaging, while the noise superposition and averaging becomes:

则输出信噪比变为:Then the output signal-to-noise ratio becomes:

式中,NO为叠加平均后输出噪声的功率。In the formula, N O is the power of the output noise after superposition and averaging.

采用本发明实施例所述的信号叠加处理后,周期信号经过m次同步叠加后,信噪比提高为原来的m倍。因此,周期性信号或可重复性信号经过多次取样积累后,其信噪比会不断提高,积累次数越多,信噪比改善越好。为了从强噪声的背景中提取出微弱的信号,优选的,所述叠加平均次数大于50次。After the signal superposition processing described in the embodiment of the present invention is adopted, after the periodic signal is synchronously superimposed m times, the signal-to-noise ratio is increased to m times of the original. Therefore, the signal-to-noise ratio of a periodic signal or repeatable signal will continue to improve after multiple sampling accumulations, and the more times of accumulation, the better the signal-to-noise ratio improvement. In order to extract a weak signal from a background of strong noise, preferably, the average number of superposition times is greater than 50.

在本发明的一个实施例中,针对50Hz工频干扰较大的场景,为了进一步提高信号信噪比,在进行步骤S110前还要对模拟信号进行一步模拟滤波。考虑到模拟滤波器精度较低,对50Hz左右频率范围内有用信号的损伤较大,因此进一步设置步骤S120对完成步骤S110之后的数字信号进行滤波,所述数字滤波器可以通过设计合理的系数在完成高精度滤波的同时对50Hz左右频率范围内有用信号损伤较小。In one embodiment of the present invention, in order to further improve the signal-to-noise ratio in a scene where 50 Hz power frequency interference is relatively large, a step of analog filtering is performed on the analog signal before step S110 is performed. Considering that the precision of the analog filter is low, and the damage to the useful signal in the frequency range of about 50 Hz is relatively large, step S120 is further set to filter the digital signal after step S110 is completed, and the digital filter can be designed with reasonable coefficients in While completing high-precision filtering, there is less damage to useful signals in the frequency range of about 50Hz.

本发明实施例中,步骤S120中,对所述经同步叠加处理过后的信号进行梳状滤波处理,包括:In the embodiment of the present invention, in step S120, comb filter processing is performed on the signal after synchronous superposition processing, including:

S121,对所述经同步叠加处理过后的信号进行延时,成为延时信号;S121. Delay the synchronously superimposed signal to become a delayed signal;

S122,将所述经同步叠加处理过后的信号与所述延时信号叠加,以产生相位抵消。S122. Superimpose the signal processed by synchronous superposition with the delayed signal to generate phase cancellation.

现有技术中,数字滤波器分为IIR和FIR两大类。其中FIR滤波器可以得到严格的线性相位,但其传递函数极点是定在原点的,只能通过改变零点位置来改变性能,其阶数相对较高。而IIR滤波器存在反馈,可以设计极点,从而能够更加高效的完成滤波。本发明通过实验发现,在相同设计指标下,FIR滤波器的阶数是IIR滤波器的5~10倍,因此本发明的实施例采用的是IIR数字梳状滤波器。In the prior art, digital filters are classified into two categories: IIR and FIR. Among them, the FIR filter can obtain a strict linear phase, but its transfer function pole is fixed at the origin, and the performance can only be changed by changing the zero position, and its order is relatively high. The IIR filter has feedback, and the pole can be designed, so that the filtering can be completed more efficiently. The present invention finds through experiments that under the same design index, the order of the FIR filter is 5-10 times that of the IIR filter, so the embodiment of the present invention adopts the IIR digital comb filter.

在本发明的一个实施例中,所述谐波梳状滤波器在时间域用如下差分方程表示:In one embodiment of the present invention, the harmonic comb filter is represented by the following difference equation in the time domain:

对方程进行Z变换,其频率域表达方式如下:Z-transform the equation, and its expression in the frequency domain is as follows:

所述谐波梳状滤波器在频率域的表达方式如下:The expression of the harmonic comb filter in the frequency domain is as follows:

频率域的另一种表达方式如下:Another expression in the frequency domain is as follows:

式中,zk和zj分别为传递函数的零点和极点,采用频率域的表达方式方式能够较为直观的设计数字滤波器。In the formula, z k and z j are the zero and pole points of the transfer function respectively, and the digital filter can be designed more intuitively by using the expression in the frequency domain.

本发明所述的数字梳状滤波器解决了50Hz的K次谐波,但是对通带内信号还是会造成损伤,为了解决此问题,优选的,在其每个零点的单位圆内加一个极点,所述梳状滤波器的频率响应由一系列规律分布的峰组成,使得梳状滤波器能够使一个信号与它的延时信号叠加,从而产生相位抵消。The digital comb filter of the present invention solves the K-th harmonic of 50 Hz, but still causes damage to the signal in the passband. In order to solve this problem, preferably, a pole is added in the unit circle of each zero point , the frequency response of the comb filter consists of a series of regularly distributed peaks, enabling the comb filter to superimpose a signal with its delayed signal, resulting in phase cancellation.

优选的,离散时间系统中的梳状滤波器满足下式:Preferably, the comb filter in the discrete-time system satisfies the following formula:

在本发明的一个实施例中,所述梳状滤波处理通过FPGA实现,包括:In one embodiment of the present invention, the comb filter processing is realized by FPGA, including:

使用带有乘加器的FPGA以实现对信号进行梳状滤波处理。优选的,设置有多个带有乘加器的FPGA级联来实现。所述实现方法如下:Use an FPGA with a multiply-adder to implement comb-filtering of the signal. Preferably, a plurality of FPGA cascades with multiplier-adders are provided for realization. The implementation method is as follows:

1)设信号采样率为fs,期望被滤除频率的基波频率为fo,则有:1) Let the signal sampling rate be fs, and the fundamental frequency of the frequency to be filtered out is fo, then:

fs=fo*Dfs=fo*D

2)上述滤波器对信号的幅度衰减率约为30dB。2) The amplitude attenuation rate of the above filter to the signal is about 30dB.

在接收机采集系统中,有连续的24kHz以及由此数据流下变频抽样得到的 2.4kHz数据流,考虑到需要滤除信号的基波频率为50Hz,本发明实施例对 24kHz数据流滤波时滤波器阶数D=96,对2.4kHz数据流D=48。In the receiver acquisition system, there are continuous 24kHz and the 2.4kHz data stream obtained by down-conversion sampling of the data stream. Considering that the fundamental frequency of the signal to be filtered is 50Hz, the embodiment of the present invention filters the 24kHz data stream. Order D=96, D=48 for 2.4kHz data stream.

由于一组滤波器的信号衰减率为30dB,不能满足实际采集中对信号滤波的需求,本发明设计了两组滤波器级联:Since the signal attenuation rate of one set of filters is 30dB, it cannot meet the requirements of signal filtering in actual collection, the present invention designs two sets of filter cascades:

使用上述滤波器,其信号衰减率约为60dB。Using the above filter, the signal attenuation rate is about 60dB.

在24kHz、2.4kHz的数据流中实现上述梳状滤波器,需要的滤波器阶数分别为192,96。系统电路中FPGA芯片的剩余资源通常难以实现192阶的滤波器,因此本发明的实施例对2.4kHz的数据流进行滤波。To implement the above comb filter in the 24kHz and 2.4kHz data streams, the required filter orders are 192 and 96 respectively. It is generally difficult to implement a 192-order filter with the remaining resources of the FPGA chip in the system circuit, so the embodiment of the present invention filters the 2.4kHz data stream.

本发明的实施例中,所述FPGA实现该滤波器需要一个乘加器和一组数据缓冲区。所述乘加器的CLK接入系统的总时钟。In the embodiment of the present invention, the FPGA needs a multiply-adder and a set of data buffers to implement the filter. The CLK of the multiplier-adder is connected to the overall clock of the system.

考虑到实际应用中,FPGA的RAM资源紧张,因此在本发明的实施例中,为了节省FPGA内部的RAM资源,使用寄存器来做缓存区,即二维寄存器数组,缓存由寄存器的移位运算实现。Considering that in practical application, the RAM resource of FPGA is tight, so in the embodiment of the present invention, in order to save the RAM resource inside FPGA, use register to do cache area, namely two-dimensional register array, cache is realized by the shift operation of register .

梳状滤波器的计算过程中,逐个通道依次进行,由于前级数据产生环节的设计,相邻通道数据产生时间间隔需要足够完成计算。每个通道的计算由一组状态机控制,所述状态机有8种状态,状态变化关系如图4所示。In the calculation process of the comb filter, it is carried out channel by channel. Due to the design of the previous data generation link, the time interval between adjacent channel data generation needs to be sufficient to complete the calculation. The calculation of each channel is controlled by a group of state machines, and the state machine has 8 states, and the state change relationship is shown in FIG. 4 .

本发明所述梳状滤波器的还包括,为了一次完成运算,提高计算速度,优选的使用多个乘加器级联,进一步优选的为5个,这样恰好可以一次完成图4 的计算,提高计算速度。The comb filter of the present invention also includes, in order to complete the operation at one time, improve the calculation speed, preferably use a plurality of multiplier-adder cascades, more preferably 5, so that the calculation of Fig. 4 can be completed at one time, and the calculation speed can be improved. Calculate speed.

本发明所述梳状滤波器的还包括,为了进一步提高阻带的衰减率,满足更复杂条件下的使用需求,本发明优选的使用片内RAM作为数据缓存区,并让多个梳状滤波器级联,提高其阶数。The comb filter of the present invention also includes that in order to further improve the attenuation rate of the stop band and meet the use requirements under more complex conditions, the present invention preferably uses on-chip RAM as a data buffer area, and allows multiple comb filters cascaded to increase its order.

完成S120的梳状滤波处理后,考虑到频率域电磁法信号受尖峰噪声干扰较大,本发明对常见尖峰噪声进行了分析研究,包括大气噪声干扰,如雷电产生的火花放电,属于脉冲宽带干扰,其覆盖频率从1Hz至百MHz以上;太阳噪声干扰(即太阳黑子的辐射噪声),在太阳黑子活动期,黑子爆发产生比平稳期高数千倍的强烈噪声;工业设备噪声干扰,如感应加热设备、高频电焊机等带来的干扰;电力设备噪声干扰,包括伺服电机、电钻以及继电器等设备通断产生电流剧变而引起的电磁干扰;汽车、内燃机点火系统产生的宽带噪声干扰。本发明将频率域电磁法常见的尖峰噪声干扰归类至脉冲型游散电流噪声,并将其模拟为三种类型:矩形脉冲、半正弦型脉冲、三角型脉冲。After completing the comb filter processing of S120, considering that the frequency-domain electromagnetic method signal is greatly disturbed by spike noise, the present invention analyzes and studies common spike noise, including atmospheric noise interference, such as spark discharge generated by lightning, which belongs to pulse broadband interference , whose coverage frequency ranges from 1 Hz to more than 100 MHz; solar noise interference (that is, the radiation noise of sunspots), during the active period of sunspots, the explosion of sunspots produces strong noise thousands of times higher than that of the stable period; industrial equipment noise interference, such as induction Interference caused by heating equipment, high-frequency electric welding machines, etc.; electrical equipment noise interference, including electromagnetic interference caused by drastic current changes caused by switching on and off of servo motors, electric drills, and relays; broadband noise interference generated by ignition systems of automobiles and internal combustion engines. The present invention classifies common peak noise interference in the frequency domain electromagnetic method into pulse-type stray current noise, and simulates it into three types: rectangular pulse, half-sine pulse and triangular pulse.

本发明通过研究发现,尖峰脉冲主要影响低频电磁信号。脉冲幅度越大,对低频信号影响越大,同时增加对高频信号的影响;脉冲幅度越小,对低频信号影响越小,同时对高频信号影响降低,甚至影响可以忽略。当脉冲宽度变宽时,对低频信号影响大,对高频信号影响也大;当脉冲宽度变窄时,影响范围扩展到高频信号。因此,通过上述尖峰噪声的时域及频域数学模型分析,本发明实施例采用α-trimmed滤波方法来对频率域电磁法的采集数据进行处理。The present invention finds through research that the peak pulse mainly affects low-frequency electromagnetic signals. The larger the pulse amplitude, the greater the impact on low-frequency signals, and at the same time increase the impact on high-frequency signals; the smaller the pulse amplitude, the smaller the impact on low-frequency signals, while the impact on high-frequency signals is reduced, or even negligible. When the pulse width becomes wider, it has a greater impact on low-frequency signals and also has a greater impact on high-frequency signals; when the pulse width becomes narrower, the range of influence extends to high-frequency signals. Therefore, through the time-domain and frequency-domain mathematical model analysis of the above-mentioned peak noise, the embodiment of the present invention uses the α-trimmed filtering method to process the data collected by the frequency-domain electromagnetic method.

本发明实施例中,步骤S130对所述经梳状滤波处理的信号进行α-Trimed 滤波处理,包括:In the embodiment of the present invention, step S130 performs α-Trimed filtering on the comb-filtered signal, including:

S131,设置窗口长度N和参数α,其中α表示被修整数据样本的百分比;S131, setting the window length N and parameter α, wherein α represents the percentage of the trimmed data sample;

S132,根据窗口长度N和参数α,对信号进行处理,所述处理步骤为:根据窗口长度N和参数α,滤除多余数据;对剩余数据做平均处理。S132. Process the signal according to the window length N and the parameter α. The processing steps include: filtering out redundant data according to the window length N and the parameter α; performing averaging processing on the remaining data.

S133,将窗口移动一个元素,对信号进行处理。S133, move the window by one element, and process the signal.

可选的,所述对信号进行处理之前,所述方法还包括:获取信号对应的数值;在信号头和信号尾分别添加附加信号,优选的,所述添加的附加信号为对称添加;对所述信号对应的数值按照从小到大排序。Optionally, before the signal is processed, the method further includes: acquiring a value corresponding to the signal; adding additional signals to the signal head and the signal tail respectively, preferably, the added additional signals are symmetrical addition; The values corresponding to the above signals are sorted from small to large.

可选的,所述在信号头和信号尾分别添加附加信号,包括:若所述窗口长度N为偶数,则在信号头和信号尾各添加N/2个数据;若所述窗口长度N为奇数,则在信号头和信号尾各添加(N-1)/2个数据。Optionally, adding additional signals at the signal header and the signal tail respectively includes: if the window length N is an even number, adding N/2 data at the signal header and the signal tail respectively; if the window length N is If the number is odd, add (N-1)/2 data to the signal head and signal tail respectively.

本发明实施例所采用的α-trimmed均值滤波器具有动态滤波器特性,其信号的任意元素值取决于该元素周边的值。选取一定大小的窗口去掉异类元素 (数据值从小到大排序之后的两端数值),用剩余元素计算平均值获得滤波结果。选择一个窗口长度N,根据相邻信号之间的相关性,产生具有目标效果的滤波器。被去掉异类元素的个数取决于参数α(α表示被修整数据样本的百分比),并且0≤α≤0.5。假设窗口长度为N的数据集合{X(j)}被重新排序,即{X(1)≤X(2)≤…≤X(i)≤…≤X(N-1)≤X(N)}。α-trimmed均值滤波器的输出如下式所示:The α-trimmed mean value filter adopted in the embodiment of the present invention has a dynamic filter characteristic, and the value of any element of the signal depends on the value around the element. Select a window of a certain size to remove heterogeneous elements (values at both ends after the data values are sorted from small to large), and use the remaining elements to calculate the average value to obtain the filtering result. Choose a window length N, according to the correlation between adjacent signals, produce a filter with the target effect. The number of removed heterogeneous elements depends on the parameter α (α represents the percentage of trimmed data samples), and 0≤α≤0.5. Assume that the data set {X(j)} with a window length of N is reordered, that is, {X(1)≤X(2)≤...≤X(i)≤...≤X(N-1)≤X(N) }. The output of the α-trimmed mean filter is shown in the following formula:

其中,[αN]表示下取整函数。当α=0时,Xα需使用窗口内全部的样本值进行平均计算获得,Xα表现为均值滤波器;当α=0.5时,、Xα等于窗口重新排序后的中间样本值(此时N为奇数),Xα表现为中值滤波器。若N 为偶数,、并且当α=0.5时,中值Xα等于中间两个样本的平均值。Among them, [αN] represents the lower integer function. When α=0, Xα needs to use all the sample values in the window for average calculation, and Xα behaves as a mean filter; when α=0.5, Xα is equal to the intermediate sample value after the window is reordered (N is an odd number at this time ), Xα behaves as a median filter. If N is an even number, and when α=0.5, the median Xα is equal to the average value of the middle two samples.

以窗口长度N取奇数为例对α-trimmed均值滤波进行说明。假设选取的滤波器窗口长度N为5,α=0.2,那么数据中有两个元素(X1、X5)被剔除。由于滤波器剔除元素方式的对称性,被剔除个数必须为偶数,且小于窗口长度 N。在实际操作中,窗口值N>1,所以在使用α-trimmed均值滤波器时,信号序列开始和结束的几个元素都无法进行符合窗口大小元素的选取。本发明采取的办法是:在运用滤波运算之前,在信号头尾添加长度约为窗口大小一半的附加信号,且以开头结尾为轴元素对称添加。Taking an odd number of window length N as an example to illustrate α-trimmed mean filtering. Assuming that the selected filter window length N is 5, and α=0.2, then two elements (X1, X5) in the data are eliminated. Due to the symmetry of the way the filter removes elements, the number of removed elements must be an even number and less than the window length N. In actual operation, the window value N>1, so when using the α-trimmed mean filter, several elements at the beginning and end of the signal sequence cannot be selected to meet the window size elements. The method adopted by the present invention is: before using the filtering operation, add an additional signal whose length is about half of the window size at the beginning and end of the signal, and add symmetrically with the beginning and the end as axis elements.

在一维信号的情况下,如图5所示,当窗口长度N为偶数时,为保证所有原数据都被执行一次滤波操作,原数据头尾各增加N/2个元素;当窗口长度N 为奇数时,原数据头尾各增加(N-1)/2个元素。In the case of a one-dimensional signal, as shown in Figure 5, when the window length N is an even number, in order to ensure that all original data are filtered once, N/2 elements are added to the head and tail of the original data; when the window length N When it is an odd number, (N-1)/2 elements are added to the head and tail of the original data.

在二维信号的情况下,信号的处理方式与一维信号基本相同,只是α -trimmed均值滤波器的窗口选取变为二维,通常为3×3的矩阵。如图6所示,二维数据的延展采取方式是在原始数据的周围增加一些轴线,然后用一维方式增加对称元素。In the case of two-dimensional signals, the signal processing method is basically the same as that of one-dimensional signals, except that the window selection of the α-trimmed mean filter becomes two-dimensional, usually a 3×3 matrix. As shown in Figure 6, the way to extend the two-dimensional data is to add some axes around the original data, and then add symmetrical elements in a one-dimensional way.

无论是MT方式还是CSAMT方式的频率域电磁法采集数据,其主要以窄脉冲干扰为主,因此需适当增加窗口长度,优选的,N=10。根据数据特点,选取适当的窗口长度N以及滤波参数α,数据滤波效果更好。No matter the frequency domain electromagnetic method of MT mode or CSAMT mode collects data, it is mainly narrow pulse interference, so it is necessary to increase the window length appropriately, preferably, N=10. According to the characteristics of the data, the appropriate window length N and the filtering parameter α are selected, and the data filtering effect is better.

本发明实施例的α-Trimed可以采用FPGA实现,本发明采用STM32嵌入式MCU芯片与XILINX Spartan6FPGA结合的方式作为频率域电磁法接收机的核心主控板处理单元。基于Verilog语言将α-trimmed滤波方法移植到FPGA 芯片程序当中。图7为本发明实施例的α-trimmed滤波器滤波的流程示意图,所述流程包括:The α-Trimed in the embodiment of the present invention can be realized by FPGA, and the present invention uses the combination of STM32 embedded MCU chip and XILINX Spartan6FPGA as the core main control board processing unit of the frequency domain electromagnetic method receiver. Based on the Verilog language, the α-trimmed filtering method is transplanted into the FPGA chip program. FIG. 7 is a schematic flow diagram of the α-trimmed filter filtering according to an embodiment of the present invention, and the flow includes:

步骤710,输入原始数据;Step 710, input raw data;

步骤720,延展数据;Step 720, extend the data;

步骤730,设置窗口长度;Step 730, setting the window length;

步骤740,数据排序;Step 740, data sorting;

步骤750,剔除异常数据;Step 750, removing abnormal data;

步骤760,剩余数据求平均;Step 760, averaging the remaining data;

步骤770,生成新的数据;Step 770, generate new data;

步骤780,判断数据处理是否完成,若是,转至步骤720,若否,则转至步骤790;Step 780, judge whether the data processing is completed, if so, go to step 720, if not, then go to step 790;

步骤790,滤波完成。Step 790, the filtering is completed.

进一步的,本发明所述用于多通道接收器的抗干扰方法在进行步骤S110 对多通道接收器接收的原始信号进行同步叠加处理之前,还包括:Further, the anti-jamming method for the multi-channel receiver according to the present invention also includes:

采集原始信号并通过读出电路,将其转换为电压信号;Collect the original signal and convert it into a voltage signal through the readout circuit;

通过模数转换器将所述电压信号转换为数字信号;converting the voltage signal into a digital signal by an analog-to-digital converter;

所述数字信号通过数据总线传输至FPGA数据缓存区;The digital signal is transmitted to the FPGA data buffer area through the data bus;

然后从FPGA数据缓存区提取信号进行叠加滤波处理,Then extract the signal from the FPGA data buffer area for superimposed filtering processing,

为了降低量化噪声,优选的,所述模数转换器为高速24位模数转换器。In order to reduce quantization noise, preferably, the analog-to-digital converter is a high-speed 24-bit analog-to-digital converter.

进一步的,本发明所述用于多通道接收器的抗干扰方法在进行步骤S130 对所述经梳状滤波处理的信号进行α-Trimed滤波处理之后,还包括:Further, the anti-interference method for a multi-channel receiver according to the present invention, after performing step S130 to perform α-Trimed filtering on the comb-filtered signal, further includes:

对所述经过α-Trimed滤波处理之后的数据进行DFT运算;Carrying out DFT operation on the data after the α-Trimed filter processing;

将数据保存至存储设备并同时传送至监控设备;Save the data to the storage device and transmit it to the monitoring device at the same time;

通过监控设备实时监控数据质量。Monitor data quality in real time through monitoring equipment.

通过读数据的实时存储和监控,当发现问题时,及时调取该段时间的相关数据进行分析和排除,能够进一步排除干扰,增加数据质量,提升多通道接收器的性能。Through the real-time storage and monitoring of read data, when a problem is found, the relevant data of this period can be retrieved in time for analysis and elimination, which can further eliminate interference, increase data quality, and improve the performance of multi-channel receivers.

本发明实施例还提供了一种用于多通道接收器的抗干扰系统,所述系统包括:An embodiment of the present invention also provides an anti-jamming system for a multi-channel receiver, the system comprising:

同步叠加滤波器,对多通道接收器接收的原始信号进行同步叠加处理;Synchronous superposition filter, which performs synchronous superposition processing on the original signal received by the multi-channel receiver;

梳状滤波器,对所述经同步叠加处理过后的信号进行梳状滤波处理;Comb filter, performing comb filter processing on the signal after synchronous superposition processing;

α-Trimed滤波器,对所述经梳状滤波处理的信号进行α-Trimed滤波处理。an α-Trimed filter, performing α-Trimed filter processing on the signal processed by the comb filter.

本发明实施例所述的多通道接收器的抗干扰方法和系统,通过引入三种滤波技术,能滤除信号中的多种干扰,实现了对宽频带弱信号采集与处理术,提高电法勘探方法及电法勘探仪器的精确度。The anti-jamming method and system of the multi-channel receiver described in the embodiment of the present invention can filter out various kinds of interference in the signal by introducing three kinds of filtering techniques, realize the acquisition and processing of wide-band weak signals, and improve the electrical method. Accuracy of prospecting methods and electrical prospecting instruments.

尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。While preferred embodiments of the invention have been described, additional changes and modifications to these embodiments can be made by those skilled in the art once the basic inventive concept is appreciated. Therefore, it is intended that the appended claims be construed to cover the preferred embodiment as well as all changes and modifications which fall within the scope of the invention.

显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalent technologies, the present invention also intends to include these modifications and variations.

Claims (10)

1. a kind of anti-interference method for multi-channel receiver, which is characterized in that the described method includes:
Superposition processing is synchronized to multi-channel receiver received signal;
Comb filtering processing is carried out to the signal after the synchronized superposition processing;
α-Trimed filtering processing is carried out to the signal after comb filtering is processed.
2. the method as described in claim 1, which is characterized in that it is described multi-channel receiver received signal is synchronized it is folded Add processing, comprising:
Multi-channel receiver received signal is divided into multiple signal segments;
Signal in each signal segment is divided into multiple signal sections;
Multiple signal sections in each signal segment are overlapped, each respective amplitude of signal segment is obtained and phase is believed Breath;
Multiple superimposed signal segments are overlapped, the amplitude and phase information of the multichannel received signal are obtained.
3. method as described in claim 1, which is characterized in that carry out pectination filter to the signal after the synchronized superposition processing Wave processing, comprising:
It is delayed to multi-channel receiver received signal, becomes time delayed signal;
The multi-channel receiver received signal is superimposed with the time delayed signal, to generate phase cancellation.
4. method as claimed in claim 3, which is characterized in that carry out pectination to the signal after the synchronized superposition processing Filtering processing, comprising:
It is realized using the field programmable gate array (FPGA) with adder and multiplier to the letter after the synchronized superposition processing Number carry out comb filtering processing.
5. method as claimed in claim 4, which is characterized in that carry out pectination to the signal after the synchronized superposition processing Filtering processing, comprising:
Pectination is carried out to the signal after the synchronized superposition processing to realize using multiple FPGA cascades with adder and multiplier Filtering processing.
6. the method as described in claim 1, which is characterized in that carry out α-Trimed to the signal handled through comb filtering Filtering processing, comprising:
Length of window N and parameter alpha are set, and wherein α indicates the percentage for being trimmed data sample;
According to length of window N and parameter alpha, the signal handled through comb filtering is handled;
It is handled by the mobile element of window, then to the signal handled through comb filtering.
7. method as claimed in claim 6, which is characterized in that it is described according to length of window N and parameter alpha, to described through pectination The signal of filtering processing is handled, comprising:
According to length of window N and parameter alpha, redundant data is filtered out;
Average treatment is done to remaining data.
8. the method for claim 7, which is characterized in that carry out α-Trimed to the signal handled through comb filtering Before filtering processing, further includes:
Obtain the corresponding numerical value of multi-channel receiver received signal;
Additional signal is added respectively in signal heads and signal tail, and the additional signal added is symmetrical addition;
To the corresponding numerical value of the multi-channel receiver received signal according to sorting from small to large.
9. method according to claim 8, which is characterized in that additional signal is added respectively in the signal heads and signal tail, Include:
If the length of window N is even number, N/2 data are respectively added in signal heads and signal tail;
If the length of window N is odd number, (N-1)/2 data are respectively added in signal heads and signal tail.
10. a kind of jamproof system for multi-channel receiver, which is characterized in that the system comprises:
Synchronous superposition filter, for synchronizing superposition processing to multi-channel receiver received signal;
Comb filter, for carrying out comb filtering processing to the signal after the synchronized superposition processing;
α-Trimed filter, for carrying out α-Trimed filtering processing to the signal handled through comb filtering.
CN201810623258.5A 2018-06-15 2018-06-15 Anti-interference method and system for multi-channel receiver Pending CN109031438A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810623258.5A CN109031438A (en) 2018-06-15 2018-06-15 Anti-interference method and system for multi-channel receiver

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810623258.5A CN109031438A (en) 2018-06-15 2018-06-15 Anti-interference method and system for multi-channel receiver

Publications (1)

Publication Number Publication Date
CN109031438A true CN109031438A (en) 2018-12-18

Family

ID=64609938

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810623258.5A Pending CN109031438A (en) 2018-06-15 2018-06-15 Anti-interference method and system for multi-channel receiver

Country Status (1)

Country Link
CN (1) CN109031438A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110542927A (en) * 2019-09-02 2019-12-06 吉林大学 Spike Noise Suppression Method of Variable Window Weighted Seismic Data
CN113687432A (en) * 2021-08-19 2021-11-23 中国地质科学院地球物理地球化学勘查研究所 Method for removing power frequency interference in time domain aviation electromagnetic data
CN114743533A (en) * 2022-03-09 2022-07-12 中科上声(苏州)电子有限公司 Vehicle noise reduction method and device for broadband noise and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100235166A1 (en) * 2006-10-19 2010-09-16 Sony Computer Entertainment Europe Limited Apparatus and method for transforming audio characteristics of an audio recording
CN103308951A (en) * 2013-06-03 2013-09-18 中国地质大学(武汉) Control system and operation method of transient electromagnetic signal acquisition instrument of handheld pseudo random
CN107329180A (en) * 2017-06-22 2017-11-07 吉林大学 Magnetic resonance underground water detection device and detection method based on simulation comb filter

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100235166A1 (en) * 2006-10-19 2010-09-16 Sony Computer Entertainment Europe Limited Apparatus and method for transforming audio characteristics of an audio recording
CN103308951A (en) * 2013-06-03 2013-09-18 中国地质大学(武汉) Control system and operation method of transient electromagnetic signal acquisition instrument of handheld pseudo random
CN107329180A (en) * 2017-06-22 2017-11-07 吉林大学 Magnetic resonance underground water detection device and detection method based on simulation comb filter

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
温金鑫 等: "Σ-ΔA/D在电机控制系统中的应用研究", 《电气传动》 *
王中兴 等: "地面核磁共振找水信号中的奇异干扰抑制", 《吉林大学学报(工学版)》 *
王文闯等: "基于α-trimmed矢量中值滤波压制同步激发邻炮干扰 ", 《石油地球物理勘探》 *
王文闯等: "基于α-trimmed矢量中值滤波压制同步激发邻炮干扰", 《石油地球物理勘探》 *
秦雪霏: "地震波多次散射波场属性分析", 《吉林大学硕士学位论文》 *
赵晓鸣等: "CSAMT资料处理中的静校正问题 ", 《工程地球物理学报》 *
赵晓鸣等: "CSAMT资料处理中的静校正问题", 《工程地球物理学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110542927A (en) * 2019-09-02 2019-12-06 吉林大学 Spike Noise Suppression Method of Variable Window Weighted Seismic Data
CN113687432A (en) * 2021-08-19 2021-11-23 中国地质科学院地球物理地球化学勘查研究所 Method for removing power frequency interference in time domain aviation electromagnetic data
CN114743533A (en) * 2022-03-09 2022-07-12 中科上声(苏州)电子有限公司 Vehicle noise reduction method and device for broadband noise and storage medium
CN114743533B (en) * 2022-03-09 2025-04-18 中科上声(苏州)电子有限公司 Vehicle noise reduction method, device and storage medium for broadband noise

Similar Documents

Publication Publication Date Title
CN110609204B (en) Single-phase-to-ground fault location method for distribution network based on morphological wavelet analysis
CN107395157B (en) Ground net potential difference filtering method based on wavelet transformation and weighted moving average
CN109031438A (en) Anti-interference method and system for multi-channel receiver
CN107976572A (en) A kind of channelizing high-precision real Time-frequency Spectrum Analysis method and its system
CN102508013A (en) Current alarm signal processing method for improved mathematical morphology filter
CN105510711A (en) Empirical mode decomposition-based improved harmonic analysis method
CN104901754B (en) Channel monitoring system based on channelization frequency spectrum perception
WO2024109031A1 (en) Hht-based voltage quality disturbance detection method
CN104111373B (en) A measurement method for digital electric energy in smart substations
CN111046791A (en) Current signal filtering and denoising method based on generalized S transform containing variable factors
CN110287853A (en) A Transient Signal Denoising Method Based on Wavelet Decomposition
CN108090270B (en) A Transient Oscillation Parameter Identification Method Based on Morphological Filtering and Blind Source Separation
CN113221615B (en) A partial discharge pulse extraction method based on noise reduction clustering
CN108776259B (en) Power analyzer based on digital bandwidth limiting technology
CN107104674A (en) Noise power estimation apparatus in electric parameter metering process
CN102721870A (en) SVD (Singular Value Decomposition) digital filter-based electric power system harmonic analysis method
CN117526980A (en) Method and system for processing signals of Roland C receiver based on interpolation FIR band-pass filtering and self-adaptive notch algorithm
CN103472376B (en) Partial discharge of transformer superfrequency positioning analysis device and method for positioning analyzing thereof
Liu et al. Clockless continuous-time neural spike sorting: Method, implementation and evaluation
CN103944572B (en) Ultra-high speed sampling rate harvester and method
CN106197480A (en) A kind of processing system of Low SNR signal
CN104836547A (en) Short group time-delay digit filtering method
Zou et al. Mathematical morphology based phase selection scheme in digital relaying
Li Design and realization of FIR digital filters based on MATLAB
CN203433082U (en) Ultrahigh frequency (UHF) positioning analyzer for local discharge of transformer

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20181218