CN116595300B - Signal processing method, system, medium and equipment based on complex evolution operation - Google Patents
Signal processing method, system, medium and equipment based on complex evolution operation Download PDFInfo
- Publication number
- CN116595300B CN116595300B CN202310720297.8A CN202310720297A CN116595300B CN 116595300 B CN116595300 B CN 116595300B CN 202310720297 A CN202310720297 A CN 202310720297A CN 116595300 B CN116595300 B CN 116595300B
- Authority
- CN
- China
- Prior art keywords
- angle
- complex
- adjusted
- signal processing
- product
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 20
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 14
- 238000012935 Averaging Methods 0.000 claims abstract description 12
- 238000003491 array Methods 0.000 claims abstract description 6
- 238000004590 computer program Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 230000014509 gene expression Effects 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 3
- 239000000047 product Substances 0.000 description 24
- 238000005259 measurement Methods 0.000 description 16
- 238000010586 diagram Methods 0.000 description 13
- 238000001914 filtration Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0602—Systems characterised by the synchronising information used
- H04J3/0617—Systems characterised by the synchronising information used the synchronising signal being characterised by the frequency or phase
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/145—Square transforms, e.g. Hadamard, Walsh, Haar, Hough, Slant transforms
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Signal Processing (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Feedback Control In General (AREA)
- Measurement Of Current Or Voltage (AREA)
Abstract
The invention belongs to the technical field of signal processing, provides a signal processing method, a system, a medium and equipment based on complex evolution operation, and provides the following scheme aiming at the evolution problem in the complex operation of signals: representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers; wherein each complex number comprises two arrays of double-precision numbers, the first number is a module, and the second number is an angle; adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted by taking one angle as a reference so that the relative position of each angle is unchanged; performing evolution operation on the product of the plurality of complex numbers after adjustment to obtain a processed signal; the square operation comprises square operation of the product of the die and angle averaging. The method realizes that the angle information is not lost by the complex prescription, so that the engineering application requiring the complex prescription is realized.
Description
Technical Field
The invention belongs to the technical field of signal processing, and particularly relates to a signal processing method, a system, a medium and equipment based on complex evolution operation.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
Since ac signals and associated digital signal processing are represented by various quantities such as voltage, current, gain, impedance (including resistance, capacitance, inductance) and admittance, etc., the angular component thereof is referred to as phase;
in the conventional processing for the evolution problem in the complex operation of the signals, the evolution result only needs to restore the data before the evolution after multiplication, for example, the partial angle information of the signals is permanently lost or the period is changed in the conventional way, so that the signal data obtained by processing cannot participate in the operation, and the obtained result cannot be interpreted and applied to the engineering.
Meanwhile, the main error source of some current phase sensitive signal measuring instruments is phase difference fluctuation caused by synchronous time jitter, and the phase difference needs to be filtered independently, but an algorithm for filtering the phase difference specially is not adopted conventionally, so that the accuracy cannot be further improved.
Disclosure of Invention
In order to solve at least one technical problem in the background art, the invention provides a signal processing method, a system, a medium and equipment based on complex evolution operation, which realize that complex evolution does not lose angle information by adjusting angles and realize engineering application needing complex evolution.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a first aspect of the present invention provides a signal processing method based on complex square operation, including the steps of:
acquiring an alternating current or related digital signal to be processed;
representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers; wherein each complex number comprises two arrays of double-precision numbers, the first number is a module, and the second number is an angle;
taking one angle as a reference, and adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted;
performing squaring operation on the product of the plurality of complex numbers subjected to angle adjustment to obtain a processed signal; the square operation comprises square operation of the product of the die and angle averaging.
Further, the adjusting the angle to be adjusted based on the difference between the reference angle and the angle to be adjusted with one of the angles as a reference includes:
at an angle of 0 to 2 pi, expressed byFor reference, pair->Make adjustments one by one, if->Then->Plus 2π, otherwise if->Then->Minus 2 pi, wherein->For reference angle, +.>For the angle to be adjusted, n= 2~N;
or (b)
0 to 360 is expressed as 0 to 2 pi angle, toFor reference, pair->Make adjustments one by one, if->Then->Plus 360 DEG, otherwise if%>Then->Minus 360 DEG, wherein->For reference angle, +.>For the angle to be adjusted, n= 2~N;
or (b)
0 to 65536 is 0 to 2 pi angle, toFor reference, pair->Make adjustments one by one, if->Then->Plus 65536 if otherwise->Then->The 65536 is subtracted; />For reference angle, +.>For the angle to be adjusted, n= 2~N.
Further, after the plurality of angles are adjusted, the plurality of angles may be ordered in order of magnitude and weighted average to achieve angle filtering.
Further, the square operation is to open the product of the N complex numbers to the power of N, where the modulus is the product of the N complex numbers modulo the power of N, and the angle is the average of the N angles.
Further, the expression of the evolution operation is:
wherein N is the number of a plurality of numbers,is the corresponding mode of the nth complex number, +.>Is the angle corresponding to the nth complex number.
Further, the alternating or related digital signal to be processed is obtained by adopting a channel exchange method.
Further, the phasors of the alternating or related digital signals are represented as products of a plurality of complex numbers using polar representations.
A second aspect of the present invention provides a signal processing system based on complex square operation, comprising:
a signal acquisition module configured to: acquiring an alternating current or related digital signal to be processed;
a complex representation module configured to: representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers; wherein each complex number comprises two arrays of double-precision numbers, the first number is a module, and the second number is an angle;
an angle adjustment module configured to: taking one angle as a reference, and adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted;
a data reduction module configured to: performing evolution operation on the product of the plurality of complex numbers subjected to angle adjustment to obtain a processed signal; the square operation comprises square operation of the product of the die and angle averaging.
A third aspect of the present invention provides a computer-readable storage medium.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, implements the steps of a signal processing method based on complex-valued evolution operation as described in the first aspect.
A fourth aspect of the invention provides a computer device.
A computer apparatus comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing steps in a signal processing method based on complex evolution operations as described in the first aspect when the program is executed.
Compared with the prior art, the invention has the beneficial effects that:
1. according to the invention, the relative positions of all angles are unchanged through the angle adjustment method, under the condition that the relative positions of all angles are unchanged, all angles participating in averaging are continuously distributed, and cycle crossing points are avoided, so that four-quadrant angle information is obtained, the situation that the angle information is not lost by complex evolution is realized, and the engineering application requiring complex evolution is realized.
2. The phase difference filter device realizes the phase difference filter function of phasors, so that the application of measuring the phase difference with high precision is realized, and the phase difference measurement precision is greatly improved.
3. The channel exchange method eliminates the error of the measurement channel and greatly improves the measurement accuracy.
Additional aspects of the invention will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the invention.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention.
Fig. 1 is a schematic flow chart of a signal processing method based on complex square operation according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a channel switching method measurement structure according to an embodiment of the present invention; in fig. 2, (a) is a schematic diagram of a structure in which channel switching measurement is not adopted, and in fig. 2, (b) is a schematic diagram of a structure in which channel switching measurement is adopted;
FIG. 3 is a schematic diagram of an angle adjustment method according to an embodiment of the present invention; wherein (a) in FIG. 3 isThe pair +.>Is (b) in FIG. 3>The pair of +.>Is adjusted;
FIG. 4 is a programming flow diagram of complex product open square provided by an embodiment of the present invention;
fig. 5 is a schematic diagram of a phase difference filtering method according to an embodiment of the present invention.
Detailed Description
The invention will be further described with reference to the drawings and examples.
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present invention. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
As mentioned in the background art, for the problem of evolution in complex traditional operation in the traditional signal processing, according to the traditional understanding, the evolution result only needs to be multiplied to restore the data before evolution, for example, according to the traditional waySquare opening is obtained->The angle information between 180 DEG and 360 DEG is permanently lost, or the result of squaring the complex number is changed from 360 DEG to 180 DEG, so that the data cannot continue to participate in the operation.
The invention provides a new idea, which is specifically as follows: a single complex square may lose angular information, so that taking a product square of multiple phasors ensures that angular information is not lost, or 360 ° periodicity of normal phasors is maintained.
Example 1
As shown in fig. 1, the present embodiment provides a signal processing method based on complex square operation, which includes the following steps:
step 1: acquiring an alternating current or related digital signal to be processed;
in this embodiment, the ac or related digital signal to be processed is obtained by a channel switching method.
As shown in fig. 2, for example, the signal to be processed is two-way voltage to be measuredThe result to be obtained is two paths of voltages to be measuredRatio->Wherein->Are phasors, FIG. 2 contains two measurement channels, wherein (a) in FIG. 2 +.>The gain is measured via the first channel to be +.>，/>The gain is measured via the second channel to be +.>(b) in FIG. 2>The gain is measured via the second channel to be +.>，/>The gain is measured via the first channel to be +.>，/>And->There is an error and only a short time to stabilize.
Before channel exchange, as shown in fig. 2 (a), the obtained measurement values are:
（1）
（2）
wherein,,to measure the voltage gain of channel a, +.>For measuring the voltage gain of channel b;
after channel exchange, as shown in fig. 2 (b), the obtained measurement values are:
（3）
（4）
the above formula is obtained according to [ (1) × (4) ]/[ (2) × (3) ]:
（5）
therefore, the voltage ratio Kv is independent of the voltage gain of the measurement channel, the above-mentioned quantities are phasors, the four phasors on the right side of (5) are expressed as the product of two complex numbers, and the square is further opened to obtainEliminates error source->And->。
The following describes the formula (5) further using polar representation;
assume that the input two voltages to be measured are respectively，V x The voltage is 3V with an initial phase of 2,V n the initial phase of the voltage 2V is 3 DEG,>i.e. +.>For measuring the voltage gain of channel a, it increases the amplitude of the signal to be measured by a factor of 1.1 and by a phase error of 1 °j->To measure the voltage gain of channel b, it increases the amplitude of the signal to be measured by a factor of 1.2 and by 356 phase error, as the main error source +.> Only a short time is kept stable, so that the error can be eliminated by adopting a channel exchange method.
According to the complex arithmetic rule, the complex arithmetic rule is obtained by the formulas (1) and (2):
channel switching is considered The reaction product is obtained by the formulas (3) and (4):
according to the formula (5), the square is obtained by adopting the angle averaging method of the invention:
this is just the ratio of the two voltages to be measured, 1.5, and the phase difference between the two voltages, namely-1 ° (359 °), the above-mentioned exchange process is continuously performed, the errors of the measurement channels are completely eliminated, and the above-mentioned square opening requires two complex numbers to participate and angular averages are taken according to rules.
The scheme has the advantages that by adopting the channel exchange method, the measurement error of the voltage amplitude ratio can be less than 1ppm (0.0001 percent), the phase difference measurement error can be less than 0.00005 degrees, and the ultra-high-precision alternating-current voltage phasor ratio measurement is realized.
Step 2: representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers;
expressing the complex number to be opened to the power of N as the product of N complex numbers, wherein N > =2, the complex number to be opened in actual engineering calculation is generally a formula with a plurality of complex variables, and the requirement can be met;
the ac signal and the associated digital signal are processed by various amounts, such as voltage, current, gain, impedance (including resistanceCapacitance->Inductance->) And admittance, etc., are called phasors, which can be represented by complex numbers, the angular portions of which are also called phases, and the calculation formulas are custom written by lettersjInstead ofi。
After calculating the final result, the result can be interpreted, for example, after calculating the impedance, the impedance can be interpreted as resistance, capacitance, inductance and combinations thereof, and after calculating the voltage and current, the active current, reactive current, active power, reactive power and the like can be obtained.
The complex number is composed of real and imaginary parts, e.g., z=a+ib, where a, b are real numbers,the complex number satisfies the four rules of operation, the addition and subtraction of the complex number are that the real part and the imaginary part are added and subtracted respectively, and the complex number is multiplied with the zero degree>The complex representation method comprises a complex component representation method and a complex polar coordinate representation method.
Wherein the component representation of the complex number is: the complex division time divides the numerator into a denominator and multiplies the conjugate of the denominator (the imaginary part of the denominator takes the negative), and can change the denominator into a real number, and then divides the real part and the imaginary part of the numerator respectively, for example:
(6)
the polar coordinate representation of the complex number is: in the X-Y plane, the X-direction represents the real part a and the Y-direction represents the imaginary part b, which is also called complex plane. Drawing a straight line from the point (a, b) to the point 0, wherein the length of the straight line is called a complex number of modes r, the included angle theta between the straight line and the X axis is called an angle, and the complex number of the polar coordinate format is expressed asThe complex number represented by the polar coordinate can be conveniently used for the multiply-divide operation and the evolution operation mentioned in the present invention.
For example:
(7)
the two expressions of the complex number can be converted from each other, and the triangle relation can be known:
(8)
both representation methods of the complex number and the respective operations and transformations between each other do not affect the final calculation result.
Each complex number is composed of an array of two double-precision numbers, the first number being the real part and the second number being the imaginary part.
In this embodiment, a complex polar coordinate representation is adopted, which specifically includes:
firstly, complex numbers are converted into a polar coordinate format, the polar coordinate format is composed of an array containing two double-precision numbers, the first number is a module, and the second number is an angle.
Is required in deriving and interpreting the results of the formulaParticipate and rely on->This important feature, however, software computation stores only the real and imaginary parts of the complex number, or the modulus and angle, and performs computation or conversion according to the imaginary rule, without giving itiThe addresses are allocated separately and are not usediParticipate in the operation.
Step 3: adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted by taking one angle as a reference so that the relative position of each angle is unchanged;
the angles after the evolution are actually the average of the angles of the respective complex numbers participating in the evolution, however complex numbers have a periodicity of 2 pi or 360, these angles cannot be directly averaged, for example, averaging 1 deg. and 359 deg. to 180 deg. is wrong, the correct answer should be 0 deg., or averaging 179 deg. and-179 deg. to 0 deg. is wrong, the correct answer is 180 deg.. The arithmetic average value of the angles is calculated, the angles are adjusted firstly, and under the condition that the relative positions of the angles are unchanged, the angles participating in the average are continuously distributed, and cycle crossing points are avoided.
By dividing the product of N complex numbers by N power, the modulus is the product of N complex numbers by N power, the angle is the average value of N angles, and since there may be a period cross breakpoint between these angles, the average value cannot be directly calculated, and therefore, adjustment of the scheduling is required.
The angle adjustment rule is as follows: taking any angle as a reference, and adjusting the rest angles one by one; the adjusting method is thatFor reference, pair->(n= 2~N) one by one, if +.>-180 °, then->Plus 360 deg., otherwise if180 DEG, then->The 360 degrees are subtracted, so that the relative positions of all angles are not changed and are continuously distributed, and the period break points of 0 degrees are not existed, thereby eliminating the period break points; the angles may then be arithmetically averaged or may be ordered to achieve weighted filtering.
The angular adjustment method schematic shown in fig. 3 explains how this adjustment method can eliminate the period break point so that all angles become continuous; in FIG. 3 (a) showsThe pair +.>The adjustment conditions of (1) at this time are 0 to->The angle between +180° is not to be adjusted, since they are identical to +.>Not exceeding + -180 DEG, but other angles are subtracted by 360 DEG to become negative angles, e.g. +.>If->Then it would be adjusted to-1 deg. to eliminate the period collapse break point and the final average result of 0 deg. is correct. In FIG. 3 (b) is shown +.>The pair of +.>Is adjusted at this time->The angle between-180 and 360 DEG is not to be adjusted, since they are identical to +.>Not exceeding + -180 DEG, but other angles will be added 360 DEG such that their value exceeds 360 DEG, e.g.>If->Then the adjustment to 361 eliminates the period collapse break point and the final average result of 360 is correct.
It should be noted that two angles differing by 180 ° are on average of no practical significance, in order toFor example, the average direction may be 90 ° or 270 °, with or without such angle adjustment, the average value obtained is still 9One of 0 ° or 270 °.
If 0-2 pi is adopted to represent the angle, the adjustment process is as follows: to be used forFor reference, pair->Make adjustments one by one, ifThen->Plus 2π, otherwise if->Then->Minus 2 pi, wherein->For reference angle, +.>For the angle to be adjusted, n= 2~N;
FIG. 4 is a complex multiplication squaring programming flow chart provided by an embodiment of the present invention, which is used for multiplying and squaring complex numbers represented by two components and outputting complex numbers represented by the components, wherein the complex numbers [ a, b ] represented by the two components are input first]And [ c, d ]]Then [ a, b ]]Converted to polar coordinate format [ r, θ ]]Wherein r represents a mode, and θ represents an angle range of 0 to 2pi; [ c, d]Converted into polar format [ s, beta ]]Wherein s represents a mode, and beta represents an angle range of 0-2 pi; further, the square value of the product of the two modes is obtained according to the following formula (9)A die that is the final result; further, the angle is adjusted, x=θ+β is the sum of the two angles, y=θ - β is the difference between the two angles, if y<-pi then x=x+2pi corresponds to θ+2pi in x, e.gFruit y>Pi then x=x-2 pi corresponds to θ -2 pi in x, thereby effecting adjustment of θ, note: here, θ may be adjusted based on β, or β may be adjusted based on θ, with the result unchanged; last step->The conversion of the modulus and angle of the polar representation into component format, in particular x/2 represents the algebraic average of the two angles.
The program written according to the flow chart is successfully applied to a high-precision digital alternating current bridge realized by adopting a channel exchange method principle, and high-precision phase difference measurement within the range of 0-360 degrees is realized.
The adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted further comprises various digital representation angles by taking one angle as a reference:
if 0-65536 is adopted to represent 0-2pi angle, the adjustment process is as follows: to be used forFor reference, pair->Make adjustments one by one, if<-32768, then->Plus 65536 if otherwise->Then->The 65536 is subtracted; />As the reference angle of the beam,for the angle to be adjusted, n= 2~N.
Fig. 5 is a schematic diagram of a phase difference filtering method provided by the embodiment of the invention, the phase difference is also called a phase angle, when the phase angle is measured in real time, the phase angle obtained by measuring each time is added into a circulation queue, the circulation queue only retains the latest 8 phase angle data, any one phase angle is used as a reference, the other 7 phase angles are adjusted, and then the algebraic average value or weighted average value of the 8 phase angles is taken, so that the phase difference filtering can be realized; in fig. 5, θ1 to θ8 are the last 8 phase angles stored in time sequence, the phase angles are represented by 16-bit unsigned fixed point numbers, 0 to 65535 represent 0 to 359.99 degrees, and the specific values of the phase angles are different, so that the phase angles are represented by 8 different phase angles in fig. 5, and the average phase angle is adjusted according to the rule and then replaced by the number average value; for simplicity, only two phase angles θ1, θ8 are taken as examples for further explanation: assuming that θ1=8192 (45 °), θ8= 62805 (345 °), and adjusting θ8 with θ1 as a reference to obtain θ8=62805-65536= -2731, the average value of both is 2730 (15 °), and the result is correct; the phase angle can be expanded into 32-bit signed numbers during specific calculation so as to realize quick operation; the phase angle measurement accuracy can be greatly improved after averaging.
Step 4: and performing evolution operation on the product of the plurality of complex numbers after adjustment, and multiplying the evolution result to obtain data with phase errors eliminated.
The square operation is to open the product of N complex numbers to the power of N, the modulus is to open the product of N complex numbers to the power of N, and the angle is the average value of N angles.
The N complex products are found to be the power N:
（9）
wherein N is the number of plural numbers,is the corresponding mode of the nth complex number, +.>Is the nth complexThe number corresponds to the angle.
After the N angles are adjusted, algebraic averaging can be adopted to obtain the angle after the squaring, and angle information is not lost in a four-quadrant range.
After the N angles are adjusted, the N angles can be sequenced according to the size and weighted and averaged, so that the angle filtering function is realized, and the phase filtering function is realized from the digital signal processing perspective.
The N angles may be represented in different units of measure or numbers.
In this embodiment, the sum of angles is adjusted according to the difference between the two angles, divided by 2 to average, and finally converted into a complex number of component expressions.
Example two
The present embodiment provides a signal processing system based on complex evolution operation, including:
a signal acquisition module configured to: acquiring an alternating current or related digital signal to be processed;
a complex representation module configured to: representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers; wherein each complex number comprises two arrays of double-precision numbers, the first number is a module, and the second number is an angle;
an angle adjustment module configured to: taking one angle as a reference, and adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted;
a data reduction module configured to: performing squaring operation on the product of the plurality of complex numbers subjected to angle adjustment to obtain a processed signal; the square operation comprises square operation of the product of the die and angle averaging.
Example III
The present embodiment provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps in a signal processing method based on complex evolution operation as described in embodiment one.
Example IV
The present embodiment provides a computer device including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the steps in a signal processing method based on complex square operations as described in embodiment one when executing the program.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of a hardware embodiment, a software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, magnetic disk storage, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Those skilled in the art will appreciate that implementing all or part of the above-described methods in accordance with the embodiments may be accomplished by way of a computer program stored on a computer readable storage medium, which when executed may comprise the steps of the embodiments of the methods described above. The storage medium may be a magnetic disk, an optical disk, a Read-Only Memory (ROM), a Random access Memory (Random AccessMemory, RAM), or the like.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (9)
1. The signal processing method based on complex evolution operation is characterized by comprising the following steps:
acquiring an alternating current or related digital signal to be processed;
representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers; wherein each complex number comprises two arrays of double-precision numbers, the first number is a module, and the second number is an angle;
taking one angle as a reference, and adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted;
the adjusting the angle to be adjusted based on the difference value between the reference angle and the angle to be adjusted includes:
if 0-2 pi is adopted to represent the angle, the adjusting process is as follows: to be used forFor reference, pair->Make adjustments one by one, if->ThenPlus 2π, otherwise if->Then->Minus 2 pi, wherein->For reference angle, +.>For the angle to be adjusted, n= 2~N;
if 0-360 degrees are adopted to represent 0-2 pi angles, the adjustment process is as follows: to be used forFor reference, pair->Make adjustments one by one, ifThen->Plus 360 DEG, otherwise if%>Then->Minus 360 DEG, wherein%>For reference angle, +.>For the angle to be adjusted, n= 2~N;
if 0-65536 is adopted to represent 0-2pi angle, the adjustment process is as follows: to be used forFor reference, pair->Make adjustments one by one, ifThen->Plus 65536 if otherwise->Then->The 65536 is subtracted; wherein (1)>For reference angle, +.>For the angle to be adjusted, n= 2~N;
and performing squaring operation on the product of the plurality of complex numbers subjected to angle adjustment to obtain a processed signal, wherein the squaring operation comprises squaring the product of the die and averaging the angle.
2. The signal processing method according to claim 1, wherein the plurality of angles are adjusted, and the plurality of angles are ordered in order of magnitude and weighted-averaged.
3. The signal processing method according to claim 1, wherein the squaring operation is a product of N complex numbers, where a modulus is a product of N complex numbers modulo N times, and an angle is an average of N angles.
4. The signal processing method according to claim 1, wherein the expression of the evolution operation is:
wherein N is the number of a plurality of numbers,is the corresponding mode of the nth complex number, +.>Is the angle corresponding to the nth complex number.
5. A signal processing method based on complex evolution operation according to claim 1, characterized in that the ac or related digital signal to be processed is obtained by means of channel switching.
6. A signal processing method based on complex evolution operation according to claim 1, characterized in that the phasors of the alternating or related digital signal are represented as products of a plurality of complex numbers using polar representation.
7. A signal processing system based on complex evolution, comprising:
a signal acquisition module configured to: acquiring an alternating current or related digital signal to be processed;
a complex representation module configured to: representing the phasors of the alternating or related digital signal as the product of a plurality of complex numbers; wherein each complex number comprises two arrays of double-precision numbers, the first number is a module, and the second number is an angle;
an angle adjustment module configured to: taking one angle as a reference, and adjusting the angle to be adjusted based on the difference value of the reference angle and the angle to be adjusted;
the adjusting the angle to be adjusted based on the difference value between the reference angle and the angle to be adjusted includes:
if 0-2 pi is adopted to represent the angle, the adjusting process is as follows: to be used forFor reference, pair->Make adjustments one by one, if->ThenPlus 2π, otherwise if->Then->Minus 2 pi, wherein->For reference angle, +.>For the angle to be adjusted, n= 2~N;
if 0 to 360 DEG is used to represent 0 to 2 pi,the adjusting process is as follows: to be used forFor reference, pair->Make adjustments one by one, ifThen->Plus 360 DEG, otherwise if%>Then->Minus 360 DEG, wherein%>For reference angle, +.>For the angle to be adjusted, n= 2~N;
if 0-65536 is adopted to represent 0-2pi angle, the adjustment process is as follows: to be used forFor reference, pair->Make adjustments one by one, ifThen->Plus 65536 if otherwise->Then->The 65536 is subtracted; wherein (1)>For reference angle, +.>For the angle to be adjusted, n= 2~N;
a data reduction module configured to: performing squaring operation on the product of the plurality of complex numbers subjected to angle adjustment to obtain a processed signal; the square operation comprises square operation of the product of the die and angle averaging.
8. A computer-readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the steps of a signal processing method based on complex evolution operation as claimed in any one of claims 1-6.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of a signal processing method based on complex evolution operations according to any one of claims 1-6 when the program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310720297.8A CN116595300B (en) | 2023-06-19 | 2023-06-19 | Signal processing method, system, medium and equipment based on complex evolution operation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310720297.8A CN116595300B (en) | 2023-06-19 | 2023-06-19 | Signal processing method, system, medium and equipment based on complex evolution operation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116595300A CN116595300A (en) | 2023-08-15 |
| CN116595300B true CN116595300B (en) | 2023-10-13 |
Family
ID=87599238
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310720297.8A Active CN116595300B (en) | 2023-06-19 | 2023-06-19 | Signal processing method, system, medium and equipment based on complex evolution operation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116595300B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017219463A (en) * | 2016-06-09 | 2017-12-14 | 三菱電機株式会社 | Signal processing apparatus and rotational frequency measuring apparatus |
| CN108196822A (en) * | 2017-12-24 | 2018-06-22 | 北京卫星信息工程研究所 | A kind of method and system of double-precision floating point extracting operation |
| CN110209373A (en) * | 2018-02-28 | 2019-09-06 | 成都鼎桥通信技术有限公司 | Realize the method and device of complex multiplication |
| JP2020076702A (en) * | 2018-11-09 | 2020-05-21 | 三菱電機株式会社 | Ac electric quantity measuring device, ac electric quantity measuring method and bus transfer system |
| CN111475773A (en) * | 2020-03-30 | 2020-07-31 | 南京国电南自电网自动化有限公司 | Rational fitting-based complex modulus and phase angle calculation method, device and system |
| CN111929491A (en) * | 2020-08-12 | 2020-11-13 | 山东大学 | Transmission line sampling synchronization method and system based on Taylor root-like power theorem |
| CN114301552A (en) * | 2022-01-06 | 2022-04-08 | 中电科思仪科技股份有限公司 | Digital modulation signal testing method and system |
| CN116011376A (en) * | 2022-10-26 | 2023-04-25 | 华南理工大学 | Modeling method and model of modularized multi-level matrix type converter based on dynamic phasor principle |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9874640B2 (en) * | 2014-08-15 | 2018-01-23 | Spirent Communications, Plc | Automated methods for suppression of spurious signals |
| CN108897716B (en) * | 2018-07-04 | 2022-07-01 | 合肥工业大学 | A data processing device and method for reducing the amount of calculation through memory read and write operations |
-
2023
- 2023-06-19 CN CN202310720297.8A patent/CN116595300B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2017219463A (en) * | 2016-06-09 | 2017-12-14 | 三菱電機株式会社 | Signal processing apparatus and rotational frequency measuring apparatus |
| CN108196822A (en) * | 2017-12-24 | 2018-06-22 | 北京卫星信息工程研究所 | A kind of method and system of double-precision floating point extracting operation |
| CN110209373A (en) * | 2018-02-28 | 2019-09-06 | 成都鼎桥通信技术有限公司 | Realize the method and device of complex multiplication |
| JP2020076702A (en) * | 2018-11-09 | 2020-05-21 | 三菱電機株式会社 | Ac electric quantity measuring device, ac electric quantity measuring method and bus transfer system |
| CN111475773A (en) * | 2020-03-30 | 2020-07-31 | 南京国电南自电网自动化有限公司 | Rational fitting-based complex modulus and phase angle calculation method, device and system |
| CN111929491A (en) * | 2020-08-12 | 2020-11-13 | 山东大学 | Transmission line sampling synchronization method and system based on Taylor root-like power theorem |
| CN114301552A (en) * | 2022-01-06 | 2022-04-08 | 中电科思仪科技股份有限公司 | Digital modulation signal testing method and system |
| CN116011376A (en) * | 2022-10-26 | 2023-04-25 | 华南理工大学 | Modeling method and model of modularized multi-level matrix type converter based on dynamic phasor principle |
Non-Patent Citations (2)
| Title |
|---|
| Better phasor calibration in turbid media using Monte Carlo simulations;Smith JT 等;《Biophysical Journal 》;第122卷(第03期);全文 * |
| 高速高精度CORDIC处理器的ASIC实现;汪润来 等;《军民两用技术与产品》(第07期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116595300A (en) | 2023-08-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yu et al. | A robust and fixed-time zeroing neural dynamics for computing time-variant nonlinear equation using a novel nonlinear activation function | |
| CN111656315A (en) | Data processing method and device based on convolutional neural network architecture | |
| Yu et al. | Small limit cycles bifurcating from fine focus points in cubic order Z2-equivariant vector fields | |
| CN116595300B (en) | Signal processing method, system, medium and equipment based on complex evolution operation | |
| CN115469259A (en) | RBF neural network-based CT error state online quantitative evaluation method and device | |
| CN114065114A (en) | Method and system for predicting metering error of capacitive voltage transformer | |
| CN114325164B (en) | Multi-fault diagnosis method for single-phase three-level rectifier | |
| CN113484818B (en) | Sliding window-based accurate positioning method for high-frequency acquisition anomaly-resistant electric energy meter | |
| CN103884910B (en) | A kind of power system phasor calculating method being applicable to frequency shift | |
| CN117040339B (en) | Vector decoupling method, device and equipment for permanent magnet synchronous motor and storage medium | |
| CN108965866B (en) | Method and system for acquiring focusing curved surface of camera movement lens | |
| CN113472228B (en) | Three-level midpoint potential balance control method and system based on voltage feedback | |
| Olver | A new method for the evaluation of zeros of Bessel functions and of other solutions of second-order differential equations | |
| CN114024320B (en) | Safety and stability control device power and phase angle measurement method, system and storage medium | |
| CN116413652A (en) | Current transformer calibration system and method | |
| JP3814834B2 (en) | Harmonic detection method by fixed sampling | |
| CN115754593A (en) | Cable fault positioning method and device and electronic equipment | |
| CN115078829A (en) | Current and voltage frequency detection method and detection system for power system | |
| CN115296574A (en) | Induction motor parameter identification method and device based on steady-state T-shaped equivalent circuit | |
| CN111950123A (en) | Gyroscope error coefficient curve fitting prediction method and system | |
| CN109327026B (en) | An Interval State Estimation Method for Low-Voltage Distribution Networks | |
| CN113485508A (en) | Voltage reference regulating circuit, processing circuit and servo driver | |
| CN105136060B (en) | The method for determining phase-shift phase main value in measuring three-dimensional morphology by plot of light intensity | |
| CN111475773A (en) | Rational fitting-based complex modulus and phase angle calculation method, device and system | |
| CN118826074B (en) | A three-phase electrical frequency estimation method, device and computer-readable storage medium |
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 |