[go: up one dir, main page]

CN115189676B - A data fusion method for distributed energy supply system - Google Patents

A data fusion method for distributed energy supply system Download PDF

Info

Publication number
CN115189676B
CN115189676B CN202210885388.2A CN202210885388A CN115189676B CN 115189676 B CN115189676 B CN 115189676B CN 202210885388 A CN202210885388 A CN 202210885388A CN 115189676 B CN115189676 B CN 115189676B
Authority
CN
China
Prior art keywords
prototype
data fusion
filter
notch
initial value
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
CN202210885388.2A
Other languages
Chinese (zh)
Other versions
CN115189676A (en
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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN202210885388.2A priority Critical patent/CN115189676B/en
Publication of CN115189676A publication Critical patent/CN115189676A/en
Application granted granted Critical
Publication of CN115189676B publication Critical patent/CN115189676B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00—Adaptive networks
    • H03H21/0012—Digital adaptive filters
    • H03H21/0043—Adaptive algorithms
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40—Arrangements for reducing harmonics

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)
  • Filters And Equalizers (AREA)

Abstract

The invention relates to the technical field of power line interference removal, in particular to a data fusion method for a distributed energy supply system, which comprises the following steps of firstly, carrying out recruitment and measurement on designated registers of each protocol device at regular time through a network port of a data fusion terminal; step two, using a high-selectivity elliptic function self-adaptive notch filter arranged in the processor to inhibit power line interference from short-delay evoked potentials through the width of a time-varying stop band and a non-zero initial value, and filtering all input sampling data; and thirdly, fusing the filtered data into corresponding messages according to preset requirements through a self-adaptive Ethernet interface of the processor, and uploading the corresponding messages to a background system. The method can simultaneously give consideration to the optimization problem of the duration and the amplitude of transient response.

Description

