CN103823996B - It is applicable to the interpolation method that power system high-speed sampling calculates - Google Patents
It is applicable to the interpolation method that power system high-speed sampling calculates Download PDFInfo
- Publication number
- CN103823996B CN103823996B CN201410098732.9A CN201410098732A CN103823996B CN 103823996 B CN103823996 B CN 103823996B CN 201410098732 A CN201410098732 A CN 201410098732A CN 103823996 B CN103823996 B CN 103823996B
- Authority
- CN
- China
- Prior art keywords
- interpolation
- power system
- frequency
- parameter table
- sampling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005070 sampling Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 238000007711 solidification Methods 0.000 abstract description 3
- 230000008023 solidification Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Complex Calculations (AREA)
Abstract
The present invention relates to a kind of interpolation method being applicable to the calculating of power system high-speed sampling.During high-speed sampling calculates, sample frequency converts according to the change of the own frequency of power system, thus meets computational accuracy requirement.In the case of hardware sample frequency is fixing, generally uses software phlase locking to complete above-mentioned requirements, but software phlase locking calculates overlong time, thus fail to be used widely.The present invention proposes solidification form, the software of off-line calculation sample interpolation desired parameters and judges that the mode of computation of table lookup solves the problems referred to above online.Institute's extracting method alleviates the operation burden of digital signal processor so that it is can carry out signal transacting efficiently;Also need not increase hardware spending, change hardware design;Meanwhile, institute's extracting method can be widely applied to other power system occasions.
Description
Technical field
The invention belongs to power system automatic field, relate to power system electric quantity signal (voltage, electric current)
High speed fixed sample rate under sample interpolation technology.
Background technology
In power system, electricity (voltage, electric current) signal high-speed sampling should at the aspect such as calculation of capacity, the quality of power supply
With more and more extensive.The sample frequency related in calculating high-speed sampling needs the own frequency according to power system
The change of rate and convert, thus reach to meet the computational accuracy of requirement.At present, two ways is had to ask to process this
Topic, one is to utilize hardware phase-locked-loop to carry out hardware phase lock high precision, is the most dynamically changing sample frequency,
Thus obtain correct sampled data;But it requires special hardware resource, need to coordinate specially before phase-locked simultaneously
The hardware frequency-survey loop of door, filter circuit ensures its correctness.It is phase-locked that another kind utilizes software mode to realize, but
Digital signal processor carries major part and calculates and control work, and the situation for supercomputing calculates due to it
Overlong time and can not get application.The present invention utilizes software interpolation method to solve the problems referred to above, and utilizes off-line meter
Calculating, the mode tabled look-up online solves this problem, has both alleviated the operation burden of digital signal processor, makes numeral letter
Number processor can carry out signal transacting efficiently;Also need not increase hardware spending, change hardware design;Meanwhile,
Can facilitate this module application in other power system occasions.
Summary of the invention
The technical problem to be solved is that proposition is a kind of is applicable to inserting of power system high-speed sampling calculating
Value method.When hardware sample frequency is fixed, in the range of the power system frequency of regulation, by this interpolation meter
Calculation method high speed computation of table lookup can go out the sampled point after interpolation, and ensures that its monocycle sampling number keeps constant.
Thus prevent the generation of lack sampling or over-sampling, it is ensured that computational accuracy.
For solving above-mentioned technical problem, the technical scheme is that
First fixing frequency range and frequency accuracy are set: frequency range is as FL≤Fc≤FH, wherein FLRefer to
The minimum frequency set, FHRefer to the peak frequency set, FCRefer to current power system frequency, generally take
FL=45Hz;FH=55Hz;Accuracy rating is set as Δ F, general Δ F≤0.01Hz.According to this setting, off-line
Calculating interpolation parameter table required during program performs, software calls according to online judged result in the process of implementation
The interpolation parameter table of curable offline.Interpolation parameter table calculation procedure is as follows:
(1) initially set: set each cycle sampling number as N, if power system standard frequency is Fs(one
As be 50Hz), the sampling period of its correspondence is Ts, have Ts=1/Fs/N;Set current power system frequency
For Fc, the sampling period of its correspondence is Tc, have Tc=1/Fc/N;Sampling period with the difference in system cycle is
Δ T=| Tc-Ts|。
(2) table number calculates: according to linear interpolation formula y (n)=(1-α) y (n-1)+α y (n+1), wherein y (n-1)
Data previous fixed time sampled point for needing interpolation, y (n+1) is for need the data later of interpolation to fix
Instance sample point, α is different with different system frequencies and different sampling instant, it is therefore desirable to offline design is inserted
Value parameter table T.Frequency accuracy is as claimed in claim 2, and the interpolation parameter table of off-line calculation amounts to
(FH-FL)/Δ F opens.
(3) single table design: set the sampled value of hardware fixed time as S0To Sn, it is judged that the interpolation after calculating
Data are D0To Dn;Wherein 0≤n < N;Sampling period and system cycle difference are Δ T,
α=Δ T/Ts=| Fc-Fs|/Fc;Interpolation method is: D0=S0;D1=(1-α) × S1+α×S2;
D2=(1-2 α) × S2+2α×S3…Dn=(1-n α) × Sn+nα×Sn+1, calculate N number of interpolation number successively
According to.Thus, single table is designed as T={0, α, 2 α, 3 α, 4 α ... (N-1) α }.With this, according to difference
Power system frequency design (FH-FL)/Δ F opens table for using at sequence of threads.
Specifically comprising the following steps that of the interpolation method being applicable to the calculating of power system high-speed sampling proposed
1) sample count is set as SAM_CNT;The sampled value of hardware fixed time is S0To Sn, it is judged that meter
Interpolated data after calculation is D0To Dn;Making all data initial values is 0;
2) sample values S is read0To Sn;
3) current system frequency value F is readc, work as Fc=Fs, enter step 6);
4) as system ongoing frequency Fc< 45Hz or Fc> 55Hz time, enter step 6);
5) as system ongoing frequency 45Hz≤FcDuring≤55Hz, according to FcValue call the interpolation of curable offline
Parameter list, uses following interpolation algorithm D0=S0;D1=(1-α) × S1+α×S2;
D2=(1-2 α) × S2+2α×S3…Dn=(1-n α) × Sn+nα×Sn+1, calculate N number of interpolated data successively, enter
Enter step 7);
6): D0=S0;D1=S1;D2=S2;……Dn=Sn;Enter step 7);
7): sample count SAM_CNT adds 1;Return step 1);
Sampling number N, system frequency range F in algorithmL≤Fc≤FHAnd frequency accuracy set Δ F can be according to reality
Border application requires to change, and first according to setting value off-line calculation interpolation parameter table T before program execution, program performs
During according to FcJudged result directly invokes the interpolation parameter table of curable offline, thus quickly obtains institute
The interpolated data needed.
Beneficial effects of the present invention:
Having the change of frequency by oneself according to power system for meeting sample frequency needs and convert, the present invention proposes profit
The problems referred to above are solved with software interpolation method.According to system frequency range and precision set, off-line calculation interpolation is joined
Number solidification table, then uses and judges that the mode tabled look-up solves the complicated calculations problem of software phlase locking online so that it is
Really applied under supercomputing occasion;Overcome the hardware spending that hardware phase lock high precision needs to increase, also simultaneously
It is not required to change hardware design.Carried interpolation method and corresponding module of software and hardware can conveniently be applied to other power systems
Occasion.
Accompanying drawing explanation
Fig. 1 is the interpolation algorithm flow chart of inventive algorithm.
Fig. 2 is the interpolation parameter solidification table calculation flow chart of inventive algorithm.
Detailed description of the invention
In conjunction with flow chart to being embodied as describing in detail.
(1) interpolation parameter table needed for utilizing interpolation method off-line calculation to go out during program performs:
As a example by every cycle fixed sample 24 under fixing power system frequency 50Hz, its sample frequency is
1200Hz.Sampling period should be Ts=1000000/ (50*24);Precision takes 0.01Hz, and off-line calculation form is altogether
Count 1000.Data in every form correspondence respectively is ongoing frequency difference multiple table.
Example: when system frequency is FcDuring=48Hz, this actual samples cycle should be Tc=1000000/ (48*24);
Ongoing frequency difference ratio i.e. cycle difference is than for α=Δ T/Ts=| Tc-Ts|/Ts=| Fc-Fs|/Fc=0.04;
Making interpolation form is T48=0,0.04,0.08,0.12,0.16.0.20.0.24 ... ..0.92};Produce it equally
His form.
(2) software judges the most online, calls the interpolation parameter of curable offline according to judged result
Table, calculating interpolated data:
Example: set sampled data as D0、D1……D23;Judge that the interpolation after calculating is S0、S1……S23;
As system ongoing frequency FcDuring=50Hz, never call interpolation table program, S0=D0,S1=D1... S23=D23;
As system ongoing frequency 45Hz≤FcDuring≤55Hz, call interpolation table routine data D0=S0;
D1=(1-α) × S1+α×S2;D2=(1-2 α) × S2+2α×S3…D23=(1-23 α) × S23+23α×S24
Calculate all 24 interpolation successively.
As system ongoing frequency Fc< 45Hz or Fc> 55Hz time, never call interpolation table program, S0=D0,
S1=D1... S23=D23;
Above-described embodiment limits the present invention, all changes made according to technical solution of the present invention, institute never in any form
When the function produced is without departing from the scope of technical solution of the present invention, belong to protection scope of the present invention.
Claims (4)
1. one kind is applicable to the interpolation method that power system high-speed sampling calculates, it is characterised in that: set fixing
Frequency range and frequency accuracy: frequency range is FL≤Fc≤FH, wherein FLRefer to the minimum frequency set,
FHRefer to the peak frequency set, FCReferring to current power system frequency, accuracy rating is set as Δ F, according to this
Setting, off-line calculation goes out interpolation parameter table required during program performs, and software is in the process of implementation according to sentencing online
The interpolation parameter table of disconnected call by result curable offline;Described interpolation parameter table calculation procedure is as follows:
Step S01: initially set: set each cycle sampling number as N, if power system standard frequency is
Fs, the sampling period of its correspondence is Ts, have Ts=1/Fs/N;Set current power system frequency as Fc, its
The corresponding sampling period is Tc, have Tc=1/Fc/N;The sampling period of current power system and power system mark
The difference in quasi-sampling period is Δ T=| Tc-Ts|;
Step S02: table number calculates: according to linear interpolation formula y (n)=(1-α) y (n-1)+α y (n+1), its
Middle y (n-1) is the data previous fixed time sampled point needing interpolation, after y (n+1) is the data needing interpolation
One fixed time sampled point, α is different with different system frequencies and different sampling instant, it is therefore desirable to from
Line design interpolation parameter table T;The interpolation parameter table of off-line calculation amounts to (FH-FL)/Δ F opens;
Step S03: single table design: set the sampled value of hardware fixed time as S0To Sn, it is judged that after calculating
Interpolated data be D0To Dn;Wherein 0≤n < N;The sampling period of current power system and power system
The difference in standard sample cycle is Δ T, α=Δ T/Ts=| Fc-Fs|/Fc;Interpolation method is: D0=S0;
D1=(1-α) × S1+α×S2;D2=(1-2 α) × S2+2α×S3…Dn=(1-n α) × Sn+nα×Sn+1
Calculate N number of interpolated data successively;Thus, single table is designed as T={0, α, 2 α, 3 α, 4 α ... (N-1) α };
With this, design (F according to different power system frequenciesH-FL)/Δ F opens table for using at sequence of threads.
The interpolation method being applicable to the calculating of power system high-speed sampling the most according to claim 1, it is special
Levy and be: described FL=45Hz;FH=55Hz;And Δ F≤0.01Hz.
The interpolation method being applicable to the calculating of power system high-speed sampling the most according to claim 1, it is special
Levy and be: described software calls the interpolation parameter table of curable offline in the process of implementation according to online judged result
Implementation is as follows:
Step S11: set sample count as SAM_CNT;The sampled value of hardware fixed time is S0To Sn,
Judge that the interpolated data after calculating is D0To Dn;Making all data initial values is 0;
Step S12: read sample values S0To Sn;
Step S13: read current system frequency value Fc, work as Fc=Fs, enter step S16;
Step S14: as system ongoing frequency Fc< 45Hz or Fc> 55Hz time, enter step S16;
Step S15: as system ongoing frequency 45Hz≤FcDuring≤55Hz, according to FcValue call curable offline
Interpolation parameter table, use following interpolation algorithm D0=S0;D1=(1-α) × S1+α×S2;
D2=(1-2 α) × S2+2α×S3…Dn=(1-n α) × Sn+nα×Sn+1, calculate N number of interpolation number successively
According to, enter step S17;
Step S16:D0=S0;D1=S1;D2=S2;……Dn=Sn;Enter step S17;
Step S17: sample count SAM_CNT adds 1;Return step S11.
The interpolation method being applicable to the calculating of power system high-speed sampling the most according to claim 1, it is special
Levy and be, described sampling number N, system frequency range FL≤Fc≤FHAnd frequency accuracy set Δ F can basis
Application request changes, and first according to setting value off-line calculation interpolation parameter table T before program execution, program is held
According to F during rowcJudged result directly invokes the interpolation parameter table of curable offline, thus quickly obtains
Required interpolated data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410098732.9A CN103823996B (en) | 2014-03-18 | 2014-03-18 | It is applicable to the interpolation method that power system high-speed sampling calculates |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410098732.9A CN103823996B (en) | 2014-03-18 | 2014-03-18 | It is applicable to the interpolation method that power system high-speed sampling calculates |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103823996A CN103823996A (en) | 2014-05-28 |
CN103823996B true CN103823996B (en) | 2016-09-07 |
Family
ID=50759052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410098732.9A Active CN103823996B (en) | 2014-03-18 | 2014-03-18 | It is applicable to the interpolation method that power system high-speed sampling calculates |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103823996B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391674B (en) * | 2014-10-22 | 2017-09-05 | 积成电子股份有限公司 | Sampled value linear interpolation arithmetic device and operation method based on FPGA |
CN104808038B (en) * | 2015-04-03 | 2017-12-26 | 王顺江 | A kind of method of soft simulation interpolation lifting AC sampling frequency |
CN107037759B (en) * | 2017-05-11 | 2019-07-23 | 国网江西省电力公司电力科学研究院 | A kind of power electronic equipment fixes points sampler with the monocycle |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2314477A (en) * | 1996-06-19 | 1997-12-24 | Quantel Ltd | Image magnification processing system employing non-linear interpolation |
JP2002237961A (en) * | 2001-02-13 | 2002-08-23 | Fuji Xerox Co Ltd | Color conversion processing unit |
CN101109718A (en) * | 2006-11-14 | 2008-01-23 | 北京国药恒瑞美联信息技术有限公司 | Virtual grid imaging method and system used for eliminating influence of scattered radiation |
CN101458329A (en) * | 2009-01-06 | 2009-06-17 | 清华大学 | Time-domain paralleling sampling rate transformation method |
CN101616249A (en) * | 2009-07-28 | 2009-12-30 | 谭洪舟 | Video scaling device based on bicubic interpolation |
CN102760281A (en) * | 2011-04-26 | 2012-10-31 | 撖龙 | Image resizing method |
-
2014
- 2014-03-18 CN CN201410098732.9A patent/CN103823996B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2314477A (en) * | 1996-06-19 | 1997-12-24 | Quantel Ltd | Image magnification processing system employing non-linear interpolation |
JP2002237961A (en) * | 2001-02-13 | 2002-08-23 | Fuji Xerox Co Ltd | Color conversion processing unit |
CN101109718A (en) * | 2006-11-14 | 2008-01-23 | 北京国药恒瑞美联信息技术有限公司 | Virtual grid imaging method and system used for eliminating influence of scattered radiation |
CN101458329A (en) * | 2009-01-06 | 2009-06-17 | 清华大学 | Time-domain paralleling sampling rate transformation method |
CN101616249A (en) * | 2009-07-28 | 2009-12-30 | 谭洪舟 | Video scaling device based on bicubic interpolation |
CN102760281A (en) * | 2011-04-26 | 2012-10-31 | 撖龙 | Image resizing method |
Also Published As
Publication number | Publication date |
---|---|
CN103823996A (en) | 2014-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
He et al. | A generalized design framework of notch filter based frequency-locked loop for three-phase grid voltage | |
Xia et al. | A complex least squares enhanced smart DFT technique for power system frequency estimation | |
CN103823996B (en) | It is applicable to the interpolation method that power system high-speed sampling calculates | |
CN105720601B (en) | The digital phase-locked loop system based on implicit PI for synchronized | |
US11939961B2 (en) | Method and system of positive and negative sequence rotor currents control for doubly-fed induction generator-based wind turbines under single dq-PI control structure | |
CN108011555B (en) | Permanent magnet synchronous motor model prediction current control method | |
CN105759119B (en) | SDFT fundamental positive sequences phase synchronization method and system | |
Ahmed et al. | Frequency adaptive parameter estimation of unbalanced and distorted power grid | |
Agustin et al. | A synchronized current difference updating technique for model-free predictive current control of PMSM drives | |
CN107124132B (en) | A kind of method for identification of rotational inertia and identifier | |
CN111625769A (en) | PMU-SCADA data time synchronization and fusion method based on Lagrange interpolation and cubic exponential smoothing | |
CN110727899A (en) | Power grid oscillation damping ratio calculation method | |
CN112985727A (en) | Control method, control system, medium and equipment of linear vibration table | |
Ko et al. | A new PLL system using full order observer and PLL system modeling in a single phase grid-connected inverter | |
CN108880357B (en) | Self-adaptive non-smooth current tracking control method of permanent magnet synchronous motor | |
Tian et al. | Robust estimation of voltage parameters via a sliding mode observer | |
JP2007288981A (en) | Phase synchronization controlling method and phase synchronization controller of power converter | |
CN104021314A (en) | Interpolation method suitable for high-speed sampling calculation of electrical power system | |
CN105337591A (en) | Circuit structure and method for realizing clock recovery on the basis of USB device | |
CN109768752B (en) | Permanent magnet synchronous motor dead-beat current prediction control method based on multipurpose disturbance observer | |
CN104484376B (en) | The method and device of displaying data in real-time | |
Zhao et al. | An enhanced transfer delay frequency locked loop method for single phase grid voltage synchronization | |
Deepak et al. | Development of a PC based multi-function recorder for the laboratory model power system on RTAI-Linux platform | |
Cacciato et al. | Multi-reference frame based PLL for single phase systems in voltage distorted grids | |
CN108871385A (en) | Encoder, motor, encoder data processing method and storage medium |
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 |