CN110046435B - Integrated digital filtering harmonic wave display method and system - Google Patents
Integrated digital filtering harmonic wave display method and system Download PDFInfo
- Publication number
- CN110046435B CN110046435B CN201910318263.XA CN201910318263A CN110046435B CN 110046435 B CN110046435 B CN 110046435B CN 201910318263 A CN201910318263 A CN 201910318263A CN 110046435 B CN110046435 B CN 110046435B
- Authority
- CN
- China
- Prior art keywords
- filtering
- module
- cursor
- data
- waveform data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
Abstract
The invention discloses an integrated digital filtering harmonic wave display method and system, wherein the method comprises the following steps: sampling the original waveform data to obtain sampling data; performing Fourier transform on the sampling data to obtain harmonic data of the original waveform data; displaying the harmonic data as a bar graph in a display area of the user interface; respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input area in the user interface; determining a filtering algorithm according to the mapping relation between the filtering parameters and a filtering algorithm classification table, and carrying out convolution operation on the original waveform data and the harmonic data according to the filtering algorithm to obtain filtered waveform data; the waveform data is displayed in a viewing area of the user interface after the raw waveform data is filtered. The scheme solves the problems of low loading efficiency and inconvenient use, improves the loading efficiency, improves the experience effect and is more convenient to use.
Description
Technical Field
The invention relates to the technical field of harmonic filtering, in particular to an integrated digital harmonic filtering display method and system.
Background
From the french mathematical fourier findings, any periodic function can be represented by an infinite series of sine and cosine functions. In many cases, the signal we are exposed to is periodic, i.e. it can be represented by the sum of waves of an infinite number of frequencies, which we call the fourier transform, and the waves of these frequencies analyzed we call the harmonics. Filtering is an operation of filtering out frequencies in a specific band in a signal, and is an important measure for suppressing and preventing interference. Is a probability theory and method for estimating another random process related to the random process according to the result of observing the random process. The filtering of digital signals is common in signal processing, and the fields of electric measurement signal processing, sound processing, earthquake monitoring and the like all need to carry out filtering analysis on the digital signals acquired by the sensors.
The filtering can be classified into high-pass filtering, low-pass filtering, band-stop filtering, etc. The high-pass filtering means that the high-frequency signal can normally pass through, and the low-frequency signal lower than the lower limit cut-off frequency is blocked and weakened. Low pass filtering means that low frequency signals can normally pass, while high frequency signals exceeding the upper cutoff frequency are blocked and attenuated. A band-pass filter refers to a filter that can pass frequency components in a certain frequency range, but attenuate frequency components in other ranges to an extremely low level, i.e., allow the passage of signals having frequencies below an upper cutoff frequency and above a lower cutoff frequency. The band reject filter attenuates (or reduces) signals of a range of frequencies, but allows the passage of signals having frequencies below a lower cutoff frequency and above an upper cutoff frequency.
The variety of filters is wide, including butterworth filters, chebyshev filters, bessel filters, elliptic filters, etc. The filtering of the digital signals generally adopts the mixed programming mode of Matlab and language platforms such as C#, C++, labview and the like, and various filtering algorithms are provided by the Matlab. However, when MatLab hybrid programming is performed, MCR needs to be loaded for the first time, namely MATLAB Compiler Runtime, and when Matlab compiles a runtime library, more than 3 seconds are needed for the first time loading of 10000 sampling points, so that the efficiency is low. The parameters of the filtering algorithm provided by Matlab are complex, require the skilled use of Matlab, and require familiarity with the methods of use of the various filtering algorithms, and even require extensive time to study the principles of each algorithm. The National Instruments control is also used for integrating the filter, and the National instruments, analysis, laboratory library provides tens of filtering algorithms, but the filter is provided in a class mode, has less data such as parameter description, and the like, is not provided with a visual filtering control, is difficult to use, and is inconvenient for a beginner or a user to use. Each algorithm of the National Instruments control is realized in a class form, and each filtering algorithm is a class which is equivalent to a component module, so that the filter is not convenient to use and needs to be understood by a relatively professional software technician. The industry now has no integrated filter control that is both comprehensive and easy to use.
Disclosure of Invention
The invention provides an integrated digital filtering harmonic wave display method and system, which are used for overcoming the defects of low loading efficiency, inconvenient use and the like in the prior art, improving the loading efficiency, integrating various filtering algorithms, realizing visual operation, being convenient and flexible to use and improving user experience.
In order to achieve the above object, the present invention provides an integrated digital filter harmonic display method, comprising:
sampling the original waveform data to obtain sampling data;
performing Fourier transform on the sampling data to obtain harmonic data of the original waveform data;
displaying the harmonic data as a bar graph in a display area of the user interface;
respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input area of the input module in the user interface;
determining a filtering algorithm according to the mapping relation between the filtering parameters and a filtering algorithm classification table, and carrying out convolution operation on the original waveform data and the harmonic data according to the filtering algorithm to obtain filtered waveform data;
the filtered waveform data is displayed in a viewing area of the user interface.
To achieve the above object, the present invention also provides an integrated digital filter harmonic display system, comprising:
the sampling module is used for sampling the original waveform data to obtain sampling data;
the transformation module is used for carrying out Fourier transformation on the sampling data to obtain harmonic data of the sampling data;
the display module is used for displaying the sampling data of the original waveform data in a display area of the user interface and displaying the harmonic data as a histogram;
the parameter module is used for respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input area of the input module in the user interface;
the convolution module is used for determining a filtering algorithm according to the mapping relation between the filtering parameters and the filtering algorithm classification table, and carrying out convolution operation on the original waveform data and the harmonic data according to the filtering algorithm to obtain filtered waveform data;
the display module is further configured to display the filtered waveform data in a viewing area of the user interface.
According to the integrated digital filtering harmonic wave display method and system, the original waveform data can be received through the data interface of the software, the original waveform data can be sampled through the sampler, namely the sampling module, the sampled data are obtained, and then the original waveform is displayed in the viewing area of the user interface according to the sampled data through the existing waveform control, so that the user can observe the original waveform conveniently; then, the original waveform sampling data is subjected to Fourier transform through a waveform control to form a Fourier series view, namely a histogram, a frequency interval is selected on the view, the waveform of the designated frequency interval can be obtained after the waveform is filtered through a filtering algorithm provided by National Instruments, and the original waveform in the checking area is changed into the filtered waveform; in the data processing process, only two steps of scheduling are needed, wherein the first step of scheduling waveform control displays a histogram, and the second step of scheduling a filtering algorithm carries out convolution operation and displays a filtering waveform; the data processing flow is realized by the following specific modes: integrating the C# control, the National Instruments control and the FFTW3 control, performing Fourier transform on the original waveform sampling data to obtain a Fourier series view, namely a histogram, through the FFTW3 control, selecting a frequency interval on the view, and filtering through the National Instruments control to obtain a waveform of a specified frequency interval. According to the invention, through the integrated control under the Winform platform, the Fourier series harmonic data display, the filter selection, the visual selection of filtering parameters such as upper and lower limit cut-off frequencies and the like, and the integrated digital filtering harmonic display control of waveform display are integrated, so that the integration degree is high, the parameters are selectable, the operation is flexible, and the use is convenient. The control C# and the control National Instruments are both netFranwork platforms, the access is not limited, the FFTW3 control is a C++ platform, the C# is not delayed in accessing the C++ platform, the efficiency performance is good, and the MCR loading efficiency is better than that of Matlab mixed programming. The control integrates the filtering algorithm provided by the National Instruments control, and the calling steps only need two steps, which is much simpler than the algorithm provided by Matlab.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to the structures shown in these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of an interface for waveform input and frequency domain analysis and filtering parameter selection in a specific implementation of a method and system for displaying digitally filtered harmonics according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of an interface of the frequency bin selection and filtering output of FIG. 1;
FIG. 3 is a second interface diagram of input waveforms and frequency domain analysis and filtering parameter selection in the first embodiment of the present invention;
FIG. 4 shows the effect of the low pass filter of FIG. 3 using a Chebyshev filter with a lower frequency limit of 25 Hz;
FIG. 5 shows the effect of the high pass filter in FIG. 3 using a Chebyshev filter with an upper frequency limit of 50 Hz;
FIG. 6 shows the effect of the band-pass filter in FIG. 3 using a Chebyshev filter with a lower frequency limit of 25Hz and an upper frequency limit of 50 Hz;
fig. 7 shows the effect of the band reject filter in fig. 3 using chebyshev filters with a lower frequency limit of 25Hz and an upper frequency limit of 50 Hz.
The achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
It should be noted that all directional indicators (such as up, down, left, right, front, and rear … …) in the embodiments of the present invention are merely used to explain the relative positional relationship, movement, etc. between the components in a particular posture (as shown in the drawings), and if the particular posture is changed, the directional indicator is changed accordingly.
Furthermore, descriptions such as those referred to as "first," "second," and the like, are provided for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implying an order of magnitude of the indicated technical features in the present disclosure. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present invention, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
In the present invention, unless specifically stated and limited otherwise, the terms "connected," "affixed," and the like are to be construed broadly, and for example, "affixed" may be a fixed connection, a removable connection, or an integral body; the device can be mechanically connected, electrically connected, physically connected or wirelessly connected; either directly or indirectly, through intermediaries, or both, may be in communication with each other or in interaction with each other, unless expressly defined otherwise. The specific meaning of the above terms in the present invention can be understood by those of ordinary skill in the art according to the specific circumstances.
In addition, the technical solutions of the embodiments of the present invention may be combined with each other, but it is necessary to be based on the fact that those skilled in the art can implement the technical solutions, and when the technical solutions are contradictory or cannot be implemented, the combination of the technical solutions should be considered as not existing, and not falling within the scope of protection claimed by the present invention.
The invention provides an integrated digital filtering harmonic wave display method and system.
Example 1
Referring to fig. 1-7, the present invention provides an integrated digital filter harmonic display method, which includes:
s1, sampling original waveform data, obtaining sampling data, and displaying the sampling data in a display area of a user interface;
the sampling points are a series of points which change with time, are graphically represented as a time-domain spectrogram, the abscissa is time, the ordinate is amplitude, and the higher the sampling frequency is, the closer the sampling data is to the original waveform;
s2, carrying out Fourier transform on the sampling data to obtain harmonic data of the original waveform data;
carrying out Fourier transform on the input original waveform through an FFTW3 control to obtain harmonic data of the original waveform; the harmonic data comprises a peak-to-peak array, a phase array and a harmonic frequency array;
s3, displaying the harmonic data in a display area of the user interface as a bar graph;
referring to fig. 1, the lower half part is a viewing area, and is displayed as a set of adopted data of original waveform data, specifically a time-domain spectrogram, the upper half part is a display area, and is displayed as harmonic data, and a fourier series view obtained by performing fourier transform on the sampled data through a waveformggraph control of a National Instruments module is obtained, wherein the view is a histogram, the histogram takes the times of harmonic as an abscissa, and the peak-to-peak value of the harmonic as an ordinate; the specific method is that the FFTSingled method of the FFTW3 control is called, the original waveform, the sampling rate, the output peak-to-peak value array, the phase array and the frequency array are input, wherein the index of the array is the harmonic frequency. The WaveformGaN control is a waveform display control, is set to be displayed as a histogram, calls a parameter transmission interface PlotY of the waveform display control, and transmits a peak-to-peak array to obtain a harmonic data histogram; taking the peak-to-peak array as an ordinate, taking zero as a starting point, calling a PlotY method of the waveformggraph control of National Instruments with a gap value of 1, and forming a histogram to show harmonic waveforms.
S4, respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input module on the input area in the user interface;
the upper limit cut-off frequency and the lower limit cut-off frequency are determined by moving two cursors in a histogram interface of harmonic data, and filtering parameters are respectively selected by a pull-down frame or a pull-down menu; the input module can be an electronic device such as a keyboard or a mouse or a touch control component;
the step of obtaining the filtering parameters according to the selection of the input module on the histogram in S4 includes:
s41a, a longitudinal cursor capable of moving along the horizontal coordinate direction along with the screen cursor on the histogram is arranged on the histogram; referring to fig. 1, two longitudinal cursors are generated on the bar graph, wherein the longitudinal cursors on the left side in the direction shown in fig. 1 are frequency lower limit cursors, the longitudinal cursors on the right side are frequency upper limit cursors, and the two longitudinal cursors can be respectively moved to a designated position under the operation of a screen cursor, and a cut-off frequency lower limit and a cut-off frequency upper limit are correspondingly selected;
s41b, controlling the screen cursor to move to the position where the longitudinal cursor is positioned in the input area according to the output signal of the input module and the output signal of the comparator; the comparator output signal is used for describing a comparison result of the current position information and the final position information of the screen cursor; the user controls the input module, the input module stores the position address or the position coordinate in the memory, and the comparator controls the screen cursor to move in the input area by calling the position coordinate of the input module and the position coordinate of the movable area and comparing the position coordinate with the position coordinate of the movable area, so that the screen cursor can be placed on one of the longitudinal cursors, and the horizontal coordinate of the screen cursor is the same as the horizontal coordinate of the longitudinal cursor at the moment; the principle is similar to the cursor movement on a display screen controlled by a mouse.
S41c, when the horizontal positions of the screen cursor and the longitudinal cursor are the same, relating the horizontal position information of the longitudinal cursor with the horizontal position information of the screen cursor according to the output signal of the input module, and controlling the longitudinal cursor to move along with the screen cursor together with the output signal of the comparator; when judging whether the horizontal positions of the screen cursor and the longitudinal cursor are the same, the comparator outputs a signal for describing a comparison result of the current position information of the screen cursor and the current position information of the longitudinal cursor; in the process of controlling the longitudinal cursor to move along with the screen cursor, the comparator output signal is used to describe the current horizontal position information of the screen cursor and the end position information in the movable region (which of the screen cursor and the longitudinal cursor has a small region boundary value), and the principle is the same as that of step S41b, and in an embodiment of the invention, for example: the user identifies through naked eyes that the water level positions of the two cursors are the same, triggers the input module to output a relevant control signal, at the moment, the comparator compares the current position information of the screen cursors with the current position information of the longitudinal cursors, when the horizontal position information (horizontal position address) of the two cursors is the same, the two horizontal position information are relevant, at the moment, the input module moves or outputs a movement control signal, and the relevant longitudinal cursors move together with the screen cursors; otherwise, when the comparator compares the current position information of the screen cursor with the current position information of the longitudinal cursor, and when the two horizontal position information (horizontal position address) are different, the two horizontal position information are not associated when the trigger input module outputs an associated control signal, and the longitudinal cursor does not move along with the movement of the screen cursor; the following is a description of one specific embodiment:
when the horizontal positions of the screen cursor and the longitudinal cursor are the same, that is, the current position information of the screen cursor and the current position information of the longitudinal cursor are compared by the comparator, and when the horizontal position information (horizontal position address) of the screen cursor and the current position information of the longitudinal cursor are the same, the left mouse button is clicked and continuously triggered, the horizontal position information of the longitudinal cursor and the horizontal position information of the screen cursor are associated, and meanwhile, when the screen cursor is moved, the longitudinal cursor is dragged to move along with the horizontal cursor in the abscissa direction, for example: the lower limit cursor of the dragging frequency moves to the position of 25Hz along with the screen cursor, the left mouse button is loosened, and the control obtains 25Hz as the lower limit of the cut-off frequency, see FIG. 2; in the operation process, when the mouse serving as an input module clicks a left button, an output signal correlates the current position coordinate of the frequency lower limit cursor with the position coordinate of the screen cursor, and at the moment, in the moving process of the screen cursor, the output signal of the comparator controls the screen cursor and the frequency lower limit cursor to move only in the horizontal direction, namely, the current horizontal position and the final horizontal position coordinate of the screen cursor are compared, and the cut-off frequency upper limit is operated similarly;
s41d, obtaining filtering parameters according to the output signals of the input module and the calibration parameters of the final position information of the longitudinal cursor mapped on the horizontal coordinates of the histogram; the filter parameters include a cut-off frequency. When the left button of the mouse is released, the output signal triggers the control to acquire a calibration parameter 25Hz of a horizontal coordinate (transverse coordinate) of the histogram, which is mapped at the position of the cursor of the upper frequency lower limit of the histogram, as a cut-off frequency lower limit; the process of obtaining the upper limit of the cut-off frequency is the same as above and will not be described again; in S4, the step of obtaining the filtering parameters according to the selection of the input module on the input area in the user interface includes:
S41A, controlling a screen cursor to move and select on each pull-down menu in an input area according to an output signal of an input module and an output signal of a comparator; the comparator output signal is used for describing a comparison result of the current position information and the final position information of the screen cursor; referring to fig. 1, an input module such as a mouse clicks a right-end hook in a window on the right side of a filter in an input area, a drop-down menu is displayed below the window, and various filter algorithms in a National Instruments control are collected in the drop-down menu, for example: butterworth, chebyshev, inverted chebyshev Butterworth, oval, bessel, etc., filter types are assembled in a drop down menu if the opposite hooks in the right window of the type are clicked, for example: high pass, low pass, band stop, etc., if the opposite hook in the right window of the order is clicked, the pull down menu displays: 1. 2, 3 … … natural numbers represent the operation orders of the filtering algorithm;
S41B, when detecting that a target option on a drop-down menu is selected, configuring the target option into a selected state; clicking one option in the drop-down menu by the mouse, and displaying the option in a corresponding window as a target option through configuration of the state;
S41C, extracting parameters corresponding to the target options in the selected state to obtain filtering parameters; the filtering parameters include: filter type, filter calculation order. The control extracts the filter parameters corresponding to the target options in the windows in the operation process, and controls the operation process;
the upper limit cut-off frequency and the lower limit cut-off frequency which are determined by moving two cursors in a histogram interface of harmonic data are respectively selected as filtering parameters by a pull-down frame; filtering the input waveform through National Instruments; the filtering parameters comprise filter types, filtering algorithms and filtering orders;
finding out frequencies of the corresponding positions of the two cursors in the harmonic frequency array to serve as a lower limit cut-off frequency and an upper limit cut-off frequency respectively; the filter type is determined according to the interface filter selection drop-down box selection, the filter algorithm is determined by the filter algorithm drop-down box selection, and the filter algorithm is determined, and the input order of the interface input box is shown at the top of fig. 1 and 2.
S5, determining a filtering algorithm according to the mapping relation between the filtering parameters and a filtering algorithm classification table, and carrying out convolution operation on the original waveform data and harmonic data according to the filtering algorithm to obtain filtered waveform data; filtering the input waveform through National Instruments; fig. 1 is an original waveform, here, a square wave signal of 10Hz is used for example, the waveform diagram of the lower part of fig. 2 is a waveform diagram after a chebyshev low-pass filter is selected, the lower frequency limit is 25Hz, and the 10-order operation is performed on the original waveform of fig. 1.
The filtering algorithm classification table is divided according to the types of the filters covered by each filter type; the filter types include high pass, low pass, band pass and band reject; the filter categories include: butterworth, chebyshev, inverted chebyshev Butterworth, ellipse, bessel; the mapping relation refers to the relation between the filtering type, the filtering method and the filtering algorithm selected by the pull-down menu. For example, if the interface selects the filter type chebyshev and the filtering method is bandpass, the mapping submodule chebyshev bandpass filter and other mapping relations are also described above. The filtering algorithm integrated in the National Instruments module is divided into the following twenty types: butterworth Wo Sigao pass filter, butterworth low pass filter, butterworth band stop filter, chebyshev high pass filter, chebyshev low pass filter, chebyshev band stop filter, inverted Chebyshev high pass filter, inverted Chebyshev Butterworth low pass filter, inverted Chebyshev Butterworth band stop filter, elliptical high pass filter, elliptical low pass filter, elliptical band stop filter, bessel high pass filter, bessel low pass filter, bessel band stop filter. The step S5 specifically comprises the following steps:
s51, determining a filtering algorithm according to the mapping relation among the filter types, the filter types and the filtering algorithm classification table; with continued reference to fig. 1, when the filter type is selected as low pass and the filter type is selected as chebyshev, the corresponding filtering algorithm is a chebyshev low pass filtering algorithm;
s52, carrying out convolution operation on the original waveform data and the harmonic data of the original waveform according to a filtering algorithm, a filtering cut-off frequency and a filtering calculation order to obtain filtered waveform data. Waveform data input waveform display control interface, displayed waveform see sine wave in the lower half of fig. 1 (viewing area);
and S6, after the original waveform data are filtered, the waveform data are displayed in a viewing area of the user interface.
In the view area, waveform data output by the waveform display control is displayed at any time, for example, in fig. 1, since only the lower limit of the cut-off frequency is selected and no filtering operation is performed yet, sampling data lower than the lower limit of the cut-off frequency in the original waveform data sampling data only displayed in the view area; after the filtering operation of step S5 is performed, the filtered waveform diagram is displayed in the viewing area, see the waveform diagram of the lower half of fig. 2. The method comprises the following specific steps:
s61, when the original waveform data is not filtered by the filtering algorithm, the sampled data of the original waveform data are displayed in a viewing area of a user interface;
and S62, when the original waveform data is filtered by a filtering algorithm, displaying the filtered waveform data in a viewing area of a user interface.
The filtered waveforms are displayed by the waveformggraph control of National Instruments. See the lower half of fig. 2; the corresponding filtering algorithm of National instruments, analysis, profestive of National Instruments is invoked, and the raw waveform data is input, upper cut-off frequency, lower cut-off frequency, and order. And outputting a filtered result waveform. The filtered result waveform is called the PlotY method of the waveformggraph control of National Instruments, and the waveform is displayed.
And determining the upper frequency limit and the lower frequency limit of the filtering by the positions of the upper frequency limit cursor and the lower frequency limit cursor of the harmonic data histogram and the frequency array output by the FFTW 3. And calling the filter algorithm FilterData of the corresponding submodule through the filter type and the filter algorithm type selected by the interface, inputting the original waveform, and providing the order, the upper frequency limit and the lower frequency limit by the interface. And outputting the filtered waveform data. And calling a PlotY method of the waveformggraph, and transmitting the output filtered waveform data to display the filtered waveform image.
Integrating Fourier series view display, a filtering algorithm provided by National Instruments according to filter types, visual selection of upper limit cut-off frequency and lower limit cut-off frequency and waveform display under a Winform platform; the method specifically comprises the following steps: the original waveform data is stored in a double-precision floating point array, and three new arrays with equal length, namely a harmonic peak value array, a harmonic frequency array and a harmonic phase array, are formed through Fourier transformation. And the harmonic frequency array is taken as an X axis, the harmonic peak value array is taken as a Y axis to form a harmonic histogram, and the upper limit cutoff frequency and the lower limit cutoff frequency of the filtering are determined by dragging the upper limit frequency and the lower limit frequency on the harmonic histogram. And (3) through selection of the filtering type, the variety and the order of the pull-down menu pair, a corresponding filtering algorithm is called to form a new double-precision floating point array with the length equal to that of the original waveform data array, and the double-precision floating point array is drawn on an interface to obtain the filtered waveform.
The common software technician needs to spend a great deal of time researching the principle of the filtering algorithm, and each of tens of filtering algorithms is researched once and then analyzed which method is the most suitable for the requirement. According to the embodiment of the invention, through dragging the filter into an application program of the filter, fourier series display, filter selection, upper and lower limit cut-off frequency visualization selection and waveform display integrated digital filter harmonic display control are integrated under a Winform platform, the use of ten filtering algorithms can be mastered, so that research and comparison are facilitated, and a novice programmer can also quickly complete the filtering function with complex principle. The choice of filtering algorithm can be determined by program developers or can be provided to the end user with comparative data. The research and development efficiency and the experience of the end user can be greatly improved. And the filtering requirement text description is very simple, the original waveform is provided, the upper limit cutoff frequency and the lower limit cutoff frequency are provided, and the effect after filtering can be shown.
The following describes in detail the implementation of the integration step of the Winform platform:
installing Visual Studio 2010, registering permissions, and creating a Visual C# language environment for Windows application development; installing NI_Measurement_studio_2010_for_VS2010, registering permission, and building an NI control library; studio 2010 creates a "Windows Window user interface application item directory":
the method comprises the steps of putting three files, namely a FilterAndFFT.dll file, a FilterExtension.dll file and Libfftw3-3.dll files in an NI control library, into a debug directory;
changing the platform to netframework4.0;
adding a reference to the NI control;
the references to the NI waveform control are correspondingly added;
clicking right key on the public control, clicking left key on the pop-up menu to select item;
clicking the left button to browse and pop up a file selection box;
finding the filterextension.dll back left key click single on;
newly adding an FFTWaveSHOW control (carried by an NI control library) in the toolbar, and dragging the FFTWaveSHOW control into a winform window for use; completing integration of each display window under the Winform platform; based on C# and National Instruments and FFTW3, the integrated digital filter harmonic display control of Fourier series display, filter selection, upper and lower limit cut-off frequency visualization selection and waveform display are integrated next in the Winform platform.
The original waveform is transmitted in, and the type, the type and the order of the filter are selected, see fig. 1;
after sliding the selection of a white cursor (lower cutoff frequency) and a red cursor (upper cutoff frequency), the left click "apply" can be filtered, see fig. 2.
The following is a description of the method of signal input:
wherein the parameters are defined as follows:
another set of raw waveforms is filtered by the method of embodiment one:
the original waveform is referred to in the lower half of fig. 3, harmonic data is obtained after fourier transformation is performed on the original waveform, and referring to the upper half of fig. 3, different filtering waveforms can be displayed by respectively selecting filtering parameters in a selection frame after the filter, the type and the order positioned at the top in fig. 3:
for example: chebyshev filter, low-pass filter type, order 10, lower frequency limit 25Hz, filtered waveform see lower half of fig. 4;
for example: chebyshev filter, high-pass filter type, order 10, upper frequency limit 50Hz, filtered waveform see lower half of fig. 5;
for example: chebyshev filter, band-pass filter type, order 10, lower frequency limit 25Hz, upper frequency limit 50Hz, filter waveform see lower half of fig. 6;
for example: chebyshev filter, band reject filter type, order 10, lower frequency limit 25Hz, upper frequency limit 50Hz, and filtered waveform see lower half of fig. 7.
Example two
On the basis of the first embodiment, the invention further provides an integrated digital filtering harmonic display system, which comprises:
the sampling module is used for sampling the original waveform data to obtain sampling data;
the transformation module is used for carrying out Fourier transformation on the sampling data to obtain harmonic data of the sampling data;
the display module is used for displaying the sampling data of the original waveform data in a display area of the user interface and displaying the harmonic data as a histogram;
the parameter module is used for respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input area in the user interface;
the convolution module is used for determining a filtering algorithm according to the mapping relation between the filtering parameters and the filtering algorithm classification table, and carrying out convolution operation on the original waveform data and the harmonic data according to the filtering algorithm to obtain filtered waveform data;
the display module is further configured to display the waveform data in a viewing area of the user interface after the original waveform data is filtered.
Preferably, the parameter module comprises a cut-off frequency module and a calculation parameter module;
the display module is also used for displaying a longitudinal cursor which can move along the horizontal coordinate direction on the bar graph along with the screen cursor on the bar graph;
the cut-off frequency module comprises:
the screen cursor control sub-module is used for controlling the screen cursor to move in the input area according to the output signal of the input module and the output signal of the comparator; the comparator output signal is used for describing a comparison result of the current position information and the final position information of the screen cursor;
the longitudinal cursor control sub-module is used for controlling the longitudinal cursor to move along with the screen cursor according to the output signal of the input module and the output signal of the comparator; the comparator output signal is used for describing a comparison result of the current position information of the screen cursor and the current position information of the longitudinal cursor;
the cut-off frequency sub-module is used for obtaining filtering parameters according to the final position information of the amplified output signal of the input module and the longitudinal cursor mapped on the calibration parameters of the horizontal coordinates of the histogram; the filter parameters include a cut-off frequency.
Preferably, the screen cursor control sub-module is further used for controlling the screen cursor to move and select on each pull-down menu in the input area according to the output signal of the input module and the output signal of the comparator; the comparator output signal is used for describing a comparison result of the current position information and the final position information of the screen cursor;
the calculation parameter module comprises:
a selecting sub-module, configured to configure a target option on a drop-down menu into a selected state when detecting that the target option is selected;
the extraction sub-module is used for extracting parameters corresponding to the target options in the selected state and obtaining filtering parameters; the filtering parameters include: filter type, filter calculation order.
Preferably, the convolution module includes:
the algorithm acquisition sub-module is used for determining a filtering algorithm according to the mapping relation among the filter types, the filter types and the filtering algorithm classification table;
and the convolution operation sub-module is used for carrying out convolution operation on the original waveform data and the harmonic data of the original waveform according to a filtering algorithm, a filtering cut-off frequency and a filtering calculation order to obtain filtered waveform data.
Preferably, the display module includes:
the first display sub-module is used for displaying the sampling data of the original waveform data in a viewing area of the user interface when the original waveform data is not filtered by the filtering algorithm;
the method is also used for displaying the filtered waveform data in a viewing area of a user interface when the original waveform data is filtered by a filtering algorithm;
and the second display sub-module is used for extracting the original waveform data, the sampling parameters and the harmonic data and displaying the histogram in the display area.
The foregoing description is only of the preferred embodiments of the present invention and is not intended to limit the scope of the invention, and all equivalent structural changes made by the description of the present invention and the accompanying drawings or direct/indirect application in other related technical fields are included in the scope of the invention.
Claims (9)
1. An integrated digital filtered harmonic display method, comprising:
sampling the original waveform data to obtain sampling data;
performing Fourier transform on the sampling data to obtain harmonic data of the original waveform data;
displaying the harmonic data as a bar graph in a display area of the user interface;
respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input area of the input module in the user interface;
determining a filtering algorithm according to the mapping relation between the filtering parameters and a filtering algorithm classification table, and carrying out convolution operation on the original waveform data and the harmonic data according to the filtering algorithm to obtain filtered waveform data;
displaying the filtered waveform data in a viewing area of a user interface;
the bar graph is provided with a longitudinal cursor which can move along the horizontal coordinate direction along with the screen cursor;
controlling the screen cursor to move to the position where the longitudinal cursor is positioned in the input area according to the output signal of the input module and the output signal of the comparator; when the horizontal positions of the screen cursor and the longitudinal cursor are the same, the longitudinal cursor horizontal position information is correlated with the screen cursor horizontal position information according to the output signal of the input module, and then the longitudinal cursor is controlled to move along with the screen cursor together with the output signal of the comparator;
obtaining filtering parameters according to calibration parameters of the final position information of the amplified output signal of the input module and the longitudinal cursor mapped on the horizontal coordinates of the histogram; the filter parameters include a cut-off frequency.
2. The integrated digital filtered harmonic display method of claim 1 wherein the step of separately obtaining the filter parameters based on selection of an input module on an input area in a user interface comprises:
controlling a screen cursor to move and select on each drop-down menu in an input area according to an output signal of the input module and an output signal of the comparator; the comparator output signal is used for describing a comparison result of the current position information and the final position information of the screen cursor;
when detecting that a target option on a drop-down menu is selected, configuring the target option into a selected state;
extracting parameters corresponding to the target options in the selected state to obtain filtering parameters; the filtering parameters include: filter type, filter calculation order.
3. The integrated digital filtering harmonic display method of claim 2, wherein the step of determining a filtering algorithm according to a mapping relation between filtering parameters and a filtering algorithm classification table, and convolving the original waveform data with the harmonic data according to the filtering algorithm, and obtaining the filtered waveform data comprises:
determining a filtering algorithm according to the mapping relation among the filter types, the filter types and the filtering algorithm classification table;
and carrying out convolution operation on the original waveform data and the harmonic data of the original waveform according to a filtering algorithm, a filtering cut-off frequency and a filtering calculation order to obtain filtered waveform data.
4. An integrated digital filtered harmonic display method as in any one of claims 1-3 wherein the step of displaying the filtered waveform data in a viewing area of a user interface comprises:
when the original waveform data is not filtered by the filtering algorithm, displaying the sampling data of the original waveform data in a viewing area of a user interface;
when the original waveform data is filtered by the filtering algorithm, the filtered waveform data is displayed in a viewing area of the user interface.
5. An integrated digitally filtered harmonic display system, comprising:
the sampling module is used for sampling the original waveform data to obtain sampling data;
the transformation module is used for carrying out Fourier transformation on the sampling data to obtain harmonic data of the sampling data;
the display module is used for displaying the sampling data of the original waveform data in a display area of the user interface and displaying the harmonic data as a histogram;
the parameter module is used for respectively acquiring filtering parameters according to the selection of the input module on the histogram and the input area of the input module in the user interface;
the convolution module is used for determining a filtering algorithm according to the mapping relation between the filtering parameters and the filtering algorithm classification table, and carrying out convolution operation on the original waveform data and the harmonic data according to the filtering algorithm to obtain filtered waveform data;
the display module is further used for displaying the filtered waveform data in a viewing area of the user interface; the bar graph is provided with a longitudinal cursor which can move along the horizontal coordinate direction along with the screen cursor; controlling the screen cursor to move to the position where the longitudinal cursor is positioned in the input area according to the output signal of the input module and the output signal of the comparator; when the horizontal positions of the screen cursor and the longitudinal cursor are the same, the longitudinal cursor horizontal position information is correlated with the screen cursor horizontal position information according to the output signal of the input module, and then the longitudinal cursor is controlled to move along with the screen cursor together with the output signal of the comparator; obtaining filtering parameters according to calibration parameters of the final position information of the amplified output signal of the input module and the longitudinal cursor mapped on the horizontal coordinates of the histogram; the filter parameters include a cut-off frequency.
6. The integrated digital filtered harmonic display system as in claim 5 wherein the parameter module comprises a cut-off frequency module and a calculate parameter module;
the display module is also used for displaying a longitudinal cursor which can move along the horizontal coordinate direction on the bar graph along with the screen cursor on the bar graph;
the cut-off frequency module comprises:
the screen cursor control sub-module is used for controlling the screen cursor to move to the position where the longitudinal cursor is positioned in the input area according to the output signal of the input module and the output signal of the comparator;
the vertical cursor control sub-module is used for controlling the vertical cursor to move along with the screen cursor together with the output signal of the comparator after the horizontal position information of the vertical cursor is correlated with the horizontal position information of the screen cursor according to the output signal of the input module when the horizontal positions of the screen cursor and the vertical cursor are the same;
the cut-off frequency sub-module is used for obtaining filtering parameters according to the final position information of the amplified output signal of the input module and the longitudinal cursor mapped on the calibration parameters of the horizontal coordinates of the histogram; the filter parameters include a cut-off frequency.
7. The integrated digital filtered harmonic display system of claim 6 wherein the screen cursor control sub-module is further configured to control the screen cursor to move and select on each of the pull-down menus in the input region based on the output signal of the input module and the comparator output signal; the comparator output signal is used for describing a comparison result of the current position information and the final position information of the screen cursor;
the calculation parameter module comprises:
a selecting sub-module, configured to configure a target option on a drop-down menu into a selected state when detecting that the target option is selected;
the extraction sub-module is used for extracting parameters corresponding to the target options in the selected state and obtaining filtering parameters; the filtering parameters include: filter type, filter calculation order.
8. The integrated digital filtered harmonic display system as in claim 7 wherein the convolution module comprises:
the algorithm acquisition sub-module is used for determining a filtering algorithm according to the mapping relation among the filter types, the filter types and the filtering algorithm classification table;
and the convolution operation sub-module is used for carrying out convolution operation on the original waveform data and the harmonic data of the original waveform according to a filtering algorithm, a filtering cut-off frequency and a filtering calculation order to obtain filtered waveform data.
9. The integrated digital filtered harmonic display system of any one of claims 5-8 wherein the display module comprises:
the first display sub-module is used for displaying the sampling data of the original waveform data in a viewing area of the user interface when the original waveform data is not filtered by the filtering algorithm;
the method is also used for displaying the filtered waveform data in a viewing area of a user interface when the original waveform data is filtered by a filtering algorithm;
and the second display sub-module is used for extracting the original waveform data, the sampling parameters and the harmonic data and displaying the histogram in the display area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910318263.XA CN110046435B (en) | 2019-04-19 | 2019-04-19 | Integrated digital filtering harmonic wave display method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910318263.XA CN110046435B (en) | 2019-04-19 | 2019-04-19 | Integrated digital filtering harmonic wave display method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110046435A CN110046435A (en) | 2019-07-23 |
| CN110046435B true CN110046435B (en) | 2023-06-09 |
Family
ID=67277906
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910318263.XA Active CN110046435B (en) | 2019-04-19 | 2019-04-19 | Integrated digital filtering harmonic wave display method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110046435B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113220256A (en) * | 2021-05-11 | 2021-08-06 | 马鞍山状元郎电子科技有限公司 | Interactive design whiteboard |
| CN115060995B (en) * | 2022-06-14 | 2023-01-24 | 费莱(浙江)科技有限公司 | Harmonic management method and system for intelligent filtering |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004184536A (en) * | 2002-11-29 | 2004-07-02 | Mitsubishi Electric Corp | Device and program for convolutional operation |
| US7649927B1 (en) * | 2006-12-29 | 2010-01-19 | Kiomars Anvari | Equalizer filter with dynamically configurable convolution filter |
| CN109448023A (en) * | 2018-10-23 | 2019-03-08 | 武汉大学 | A kind of satellite video Small object method for real time tracking of combination space confidence map and track estimation |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5053983A (en) * | 1971-04-19 | 1991-10-01 | Hyatt Gilbert P | Filter system having an adaptive control for updating filter samples |
| CA2324832A1 (en) * | 2000-10-31 | 2002-04-30 | Roger Lecomte | Real-time image reconstruction for computed tomography systems |
| CN101860344A (en) * | 2010-06-12 | 2010-10-13 | 刘海成 | Construction method of frequency-selecting filter and construction method for realizing FIR-type and IIR-type filters by adopting same |
| CN104572388B (en) * | 2013-10-15 | 2017-01-25 | 鸿富锦精密工业（深圳）有限公司 | Data filtering system, data filtering device and method |
| CN106059529A (en) * | 2016-05-27 | 2016-10-26 | 大连楼兰科技股份有限公司 | Data filtering method integrating multiple kinds of adjustable filters and system integrating multiple kinds of adjustable filters |
| CN106487357B (en) * | 2016-07-08 | 2019-02-19 | 青岛歌尔声学科技有限公司 | A kind of digital filter design device and method |
| JP6724828B2 (en) * | 2017-03-15 | 2020-07-15 | カシオ計算機株式会社 | Filter calculation processing device, filter calculation method, and effect imparting device |
-
2019
- 2019-04-19 CN CN201910318263.XA patent/CN110046435B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004184536A (en) * | 2002-11-29 | 2004-07-02 | Mitsubishi Electric Corp | Device and program for convolutional operation |
| US7649927B1 (en) * | 2006-12-29 | 2010-01-19 | Kiomars Anvari | Equalizer filter with dynamically configurable convolution filter |
| CN109448023A (en) * | 2018-10-23 | 2019-03-08 | 武汉大学 | A kind of satellite video Small object method for real time tracking of combination space confidence map and track estimation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110046435A (en) | 2019-07-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110046435B (en) | Integrated digital filtering harmonic wave display method and system | |
| US6477472B2 (en) | Analyzing signals generated by rotating machines using an order mask to select desired order components of the signals | |
| CN1725024A (en) | Realtime power mask trigger | |
| US10945602B2 (en) | Monitoring parameter threshold value setting method and monitoring system | |
| EP3827226A1 (en) | Systems and methods for tagging and linking acoustic images | |
| CN106679741A (en) | Anti-interference signal processing method and system based on vortex shedding flowmeter | |
| US20070176931A1 (en) | Method of automatically displaying medical measurement data | |
| US8379703B2 (en) | Signal analyzer and method for displaying frequency domain data | |
| US20020183948A1 (en) | Time varying harmonic analysis including determination of order components | |
| Abed et al. | Performance analysis and computational cost evaluation of high-resolution time-frequency distributions derived from compact support time-lag kernels | |
| CN104378115A (en) | AD sampling software filtering method | |
| CN109350031A (en) | A kind of photoplethysmographic signal feature point detecting method and device | |
| US20040148116A1 (en) | Estimating a plurality of tones in an input signal | |
| US8108166B2 (en) | Analysis of chirp frequency response using arbitrary resampling filters | |
| CN108344795A (en) | Oil-gas pipeline defect identification method, device and electronic equipment | |
| CN109949897A (en) | Display methods, device and the terminal device of monitoring curve | |
| CN106725415A (en) | The treating method and apparatus of electricity physiological signal | |
| CN107843767B (en) | Signal processing method and system for spectrum analysis | |
| CN105738697A (en) | PXI-based virtual frequency spectrum analyzer | |
| CN110441822A (en) | A kind of seismic data noise attenuation method, apparatus, electronic equipment and storage medium | |
| CN107703541B (en) | The determination method and apparatus of stratigraphic dip | |
| US20150301085A1 (en) | Automatically Capturing Data Sets of Interest from a Data Acquisition Data Stream | |
| Mendes et al. | A particle counting system for calculation of bedload fluxes | |
| JP2730876B2 (en) | Novel leak detection method in a small leak gas atmosphere in a plant | |
| US20190200887A1 (en) | Information processing device, biomedical-signal measuring system, and recording medium storing program code |
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 |