Data fusion method for distributed energy supply system
Technical Field
The invention relates to the technical field of power line interference removal, in particular to a data fusion method for a distributed energy supply system.
Background
Power line (electric field line) interference is one of the common interferences affecting data fusion terminals for distributed energy supply systems, and its effective removal remains a challenge in distributed energy supply systems. Power line interference is generally characterized by a fixed frequency sine wave (50/60 Hz, depending on country/region) with random phase and amplitude and higher harmonics.
Many signal processing methods and algorithms have been proposed to solve the power line interference removal problem, with notch filtering being the most common method of removing power line noise, and a series of multiple notch filters (fixed notch filters) centered around the power line fundamental frequency and higher harmonics. The multi-notch filter may be implemented as an FIR or IIR structure, and the order of the IIR filter may be designed to be much smaller than an FIR filter that satisfies the same frequency assumption. The cascading method is a direct design method of the IIR multi-notch filter, the multi-notch filter is designed by cascading a required number of second-order single-notch filters, and other methods are usually based on optimization of pole positions, and advanced iterative algorithm or complex searching technology is needed. The multi-notch filter is capable of suppressing selected frequencies in the input signal and if the notch (selected frequency) is sufficiently narrow, attenuation of other frequencies can be minimized, the bandwidth of the notch is directly dependent on the pole radius of the filter, which becomes more selective as the pole radius approaches a unit circle, however, shrinking the notch of the filter results in an increase in the transient response duration of the filter. Thus, the interference rejection filter should on the one hand be very selective and narrowband, while it should have good time domain properties, choosing between a very sharp amplitude response, a narrow notch and a short transient response duration is one of the important issues in designing an optimal multi-notch filter.
In summary, how to simultaneously consider the optimization problem of the duration and the amplitude of transient response becomes the problem to be solved in the data fusion of the power line to the distributed energy supply system at present.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a data fusion method for a distributed energy supply system, which can simultaneously solve the problem of optimizing the duration and the amplitude of transient response.
In order to solve the technical problems, the invention adopts the following technical scheme:
A data fusion method for a distributed energy supply system comprises the following steps:
step one, through a data fusion terminal network port, carrying out recruitment and measurement on each protocol equipment designated register at regular time;
step two, using a high-selectivity elliptic function self-adaptive notch filter arranged in the processor to inhibit power line interference from short-delay evoked potentials through the width of a time-varying stop band and a non-zero initial value, and filtering all input sampling data;
step three, fusing the filtered data into corresponding messages according to preset requirements through a self-adaptive Ethernet interface of the processor, and uploading the corresponding messages to a background system;
In the second step, the design process of the high-selectivity elliptic function adaptive notch filter includes:
S1, synthesizing a plurality of prototype notch filters into prototype multi-notch filters, wherein the prototype notch filters are elliptical filters;
S2, introducing a non-zero initial value to the prototype multi-notch filter, and inhibiting transient response of the prototype multi-notch filter;
And S3, introducing a time-varying stop band width into the prototype multi-notch filter to obtain the high-selectivity elliptic function self-adaptive notch filter.
Preferably, in S2, the calculating process of the non-zero initial value includes:
S201, calculating a non-zero initial value vector
Wherein, I is an identity matrix, P is a projection operator, and X is a vector of input samples;
S202, analyzing the filtering quality and initial value vector measured by the MSE The correlation between the lengths k of the filter is obtained, and the k value which enables the initial value of the selected filter to have the optimal interference suppression effect is obtained;
S203, according to the optimal k value and initial value vector A non-zero initial value is obtained.
Preferably, in S201, the calculation formula of the projection operator P is:
P=B(BTB)-1BT;
Where ω represents the power supply disturbance frequency and h is the number of highest harmonics.
Preferably, the calculation formula of the vector X of the input samples is:
X=U+D=[x(0)x(1)...x(k-1)]T;
x(n)=u(n)+d(n);
Where x (n) is the input signal, u (n) is the ideal signal, d (n) is the sinusoidal disturbance of known frequency, i is the index of the harmonic, h is the number of the highest harmonics, A i is the amplitude of the harmonic, ω is the power supply disturbance frequency, and φ i is the phase of the harmonic.
Preferably, in S202, the mean square error MSE is calculated as:
where N is the window size and y (N) is the filtered output signal.
Preferably, in S3, after the time-varying stop-band width is introduced, the description formula of the adaptive notch filter with high selectivity elliptic function is:
Where y (n) is the signal after filtering, l is the length of the variation range in the sample, and b 0m(n)、…、bNm (n) and a 1m(n)、…、aNm (n) are coefficients in the corresponding time domain equation for the prototype multi-notch filter transfer function.
Preferably, in S1, each of the plurality of prototype notch filters is centered on a fundamental frequency of the power line at 60Hz or higher.
Preferably, in S1, the number of prototype notch filters is five, and the center frequencies of the five prototype notch filters are ω 1＝60Hz、ω2＝180Hz、ω3＝540Hz、ω4＝900Hz、ω5 =1260 Hz.
Preferably, the maximum passband ripple, the minimum stopband attenuation, and the stopband width of the five prototype notch filters are all the same, and the maximum passband ripple rp=0.1 dB, the minimum stopband attenuation rs=40 dB, and the stopband width bw=4 Hz.
Preferably, the transfer function of the prototype multi-notch filter is:
Wherein, H 1(z)、H2(z)、H3(z)、H4(z)、H5 (z) is the transfer function of the five prototype notch filters, z -1、z-2、…、z-30 is the unit impulse response sequence of the prototype notch filter, and a 1、a2、…、a30 and b 0、b1、b2、…、b30 are the coefficients before z -1、z-2、…、z-30 are calculated after multiplying H 1(z)、H2(z)、H3(z)、H4(z)、H5 (z).
Compared with the prior art, the invention has the following beneficial effects:
1. In the invention, the acquired data of the protocol equipment is filtered by the high-selectivity elliptic function self-adaptive notch filter in the processor and then is uploaded.
The high-selectivity elliptic function self-adaptive notch filter is composed of a plurality of elliptic filters, the elliptic filters provide a sharp transition zone, so that narrow-band power line noise in an evoked potential is effectively eliminated, the noise is separated from useful signals, the elliptic filters have high-quality factors, sharp transition zones and enough selectivity, the narrow-band power line noise in the evoked potential can be effectively eliminated, and compared with a Butterworth filter or a Chebyshev filter, the elliptic filters can provide better suppression. The high-selectivity elliptic function adaptive notch filter formed by a plurality of elliptic filters, namely the digital IIR multi-notch filter based on the elliptic prototype analog filter, has the transient response duration and amplitude which are obviously lower than those of the traditional filter, and can eliminate the power line interference more effectively. Besides, in the design process of the high-selectivity elliptic function self-adaptive notch filter, a non-zero initial value is introduced, and compared with a filter adopting the zero initial value, the power line interference can be further eliminated, and a time-varying stop band width is introduced into a structure corresponding to the second stage of the transient suppression process, and as the transient response duration of the filter depends on the bandwidth of the notch, the transient response duration of the filter can be increased by reducing the notch of the filter, and the transient response duration of the filter can be effectively reduced after the time-varying stop band width is introduced.
In conclusion, the transient response duration and amplitude of the invention are obviously lower than those of the traditional filter, and the optimization problem of the transient response duration and amplitude can be simultaneously considered.
2. The invention provides a design process of a high-selectivity elliptic function self-adaptive notch filter, and the corresponding non-zero initial value can be accurately, rapidly and boundaiy obtained through the calculation process of the non-zero initial value provided by the invention.
3. The invention also provides a selection mode of the superior prototype notch filter, which comprises the number of the selected prototype notch filters, the center frequency, the maximum passband ripple, the minimum stopband attenuation and the stopband width, and the prototype notch filter reference provided by the invention, by combining the processing mode of the non-zero initial value and the time-varying stop band width, the high-selectivity elliptic function self-adaptive notch filter with the transient response duration and amplitude being obviously lower than those of the traditional filter can be obtained.
Drawings
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings, in which:
FIG. 1 is a block diagram of a data fusion terminal for a distributed energy supply system in an embodiment;
FIG. 2 is a flow chart of a data fusion method for a distributed energy supply system in an embodiment;
FIG. 3 is a flow chart of a design process of a high selectivity elliptic function adaptive notch filter according to an embodiment;
FIG. 4 is a schematic diagram showing the amplitude and time response of a raw multi-notch filter according to an embodiment;
FIG. 5 is a schematic diagram showing the coefficient values of some of the embodiments;
FIG. 6 is a mean square error and initial value vector in the embodiment A correlation diagram between lengths k;
Fig. 7 is a graph showing the time response of the 4Hz and 30Hz bandwidths to a 60Hz sinusoidal input signal in an exemplary embodiment of a raw multi-notch filter.
Detailed Description
The following is a further detailed description of the embodiments:
examples:
It should be noted that, for ease of understanding, a block diagram of the data fusion terminal for a distributed energy supply system in this embodiment is shown in fig. 1, and the data fusion terminal for a distributed energy supply system includes a Cortex-a7 processor, 4 sets of RS485 ports, 20 MODBUS protocol devices, a 10/100M adaptive ethernet interface, an ethernet and a power supply.
The specific application scene in the embodiment is that a background system regularly carries out call testing on 20 MODBUS protocol equipment designated registers accessed to 4 groups of RS485 ports through a data fusion terminal network port, a Cortex-a7 processor filters all input sampling data, and finally, the acquired data are fused into specific messages as required through a 10/100M self-adaptive Ethernet interface of the processor, and then the specific messages are actively and regularly uploaded to the background system. In specific implementation, the main frequency of the Cortex-a7 processor is 1.2GHz, each path of RS485 independent thread and each network independent thread, each group of RS485 interface is independently provided with communication parameters, the baud rate supports 1200-115200 bp/s, the shortest time interval between the two detection is 1 minute, 4 groups of RS485 buses are independently provided with automatic acquisition registers, each register is independently provided with a data conversion mode (fixed point to floating point and multi-byte arrangement sequence), and the processor supports network port upgrading.
As shown in fig. 2, the embodiment discloses a data fusion method for a distributed energy supply system, which includes the following steps:
step one, through a data fusion terminal network port, carrying out recruitment and measurement on each protocol equipment designated register at regular time;
step two, using a high-selectivity elliptic function self-adaptive notch filter arranged in the processor to inhibit power line interference from short-delay evoked potentials through the width of a time-varying stop band and a non-zero initial value, and filtering all input sampling data;
step three, fusing the filtered data into corresponding messages according to preset requirements through a self-adaptive Ethernet interface of the processor, and uploading the corresponding messages to a background system;
In the second step, as shown in fig. 3, the design process of the adaptive notch filter for a high-selectivity elliptic function includes:
S1, synthesizing a plurality of prototype notch filters into a prototype multi-notch filter, wherein the prototype notch filters are elliptical filters.
As shown in fig. 4, in the present embodiment, the number of prototype notch filters is five, the center frequencies of the five prototype notch filters are ω 1＝60Hz、ω2＝180Hz、ω3＝540Hz、ω4＝900Hz、ω5 =1260 Hz, the maximum passband ripple, the minimum stopband attenuation and the stopband width of the five prototype notch filters are the same, the maximum passband ripple rp=0.1 dB, the minimum stopband attenuation rs=40 dB, and the stopband width bw=4 Hz.
The transfer function of the prototype multi-notch filter is:
Wherein, H 1(z)、H2(z)、H3(z)、H4(z)、H5 (z) is the transfer function of the five prototype notch filters, z -1、z-2、…、z-30 is the unit impulse response sequence of the prototype notch filter, and a 1、a2、…、a30 and b 0、b1、b2、…、b30 are the coefficients before z -1、z-2、…、z-30 are calculated after multiplying H 1(z)、H2(z)、H3(z)、H4(z)、H5 (z).
In this example, values of a1, a2, a30, b0, b1, b2, and b30 are shown in fig. 5.
S2, introducing a non-zero initial value to the prototype multi-notch filter, and inhibiting transient response of the prototype multi-notch filter.
The calculation process of the non-zero initial value comprises the following steps:
S201, calculating a non-zero initial value vector
Wherein, I is an identity matrix, P is a projection operator, and X is a vector of input samples;
the calculation formula of the projection operator P is:
P=B(BTB)-1BT;
Where ω represents the power supply disturbance frequency and h is the number of highest harmonics.
The calculation formula of the vector X of the input samples is:
X=U+D=[x(0)x(1)...x(k-1)]T;
x(n)=u(n)+d(n);
Where x (n) is the input signal, u (n) is the ideal signal, d (n) is the sinusoidal disturbance of known frequency, i is the index of the harmonic, h is the number of the highest harmonics, A i is the amplitude of the harmonic, ω is the power supply disturbance frequency, and φ i is the phase of the harmonic.
S202, analyzing the filtering quality and initial value vector measured by the MSEThe correlation between the lengths k of the filter is obtained, and the k value which enables the initial value of the selected filter to have the optimal interference suppression effect is obtained;
Wherein, the calculation formula of the MSE is:
where N is the window size and y (N) is the filtered output signal.
In this embodiment, as shown in FIG. 6, when k is an integer multiple of f s/f0, the MSE reaches a local minimum, and in summary, the local minimum k min of the function MSE (k) can be determined by the following relationship:
Where f s is the sampling rate, f 0 is the power supply interference frequency, j is an integer, an initial condition vector The choice of the length k is a compromise between the filtering quality and the delay due to the need to record k samples in the signal x (n).
S203, according to the optimal k value and initial value vectorA non-zero initial value is obtained.
And S3, introducing a time-varying stop band width into the prototype multi-notch filter to obtain the high-selectivity elliptic function self-adaptive notch filter.
As shown in fig. 7, the transient response duration of the filter depends on the bandwidth of the notch, and shrinking the notch increases the transient response duration of the filter, thus introducing a time-varying stop-band width into the filter structure during the second phase of the considered multi-notch filter transient suppression process.
After the time-varying stop band width is introduced, the description formula of the high-selectivity elliptic function self-adaptive notch filter is as follows:
Where y (n) is the signal after filtering, l is the length of the variation range in the sample, and b 0m(n)、…、bNm (n) and a 1m(n)、…、aNm (n) are coefficients in the corresponding time domain equation for the prototype multi-notch filter transfer function.
The linear function of the filter stopband width length is the simplest form of stopband width variation, in this embodiment, to ensure that the filter is operating at a sufficient speed at its beginning, the initial value of stopband width BW start = 30Hz, and the final stopband width BW final = 4Hz is selectively determined by the established filter.
In the invention, the acquired data of the protocol equipment is filtered by the high-selectivity elliptic function self-adaptive notch filter in the processor and then is uploaded. The high-selectivity elliptic function self-adaptive notch filter is composed of a plurality of elliptic filters, the elliptic filters provide a sharp transition zone, so that narrow-band power line noise in an evoked potential is effectively eliminated, the noise is separated from useful signals, the elliptic filters have high-quality factors, sharp transition zones and enough selectivity, the narrow-band power line noise in the evoked potential can be effectively eliminated, and compared with a Butterworth filter or a Chebyshev filter, the elliptic filters can provide better suppression. The high-selectivity elliptic function adaptive notch filter formed by a plurality of elliptic filters, namely the digital IIR multi-notch filter based on the elliptic prototype analog filter, has the transient response duration and amplitude which are obviously lower than those of the traditional filter, and can eliminate the power line interference more effectively. Besides, in the design process of the high-selectivity elliptic function self-adaptive notch filter, a non-zero initial value is introduced, and compared with a filter adopting the zero initial value, the power line interference can be further eliminated, and a time-varying stop band width is introduced into a structure corresponding to the second stage of the transient suppression process, and as the transient response duration of the filter depends on the bandwidth of the notch, the transient response duration of the filter can be increased by reducing the notch of the filter, and the transient response duration of the filter can be effectively reduced after the time-varying stop band width is introduced.
In addition, in the invention, a design process of a high-selectivity elliptic function self-adaptive notch filter is provided, and the corresponding non-zero initial value can be accurately, rapidly and boundry obtained through the calculation process of the non-zero initial value provided by the invention. The invention also provides a selection mode of the superior prototype notch filter, which comprises the number of the selected prototype notch filters, the center frequency, the maximum passband ripple, the minimum stopband attenuation and the stopband width, and the prototype notch filter reference provided by the invention, by combining the processing mode of the non-zero initial value and the time-varying stop band width, the high-selectivity elliptic function self-adaptive notch filter with the transient response duration and amplitude being obviously lower than those of the traditional filter can be obtained.
In conclusion, the transient response duration and amplitude of the invention are obviously lower than those of the traditional filter, and the optimization problem of the transient response duration and amplitude can be simultaneously considered.
Finally, it should be noted that the above embodiments are only for illustrating the technical solution of the present invention and not for limiting the technical solution, and those skilled in the art should understand that modifications and equivalents may be made to the technical solution of the present invention without departing from the spirit and scope of the present invention, and all such modifications and equivalents are included in the scope of the claims.

