CN109991844A - A kind of d-q Decoupling Controller Design method using embedded decoupling synchronous reference coordinate transform - Google Patents

A kind of d-q Decoupling Controller Design method using embedded decoupling synchronous reference coordinate transform Download PDF

Info

Publication number
CN109991844A
CN109991844A CN201910325121.6A CN201910325121A CN109991844A CN 109991844 A CN109991844 A CN 109991844A CN 201910325121 A CN201910325121 A CN 201910325121A CN 109991844 A CN109991844 A CN 109991844A
Authority
CN
China
Prior art keywords
coordinate system
transformation
axis component
outer layer
component
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910325121.6A
Other languages
Chinese (zh)
Other versions
CN109991844B (en
Inventor
邵振国
杨少华
郑文迪
陈飞雄
周琪琪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN201910325121.6A priority Critical patent/CN109991844B/en
Publication of CN109991844A publication Critical patent/CN109991844A/en
Application granted granted Critical
Publication of CN109991844B publication Critical patent/CN109991844B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0205Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/024Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Inverter Devices (AREA)

Abstract

The present invention relates to a kind of d-q Decoupling Controller Design methods using embedded decoupling synchronous reference coordinate transform, comprising the following steps: step S1: extracts the of ac under three-phase static coordinate system;Step S2: outer layer synchronous reference coordinate system transformation is carried out;Step S3: direct current aim parameter under d-q reference frame and harmonic wave interference amount are subjected to phase shift to three phasors;Step S4: embedded synchronous reference coordinate system transformation is carried out;Step S5: decoupling direct current aim parameter and harmonic wave interference amount, obtains direct current aim parameter;Step S6: judging whether ripple content meets required precision, if so, enter step S7, otherwise return step S3;Step S7: will be in direct current aim parameter input value d-q decoupling controller.The present invention can be in the case where having powerful connections harmonic wave, accurate easy realization synchronous reference coordinate system transformation.

Description

It is a kind of to be set using the embedded d-q decoupling controller for decoupling synchronous reference coordinate transform Meter method
Technical field
The present invention relates to field of power systems to mention especially suitable for needing the controller design field using coordinate transform A kind of d-q Decoupling Controller Design method using embedded decoupling synchronous reference coordinate transform is gone out.
Background technique
D-q decoupling controller based on synchronous reference coordinate system transformation is very widely used in the power system, such as Phaselocked loop (Phase-Locked Loop, PLL).Due to the presence of background harmonics, the result of synchronous reference coordinate system transformation can be made It is mingled with ripple component, this causes very big influence to subsequent d-q decoupling controller.It is accurate to realize that synchronous reference coordinate system becomes It changes, the d- axis, q- axis component after making decoupling uniquely correspond to positive sequence fundametal compoment, this is the weight of d-q Decoupling Controller Design, operation It will foundation.
The synchronous reference coordinate transformation method that traditional d-q decoupling controller uses mainly has single synchronous reference coordinate system method (Single Synchronous Reference Frame, SSRF) and decouple double synchronous rotating frame method (Decoupled Double Synchronous Reference Frame,DDSRF).Both methods can not all be completely free of background harmonics pair Coordinate transform bring adverse effect.SSRF can make background harmonics according to the progress of disease rule of Park Transformation to d- axis and q- axis component Transmitting.DDSRF has carried out simple processing to background harmonics using low-pass filter (Low Pass Filter, LPF), still Again without avoiding background harmonics from being transmitted from three-phase static coordinate system to transformed synchronous reference coordinate system.
SSRF is by electric signal using Park Transformation (abc-to-dq0) from static three-phase coordinate system transformation to rotation In synchronous coordinate system, Principle representation formula is as follows:
udq=T3s-dq(θ)uabc
Wherein,
It can be obtained from expression parsing: after background harmonics can be by the harmonic wave progress of disease rule feed-in coordinate transform of Park Transformation D- axis and q- axis.SSRF can not be handled background harmonics, therefore there are wretched insufficiencies for the d-q decoupling controller based on SSRF.
The basic principle of DDSRF technology is to eliminate negative phase-sequence fundametal compoment with double synchronous rotating angles, uses low-pass filtering The method of device removes harmonic component, to achieve the purpose that output signal does not include negative phase-sequence fundametal compoment completely.
To offset negative phase-sequence fundametal compoment, this method is to usabcDq is obtained by forward and reverse synchronous rotating angle to sit Mark system and d-1q-1Component in coordinate systemWithThe two components can be denoted as respectively as positive sequence fundametal compoment, negative phase-sequence base The superposition of wave component and harmonic component, is expressed asWithWhereinIt is right Positive sequence fundametal compoment is answered,Corresponding negative phase-sequence fundametal compoment,Corresponding harmonic component.
Wherein,WithThere are certain relationships:
It can be designed the synchronous reference coordinate transformation method DDSRF that can offset negative phase-sequence fundametal compoment according to relation above. Its functional block diagram is as shown in Figure 1.
By analyzing above: DDSRF technology has very strong superiority compared to SRF, can be completely counterbalanced by negative phase-sequence base Wave component.Suitable for power grid, there are the non-equilibrium states of negative phase-sequence fundametal compoment and extremely a small amount of harmonic component.However DDSRF does not have It solves the problems, such as harmonic wave, is filtered simply by LPF.
To sum up, the d-q decoupling controller based on DDSRF technology has the following disadvantages: for different background harmonics Situation needs constantly to carry out parameter adjustment;The change for being likely to result in poles and zeros assignment is combined with subsequent controllers by LPF, Controller stability is in turn resulted in go wrong;It cannot achieve when background harmonics pollution condition complexity and filter out completely.
Summary of the invention
In view of this, the purpose of the present invention is to propose to a kind of d-q solutions using embedded decoupling synchronous reference coordinate transform Coupling controller design method makes d-q decoupling controller be suitable for the situation of background harmonics complexity, can be in the feelings of harmonic wave that have powerful connections Under condition, accurate easy realization synchronous reference coordinate system transformation.
The present invention is realized using following scheme: a kind of d-q decoupling control using embedded decoupling synchronous reference coordinate transform Device design method processed, specifically includes the following steps:
Step S1: the of ac under three-phase static coordinate system is extracted;
Step S2: outer layer synchronous reference coordinate system transformation is carried out;
Step S3: direct current aim parameter under d-q reference frame and harmonic wave interference amount are subjected to phase shift to three phasors;
Step S4: embedded synchronous reference coordinate system transformation is carried out;
Step S5: decoupling direct current aim parameter and harmonic wave interference amount, obtains direct current aim parameter;
Step S6: judging whether ripple content meets required precision, if so, enter step S7, otherwise return step S3;
Step S7: will be in direct current aim parameter input value d-q decoupling controller.
Further, step S3 specifically includes the following steps:
Step S31: three phase harmonic is compared is defined as:
To ua(t)=h (ω t+ θ) function is translated, and translates in different directions t=1/3f respectivelyn, wherein fnFundamental wave Frequency is 50Hz, obtains following formula:
Step S32: two formula in simultaneous step S31 obtain:
Further, step S4 specifically: embedded synchronous reference is carried out to the three-phase electrical quantity containing harmonic wave using following formula Coordinate system transformation:
udq0_new=T3s-dq(θ)u3s_new
In formula, udq0_newIt indicates through the transformed d-q reference frame component of internal layer, T3s-dq(θ) indicates Park Transformation square Battle array, u3s_newIndicate that the result of outer layer transformation carries out the three phase static reference frame component after phase shift.
Further, step S5 specifically includes the following steps:
Step S51: the d axis of internal layer rotating coordinate system, q axis component correspond to the ripple interference of outer layer rotating coordinate system Amount, such as following formula:
In formula,Indicate that outer layer transformation d- axis component carries out the d- in the new rotating coordinate system obtained after internal layer transformation The DC component of axis component,Indicate that outer layer transformation d- axis component carries out in the new rotating coordinate system obtained after internal layer transformation D- axis component positive sequence fundametal compoment,Indicate that outer layer transformation d- axis component carries out the new rotation obtained after internal layer transformation The AC compounent of d- axis component in coordinate system,It is obtained after expression outer layer transformation d- axis component progress internal layer transformation new The positive sequence harmonic component of d- axis component in rotating coordinate system,After indicating that outer layer transformation d- axis component carries out internal layer transformation The Negative sequence harmonic component of the obtained d- axis component in new rotating coordinate system;In indicating that outer layer transformation d- axis component carries out The DC component of the q- axis component in new rotating coordinate system obtained after layer transformation,Indicate outer layer convert d- axis component into The positive sequence fundametal compoment of the q- axis component in new rotating coordinate system obtained after the transformation of row internal layer,Indicate that outer layer converts d- Axis component carries out the AC compounent of the q- axis component in the new rotating coordinate system obtained after internal layer transformation,Indicate that outer layer becomes The positive sequence harmonic component for the q- axis component that d- axis component carries out in the new rotating coordinate system obtained after internal layer transformation is changed,Table Show that outer layer transformation d- axis component carries out the negative phase-sequence fundametal compoment of the q- axis component in the new rotating coordinate system obtained after internal layer transformation;Indicate that outer layer transformation q- axis component carries out the direct current point of the d- axis component in the new rotating coordinate system obtained after internal layer transformation Amount,Indicate that outer layer transformation q- axis component is carrying out the d- axis component in the new rotating coordinate system obtained after internal layer transformation just Sequence fundametal compoment,Indicate that outer layer transformation q- axis component carries out the d- axis in the new rotating coordinate system obtained after internal layer transformation The AC compounent of component,Indicate that outer layer transformation q- axis component carries out in the new rotating coordinate system obtained after internal layer transformation The positive sequence harmonic component of d- axis component,Indicate that outer layer transformation q- axis component carries out the new rotation obtained after internal layer transformation and sits The Negative sequence harmonic component of d- axis component in mark system;It is obtained after expression outer layer transformation q- axis component progress internal layer transformation The DC component of q- axis component in new rotating coordinate system,It indicates to obtain after outer layer transformation q- axis component carries out internal layer transformation To new rotating coordinate system in q- axis component positive sequence fundametal compoment,Indicate that outer layer transformation q- axis component carries out internal layer The AC compounent of the q- axis component in new rotating coordinate system obtained after transformation,Indicate that outer layer transformation q- axis component carries out The positive sequence harmonic component of the q- axis component in new rotating coordinate system obtained after internal layer transformation,Indicate that outer layer converts q- axis Component carries out the Negative sequence harmonic component of the q- axis component in the new rotating coordinate system obtained after internal layer transformation;
Step S52:0 axis component corresponds to the direct current aim parameter of outer layer rotating coordinate system, such as following formula:
In formula,Indicate that outer layer transformation d- axis component carries out the 0- in the new rotating coordinate system obtained after internal layer transformation The DC component of axis component,Indicate that the of ac under three-phase static coordinate system is obtained by outer layer synchronous reference coordinate system transformation D- axis direct current aim parameter,Indicate that outer layer transformation q- axis component carries out in the new rotating coordinate system obtained after internal layer transformation 0- axis component DC component,Indicate that the of ac under three-phase static coordinate system passes through outer layer synchronous reference coordinate system transformation Obtained q- axis direct current aim parameter;
Step S53: being divided into one kind according to d axis and q axis for the result of embedded coordinate transform, and 0 axis is divided into another kind of, is come with this Realize the decoupling of aim parameter and interference volume.
Method proposed by the present invention is a kind of synchronous reference coordinate system transformation method of multilayer.The electrical quantity for the harmonic wave that has powerful connections After undergoing first time synchronous reference coordinate system transformation (outer layer synchronous reference coordinate system transformation), obtained d axis and q axis component Understand the reason due to background harmonics and includes ripple quantity.Accurate in order to obtain and unique rotational coordinates corresponding with positive sequence fundamental wave It is component, needs to decouple DC quantity and ripple quantity.By phase shift, the direct current tolerance containing ripple is integrated into Three-phase electrical quantity containing harmonic wave then carries out second of synchronous reference coordinate system transformation to the three-phase electrical quantity containing harmonic wave. The result of internal layer coordinate transform is divided into one kind according to d axis and q axis, 0 axis is divided into another kind of, realizes aim parameter and interference with this The decoupling of amount.If the ripple content in aim parameter is unsatisfactory for required precision, after carrying out phase shift, continue embedded synchronization Reference coordinate transformation, further decouples aim parameter and interference volume.
Compared with prior art, the invention has the following beneficial effects: the present invention can be quasi- in the case where having powerful connections harmonic wave Really easy realization synchronous reference coordinate system transformation.The present invention can substantially reduce background harmonics to d- and q- axis DC component not Good transmitting.Huge advantage has been compared compared to using the d-q decoupling controller of SRF with DDSRF technology.
Detailed description of the invention
Fig. 1 is DDSRF functional block diagram in the background technique of the embodiment of the present invention.
Fig. 2 compares for the ripple amplitude of DC component under three kinds of coordinate transformating technologies of the embodiment of the present invention.Wherein (a) is The DC ripple amplitude of corresponding 2 background harmonics;It (b) is the DC ripple amplitude of corresponding 5 background harmonics;It (c) is 7 times corresponding The DC ripple amplitude of background harmonics.
Fig. 3 is the flow diagram of the embodiment of the present invention.
Specific embodiment
The present invention will be further described with reference to the accompanying drawings and embodiments.
It is noted that described further below be all exemplary, it is intended to provide further instruction to the application.Unless another It indicates, all technical and scientific terms used herein has usual with the application person of an ordinary skill in the technical field The identical meanings of understanding.
It should be noted that term used herein above is merely to describe specific embodiment, and be not intended to restricted root According to the illustrative embodiments of the application.As used herein, unless the context clearly indicates otherwise, otherwise singular Also it is intended to include plural form, additionally, it should be understood that, when in the present specification using term "comprising" and/or " packet Include " when, indicate existing characteristics, step, operation, device, component and/or their combination.
As shown in figure 3, present embodiments providing a kind of d-q decoupling control using embedded decoupling synchronous reference coordinate transform Device design method processed, specifically includes the following steps:
Step S1: the of ac under three-phase static coordinate system is extracted;
Step S2: outer layer synchronous reference coordinate system transformation is carried out;
Step S3: direct current aim parameter under d-q reference frame and harmonic wave interference amount are subjected to phase shift to three phasors;
Step S4: embedded synchronous reference coordinate system transformation is carried out;
Step S5: decoupling direct current aim parameter and harmonic wave interference amount, obtains direct current aim parameter;
Step S6: judging whether ripple content meets required precision, if so, enter step S7, otherwise return step S3;
Step S7: will be in direct current aim parameter input value d-q decoupling controller.
In the present embodiment, step S3 specifically includes the following steps:
Step S31: three phase harmonic is compared is defined as:
To ua(t)=h (ω t+ θ) function is translated, and translates in different directions t=1/3f respectivelyn, wherein fnFundamental wave Frequency is 50Hz, obtains following formula:
Step S32: two formula in simultaneous step S31 obtain:
In the present embodiment, step S4 specifically: embedded synchronization is carried out to the three-phase electrical quantity containing harmonic wave using following formula Reference frame transformation:
udq0_new=T3s-dq(θ)u3s_new
In formula, udq0_newIt indicates through the transformed d-q reference frame component of internal layer, T3s-dq(θ) indicates Park Transformation square Battle array, u3s_newIndicate that the result of outer layer transformation carries out the three phase static reference frame component after phase shift.
In the present embodiment, step S5 specifically includes the following steps:
Step S51: the d axis of internal layer rotating coordinate system, q axis component correspond to the ripple interference of outer layer rotating coordinate system Amount, such as following formula:
In formula,Indicate that outer layer transformation d- axis component carries out the d- in the new rotating coordinate system obtained after internal layer transformation The DC component of axis component,Indicate that outer layer transformation d- axis component carries out in the new rotating coordinate system obtained after internal layer transformation D- axis component positive sequence fundametal compoment,Indicate that outer layer transformation d- axis component carries out the new rotation obtained after internal layer transformation The AC compounent of d- axis component in coordinate system,It is obtained after expression outer layer transformation d- axis component progress internal layer transformation new The positive sequence harmonic component of d- axis component in rotating coordinate system,After indicating that outer layer transformation d- axis component carries out internal layer transformation The Negative sequence harmonic component of the obtained d- axis component in new rotating coordinate system;In indicating that outer layer transformation d- axis component carries out The DC component of the q- axis component in new rotating coordinate system obtained after layer transformation,Indicate outer layer convert d- axis component into The positive sequence fundametal compoment of the q- axis component in new rotating coordinate system obtained after the transformation of row internal layer,Indicate that outer layer converts d- Axis component carries out the AC compounent of the q- axis component in the new rotating coordinate system obtained after internal layer transformation,Indicate that outer layer becomes The positive sequence harmonic component for the q- axis component that d- axis component carries out in the new rotating coordinate system obtained after internal layer transformation is changed,Table Show that outer layer transformation d- axis component carries out the negative phase-sequence fundametal compoment of the q- axis component in the new rotating coordinate system obtained after internal layer transformation;Indicate that outer layer transformation q- axis component carries out the direct current point of the d- axis component in the new rotating coordinate system obtained after internal layer transformation Amount,Indicate that outer layer transformation q- axis component is carrying out the d- axis component in the new rotating coordinate system obtained after internal layer transformation just Sequence fundametal compoment,Indicate that outer layer transformation q- axis component carries out the d- axis in the new rotating coordinate system obtained after internal layer transformation The AC compounent of component,Indicate that outer layer transformation q- axis component carries out in the new rotating coordinate system obtained after internal layer transformation The positive sequence harmonic component of d- axis component,Indicate that outer layer transformation q- axis component carries out the new rotation obtained after internal layer transformation and sits The Negative sequence harmonic component of d- axis component in mark system;It is obtained after expression outer layer transformation q- axis component progress internal layer transformation The DC component of q- axis component in new rotating coordinate system,It indicates to obtain after outer layer transformation q- axis component carries out internal layer transformation To new rotating coordinate system in q- axis component positive sequence fundametal compoment,Indicate that outer layer transformation q- axis component carries out internal layer The AC compounent of the q- axis component in new rotating coordinate system obtained after transformation,Indicate that outer layer transformation q- axis component carries out The positive sequence harmonic component of the q- axis component in new rotating coordinate system obtained after internal layer transformation,Indicate that outer layer converts q- axis Component carries out the Negative sequence harmonic component of the q- axis component in the new rotating coordinate system obtained after internal layer transformation;
Step S52:0 axis component corresponds to the direct current aim parameter of outer layer rotating coordinate system, such as following formula:
In formula,Indicate that outer layer transformation d- axis component carries out the 0- in the new rotating coordinate system obtained after internal layer transformation The DC component of axis component,Indicate that the of ac under three-phase static coordinate system is obtained by outer layer synchronous reference coordinate system transformation D- axis direct current aim parameter,Indicate that outer layer transformation q- axis component carries out in the new rotating coordinate system obtained after internal layer transformation 0- axis component DC component,Indicate that the of ac under three-phase static coordinate system passes through outer layer synchronous reference coordinate system transformation Obtained q- axis direct current aim parameter;
Step S53: being divided into one kind according to d axis and q axis for the result of embedded coordinate transform, and 0 axis is divided into another kind of, is come with this Realize the decoupling of aim parameter and interference volume.
What the present embodiment proposed influences decoupling d, q axis component after coordinate transform by background harmonics completely.For The effect for embodying the present embodiment, by the background harmonics of three-phase static coordinate system sets 2 subharmonic of negative phase-sequence, 5 subharmonic of negative phase-sequence With 7 subharmonic of positive sequence.Their amplitude is the 10% of fundamental wave.
The signal containing background harmonics is coordinately transformed with tri- kinds of technologies of SRF, DDSRF, NDSRF respectively again, is obtained The ripple amplitude (per unit value) of DC component after transformation, as shown in Figure 2.
By observe Fig. 2 it is found that using the d-q decoupling controller of NDSRF technology can substantially reduce background harmonics to d- and The bad transmitting of q- axis DC component.Huge advantage has been compared compared to using the d-q decoupling controller of SRF with DDSRF technology.
The method that the present embodiment proposes is a kind of synchronous reference coordinate system transformation method of multilayer.The harmonic wave that has powerful connections it is electrical Amount is after undergoing first time synchronous reference coordinate system transformation (outer layer synchronous reference coordinate system transformation), obtained d axis and q axis point Amount can due to background harmonics reason and include ripple quantity.Accurate in order to obtain and unique rotation seat corresponding with positive sequence fundamental wave Mark system component, needs to decouple DC quantity and ripple quantity.By phase shift, the direct current tolerance containing ripple is integrated At the three-phase electrical quantity containing harmonic wave, second of synchronous reference coordinate system then is carried out to the three-phase electrical quantity containing harmonic wave and is become It changes.The result of internal layer coordinate transform is divided into one kind according to d axis and q axis, 0 axis be divided into it is another kind of, with this come realize aim parameter and The decoupling of interference volume.If the ripple content in aim parameter is unsatisfactory for required precision, after carrying out phase shift, continue to embed Synchronous reference coordinate transform further decouples aim parameter and interference volume.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more, The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates, Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one The step of function of being specified in a box or multiple boxes.
The above described is only a preferred embodiment of the present invention, being not that the invention has other forms of limitations, appoint What those skilled in the art changed or be modified as possibly also with the technology contents of the disclosure above equivalent variations etc. Imitate embodiment.But without departing from the technical solutions of the present invention, according to the technical essence of the invention to above embodiments institute Any simple modification, equivalent variations and the remodeling made, still fall within the protection scope of technical solution of the present invention.

Claims (4)

1. a kind of d-q Decoupling Controller Design method using embedded decoupling synchronous reference coordinate transform, it is characterised in that: packet Include following steps:
Step S1: the of ac under three-phase static coordinate system is extracted;
Step S2: outer layer synchronous reference coordinate system transformation is carried out;
Step S3: direct current aim parameter under d-q reference frame and harmonic wave interference amount are subjected to phase shift to three phasors;
Step S4: embedded synchronous reference coordinate system transformation is carried out;
Step S5: decoupling direct current aim parameter and harmonic wave interference amount, obtains direct current aim parameter;
Step S6: judging whether ripple content meets required precision, if so, enter step S7, otherwise return step S3;
Step S7: will be in direct current aim parameter input value d-q decoupling controller.
2. a kind of d-q decoupling controller using embedded decoupling synchronous reference coordinate transform according to claim 1 is set Meter method, it is characterised in that: step S3 specifically includes the following steps:
Step S31: three phase harmonic is compared is defined as:
To ua(t)=h (ω t+ θ) function is translated, and translates in different directions t=1/3f respectivelyn, wherein fnFundamental frequency For 50Hz, following formula is obtained:
Step S32: two formula in simultaneous step S31 obtain:
3. a kind of d-q decoupling controller using embedded decoupling synchronous reference coordinate transform according to claim 1 is set Meter method, it is characterised in that: step S4 specifically: embedded synchronous reference is carried out to the three-phase electrical quantity containing harmonic wave using following formula Coordinate system transformation:
udq0_new=T3s-dq(θ)u3s_new
In formula, udq0_newIt indicates through the transformed d-q reference frame component of internal layer, T3s-dq(θ) indicates Park Transformation matrix, u3s_newIndicate that the result of outer layer transformation carries out the three phase static reference frame component after phase shift.
4. a kind of d-q decoupling controller using embedded decoupling synchronous reference coordinate transform according to claim 1 is set Meter method, it is characterised in that: step S5 specifically includes the following steps:
Step S51: the d axis of internal layer rotating coordinate system, q axis component correspond to the ripple interference volume of outer layer rotating coordinate system, such as Following formula:
In formula,Indicate that outer layer transformation d- axis component carries out the d- axis component in the new rotating coordinate system obtained after internal layer transformation DC component,Indicate that outer layer transformation d- axis component carries out the d- axis in the new rotating coordinate system obtained after internal layer transformation The positive sequence fundametal compoment of component,Indicate that outer layer transformation d- axis component carries out the new rotating coordinate system obtained after internal layer transformation In d- axis component AC compounent,Indicate that outer layer transformation d- axis component carries out the new rotation obtained after internal layer transformation and sits The positive sequence harmonic component of d- axis component in mark system,It is obtained after expression outer layer transformation d- axis component progress internal layer transformation The Negative sequence harmonic component of d- axis component in new rotating coordinate system;Indicate that outer layer transformation d- axis component carries out internal layer transformation The DC component of the q- axis component in new rotating coordinate system obtained afterwards,Indicate that outer layer transformation d- axis component carries out internal layer The positive sequence fundametal compoment of the q- axis component in new rotating coordinate system obtained after transformation,Indicate that outer layer converts d- axis component The AC compounent of the q- axis component in the new rotating coordinate system obtained after internal layer transformation is carried out,Indicate that outer layer converts d- axis Component carries out the positive sequence harmonic component of the q- axis component in the new rotating coordinate system obtained after internal layer transformation,Indicate outer layer Transformation d- axis component carries out the negative phase-sequence fundametal compoment of the q- axis component in the new rotating coordinate system obtained after internal layer transformation;Table Show that outer layer transformation q- axis component carries out the DC component of the d- axis component in the new rotating coordinate system obtained after internal layer transformation,Indicate that outer layer transformation q- axis component carries out the positive sequence base of the d- axis component in the new rotating coordinate system obtained after internal layer transformation Wave component,Indicate that outer layer transformation q- axis component carries out the d- axis component in the new rotating coordinate system obtained after internal layer transformation AC compounent,Indicate that outer layer transformation q- axis component carries out the d- axis in the new rotating coordinate system obtained after internal layer transformation The positive sequence harmonic component of component,Indicate that outer layer transformation q- axis component carries out the new rotating coordinate system obtained after internal layer transformation In d- axis component Negative sequence harmonic component;Indicate that outer layer transformation q- axis component carries out the new rotation obtained after internal layer transformation Turn the DC component of the q- axis component in coordinate system,It is obtained after expression outer layer transformation q- axis component progress internal layer transformation The positive sequence fundametal compoment of q- axis component in new rotating coordinate system,Indicate that outer layer transformation q- axis component carries out internal layer transformation The AC compounent of the q- axis component in new rotating coordinate system obtained afterwards,Indicate that outer layer transformation q- axis component carries out internal layer The positive sequence harmonic component of the q- axis component in new rotating coordinate system obtained after transformation,Indicate that outer layer converts q- axis component Carry out the Negative sequence harmonic component of the q- axis component in the new rotating coordinate system obtained after internal layer transformation;
Step S52:0 axis component corresponds to the direct current aim parameter of outer layer rotating coordinate system, such as following formula:
In formula,Indicate that outer layer transformation d- axis component carries out the 0- axis component in the new rotating coordinate system obtained after internal layer transformation DC component,Indicate that the of ac under three-phase static coordinate system passes through the d- axis that outer layer synchronous reference coordinate system transformation obtains Direct current aim parameter,Indicate that outer layer transformation q- axis component carries out the 0- axis in the new rotating coordinate system obtained after internal layer transformation The DC component of component,Indicate what the of ac under three-phase static coordinate system was obtained by outer layer synchronous reference coordinate system transformation Q- axis direct current aim parameter;
Step S53: being divided into one kind according to d axis and q axis for the result of embedded coordinate transform, and 0 axis is divided into another kind of, is realized with this The decoupling of aim parameter and interference volume.
CN201910325121.6A 2019-04-22 2019-04-22 Design method of d-q decoupling controller adopting embedded decoupling synchronous reference coordinate transformation Active CN109991844B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910325121.6A CN109991844B (en) 2019-04-22 2019-04-22 Design method of d-q decoupling controller adopting embedded decoupling synchronous reference coordinate transformation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910325121.6A CN109991844B (en) 2019-04-22 2019-04-22 Design method of d-q decoupling controller adopting embedded decoupling synchronous reference coordinate transformation

Publications (2)

Publication Number Publication Date
CN109991844A true CN109991844A (en) 2019-07-09
CN109991844B CN109991844B (en) 2021-08-31

Family

ID=67132923

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910325121.6A Active CN109991844B (en) 2019-04-22 2019-04-22 Design method of d-q decoupling controller adopting embedded decoupling synchronous reference coordinate transformation

Country Status (1)

Country Link
CN (1) CN109991844B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6956352B2 (en) * 2001-03-19 2005-10-18 Abb Oy Vector control of an induction motor
CN101673952A (en) * 2009-08-14 2010-03-17 燕山大学 Precise phase locking method based on cross decoupling self-adaptive complex filter
CN102095915A (en) * 2010-11-24 2011-06-15 重庆大学 Voltage signal detecting device using multiple synchronous reference coordinate system transformation
CN103391095A (en) * 2013-07-18 2013-11-13 深圳市晶福源电子技术有限公司 Three-phase voltage unbalance phase-locked loop based on decoupling control
CN106410858A (en) * 2016-08-30 2017-02-15 成都麦隆电气有限公司 Software digital phase-locking method based on dual dq coordination conversion
US10215784B1 (en) * 2017-12-05 2019-02-26 Industrial Technology Research Institute Measuring apparatus including phase locked loop and measuring method thereof

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6956352B2 (en) * 2001-03-19 2005-10-18 Abb Oy Vector control of an induction motor
CN101673952A (en) * 2009-08-14 2010-03-17 燕山大学 Precise phase locking method based on cross decoupling self-adaptive complex filter
CN102095915A (en) * 2010-11-24 2011-06-15 重庆大学 Voltage signal detecting device using multiple synchronous reference coordinate system transformation
CN103391095A (en) * 2013-07-18 2013-11-13 深圳市晶福源电子技术有限公司 Three-phase voltage unbalance phase-locked loop based on decoupling control
CN106410858A (en) * 2016-08-30 2017-02-15 成都麦隆电气有限公司 Software digital phase-locking method based on dual dq coordination conversion
US10215784B1 (en) * 2017-12-05 2019-02-26 Industrial Technology Research Institute Measuring apparatus including phase locked loop and measuring method thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李珊瑚: "解耦多同步参考坐标系电网电压同步信号检测方法", 《电工技术学报》 *

Also Published As

Publication number Publication date
CN109991844B (en) 2021-08-31

Similar Documents

Publication Publication Date Title
CN108462206B (en) Method and device for determining selectable range of virtual inertia and damping coefficient of VSG
Xia et al. A novel direct torque and flux control method of matrix converter-fed PMSM drives
CN106602895B (en) The detection method and system of high voltage direct current transmission device commutation parameter
CN110417056B (en) Equipment-level oscillation source positioning method and device of doubly-fed fan grid-connected system
CN107834830A (en) The control method and control system that a kind of mixed type MMC runs without interruption
CN106410858A (en) Software digital phase-locking method based on dual dq coordination conversion
CN101771361A (en) Method for controlling direct power of grid-connected inverter without non-AC voltage sensor
CN103267897A (en) Three-phase-locked loop based on reversed Park conversion
CN106655277A (en) Improved permanent magnet synchronous motor phase-locked loop method
Martí et al. Shifted Frequency Analysis (SFA) concepts for EMTP modelling and simulation of Power System Dynamics
WO2022257405A1 (en) Maximum torque per ampere control method and apparatus, terminal device, and storage medium
CN104201943B (en) Determination method and device for number of pole pairs
WO2020220186A1 (en) Phase sequence detection method for three-phase power grid, and related apparatus
CN109991844A (en) A kind of d-q Decoupling Controller Design method using embedded decoupling synchronous reference coordinate transform
CN110596455B (en) Power frequency electrical parameter extraction method, system and computer readable storage medium
Yang et al. Mitigation of background harmonics effect on MMC controller based on a novel coordinate transformation technique
CN115000982A (en) Doubly-fed wind power grid-connected system subsynchronous oscillation inducement determination method
CN107994565A (en) The emulation mode and its system of a kind of THE UPFC
CN107085136A (en) A kind of three-phase alternating voltage falls the determination methods and system of state
CN111106618B (en) Harmonic analysis method and device for access of new energy power generation equipment to power system
CN207967943U (en) A kind of Active Power Filter-APF
CN117614020B (en) Method, device and medium for software phase-locked loop
CN113625066B (en) Distribution transformer phase unbalance detection method, system, device and storage medium
CN108879722A (en) A kind of virtual inertia generation method of virtual synchronous generator
Wei et al. Cooperative Control to Suppress Unbalanced and Harmonic Distortion in the Distribution Network With Inverter-Based Distributed Generators

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