CN102393866A - Transfer function determination method of system object parameter model and device - Google Patents
Transfer function determination method of system object parameter model and device Download PDFInfo
- Publication number
- CN102393866A CN102393866A CN201110282741XA CN201110282741A CN102393866A CN 102393866 A CN102393866 A CN 102393866A CN 201110282741X A CN201110282741X A CN 201110282741XA CN 201110282741 A CN201110282741 A CN 201110282741A CN 102393866 A CN102393866 A CN 102393866A
- Authority
- CN
- China
- Prior art keywords
- inertia
- exponent number
- constant
- final
- unit
- 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
Links
Images
Landscapes
- Feedback Control In General (AREA)
Abstract
The invention discloses a transfer function determination method of a system object parameter model and a device. The method comprises the following steps of: presetting the transfer function of the system object parameter model to comprise n orders of isochoric time-lag links; selecting a plurality of preset identification frequency points, obtaining amplitude value transfer coefficients of the identification frequency points and phase transfer values of the identification frequency points; according to the amplitude transfer coefficients of the identification frequency points, computing inertia order n and inertia constants of the isochoric links; and according to the current inertia order and the current constants, computing the proportional gain and time-lag constant of the transfer function. By adopting the method and the device, disclosed by the invention, through obtaining the frequency characteristic data of a limited number of frequency points, the transfer function of an object model is computed. The method and the device, disclosed by the invention, have the characteristics that the algorithm is simple and the implementation on engineering application is easy, and have higher identification accuracy and stronger anti-interference property, smaller influence to the process. The effectiveness of the method is verified by a simulation experimental study.
Description
Technical field
The present invention relates to thermal technology's control technology field, the transport function that relates in particular to the system object parameter model is confirmed method and device.
Background technology
Thermal power plant's heat power engineering system image parameter model has characteristics such as big inertia, pure time-delay and time variation, the particularly characteristic of system object parameter model and is difficult to deeply factors such as accurately understanding, and the system that makes is difficult to control.How effectively to overcome various adverse effects, thereby improve the regulation quality of thermal technology's control loop, be that the thermal technology controls one of target of constantly pursuing in the field always.
For improving the controlling performance of thermal object, it is very important that the transport function of system object parameter model is carried out identification.Existing object model identifying approach mainly is based on system's time domain specification the process object model is debated knowledge; There are complex algorithm and the problem that is difficult in existing DCS, use; In addition process is exported the external ornamental of debating the knowledge signal higher requirement is also arranged; Need the pumping signal of strong amplitude, this is also bigger to the process influence.From the frequency domain angle, the frequency spectrum characteristic of thermal object is: be limit with the zero frequency, from high frequency that limit for height is arranged to unlimited low low frequency and continuous frequency spectrum.Model in Time Domain is debated and is known strategy and need use all continuous frequency spectrum resources of thermal object, and spectrum environment debates to Model in Time Domain that to know result's influence very big, and in simulation process, owing to there not being interference, frequency spectrum is very clean, can obtain good effect; But in the process of reality, owing to exist to disturb, therefore Model in Time Domain is debated that to know the strategy influence very big, even poor.
Summary of the invention
In order to solve the problems of the technologies described above, the transport function that the invention provides the system object parameter model is confirmed method and device, can carry out identification completely to system object parameter model Model Transfer function through the frequency characteristic data of preset several Frequency points.
The transport function that the invention provides the system object parameter model is confirmed method, comprising:
The transport function of predetermined system image parameter model comprises appearance time lag links such as n rank;
Select preset several identification Frequency points, obtain the amplitude transfer coefficient of each identification Frequency point and the phase transfer value of said each identification Frequency point;
Amplitude transfer coefficient according to the identification Frequency point calculates said inertia exponent number n and the inertia constant of holding link that wait;
Calculate the proportional gain and the time lag constant of said transport function according to current inertia exponent number and current inertia constant.
Correspondingly, the present invention also provides the transport function of system object parameter model to confirm device, comprising:
The parameter model tectonic element is used for the transport function of predetermined system image parameter model, and wherein said transport function comprises that n rank etc. hold the time lag link;
Amplitude transfer coefficient acquiring unit is used to obtain the amplitude transfer coefficient of several preset identification Frequency points;
The phase transfer value acquiring unit is used to obtain the phase transfer value of several preset identification Frequency points;
With the parameter calculation unit that said amplitude transfer coefficient acquiring unit, said phase transfer value acquiring unit link to each other respectively, be used for calculating inertia exponent number, inertia constant, proportional gain, the time lag constant of said transport function according to the phase transfer value of the amplitude transfer coefficient of identification Frequency point and any identification Frequency point;
With the model acquiring unit that said parameter model tectonic element, parameter calculation unit link to each other respectively, be used for according to the get parms transport function of model of said parameter calculation unit parameters calculated.
Embodiment of the present invention has following beneficial effect:
Different with the time domain identifying approach; The frequency domain identifying approach is the conserve on spectrum resource very; Need not use all continuous frequency spectrum resources of thermal object, only use the discrete unifrequency point frequency spectrum of limited quantity, quantitatively only account for the minimum part of the whole continuous frequency spectrum resources of thermal object.Can use BPF. to leach these unifrequency point spectrum signals, the undesired signal of the overwhelming majority is foreclosed.The present invention is through the frequency characteristic data of several Frequency points of acquisition limited quantity, and the iterative algorithm on the applied mathematics calculates the object model transport function.It is simple that the present invention has an algorithm, and the characteristics that are easy on practical applications, realize, also has higher debating and know precision and stronger anti-interference, and less to the process influence, Simulating Test Study has been verified the validity of this method.
Description of drawings
Fig. 1 is the process flow diagram that the transport function of system object parameter model of the present invention is confirmed method;
Fig. 2 is the first embodiment process flow diagram that the transport function of system object parameter model of the present invention is confirmed method;
Fig. 3 is the first embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed method;
Fig. 4 is the second embodiment process flow diagram that the transport function of system object parameter model of the present invention is confirmed method;
Fig. 5 is the second embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed method;
Fig. 6 is the first embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed device;
Fig. 7 is the second embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed device.
Embodiment
For making the object of the invention, technical scheme and advantage clearer, will combine accompanying drawing that the present invention is done to describe in detail further below.
Fig. 1 is the process flow diagram that the transport function of system object parameter model of the present invention is confirmed method, comprising:
S101: the transport function of predetermined system image parameter model comprises appearance time lag links such as n rank;
S102: select preset several identification Frequency points, obtain the amplitude transfer coefficient of each identification Frequency point and the phase transfer value of said each identification Frequency point;
S103: the amplitude transfer coefficient according to the identification Frequency point calculates said inertia exponent number n and the inertia constant of holding link that wait;
S104: the proportional gain and the time lag constant that calculate said transport function according to current inertia exponent number and current inertia constant.
Because thermal power plant's heat power engineering system image parameter model also has pure time-delay characteristics except having characteristics such as big inertial time sex change.So, for satisfying application of practical project, need to confirm that system object parameter model model is appearance time lag links such as n rank, be specially:
In the middle of prior art, the time domain of employing is debated knowledge need use all continuous frequency spectrum resources of thermal object, and result's influence that spectrum environment is debated knowledge to time domain is very big, disturbs owing to exist, and therefore time domain is debated the poor effect that knowledge obtains.If can identify the continuous full spectrum of frequencies characteristic of object, also can identify the object model transport function according to continuous full spectrum of frequencies characteristic fully.But see then very difficulty and unrealistic from engineering practice.The present invention selects preset several identification Frequency points; Only be that this limited several frequency characteristic of unifrequency point is debated knowledge; The frequency characteristic data that obtains these unifrequency points is very easy to; And use BPF. that these unifrequency points are debated the knowledge frequency signal and carry out Filtering Processing, therefore have stronger anti-interference again.This also is to be made up of various frequency contents because of undesired signal, and " bandpass filtering " maximum characteristics are can most unwanted frequency compositions be discharged outside, so the Points Identification for Process Frequency Response sharpest edges of unifrequency point also are stronger anti-interference.Preferably, the present invention only need select the frequency characteristic data of 3 limited identification Frequency points, just can debate knowledge fully to appearance Object with Time Delay Model Transfer functions such as high-orders.Preferably, choosing 3 identification Frequency points is respectively:
,
,
; Obtain the amplitude transfer coefficient of each identification Frequency point:;
,
; The phase transfer value that obtains said each identification Frequency point is:
;
,
.
Fig. 2 is the first embodiment process flow diagram that the transport function of system object parameter model of the present invention is confirmed method, and compared to Figure 1, Fig. 2 is the concrete refinement step of step S103.
S201: confirm first functional relation between inertia exponent number n and the inertia constant T according to first group of amplitude transfer coefficient;
S202: confirm second functional relation between inertia exponent number n and the inertia constant T according to second group of amplitude transfer coefficient;
S203: calculate inertia exponent number n and inertia constant T according to said first functional relation and said second functional relation, wherein, the said first functional relation non-equivalence is in second functional relation.
The object model transport function is (1) formula:
(1)
Ask for the frequency domain amplitude characteristic function and obtain (2) formula:
Amplitude transfer coefficient
and frequency
thereof for people (2) formula, are obtained (3) formula:
Amplitude transfer coefficient
and frequency
thereof for people (2) formula, are obtained (4) formula:
Amplitude transfer coefficient
and frequency
thereof for people (2) formula, are obtained (5) formula:
Arrangement (4) formula obtains (6) formula with (5) formula:
（6）
(6) formula is obtained first functional relation between inertia exponent number n and the inertia constant T:
（7）
Arrangement (3) formula and (4) formula obtain second functional relation between inertia exponent number n and the inertia constant T:
Can learn according to said first functional relation (7) and said second functional relation (8); Have only n and two unknown numbers of T in these two function formulas; And first functional relation (7) non-equivalence is in said second functional relation (8); So, be the value that can calculate n and T according to the equation principle of mathematics.Particularly, be the first embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed method like Fig. 3.Wherein, horizontal ordinate is represented inertia exponent number n; Ordinate is represented inertia time system T; Curve 1 expression first functional relation; Curve 2 expressions second functional relation.So, the point of crossing Sn of two curves is just by being asked.
In addition, said first functional relation (7) is one embodiment of the present invention with said second functional relation (8), and the present invention does not get rid of and can also obtain have the function formula that realizes identical function through other similar mapping mode.For example, above-mentioned first functional relation (7) is that amplitude transfer coefficient through
and
obtains.Certainly; Under the prerequisite that guarantees first functional relation and the second functional relation non-equivalence, can obtain through the amplitude transfer coefficient of
and
.
The present invention can only select the frequency characteristic data of 3 limited identification Frequency points, can choose the Frequency point more than 3, selects several Frequency point data wherein to obtain first funtcional relationship and second functional relation.For example, can choose 4 identification Frequency points is respectively:
,
,
,
.Then, using the
,
to obtain the amplitude transmission coefficient of the first function formula; using
,
amplitude transfer coefficient for the second function formula.
The proportional gain of calculating said transport function according to the amplitude transfer coefficient of identification Frequency point, current inertia exponent number and current inertia constant.Preferably, arrangement (4) formula obtains gain coefficient formula (9) formula:
Certainly, also can ask for proportional gain K through the formula that concerns of (3) formula or (5) or other Frequency point.
The time lag constant that calculates said transport function according to the phase transfer value of identification Frequency point, current inertia exponent number and current inertia constant.
Object is worth theoretical formula in
frequency plot transmission and is (10) formula:
（10）
For people (10) formula, and arrangement obtains the calculating formula of time lag constant with phase transfer value
and frequency
thereof:
Certainly, the above-mentioned formula of phase transfer value substitution (11) that also can be through
or
or other Frequency point is to calculate the time lag constant.
Fig. 4 is the second embodiment process flow diagram that the transport function of system object parameter model of the present invention is confirmed method; Compare with first embodiment of Fig. 2, Fig. 3; First embodiment illustrates on principle can obtain inertia exponent number n and inertia constant T comparatively exactly through calculation procedure; Second embodiment is then from practical application; Obtain inertia exponent number n and inertia constant T comparatively apace through simple iterative algorithm, be convenient on hardware view, realize.
S301: when calculating for the first time, said n is set to the initial value exponent number, promptly current exponent number n
NBe made as the initial value exponent number, according to second functional relation and current exponent number n
NCalculate the N time current inertia constant T
N
S302: confirm exponent number again according to first functional relation and the N time current inertia constant, obtain (N+1) inferior inertia exponent number n
N+1
S303: according to said second functional relation and said (N+1) inferior inertia exponent number n
N+1Again confirm inertia constant, obtain (N+1) inferior inertia constant T
N+1
S304: when the N time exponent number and said (N+1) inferior exponent number of confirming again satisfy pre-conditioned; Then (N+1) inferior exponent number is confirmed as current exponent number, and will confirm as current inertia constant according to (N+1) inferior inertia constant that (N+1) inferior exponent number obtains.
At first; First functional relation (7) is done suitable distortion; Make it be applicable to iterative algorithm; Replace
with
, replace
with
; And preferably, said n is set to initial value exponent number 2:
Above result of calculation is: the exponent number value of
inferior appearance inertia.
is the cycle calculations number of times;
is last calculation times;
is the exponent number of
inferior appearance inertia of confirming again, and
is current inertia constant.Wherein
is for calculating for the first time; Because current inertia constant is not calculated as yet, so its default value is represented with 0.
Secondly; Second functional relation (8) is done suitable distortion; Make it be applicable to iterative algorithm; Replace
with
, replace
, obtain (13) formula with
:
（13）
Above result of calculation is:
inferior appearance inertia constant calculated value.The inertia exponent number n that promptly utilizes (12) formula to obtain
NAgain confirm inertia constant T
N
Do the gain coefficient of suitable The deformation calculation system object parameter model with (9) formula; Replace
with
, replace
, replace
, obtain (14) formula with
with
:
Above result of calculation is:
inferior gain coefficient calculated value.
With (11) formula computing system image parameter model time lag constant; Replace
with
, replace
with
; Replace
with
, obtain (15) formula:
Above result of calculation is:
inferior time lag constant calculations value.
Fig. 5 is the second embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed method.In order specifically to illustrate the iterative algorithm of present embodiment, enumerate an instantiation below in conjunction with Fig. 5.
Supposing the system image parameter model is appearance time lag links such as one 8 rank, and the theoretical transport function of this object is:
Below the utilization computing method of the present invention to its calculate, emulation.
Select 3 identification Frequency points to be:
The amplitude transfer coefficient that obtains each identification Frequency point is:
Obtain the phase transfer value of said each identification Frequency point:
N is preset as initial value 2 with the inertia exponent number, and n1=2 substitution (13) formula is obtained current inertia constant T1; This inertia constant T1 substitution (12) formula is confirmed exponent number again, obtain inertia exponent number n2; This inertia exponent number n2 substitution (13) formula is confirmed inertia constant again, obtain inertia constant T2; So analogize, constantly iteration., calculate inertia constant T1 according to inertia exponent number n1 for the first time and be some S1 as description of drawings with shown in Figure 5, the inertia exponent number n2 that calculates according to inertia constant T1 for the second time is some S2; The inertia constant T2 that calculates according to inertia exponent number n2 for the second time is some S3 ... So analogize; As shown in Figure 5, ordered series of numbers S1, S2; S3 ... Constantly converge on Sn, be institute and ask.Simulation calculation is as a result shown in the following table:
Can find out that from last table after cycle index was greater than 49 times, each calculation of parameter value stabilization was constant, but also have certain error that reason is to have due to the error in the identified parameters that each single-point identification frequency obtains with target transfer function theoretical parameter value.The object model transport function final calculation result that obtains is:
In the scope that error allows; Can obtain inertia exponent number n and inertia constant T fast through the method shown in second embodiment; When the N time inertia exponent number and said (N+1) inferior inertia exponent number of confirming again satisfy pre-conditioned; Then that (N+1) is inferior inertia exponent number is confirmed as the inertia exponent number of final value, and will confirm as the final value inertia constant according to (N+1) inferior inertia constant that (N+1) inferior inertia exponent number obtains.Preferably, for present embodiment determined (12) formula and (13) formula, saidly pre-conditionedly can comprise (14) formula:
As current inertia exponent number n
N+1Inertia exponent number n with the last time
NSatisfy (14) formula, can utilize current inertia exponent number n
N+1Calculate inertia constant T
N+1After, and stop circulation.As above shown in the table; After cycle index is greater than 8 times;
promptly satisfies condition; So, need not in the calculating of carrying out N9.
When current inertia exponent number is integer, then this inertia exponent number is confirmed as final inertia exponent number, current inertia constant is confirmed as final inertia constant;
When current inertia exponent number is not integer, to its processing that rounds up, obtain final inertia exponent number, and confirm inertia constant again according to said second functional relation and said final inertia exponent number, obtain final inertia constant;
According to the amplitude transfer coefficient of any identification Frequency point, final inertia exponent number and final inertia constant, calculate said the grade and hold the final proportional gain of link, obtain final proportional gain;
According to the phase transfer value of any identification Frequency point, final inertia exponent number and final inertia constant, calculate the time lag constant, obtain final time lag constant;
According to said final inertia exponent number, said final inertia constant, said final proportional gain, said final time lag constant, obtain the transport function of system object parameter model.
Since heat power engineering system image parameter model etc. the appearance exponent number generally can be preset as integer, so, the processing that can the final calculation result of inertia exponent number be rounded up.Confirm parameters such as final inertia constant, proportional gain, time lag constant again according to final inertia exponent number.In sum, available object model transport function final calculation result is:
.Ining contrast to original system resembles the theoretical transport function of parameter model and can know that this result of calculation is more accurate.
Need to prove;
carried out the initial value setting;
initial value scope is so long as positive integer gets final product, and reason is can finally converge on the Sn point through this iterative algorithm.Preferably,
initial value is made as 2.
In sum; The transport function of the system object parameter model that the embodiment of the invention provides is confirmed method; Have calculated amount less higher with precision, advantage such as be easy in active computer control (DCS) system, use, significant to the characteristic of accurate understanding and recognition system image parameter model.
Fig. 6 is the first embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed device, comprising:
The parameter model tectonic element is used for the transport function of predetermined system image parameter model, and wherein said transport function comprises that n rank etc. hold the time lag link;
Amplitude transfer coefficient acquiring unit is used to obtain the amplitude transfer coefficient of several preset identification Frequency points;
The phase transfer value acquiring unit is used to obtain the phase transfer value of several preset identification Frequency points;
With the parameter calculation unit that said amplitude transfer coefficient acquiring unit, said phase transfer value acquiring unit link to each other respectively, be used for calculating inertia exponent number, inertia constant, proportional gain, the time lag constant of said transport function according to the phase transfer value of the amplitude transfer coefficient of identification Frequency point and any identification Frequency point;
With the model acquiring unit that said parameter model tectonic element, parameter calculation unit link to each other respectively, be used for according to the get parms transport function of model of said parameter calculation unit parameters calculated.
Because thermal power plant's heat power engineering system image parameter model also has pure time-delay characteristics except having characteristics such as big inertial time sex change.So, need to confirm that system object parameter model model is appearance time lag links such as n rank, be specially:
As shown in Figure 6, also comprise first algorithm unit and second algorithm unit that link to each other with said parameter calculation unit,
Said first algorithm unit comprises: first function is confirmed the unit, is used for confirming first functional relation between inertia exponent number n and the inertia constant T; Second function is confirmed the unit, is used for confirming second functional relation between inertia exponent number n and the inertia constant T;
Said second algorithm unit comprises: the scale-up factor function is confirmed the unit, is used for confirming to calculate the functional relation of scale-up factor; The time lag constant function is confirmed the unit, is used for confirming to calculate the functional relation of time lag constant.
Fig. 7 is the second embodiment synoptic diagram that the transport function of system object parameter model of the present invention is confirmed device, and said parameter calculation unit comprises:
Initial value is provided with the unit, is used for said n is carried out the initial value setting;
With said initial value unit, said amplitude transfer coefficient acquiring unit, said second function are set and confirm the inertia constant computing unit that the unit links to each other respectively, be used for confirming inertia constant again according to the initial setting up value of said second functional relation and inertia exponent number;
Confirm the inertia exponent number computing unit that the unit links to each other respectively with said inertia constant computing unit, said amplitude transfer coefficient acquiring unit, said first function, be used for confirming exponent number again according to first functional relation and current inertia constant;
Said inertia constant computing unit also is used for confirming current inertia constant again according to said second functional relation and current inertia exponent number.
Said parameter calculation unit also comprises:
Be connected the circulation stop element between said inertia exponent number computing unit and the said inertia constant computing unit; Be used for satisfying pre-conditioned when the N time exponent number and said (N+1) inferior exponent number of confirming again; Then (N+1) inferior exponent number is confirmed as current exponent number, and will confirm as current inertia constant according to (N+1) inferior inertia constant that (N+1) inferior exponent number obtains.
Confirm the scale-up factor computing unit that unit, said final argument arrangement unit links to each other respectively with said amplitude transfer coefficient acquiring unit, said scale-up factor function, be used for the current scale-up factors of appearance link such as amplitude transfer coefficient, scale-up factor function calculation according to any identification Frequency point be said;
Confirm the time lag constant computation unit that unit, said final argument arrangement unit link to each other respectively with said phase transfer value acquiring unit, said time lag constant function, the phase transfer value, the time lag constant function that are used for according to any identification Frequency point calculate the current time lag constants of appearance link such as said;
Put the unit in order with the final argument that said inertia constant computing unit, said inertia exponent number computing unit, said circulation stop element, said scale-up factor computing unit, said time lag constant computation unit link to each other respectively; Be used for when current inertia exponent number is integer; Then this inertia exponent number is confirmed as final inertia exponent number, current inertia constant is confirmed as final inertia constant; When current inertia exponent number is not integer, to its processing that rounds up, obtain final inertia exponent number, and confirm inertia constant again according to said second functional relation and said final inertia exponent number, obtain final inertia constant;
Said scale-up factor computing unit also is used for the amplitude transfer coefficient according to any identification Frequency point, final inertia exponent number and final inertia constant, confirms the said scale-up factor that holds link that waits again, obtains to wait the final scale-up factor that holds link;
Said time lag constant computation unit also is used for the phase place transfer coefficient according to any identification Frequency point, final inertia exponent number and final inertia constant, confirms again to obtain final time lag constant by the time lag constant;
Said model acquiring unit also is used for obtaining the transport function of system object parameter model according to said final inertia exponent number, said final inertia constant, said final proportional gain, said final time lag constant.
Step in the apparatus and method of Fig. 7 is corresponding, identical in the embodiment of each unit and the method.
This paper has proposed a kind of strategy of knowing appearance Object with Time Delay Model Transfer functions such as high-order of debating based on the finite frequency performance data; And this strategy conducted in-depth analysis, provide frequency domain and debated the principle of knowledge, have the simple and stronger anti-interference of algorithm; Also has the higher knowledge precision of debating; Be easy in existing DCS, realize, to characteristics such as process influence are less, simulation study verified should strategy validity.
Through the description of above embodiment, those skilled in the art can be well understood to the present invention and can realize by the mode that software adds essential hardware platform, can certainly all implement through hardware.Based on such understanding; All or part of can the coming out that technical scheme of the present invention contributes to background technology with the embodied of software product; This computer software product can be stored in the storage medium, like ROM/RAM, magnetic disc, CD etc., comprises that some instructions are with so that a computer equipment (can be a personal computer; Server, the perhaps network equipment etc.) carry out the described method of some part of each embodiment of the present invention or embodiment.
Above-described embodiment of the present invention does not constitute the qualification to protection domain of the present invention.Any modification of within spirit of the present invention and principle, being done, be equal to replacement and improvement etc., all should be included within the claim protection domain of the present invention.
Claims (10)
1. the transport function of a system object parameter model is confirmed method, it is characterized in that, comprising:
The transport function of predetermined system image parameter model comprises appearance time lag links such as n rank;
Select preset several identification Frequency points, obtain the amplitude transfer coefficient of each identification Frequency point and the phase transfer value of said each identification Frequency point;
Amplitude transfer coefficient according to the identification Frequency point calculates said inertia exponent number n and the inertia constant of holding link that wait;
Calculate the proportional gain and the time lag constant of said transport function according to current inertia exponent number and current inertia constant.
2. the transport function of system object parameter model according to claim 1 is confirmed method, it is characterized in that, according to the said inertia exponent number n of appearance link and the step of inertia constant of waiting of the amplitude transfer coefficient calculating of identification Frequency point, comprising:
Confirm first functional relation between inertia exponent number n and the inertia constant T according to first group of amplitude transfer coefficient;
Confirm second functional relation between inertia exponent number n and the inertia constant T according to second group of amplitude transfer coefficient;
Calculate inertia exponent number n and inertia constant T according to said first functional relation and said second functional relation, wherein, the said first functional relation non-equivalence is in second functional relation.
3. the transport function of system object parameter model according to claim 2 is confirmed method, it is characterized in that, the step according to said first functional relation and said second functional relation calculating inertia exponent number n and inertia constant T comprises:
During first calculating, said n is set to the initial value exponent number, promptly current exponent number n
NBe made as the initial value exponent number, according to second functional relation and current exponent number n
NCalculate the N time current inertia constant T
N
Again confirm exponent number according to first functional relation and the N time current inertia constant, obtain (N+1) inferior inertia exponent number n
N+1
According to said second functional relation and said (N+1) inferior inertia exponent number n
N+1Again confirm inertia constant, obtain (N+1) inferior inertia constant T
N+1
4. the transport function of system object parameter model according to claim 3 is confirmed method, it is characterized in that:
When the N time exponent number and said (N+1) inferior exponent number of confirming again satisfy pre-conditioned; Then (N+1) inferior exponent number is confirmed as current exponent number, and will confirm as current inertia constant according to (N+1) inferior inertia constant that (N+1) inferior exponent number obtains.
5. confirm method according to the transport function of each described system object parameter model of claim 2 to 4, it is characterized in that:
When current inertia exponent number is integer, then this inertia exponent number is confirmed as final inertia exponent number, current inertia constant is confirmed as final inertia constant;
When current inertia exponent number is not integer, to its processing that rounds up, obtain final inertia exponent number, and confirm inertia constant again according to said second functional relation and said final inertia exponent number, obtain final inertia constant;
According to the amplitude transfer coefficient of any identification Frequency point, final inertia exponent number and final inertia constant, calculate said the grade and hold the final proportional gain of link, obtain final proportional gain;
According to the phase transfer value of any identification Frequency point, final inertia exponent number and final inertia constant, calculate the time lag constant, obtain final time lag constant;
According to said final inertia exponent number, said final inertia constant, said final proportional gain, said final time lag constant, obtain the transport function of system object parameter model.
6. the transport function of a system object parameter model is confirmed device, it is characterized in that, comprising:
The parameter model tectonic element is used for the transport function of predetermined system image parameter model, and wherein said transport function comprises that n rank etc. hold the time lag link;
Amplitude transfer coefficient acquiring unit is used to obtain the amplitude transfer coefficient of several preset identification Frequency points;
The phase transfer value acquiring unit is used to obtain the phase transfer value of several preset identification Frequency points;
With the parameter calculation unit that said amplitude transfer coefficient acquiring unit, said phase transfer value acquiring unit link to each other respectively, be used for calculating inertia exponent number, inertia constant, proportional gain, the time lag constant of said transport function according to the phase transfer value of the amplitude transfer coefficient of identification Frequency point and any identification Frequency point;
With the model acquiring unit that said parameter model tectonic element, parameter calculation unit link to each other respectively, be used for according to the get parms transport function of model of said parameter calculation unit parameters calculated.
7. the transport function of system object parameter model according to claim 6 is confirmed device, it is characterized in that, comprises first algorithm unit and second algorithm unit that link to each other with said parameter calculation unit,
Said first algorithm unit comprises: first function is confirmed the unit, is used for confirming first functional relation between inertia exponent number n and the inertia constant T; Second function is confirmed the unit, is used for confirming second functional relation between inertia exponent number n and the inertia constant T;
Said second algorithm unit comprises: the scale-up factor function is confirmed the unit, is used for confirming to calculate the functional relation of scale-up factor; The time lag constant function is confirmed the unit, is used for confirming to calculate the functional relation of time lag constant.
8. the transport function of system object parameter model according to claim 7 is confirmed device, it is characterized in that, said parameter calculation unit comprises:
Initial value is provided with the unit, is used for said n is carried out the initial value setting;
With said initial value unit, said amplitude transfer coefficient acquiring unit, said second function are set and confirm the inertia constant computing unit that the unit links to each other respectively, be used for confirming inertia constant again according to the initial setting up value of said second functional relation and inertia exponent number;
Confirm the inertia exponent number computing unit that the unit links to each other respectively with said inertia constant computing unit, said amplitude transfer coefficient acquiring unit, said first function, be used for confirming exponent number again according to first functional relation and current inertia constant;
Said inertia constant computing unit also is used for confirming current inertia constant again according to said second functional relation and current inertia exponent number.
9. the transport function of system object parameter model according to claim 8 is confirmed device, it is characterized in that, said parameter calculation unit also comprises:
Be connected the circulation stop element between said inertia exponent number computing unit and the said inertia constant computing unit; Be used for satisfying pre-conditioned when the N time exponent number and said (N+1) inferior exponent number of confirming again; Then (N+1) inferior exponent number is confirmed as current exponent number, and will confirm as current inertia constant according to (N+1) inferior inertia constant that (N+1) inferior exponent number obtains.
10. confirm device according to the transport function of each described system object parameter model of claim 7 to 9, it is characterized in that said parameter calculation unit also comprises:
Confirm the scale-up factor computing unit that unit, said final argument arrangement unit links to each other respectively with said amplitude transfer coefficient acquiring unit, said scale-up factor function, be used for the current scale-up factors of appearance link such as amplitude transfer coefficient, scale-up factor function calculation according to any identification Frequency point be said;
Confirm the time lag constant computation unit that unit, said final argument arrangement unit link to each other respectively with said phase transfer value acquiring unit, said time lag constant function, the phase transfer value, the time lag constant function that are used for according to any identification Frequency point calculate the current time lag constants of appearance link such as said;
Put the unit in order with the final argument that said inertia constant computing unit, said inertia exponent number computing unit, said circulation stop element, said scale-up factor computing unit, said time lag constant computation unit link to each other respectively; Be used for when current inertia exponent number is integer; Then this inertia exponent number is confirmed as final inertia exponent number, current inertia constant is confirmed as final inertia constant; When current inertia exponent number is not integer, to its processing that rounds up, obtain final inertia exponent number, and confirm inertia constant again according to said second functional relation and said final inertia exponent number, obtain final inertia constant;
Said scale-up factor computing unit also is used for the amplitude transfer coefficient according to any identification Frequency point, final inertia exponent number and final inertia constant, confirms the said scale-up factor that holds link that waits again, obtains to wait the final scale-up factor that holds link;
Said time lag constant computation unit also is used for the phase place transfer coefficient according to any identification Frequency point, final inertia exponent number and final inertia constant, confirms again to obtain final time lag constant by the time lag constant;
Said model acquiring unit also is used for obtaining the transport function of system object parameter model according to said final inertia exponent number, said final inertia constant, said final proportional gain, said final time lag constant.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110282741 CN102393866B (en) | 2011-09-22 | 2011-09-22 | Transfer function determination method of system object parameter model and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201110282741 CN102393866B (en) | 2011-09-22 | 2011-09-22 | Transfer function determination method of system object parameter model and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102393866A true CN102393866A (en) | 2012-03-28 |
| CN102393866B CN102393866B (en) | 2013-08-14 |
Family
ID=45861189
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201110282741 Active CN102393866B (en) | 2011-09-22 | 2011-09-22 | Transfer function determination method of system object parameter model and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102393866B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105242572A (en) * | 2015-10-23 | 2016-01-13 | 广东电网有限责任公司电力科学研究院 | Mixing identification method and system for thermal process of thermal power plant |
| CN106094509A (en) * | 2016-06-22 | 2016-11-09 | 广东电网有限责任公司电力科学研究院 | A kind of LCR filtering control method and device |
| CN106200373A (en) * | 2016-06-22 | 2016-12-07 | 广东电网有限责任公司电力科学研究院 | A kind of method and device of the approximate differential signal of acquisition process signal |
| CN106200372A (en) * | 2016-06-22 | 2016-12-07 | 广东电网有限责任公司电力科学研究院 | A kind of method and device of the approximate differential signal obtaining low order Object Process signal |
| CN106406123A (en) * | 2016-11-30 | 2017-02-15 | 中国兵器装备集团自动化研究所 | Transfer function-based rotary table index decomposition method |
| CN107748497A (en) * | 2017-09-28 | 2018-03-02 | 中国科学院长春光学精密机械与物理研究所 | The Model Distinguish of kinetic control system and Parameters design and system |
| CN112882384A (en) * | 2021-01-11 | 2021-06-01 | 清华大学 | Current-torque transfer model identification method and device for servo motor on robot |
| CN114184192A (en) * | 2021-12-27 | 2022-03-15 | 北京计算机技术及应用研究所 | Method for obtaining angular velocity measurement channel transfer function of inertial measurement unit |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100358235C (en) * | 2004-12-03 | 2007-12-26 | 三菱电机株式会社 | Motor controller |
| CN101403924A (en) * | 2008-11-03 | 2009-04-08 | 东华大学 | Method for on-line delay detection of control loop |
-
2011
- 2011-09-22 CN CN 201110282741 patent/CN102393866B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100358235C (en) * | 2004-12-03 | 2007-12-26 | 三菱电机株式会社 | Motor controller |
| CN101403924A (en) * | 2008-11-03 | 2009-04-08 | 东华大学 | Method for on-line delay detection of control loop |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105242572A (en) * | 2015-10-23 | 2016-01-13 | 广东电网有限责任公司电力科学研究院 | Mixing identification method and system for thermal process of thermal power plant |
| CN106094509A (en) * | 2016-06-22 | 2016-11-09 | 广东电网有限责任公司电力科学研究院 | A kind of LCR filtering control method and device |
| CN106200373A (en) * | 2016-06-22 | 2016-12-07 | 广东电网有限责任公司电力科学研究院 | A kind of method and device of the approximate differential signal of acquisition process signal |
| CN106200372A (en) * | 2016-06-22 | 2016-12-07 | 广东电网有限责任公司电力科学研究院 | A kind of method and device of the approximate differential signal obtaining low order Object Process signal |
| CN106200372B (en) * | 2016-06-22 | 2019-03-15 | 广东电网有限责任公司电力科学研究院 | A kind of method and device for the approximate differential signal obtaining low order Object Process signal |
| CN106406123A (en) * | 2016-11-30 | 2017-02-15 | 中国兵器装备集团自动化研究所 | Transfer function-based rotary table index decomposition method |
| CN107748497A (en) * | 2017-09-28 | 2018-03-02 | 中国科学院长春光学精密机械与物理研究所 | The Model Distinguish of kinetic control system and Parameters design and system |
| CN107748497B (en) * | 2017-09-28 | 2019-10-15 | 中国科学院长春光学精密机械与物理研究所 | The Model Distinguish of kinetic control system and Parameters design and system |
| CN112882384A (en) * | 2021-01-11 | 2021-06-01 | 清华大学 | Current-torque transfer model identification method and device for servo motor on robot |
| CN114184192A (en) * | 2021-12-27 | 2022-03-15 | 北京计算机技术及应用研究所 | Method for obtaining angular velocity measurement channel transfer function of inertial measurement unit |
| CN114184192B (en) * | 2021-12-27 | 2023-09-26 | 北京计算机技术及应用研究所 | Method for acquiring angular velocity measurement channel transfer function of inertial measurement device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102393866B (en) | 2013-08-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102393866B (en) | Transfer function determination method of system object parameter model and device | |
| TWI638272B (en) | System and method for performing neural network computations for a neural network and related normalization circuitry | |
| CN111639429B (en) | Underwater sound field numerical simulation method, system and medium based on Chebyshev polynomial spectrum | |
| Xu et al. | Feasibility study on the least square method for fitting non-Gaussian noise data | |
| CN108008632B (en) | State estimation method and system of time-lag Markov system based on protocol | |
| CN112737854B (en) | Service chain migration method and device based on energy consumption and service quality | |
| Bozorg et al. | Control of time delay processes with uncertain delays: Time delay stability margins | |
| Pakzad et al. | Stability analysis of time-delayed linear fractional-order systems | |
| CN107766291A (en) | A kind of method of remaining time delay in acquisition very long baseline interferometry(VLBI | |
| CN109218769A (en) | A kind of recommended method and relevant device of direct broadcasting room | |
| Oprzędkiewicz et al. | An estimation of accuracy of Oustaloup approximation | |
| Wirtz et al. | Efficient a-posteriori error estimation for nonlinear kernel-based reduced systems | |
| CN111476837A (en) | Adaptive stereo matching optimization method and device, equipment and storage medium thereof | |
| CN103473437A (en) | Cloud service platform selection method based on fuzzy ideal target analysis | |
| CN110389365A (en) | A kind of satellite navigation locating method and device, terminal, storage medium | |
| CN110493029A (en) | A kind of transmission method, transmitting device, electronic equipment and storage medium | |
| CN108574649B (en) | Method and device for determining digital predistortion coefficient | |
| CN110286275A (en) | The method and device of system stability when a kind of determining equipment is grid-connected | |
| CN109358609B (en) | A kind of adjustment method and system of model predictive controller | |
| CN103926561B (en) | A kind of parameter estimation weights method for designing eliminated for the singular value of ultra-short baseline alignment error calibration | |
| Kim | Two simple numerical methods for the free boundary in one‐phase Stefan problem | |
| Zhao et al. | Robust control against uncertainty quartet: A polynomial approach | |
| CN103513618B (en) | The control method of industrial process and equipment | |
| WO2013027332A1 (en) | Information processing device, information processing method, and program | |
| Roy et al. | Designing fractional-order PIλDμ controller using differential harmony search algorithm |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8 Patentee after: ELECTRIC POWER RESEARCH INSTITUTE OF GUANGDONG POWER GRID Co.,Ltd. Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8 Patentee before: ELECTRIC POWER RESEARCH INSTITUTE OF GUANGDONG POWER GRID Corp. |
|
| CP01 | Change in the name or title of a patent holder | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20171117 Address after: 510655, No. fifth, No. 146-150, Whampoa Avenue, Tianhe District, Guangdong, Guangzhou Patentee after: GUANGDONG DIANKEYUAN ENERGY TECHNOLOGY Co.,Ltd. Address before: 510080 Dongfeng East Road, Dongfeng, Guangdong, Guangzhou, Zhejiang Province, No. 8 Patentee before: ELECTRIC POWER RESEARCH INSTITUTE OF GUANGDONG POWER GRID Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| CP03 | Change of name, title or address |
Address after: Room 501-503, annex building, Huaye building, No.1-3 Chuimao new street, Xihua Road, Yuexiu District, Guangzhou City, Guangdong Province 510000 Patentee after: China Southern Power Grid Power Technology Co.,Ltd. Address before: 510655 5th floor, no.146-150, Huangpu Avenue Middle Road, Tianhe District, Guangzhou City, Guangdong Province Patentee before: GUANGDONG DIANKEYUAN ENERGY TECHNOLOGY Co.,Ltd. |
|
| CP03 | Change of name, title or address |