Claims (8)

1.一种分布式供能系统用数据融合方法，其特征在于，包括以下步骤：1. A data fusion method for a distributed energy supply system, characterized in that it comprises the following steps: 步骤一、通过数据融合终端网口，定时对各协议设备指定寄存器进行招测；Step 1: Through the network port of the data fusion terminal, the designated registers of each protocol device are tested regularly; 步骤二、用处理器中设置的高选择性椭圆函数自适应陷波滤波器，通过时变阻带宽度和非零初始值抑制来自短延迟诱发电位的电力线干扰，对所有输入的采样数据进行滤波；Step 2: Using a highly selective elliptic function adaptive notch filter set in the processor to suppress power line interference from short-delay evoked potentials through a time-varying stopband width and a non-zero initial value, all input sampled data are filtered; 步骤三、通过处理器的自适应以太网接口，将滤波后的数据根据预设需求融合为对应报文后，上传后台系统；Step 3: Through the adaptive Ethernet interface of the processor, the filtered data is merged into corresponding messages according to preset requirements and uploaded to the background system; 其中，步骤二中，所述高选择性椭圆函数自适应陷波滤波器的设计过程包括：Among them, in step 2, the design process of the highly selective elliptic function adaptive notch filter includes: S1、将多个原型陷波滤波器合成为原型多陷波滤波器；所述多个原型陷波滤波器均为椭圆滤波器；S1, synthesizing a plurality of prototype notch filters into a prototype multi-notch filter; the plurality of prototype notch filters are all elliptic filters; S2、为原型多陷波滤波器引入非零初始值，抑制原型多陷波滤波器的瞬态响应；S2, introducing a non-zero initial value for the prototype multi-notch filter to suppress the transient response of the prototype multi-notch filter; S3、在原型多陷波滤波器中引入时变阻带宽度，得到高选择性椭圆函数自适应陷波滤波器；S3, introducing a time-varying stopband width into the prototype multi-notch filter to obtain a highly selective elliptic function adaptive notch filter; 其中，S2中，所述非零初始值的计算过程包括：Wherein, in S2, the calculation process of the non-zero initial value includes: S201、计算非零初始值向量 S201, calculate non-zero initial value vector 其中，I是单位矩阵、P为投影算子、X为输入样本的向量；Where I is the identity matrix, P is the projection operator, and X is the vector of input samples; S202、分析由均方误差MSE测量的滤波质量与初始值向量的长度k之间的相关性，得到令所选滤波器初始值抑制干扰效果最优的k值；S202, analyzing the filtering quality measured by the mean square error MSE and the initial value vector The correlation between the length k of the filter is used to obtain the k value that makes the selected filter initial value have the best interference suppression effect; S203、根据最优的k值及初始值向量得到非零初始值；S203, according to the optimal k value and initial value vector Get a non-zero initial value; S3中，引入时变阻带宽度后，所述高选择性椭圆函数自适应陷波滤波器的描述式为：In S3, after the time-varying stopband width is introduced, the description formula of the highly selective elliptic function adaptive notch filter is: 式中，x(n)为输入信号，y(n)是表示滤波后的信号，l是样本中变化范围的长度，b0m(n)、…、bNm(n)和a1m(n)、…、aNm(n)为所述原型多陷波滤波器传递函数对应时域方程中的系数。In the formula, x(n) is the input signal, y(n) is the filtered signal, l is the length of the variation range in the sample, b0m (n),…, bNm (n) and a1m (n),…, aNm (n) are the coefficients in the time domain equation corresponding to the transfer function of the prototype multi-notch filter. 2.如权利要求1所述的分布式供能系统用数据融合方法，其特征在于：S201中，投影算子P的计算式为:2. The data fusion method for a distributed energy supply system according to claim 1, characterized in that: in S201, the calculation formula of the projection operator P is: P＝B(BTB)-1BT；P = B (B T B) -1 B T ; 其中，ω表示电源干扰频率，h为最高谐波的数量。Where ω is the power disturbance frequency and h is the number of the highest harmonic. 3.如权利要求2所述的分布式供能系统用数据融合方法，其特征在于：输入样本的向量X的计算式为：3. The data fusion method for a distributed energy supply system according to claim 2, wherein the calculation formula of the vector X of the input sample is: X＝U+D＝[x(0)x(1)…x(k-1)]T；X=U+D=[x(0)x(1)...x(k-1)] T ; U＝[u(0)u(1)…u(k-1)]T U＝[u(0)u(1)…u(k-1)] T D＝[d(0)d(1)…d(k-1)]T；D = [d(0)d(1)…d(k-1)] T ; x(n)＝u(n)+d(n)；x(n)＝u(n)+d(n); 其中，x(n)为输入信号、u(n)为理想信号、d(n)为已知频率的正弦干扰、i为谐波的指数、h为最高谐波的数量、Ai为谐波的幅度、ω为电源干扰频率、φi为谐波的相位。Where x(n) is the input signal, u(n) is the ideal signal, d(n) is a sinusoidal disturbance of known frequency, i is the index of the harmonic, h is the number of the highest harmonic, Ai is the amplitude of the harmonic, ω is the power supply disturbance frequency, and φi is the phase of the harmonic. 4.如权利要求3所述的分布式供能系统用数据融合方法，其特征在于：S202中，均方误差MSE的计算式为：4. The data fusion method for a distributed energy supply system according to claim 3, characterized in that: in S202, the calculation formula of the mean square error MSE is: 式中，N为窗口大小，y(n)是滤波后的输出信号。Where N is the window size and y(n) is the output signal after filtering. 5.如权利要求4所述的分布式供能系统用数据融合方法，其特征在于：S1中，多个原型陷波滤波器均以60Hz或更高谐波的电力线基频为中心。5. The data fusion method for a distributed energy supply system as described in claim 4 is characterized in that: in S1, multiple prototype notch filters are centered on the power line fundamental frequency of 60 Hz or higher harmonics. 6.如权利要求5所述的分布式供能系统用数据融合方法，其特征在于：S1中，原型陷波滤波器的数量为五个；五个原型陷波滤波器的中心频率分别为：ω1＝60Hz、ω2＝180Hz、ω3＝540Hz、ω4＝900Hz、ω5＝1260Hz。6. The data fusion method for a distributed energy supply system as claimed in claim 5, characterized in that: in S1, the number of prototype notch filters is five; the center frequencies of the five prototype notch filters are respectively: ω 1 =60 Hz, ω 2 =180 Hz, ω 3 =540 Hz, ω 4 =900 Hz, ω 5 =1260 Hz. 7.如权利要求6所述的分布式供能系统用数据融合方法，其特征在于：五个原型陷波滤波器的最大通带纹波、最小阻带衰减和阻带宽度均相同，且最大通带纹波Rp＝0.1dB，最小阻带衰减Rs＝40dB，阻带宽度BW＝4Hz。7. The data fusion method for a distributed energy supply system as described in claim 6 is characterized in that the maximum passband ripple, minimum stopband attenuation and stopband width of the five prototype notch filters are the same, and the maximum passband ripple Rp = 0.1dB, the minimum stopband attenuation Rs = 40dB, and the stopband width BW = 4Hz. 8.如权利要求7所述的分布式供能系统用数据融合方法，其特征在于：所述原型多陷波滤波器的传递函数为：8. The data fusion method for a distributed energy supply system according to claim 7, wherein the transfer function of the prototype multi-notch filter is: 其中，H1(z)、H2(z)、H3(z)、H4(z)、H5(z)分别为五个原型陷波滤波器的传递函数，z-1、z-2、…、z-30为原型陷波滤波器的单位冲激响应序列，a1、a2、…、a30和b0、b1、b2、…、b30为H1(z)、H2(z)、H3(z)、H4(z)、H5(z)相乘后计算得到对应z-1、z-2、…、z-30前的系数。Among them, H1 (z), H2 (z), H3 (z), H4 (z), and H5 (z) are the transfer functions of the five prototype notch filters respectively, z -1 , z -2 , …, z -30 are the unit impulse response sequences of the prototype notch filters, a1 , a2 , … , a30 and b0, b1 , b2 , …, b30 are the coefficients before z -1 , z- 2 , … , z -30 calculated by multiplying H1(z), H2 (z), H3 (z), H4(z), and H5 (z).
CN202210885388.2A 2022-07-26 2022-07-26 A data fusion method for distributed energy supply system Active CN115189676B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210885388.2A CN115189676B (en) 2022-07-26 2022-07-26 A data fusion method for distributed energy supply system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210885388.2A CN115189676B (en) 2022-07-26 2022-07-26 A data fusion method for distributed energy supply system

