CN101701984A - Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT - Google Patents
Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT Download PDFInfo
- Publication number
- CN101701984A CN101701984A CN200910154681A CN200910154681A CN101701984A CN 101701984 A CN101701984 A CN 101701984A CN 200910154681 A CN200910154681 A CN 200910154681A CN 200910154681 A CN200910154681 A CN 200910154681A CN 101701984 A CN101701984 A CN 101701984A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- coefficient
- interpolation
- fft
- amplitude
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000003595 spectral effect Effects 0.000 claims abstract description 27
- 238000001228 spectrum Methods 0.000 claims description 5
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 claims description 3
- 238000005070 sampling Methods 0.000 abstract description 3
- 238000009472 formulation Methods 0.000 abstract 2
- 239000000203 mixture Substances 0.000 abstract 2
- 238000001514 detection method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
The invention discloses a fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT. The method comprises the following steps of: obtaining a sampling data of a detected signal by an analog-to-digital converter; carrying out a three-coefficient Nuttall windowed FFT calculation on the sampling data; searching for a largest spectral line and an adjacent second largest spectral line of an amplitude spectral line for the fundamental wave and each harmonic in an FFT spectral line, according to the amplitude ratio of adjacent spectral peaks and an explicit formulation, directly calculating interpolation coefficients of the fundamental wave and each harmonic; and finally, obtaining the frequency, the amplitude and the phase of the fundamental wave and each harmonic by the interpolation calculation. The invention provides the three-coefficient Nuttall windowed interpolation FFT fundamental wave and harmonic detecting method for directly calculating the interpolation coefficients based on the explicit formulation. The invention has advantages of providing the windowed interpolation FFT fundamental wave and harmonic wave detecting method thereof with small calculated amount and high analysis precision.
Description
Technical field
The present invention relates to first-harmonic and harmonic wave detection range in a kind of signal, relate in particular to a kind of first-harmonic and harmonic detecting method based on three coefficient Nuttall window interpolation FFTs.
Background technology
With the electric power signal detection is example, and Harmonious Waves in Power Systems influences the normal operation of power equipment, and it is the basic demand of electrical network and power equipment safety stable operation that first-harmonic and harmonic wave are accurately detected.Windowed interpolation FFT (Fast Fourier Transform (FFT)) method is the common method of Measurement of Harmonics in Power System.When non-synchronous sampling, Fourier transform exists spectrum leakage and fence effect.Windowed interpolation FFT suppresses spectrum leakage by the windowing computing, eliminates the influence of fence effect by interpolation arithmetic.Windowed interpolation FFT method Cosine Window function commonly used has Hanning window, Blackman window, Blackman-Harris window, Rife-Vincent window, optimum Cosine Window etc.
Existing patent document " based on the first-harmonic and the harmonic detecting method of Nuttall window double peak interpolation FFT " (200710035653.3), the purpose of its invention is to provide the double peak interpolation FFT method based on the Nuttall window, this method is used four three rank Nuttall windows, adopt the double peak interpolation algorithm, do not have the interpolation coefficient computing formula of explicit direct calculating, calculated amount is big.
The weak point of prior art is, and is little based on the windowed interpolation FFT method calculated amount of the few window function of item number, have explicit interpolation coefficient computing formula, but analysis precision is low; Based on the windowed interpolation FFT methods analyst precision height of the many window functions of item number, but need find the solution repeatedly equation, can't directly calculate interpolation coefficient by explicit expression, calculated amount is big.For example, the Hanning window is two coefficient Cosine Window, and its interpolation formula is explicit computing formula, and is simple and clear, can directly calculate, and calculated amount is little, and computational stability is good, but analysis precision is lower; The item number of Blackman-Harris window and four three rank Nuttall window functions is all more than the Hanning window, analysis precision is higher, but when calculating, interpolation coefficient need find the solution repeatedly equation or fitting of a polynomial approaches, can't directly calculate interpolation coefficient, calculated amount is big, may cause effectively separating of equation not exist at noise with under disturbing.
Summary of the invention
Purpose of the present invention is at the deficiency of above-mentioned technology, a kind of first-harmonic and harmonic detecting method based on three coefficient Nuttall window interpolation FFTs is provided, and it has explicit interpolation coefficient computing formula, and calculated amount is little, computational stability is good, the accuracy of detection height of first-harmonic and harmonic wave.
First-harmonic and harmonic detecting method based on three coefficient Nuttall window interpolation FFTs comprise the steps:
1) by analog to digital converter with sample frequency f
sTested voltage and current signal is converted into digital signal from simulating signal, obtains the sampled data of N point length;
2) three coefficient Nuttall window functions of structure N point length add three coefficient Nuttall window FFT computing to the sampled data of N point length, obtain FFT spectral line X (k), k=0, and 1 ..., N, wherein the data length of FFT computing is N;
3) first-harmonic and each harmonic are searched for maximum and adjacent big spectral line of amplitude spectral line in the FFT spectral line, directly calculated the interpolation coefficient of first-harmonic and each harmonic according to the ratio of the amplitude of adjacent spectral peaks by explicit expression;
4) obtain frequency, amplitude and the phase place of first-harmonic and each harmonic by interpolation arithmetic.
Above-mentioned steps 2) in, the building method of three coefficient Nuttall window functions is:
A wherein
0=0.375, a
1=0.5, a
2=0.125.
Above-mentioned steps 3) in, the explicit expression that calculates first-harmonic and each harmonic interpolation coefficient is:
In the formula, δ
mBe the interpolation coefficient of m subharmonic,
It is the ratio of the amplitude of adjacent maximum of m subharmonic and time big spectral line.
Above-mentioned steps 4) in, the interpolation arithmetic formula that calculates frequency, amplitude and the phase place of first-harmonic and each harmonic is:
f
m=(k
m+δ
m)f
s/N (3)
In the formula, f
m, A
m,
Be respectively frequency, amplitude and the phase place of m subharmonic, f
sBe sample frequency, the phase place that frequency spectrum is got in arg () expression.
The present invention proposes directly to calculate based on explicit expression three the coefficient Nuttall window FFT first-harmonics and the harmonic detecting method of interpolation coefficient first, has reduced calculated amount, has improved computing stability, has improved the analysis precision of windowed interpolation FFT.Advantage of the present invention is: 1, calculated amount of the present invention is little, and the computing good stability owing to adopted explicit interpolation coefficient to calculate, need not to find the solution repeatedly equation, at noise with under disturbing good applicability is arranged also, is easy to single-chip microcomputer or Implementation of Embedded System; 2, accuracy of detection height of the present invention is because the maximum secondary lobe of three coefficient Nuttall windows is-47dB that the rate of decay of each secondary lobe is 30dB, so the accuracy of detection height of first-harmonic and harmonic wave.
Description of drawings
Fig. 1 is based on the first-harmonic of three coefficient Nuttall window interpolation FFTs and the block diagram of harmonic detecting method.
Embodiment
First-harmonic and harmonic detecting method based on three coefficient Nuttall window interpolation FFTs comprise the steps:
1) by analog to digital converter with sample frequency f
sTested voltage and current signal is converted into digital signal from simulating signal, obtains the sampled data of N point length;
2) three coefficient Nuttall window functions of structure N point length add three coefficient Nuttall window FFT computing to the sampled data of N point length, obtain FFT spectral line X (k), k=0, and 1 ..., N, wherein the data length of FFT computing is N;
3) first-harmonic and each harmonic are searched for maximum and adjacent big spectral line of amplitude spectral line in the FFT spectral line, directly calculated the interpolation coefficient of first-harmonic and each harmonic according to the ratio of the amplitude of adjacent spectral peaks by explicit expression;
4) obtain frequency, amplitude and the phase place of first-harmonic and each harmonic by interpolation arithmetic.
Above-mentioned steps 2) in, the building method of three coefficient Nuttall window functions is:
A wherein
0=0.375, a
1=0.5, a
2=0.125.
Above-mentioned steps 3) in, the explicit expression that calculates first-harmonic and each harmonic interpolation coefficient is:
In the formula, δ
mBe the interpolation coefficient of m subharmonic,
It is the ratio of the amplitude of adjacent maximum of m subharmonic and time big spectral line.
Above-mentioned steps 4) in, the interpolation arithmetic formula that calculates frequency, amplitude and the phase place of first-harmonic and each harmonic is:
f
m=(k
m+δ
m)f
s/N (3)
In the formula, f
m, A
m,
Be respectively frequency, amplitude and the phase place of m subharmonic, f
sBe sample frequency, the phase place that frequency spectrum is got in arg () expression.
Embodiment 1
It is example that current harmonics during with certain electrical work detects, and the current expression of establishing this electrical equipment is
Its setting value is as shown in table 1.Application the present invention is based on the first-harmonic and the harmonic detecting method of three coefficient Nuttall window interpolation FFTs and measures its first-harmonic and 2~9 subharmonic (but the present invention is not limited to 2~9 subharmonic), by the Matlab simulation software process of executing in fact is described in the present embodiment.
(1) obtain the sampled data of this electric current by analog to digital converter, wherein analog to digital converter adopts U.S. letter MAX125CEAX integrated circuit (IC) chip, sample frequency f
s=10kHz, data length N=2048.Matlab software increases the white Gaussian noise of 80dB and measures noise to represent it in original signal.
(2) structure 2048 3 coefficient Nuttall windows, to the sampled data windowing, and carry out 2048 FFT conversion, obtain 1024 FFT spectral lines, be designated as X (0), X (1) ..., X (1023).
(3) first-harmonic and each harmonic are searched for maximum and adjacent big spectral line of amplitude spectral line in spectral line FFT, obtained: k
1=10, k
2=20, k
3=30, k
4=40, k
5=51, k
6=61, k
7=71, k
8=81, k
9=91; Interpolation coefficient according to formula (2) calculating first-harmonic and each harmonic obtains: δ
1=0.2195, δ
2=0.4398, δ
3=0.6585, δ
4=0.8777, δ
5=0.0977, δ
6=0.3167, δ
7=0.5363, δ
8=0.7565, δ
9=0.9757.
(4) obtain frequency, amplitude and the phase place of first-harmonic and each harmonic at last by interpolation arithmetic formula (3), (4) and (5), the result is as shown in table 1.
Table 1 embodiment testing result
In sum, the present invention is based on the first-harmonic of three coefficient Nuttall window interpolation FFTs and the interpolation coefficient that harmonic detecting method can directly calculate first-harmonic and each harmonic, calculated amount is little, computing good stability, the accuracy of detection height of first-harmonic and harmonic wave.
Claims (4)
1. first-harmonic and harmonic detecting method based on three coefficient Nuttall window interpolation FFTs is characterized in that comprising the steps:
1) by analog to digital converter with sample frequency f
sTested voltage and current signal is converted into digital signal from simulating signal, obtains the sampled data of N point length;
2) three coefficient Nuttall window functions of structure N point length add three coefficient Nuttall window FFT computing to the sampled data of N point length, obtain FFT spectral line X (k), k=0, and 1 ..., N, wherein the data length of FFT computing is N;
3) first-harmonic and each harmonic are searched for maximum and adjacent big spectral line of amplitude spectral line in the FFT spectral line, directly calculated the interpolation coefficient of first-harmonic and each harmonic according to the ratio of the amplitude of adjacent spectral peaks by explicit expression;
4) obtain frequency, amplitude and the phase place of first-harmonic and each harmonic by interpolation arithmetic.
2. method according to claim 1 is characterized in that: step 2) in, the building method of described three coefficient Nuttall window functions is:
A wherein
0=0.375, a
1=0.5, a
2=0.125.
3. method according to claim 1 is characterized in that: in the step 3), the explicit expression of described calculating first-harmonic and each harmonic interpolation coefficient is:
4. method according to claim 1 is characterized in that: in the step 4), the interpolation arithmetic formula of the frequency of described first-harmonic and each harmonic, amplitude and phase place is:
f
m=(k
m+δ
m)f
s/N (3)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101546816A CN101701984B (en) | 2009-11-23 | 2009-11-23 | Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101546816A CN101701984B (en) | 2009-11-23 | 2009-11-23 | Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101701984A true CN101701984A (en) | 2010-05-05 |
CN101701984B CN101701984B (en) | 2011-05-18 |
Family
ID=42156904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101546816A Expired - Fee Related CN101701984B (en) | 2009-11-23 | 2009-11-23 | Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101701984B (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243272A (en) * | 2011-04-01 | 2011-11-16 | 重庆大学 | High-precision harmonic analysis method of sampling data synchronization |
CN102253282A (en) * | 2011-06-10 | 2011-11-23 | 中国矿业大学 | Method for obtaining continuous frequency spectrum interpolation power harmonic parameter of Nuttall window function |
CN102331526A (en) * | 2011-06-10 | 2012-01-25 | 中国矿业大学 | Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation |
CN104062528A (en) * | 2014-07-04 | 2014-09-24 | 武汉大学 | Signal harmonic analysis method and system based on Hanning product window |
CN104391178A (en) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | Time shift phase difference steady harmonic signal correction method based on Nuttall window |
CN104749433A (en) * | 2015-01-04 | 2015-07-01 | 云南电网公司电力科学研究院 | Method for accurately calculating and analyzing harmonic and electric energy quality |
CN104833937A (en) * | 2015-05-21 | 2015-08-12 | 湖南大学 | Harmonic measurement channel calibration method based on MIR-RSD high-precision cosine window interpolation FFT algorithm |
CN105137181A (en) * | 2015-09-15 | 2015-12-09 | 福州大学 | Double-spectrum-line interpolation harmonic analysis algorithm based on Nuttall-Kaiser composite window |
CN105137183A (en) * | 2015-09-15 | 2015-12-09 | 湖北工业大学 | Analysis method and analysis system of harmonious waves in power systems |
CN105353215A (en) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | Harmonic detection method based on Nuttall window four-spectral-line interpolation FFT (fast Fourier transform) |
CN105717359A (en) * | 2016-02-19 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | Harmonic analysis algorithm |
CN107271002A (en) * | 2017-06-19 | 2017-10-20 | 重庆邮电大学 | A kind of Spectrum Correction interpolation algorithm of quick high accuracy |
CN109782063A (en) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | A kind of dynamic m-Acetyl chlorophosphonazo analysis method based on three spectral line interpolation FFT of Nuttall self-convolution window |
CN110007148A (en) * | 2019-03-28 | 2019-07-12 | 东南大学 | A kind of simple signal frequency estimating methods based on the comprehensive interpolation of discrete spectrum phase and amplitude |
CN114659791A (en) * | 2022-02-28 | 2022-06-24 | 广东机电职业技术学院 | Steam turbine fault detection method, system, device and storage medium |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SE513892C2 (en) * | 1995-06-21 | 2000-11-20 | Ericsson Telefon Ab L M | Spectral power density estimation of speech signal Method and device with LPC analysis |
CN101113995A (en) * | 2007-08-29 | 2008-01-30 | 湖南大学 | Base wave and harmonic detecting method based on Nuttall window double peak interpolation FFT |
CN101261618A (en) * | 2008-04-14 | 2008-09-10 | 湖南大学 | Triangular self-convolution window weighting algorithm |
-
2009
- 2009-11-23 CN CN2009101546816A patent/CN101701984B/en not_active Expired - Fee Related
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102243272A (en) * | 2011-04-01 | 2011-11-16 | 重庆大学 | High-precision harmonic analysis method of sampling data synchronization |
CN102253282A (en) * | 2011-06-10 | 2011-11-23 | 中国矿业大学 | Method for obtaining continuous frequency spectrum interpolation power harmonic parameter of Nuttall window function |
CN102331526A (en) * | 2011-06-10 | 2012-01-25 | 中国矿业大学 | Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation |
CN102331526B (en) * | 2011-06-10 | 2013-07-17 | 中国矿业大学 | Method for acquiring parameters of electric power harmonic waves by using Hanniing window function continuous frequency spectrum interpolation |
CN102253282B (en) * | 2011-06-10 | 2014-07-16 | 中国矿业大学 | Method for obtaining continuous frequency spectrum interpolation power harmonic parameter of Nuttall window function |
CN104062528A (en) * | 2014-07-04 | 2014-09-24 | 武汉大学 | Signal harmonic analysis method and system based on Hanning product window |
CN104391178B (en) * | 2014-12-05 | 2017-12-15 | 国家电网公司 | A kind of time shift phase difference stable state harmonic signal bearing calibration based on Nuttall windows |
CN104391178A (en) * | 2014-12-05 | 2015-03-04 | 国家电网公司 | Time shift phase difference steady harmonic signal correction method based on Nuttall window |
CN104749433A (en) * | 2015-01-04 | 2015-07-01 | 云南电网公司电力科学研究院 | Method for accurately calculating and analyzing harmonic and electric energy quality |
CN104833937B (en) * | 2015-05-21 | 2017-08-11 | 湖南大学 | A kind of harmonic measure channel calibration method based on the high-precision Cosine Window interpolated FFT algorithms of MIR RSD |
CN104833937A (en) * | 2015-05-21 | 2015-08-12 | 湖南大学 | Harmonic measurement channel calibration method based on MIR-RSD high-precision cosine window interpolation FFT algorithm |
CN105137183A (en) * | 2015-09-15 | 2015-12-09 | 湖北工业大学 | Analysis method and analysis system of harmonious waves in power systems |
CN105137181A (en) * | 2015-09-15 | 2015-12-09 | 福州大学 | Double-spectrum-line interpolation harmonic analysis algorithm based on Nuttall-Kaiser composite window |
CN105353215A (en) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | Harmonic detection method based on Nuttall window four-spectral-line interpolation FFT (fast Fourier transform) |
CN105717359A (en) * | 2016-02-19 | 2016-06-29 | 云南电网有限责任公司电力科学研究院 | Harmonic analysis algorithm |
CN107271002A (en) * | 2017-06-19 | 2017-10-20 | 重庆邮电大学 | A kind of Spectrum Correction interpolation algorithm of quick high accuracy |
CN109782063A (en) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | A kind of dynamic m-Acetyl chlorophosphonazo analysis method based on three spectral line interpolation FFT of Nuttall self-convolution window |
CN110007148A (en) * | 2019-03-28 | 2019-07-12 | 东南大学 | A kind of simple signal frequency estimating methods based on the comprehensive interpolation of discrete spectrum phase and amplitude |
CN114659791A (en) * | 2022-02-28 | 2022-06-24 | 广东机电职业技术学院 | Steam turbine fault detection method, system, device and storage medium |
CN114659791B (en) * | 2022-02-28 | 2023-07-04 | 广东机电职业技术学院 | Turbine fault detection method, system, device and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN101701984B (en) | 2011-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101701984B (en) | Fundamental wave and harmonic wave detecting method based on three-coefficient Nuttall windowed interpolation FFT | |
CN101701982B (en) | Method for detecting harmonic waves of electric system based on window and interpolated FFT | |
CN101603985B (en) | Method for measuring sine signal with high accuracy | |
CN103454497B (en) | Based on the method for measuring phase difference improving windowed DFT | |
CN106483374B (en) | A kind of harmonic wave harmonic detection method based on Nuttall double window whole phase FFT | |
CN102288807B (en) | Method for measuring electric network voltage flicker | |
CN103245832B (en) | Based on harmonic wave time-frequency characteristic method for parameter estimation and the analyser of quick S-transformation | |
CN101718816B (en) | Fundamental wave and harmonic wave detection method based on four-item coefficient Nuttall window interpolation FFT | |
CN104897960B (en) | Harmonic wave rapid analysis method and system based on the spectral line interpolation FFT of adding window four | |
CN103308766A (en) | Harmonic analysis method based on Kaiser self-convolution window dual-spectrum line interpolation FFT (Fast Fourier Transform) and device thereof | |
CN103575984A (en) | Harmonic analysis method based on Kaiser window double-spectral-line interpolation FFT | |
CN101261292A (en) | Base wave and harmonic detection method based on fiver item Rife-Vincent(1)window double spectral line interpolation FFT | |
CN109782063A (en) | A kind of dynamic m-Acetyl chlorophosphonazo analysis method based on three spectral line interpolation FFT of Nuttall self-convolution window | |
CN103197141A (en) | Method of measuring electrical power system signal frequency and harmonic wave parameters | |
CN109946512B (en) | A kind of dynamic power analysis method for improving frequency domain interpolation | |
CN109541312A (en) | A kind of new energy collects the subsynchronous harmonic detecting method in area | |
CN103399204A (en) | Rife-Vincent (II) window interpolation FFT (Fast Fourier Transform)-based harmonic and inter-harmonic detection method | |
CN108535613A (en) | A kind of voltage flicker parameter detection method based on combination window function | |
CN105137181A (en) | Double-spectrum-line interpolation harmonic analysis algorithm based on Nuttall-Kaiser composite window | |
CN110133564A (en) | A kind of synchronous detecting method and device of current sensor phase characteristic | |
CN102590598B (en) | Method for predicting zero crossing point of ultrasonic signal based on multi-threshold comparison | |
CN103795411A (en) | SFDR testing method based on five-maximum-sidelobe-damping-window three-spectral-line interpolation | |
CN103941090A (en) | Harmonic wave measurement method based on spectral energy interpolation | |
CN108776263A (en) | Harmonic detecting method based on the peaceful self-convolution window of the high-order Chinese and improvement interpolation algorithm | |
CN109581045B (en) | Inter-harmonic power metering method meeting IEC standard framework |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110518 Termination date: 20111123 |