Publications (2)

Publication Number Publication Date
CN115189676A CN115189676A (en) 2022-10-14
CN115189676B true CN115189676B (en) 2025-06-13

Family

ID=83521457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210885388.2A Active CN115189676B (en) 2022-07-26 2022-07-26 A data fusion method for distributed energy supply system

Country Status (1)

Country Link
CN (1) CN115189676B (en)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7492675B2 (en) * 2001-01-25 2009-02-17 Dphi Acquisitions, Inc. Digital servo system with calibrated notch filters
CN102279327B (en) * 2011-05-06 2014-06-18 重庆大学 On-line monitoring and state evaluation system for photovoltaic grid-connected power generation
CN102723601B (en) * 2012-06-19 2015-01-07 北京航空航天大学 Ultra-wide-band dual-notch paster antenna adopting wide-attenuation-band electromagnetic band gap structure
CN110828959B (en) * 2019-10-31 2020-11-10 西安交通大学 Deformed ellipsoid resonant cavity and dual-mode waveguide filter based on same and without tuning

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
配网用模块化多电平换流器损耗分析与降损控制策略研究;李奇芬;重庆大学硕士学位论文;20230601;1-98 *

Also Published As

Publication number Publication date
CN115189676A (en) 2022-10-14

Similar Documents

Publication Publication Date Title
CN108763720B (en) DDC implementation method with sampling rate capable of being adjusted down at will
CN111510110A (en) A parallel processing interpolation matched filtering method and filter
CN115189676B (en) A data fusion method for distributed energy supply system
Van Gerwen et al. A new type of digital filter for data transmission
Kumar et al. Design and comparison between IIR Butterworth and Chebyshev digital filters using MATLAB
US4689759A (en) Process and installation for the analysis and retrieval of a sampling and interpolation signal
Dhillon et al. Power line interference removal from electrocardiogram using a simplified lattice based adaptive IIR notch filter
CN105915193B (en) A kind of improvement generation method for multiphase filter
Kurosu et al. A technique to truncate IIR filter impulse response and its application to real-time implementation of linear-phase IIR filters
CN113037248B (en) A Fractional Delay Filter Design Method Based on Fractional Frequency Domain Optimization
Yunlong et al. An extreme simple method for digital FIR filter design
CN115882820A (en) Filter circuit and analog-to-digital converter
CN115308486A (en) Electric energy harmonic wave metering method and system
CN106953615A (en) A kind of FIR digital audio frequency dividers being adjusted flexibly
Razee et al. Comparative performance analysis of IIR and FIR filters for 5G networks
Rathore et al. Design and Implementation of Time Varying pole Radius IIR Notch filter with improved performance
CN110212889B (en) A digital signal sampling device and method
Chandel et al. Digital Filter Designing: Minimum-Order IIR Filter Design
CN118114026B (en) Fuel cell impedance spectrum test signal generating and processing system and method
CN2326965Y (en) Filter for generator falt discharge on-line monitoring
Wang et al. Discussion on the-3 dB rejection bandwidth of IIR notch filters
CN102082559B (en) Method for realizing linear phase IIR (infinite impulse response) filter
Niyama et al. Implementation of IIR and FIR filters in Simulink MATLAB and its application in ECG
CN119652710B (en) Baseline drift correction circuit, baseline drift correction method and communication system
Okoniewski et al. Allpass based multi-notch IIR filter with equalized group delay

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