CN103983852A - Harmonic analysis method of power quality harmonic analyzer - Google Patents
Harmonic analysis method of power quality harmonic analyzer Download PDFInfo
- Publication number
- CN103983852A CN103983852A CN201410233481.0A CN201410233481A CN103983852A CN 103983852 A CN103983852 A CN 103983852A CN 201410233481 A CN201410233481 A CN 201410233481A CN 103983852 A CN103983852 A CN 103983852A
- Authority
- CN
- China
- Prior art keywords
- harmonic
- frequency
- amplitude
- electric signal
- described electric
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 230000001360 synchronised effect Effects 0.000 claims description 19
- 230000006872 improvement Effects 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 25
- 230000010354 integration Effects 0.000 description 21
- 238000005516 engineering process Methods 0.000 description 12
- 230000003750 conditioning effect Effects 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241001515806 Stictis Species 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Abstract
The invention discloses a harmonic analysis method of a power quality harmonic analyzer, wherein the harmonic analysis method is achieved based on the improvement on the quasi-synchronization DFT. The harmonic analysis method comprises the steps that an electrical signal X on which the harmonic analysis needs to be carried out is sampled in a power grid, and data of W+2 sampling points are sampled at equal intervals; when the harmonic analysis is carried out according to the quasi-synchronization DFT, the position of frequency domain sampling is changed according to the drifting value, seen in the specification, of the signal frequency and then the frequency domain sampling position confirms to the value in the specification; the amplitude, seen in the specification, of all harmonics of the electrical signal X and the initial phase angles, seen in the specification, of the harmonics of the electrical signals X are worked out; the amplitude percentage, seen in the specification, of the amplitude Pk of the higher harmonics of the voltage signal V and the current signal I to the amplitude of the fundamental harmonic P1 is worked out, and the phase angle difference seen in the specification is worked out, and the amplitude percentage and the phase angle difference are output and displayed.
Description
The application is that application number is: 201110322629.4, and invention and created name is " a kind of harmonic analysis method of electric energy quality harmonic analyser ", the applying date is: the divisional application of the application for a patent for invention on October 21st, 2011.
Technical field
The present invention relates to the technical field that electric energy quality harmonic is analyzed, specifically a kind of high-precision harmonic analysis method.
Background technology
Frequency analysis technology is widely used in various fields such as electric energy quality monitoring, electronic product production testing, electric equipment monitoring, is the important technical of carrying out power system monitor, quality inspection, monitoring of tools.
In power supply grid, the generation of harmonic phenomena is mainly due to large-capacity power equipment, cause with electric rectification or converter and other nonlinear-loads.Along with the develop rapidly of Power Electronic Technique, many industrial enterprises introduce a large amount of impact load, asymmetric load and nonlinear load in recent years, and they inevitably occur harmonic pollution phenomenon in the time of access electrical network.Therefore, voltage, current signal to power supply grid carry out frequency analysis, can observe better electric network pollution situation, and then provide theoretical foundation for power network compensation and purification.
The most widely used technology of frequency analysis is discrete Fourier transformation (DFT) and Fast Fourier Transform (FFT) (FFT) at present.The frequency analysis technology (that is: accurate synchronous DFT frequency analysis technology) that accurate synchronous sampling technique and DFT technology combine can improve the precision of frequency analysis, and its formula is:
，
In formula: k is the number of times (as first-harmonic k=1,3 subharmonic k=3) that needs the harmonic wave obtaining; Sin and cos are respectively sine and cosine functions; And a
kand b
kbe respectively real part and the imaginary part of k subharmonic; N is iterations; W determines by integration method, while adopting complexification trapezoidal integration method, and W=nN;
it is a weighting coefficient;
, be all weighting coefficient sums; I the sampled value that f (i) is analysis waveform; N is sampling number in the cycle.
In engineering application, frequency analysis is always carried out the sampling of finite point and is difficult to accomplish the synchronized sampling of stricti jurise.Like this, in the time that the accurate synchronous DFT of application carries out frequency analysis, the short scope leakage that will exist the long scope causing due to truncation effect to leak and cause due to fence effect, makes analysis result precision not high, even not credible.
Summary of the invention
The technical problem to be solved in the present invention is to provide the harmonic analysis method of the electric energy quality harmonic analyser that a kind of precision is higher, effectively to improve the analytical error of accurate synchronous DFT frequency analysis technology, obtain high-precision frequency analysis result, thereby improve the reliability of the electric energy quality monitoring equipment based on frequency analysis theory and the validity of state judgement.
For solving the problems of the technologies described above, the invention provides a kind of harmonic analysis method of electric energy quality harmonic analyser, described electric energy quality harmonic analyser comprises successively and to be connected: sensor, signal conditioning circuit, data acquisition circuit and CPU element; Described harmonic analysis method comprises the steps:
(1) in described CPU element sampling electrical network, need the electric signal X of frequency analysis, and equal interval sampling W+2 sampling number is according to { f (i), i=0,1, w+1}(W is determined by selected integration method, the present invention does not specify a certain integration method, conventional integration method has complexification trapezoidal integration method W=nN, complexification rectangular integration method W=n (N-1), iterative Simpson integration method W=n (N-1)/2 etc., and the actual conditions that can apply according to the present invention are selected suitable integration method.Generally more satisfactory with complexification trapezoidal integration method effect);
(2) CPU element starts the accurate synchronous DFT formula of application from the sampled point i=0 of described electric signal X:
Analyze W+1 data acquisition first-harmonic information
with
;
Then, CPU element starts the accurate synchronous DFT formula of application from the sampled point i=1 of described electric signal X:
Analyze W+1 data acquisition first-harmonic information
;
Application of formula
calculate the frequency drift of described electric signal X
;
Application of formula
, draw respectively the first-harmonic of described electric signal X and higher hamonic wave
;
(3), CPU element calculates respectively the first-harmonic of described electric signal X and the amplitude of higher hamonic wave
and initial phase angle
;
(5), CPU element calculates respectively the amplitude P of the higher hamonic wave of described electric signal X
kwith first-harmonic P
1amplitude number percent
and phase angle difference
, and output display (the general LCD unit being connected with this CPU element that adopts shows).
Described electric signal X is generally grid ac voltage signal or grid alternating current stream signal.
Accurate synchronous DFT frequency analysis can effectively suppress long scope and leak, the main cause of its spectrum leakage is that the short scope that signal frequency drift causes is leaked, and signal frequency drift causes principal character that short scope leaks is spectrum peak-to-peak value occurs position along with signal frequency drift synchronous change, so variable fence frequency domain sample can effectively catch according to signal drift the position that spectrum peak-to-peak value occurs, thereby obtains high-precision harmonic information.
Equal interval sampling is cycle T and the frequency f (if power frequency component frequency f is 50Hz, the cycle is 20mS) that basis is carried out the ideal signal of frequency analysis, the N point of sampling in one-period, and sample frequency is f
s=Nf, and N>=64.
Described sampling W+2 sampling number is according to being to do corresponding selection according to selected integration method, if adopt complexification trapezoidal integration method, W=nN; If adopt complexification rectangular integration method, W=n (N-1); If adopt iterative Simpson integration method, W=n (N-1)/2.Then according to sample frequency f
s=Nf, acquisition sampled point data sequence f (i), and i=0,1 ..., w+1}, n>=3, finally carry out frequency analysis to this data sequence.
An iteration coefficient
determined by integration method, ideal period sampled point N and iterations n, concrete derivation referring to document [Dai Xianzhong. the some problems [J] in accurate synchronized sampling application. electrical measurement and instrument, 1988, (2): 2-7.].
for all weighting coefficient sums.
for imaginary part and the real part of k subharmonic, according to
just can obtain harmonic amplitude and initial phase angle.
The drift of signal frequency
to obtain according to the fixed relationship of sampling number N in neighbouring sample point first-harmonic phase angle difference and ideal period, the drift of signal frequency
also can be used for revising the frequency f of first-harmonic and higher hamonic wave
1frequency f with higher hamonic wave
k.
The present invention is mainly used in the analysis software of electric energy quality harmonic analyser, for improvement of with improve the quality of frequency analysis, obtain amplitude, phase angle and the frequency of high-precision first-harmonic and higher hamonic wave.Analyzed signal is carried out sampling and analyzing by the present invention, and analysis result (amplitude, phase angle and frequency etc.) is exported according to the requirement of electric energy quality harmonic analyser.The harmonic analysis method that electric energy quality harmonic analyser of the present invention adopts, is the frequency analysis technology based on variable fence thought, has following technical advantage:
(1) high-precision frequency analysis result.No matter the analysis result that frequency analysis technology of the present invention obtains is that amplitude or phase angle error improve 4 more than the order of magnitude.
(2) frequency analysis technology of the present invention has fundamentally solved the low problem of accurate synchronous DFT analysis precision, and without carrying out complicated inverting and correction, algorithm is simple.
(3) with respect to the synchronous DFT of standard, frequency analysis technology of the present invention only need to increase a sampled point and just solve the large problem of accurate synchronous DFT analytical error, is easy to realize.
(4) application the present invention improves existing instrument and equipment, and technical is feasible, and does not need to increase any hardware spending and just can make analysis result can improve 4 more than the order of magnitude.
(5) variable fence thought is also applicable to carry out repeatedly iteration and the frequency analysis process of non-once iteration too, now only an iteration need to be resolved into repeatedly to iteration and realize just passable.Iteration is the same with iteration repeatedly in essence, and just in the time calculating, repeatedly iteration is carried out substep calculating, and iteration is that the process of iteration is repeatedly merged to iteration coefficient
in once calculated, so the present invention is equally applicable to repeatedly iterative process.
(6) the present invention applies the position that accurate synchronous DFT carries out the sampling of frequency analysis time-frequency domain and changes according to the drift of signal frequency, and described frequency domain sampling position is μ 2 π/N, wherein:
for the drift of signal frequency, during without drift
be 1.The thought of harmonic analysis method of the present invention based on variable fence, realizes by 5 analytical procedures.The thought of variable fence: the main cause of accurate synchronous DFT analytical error is position and the ideal position generation deviation that the drift of signal frequency causes spectrum peak to occur, if the analysis result still obtaining to sample in frequency domain according to 2 π/N is extremely incorrect.Variable fence refers to: the position of frequency domain sampling be not the 2 π/N fixing, but change according to the drift of signal frequency, frequency sampling position is the drift that μ 2 π/N(μ is signal frequency).Frequency domain sampling fence changes along with the drift of signal frequency can accurately estimate the position that higher hamonic wave peak value occurs, and then obtains high-precision amplitude and phase angle information.
Embodiment
Electric energy quality harmonic analyser is mainly made up of sensor, signal conditioning circuit, data acquisition circuit, CPU, LCD unit and corresponding analysis software.Sensor generally adopts voltage transformer (VT) and current transformer, and the high voltage of tested electrical network and large current signal are changed by a certain percentage; Voltage, the current signal that signal conditioning circuit is sent sensor here is transformed to the voltage signal that is applicable to data acquisition circuit sampling; By data acquisition circuit, the signal after conditioning is sampled, then complete analysis and the detection to harmonic wave by CPU by analysis software.
The harmonic analysis method of electric energy quality harmonic analyser of the present invention, comprises the following steps:
First, in CPU element sampling electrical network, need the electric signal X of frequency analysis, and an equal interval sampling W+2 sampled point, with obtain described electric signal X discrete series f (k), k=0,1 ..., w+1}.W is determined jointly by sampling number N in integration method, iterations n and ideal period.Equal interval sampling refers to according to the frequency f (if power frequency component frequency is 50Hz, the cycle is 20mS) of the ideal signal that carries out frequency analysis determines sample frequency f
s=Nf, at sample frequency f
seffect under the N point of sampling equably in one-period.Usually, periodic sampling point N=64 or just can obtain good frequency analysis result above, and iterations n=3-5 just can obtain comparatively ideal frequency analysis result.Integration method has complexification trapezoidal integration method W=nN, complexification rectangular integration method W=n (N-1), Simpson's integration method W=n (N-1)/2 etc. multiple, can select according to actual conditions.
Secondly, CPU element starts the accurate synchronous DFT formula of application from the sampled point i=0 of described electric signal X:
analyze W+1 data acquisition first-harmonic information
with
.Wherein, an iteration coefficient
determined by integration method, ideal period sampled point N and iterations n, and
for all weighting coefficient sums.
Then, CPU element starts the accurate synchronous DFT formula of application from the sampled point i=1 of described electric signal X:
analyze W+1 data acquisition first-harmonic information
.
Then, application of formula
calculate the frequency drift of described electric signal X
.Obtain frequency drift
after, can be according to sample frequency f
scalculate and obtain the first-harmonic of analyzed signal and the frequency f of higher hamonic wave with sampling number N in ideal period.
Then, application of formula
, calculate respectively the first-harmonic of described electric signal X and the real part of higher hamonic wave
and imaginary part
;
CPU element and then according to formula:
calculate respectively the first-harmonic of described electric signal X and the amplitude of higher hamonic wave
and initial phase angle:
.
Finally, CPU element calculates respectively the amplitude P of the higher hamonic wave of described electric signal X
kwith first-harmonic P
1amplitude number percent
and phase angle difference
, and output display (the general LCD unit being connected with this CPU element that adopts shows).
Those skilled in the art will be appreciated that, above embodiment is only for the present invention is described, and not as limitation of the invention, the present invention can also be varied to more mode, as long as within the scope of connotation of the present invention, variation, modification to the above embodiment all will drop within the scope of claims of the present invention.
Claims (1)
1. a harmonic analysis method for electric energy quality harmonic analyser, is characterized in that comprising the following steps:
(1), need the electric signal X of frequency analysis equal interval sampling W+2 sampling number certificate in the CPU element of electric energy quality harmonic analyser sampling electrical network: f (i), i=0,1 ..., w+1};
(2), CPU element starts the accurate synchronous DFT formula of application from the sampled point i=0 of described electric signal X:
analyze W+1 data acquisition first-harmonic information
with
;
Then, CPU element starts the accurate synchronous DFT formula of application from the sampled point i=1 of described electric signal X:
, analyze W+1 data acquisition first-harmonic information
;
Application of formula:
,
Calculate the frequency drift of described electric signal X
, frequency drift
to obtain according to the fixed relationship of sampling number N in neighbouring sample point first-harmonic phase angle difference and ideal period, frequency drift
for revising the frequency f of first-harmonic
1frequency f with higher hamonic wave
k;
Application of formula:
, draw respectively the first-harmonic of described electric signal X and higher hamonic wave
;
(3), CPU element calculates respectively the first-harmonic of described electric signal X and the amplitude of higher hamonic wave
and initial phase angle
(4), CPU element calculates respectively the amplitude P of the higher hamonic wave of described electric signal X
kwith first-harmonic P
1amplitude number percent
and phase angle difference
, and output display.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410233481.0A CN103983852B (en) | 2011-10-21 | 2011-10-21 | Harmonic analysis method of power quality harmonic analyzer |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410233481.0A CN103983852B (en) | 2011-10-21 | 2011-10-21 | Harmonic analysis method of power quality harmonic analyzer |
| CN201110322629.4A CN102508026B (en) | 2011-10-21 | 2011-10-21 | Harmonic wave analysis method for electric energy quality harmonic wave analyzer |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201110322629.4A Division CN102508026B (en) | 2011-10-21 | 2011-10-21 | Harmonic wave analysis method for electric energy quality harmonic wave analyzer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103983852A true CN103983852A (en) | 2014-08-13 |
| CN103983852B CN103983852B (en) | 2017-04-12 |
Family
ID=51275898
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410233481.0A Active CN103983852B (en) | 2011-10-21 | 2011-10-21 | Harmonic analysis method of power quality harmonic analyzer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103983852B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104898041A (en) * | 2015-06-30 | 2015-09-09 | 绵阳市维博电子有限责任公司 | Phase-sensitive rail signal detecting method and system |
| CN109061345A (en) * | 2018-07-13 | 2018-12-21 | 许继集团有限公司 | Effective value measurement method and system suitable for electric system |
| CN109239463A (en) * | 2015-05-19 | 2019-01-18 | 江苏理工学院 | A kind of dielectric loss measurement method based on linear correction algorithm |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006276006A (en) * | 2005-03-01 | 2006-10-12 | Nagoya Institute Of Technology | Harmonic analysis method in power system |
| CN101900761A (en) * | 2009-11-05 | 2010-12-01 | 中国航天科技集团公司第五研究院第五一四研究所 | High-accuracy non-integer-period sampled harmonic analysis and measurement method |
| CN102128982A (en) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | Harmonic analysis method based on windowed interpolation FFT (Fast Fourier Transform) base frequency tracking technology |
| CN102508026B (en) * | 2011-10-21 | 2014-08-06 | 江苏理工学院 | Harmonic wave analysis method for electric energy quality harmonic wave analyzer |
-
2011
- 2011-10-21 CN CN201410233481.0A patent/CN103983852B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006276006A (en) * | 2005-03-01 | 2006-10-12 | Nagoya Institute Of Technology | Harmonic analysis method in power system |
| CN101900761A (en) * | 2009-11-05 | 2010-12-01 | 中国航天科技集团公司第五研究院第五一四研究所 | High-accuracy non-integer-period sampled harmonic analysis and measurement method |
| CN102128982A (en) * | 2010-12-23 | 2011-07-20 | 东莞市创锐电子技术有限公司 | Harmonic analysis method based on windowed interpolation FFT (Fast Fourier Transform) base frequency tracking technology |
| CN102508026B (en) * | 2011-10-21 | 2014-08-06 | 江苏理工学院 | Harmonic wave analysis method for electric energy quality harmonic wave analyzer |
Non-Patent Citations (3)
| Title |
|---|
| 张建秋等: "非整周期采样应用于周期信号的谐波分析", 《仪器仪表学报》, vol. 16, no. 1, 28 February 1995 (1995-02-28), pages 1 - 7 * |
| 李芙英等: "用准同步离散Fourier变换实现高准确度谐波分析", 《清华大学学报》, vol. 39, no. 5, 31 May 1999 (1999-05-31), pages 47 - 50 * |
| 沈国峰; 王祁: "进一步提高准同步采样谐波分析法准确度的方案", 《仪器仪表学报》, vol. 22, no. 5, 30 October 2001 (2001-10-30), pages 456 - 465 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109239463A (en) * | 2015-05-19 | 2019-01-18 | 江苏理工学院 | A kind of dielectric loss measurement method based on linear correction algorithm |
| CN109239463B (en) * | 2015-05-19 | 2020-07-31 | 江苏理工学院 | Dielectric loss measurement method based on linear correction algorithm |
| CN104898041A (en) * | 2015-06-30 | 2015-09-09 | 绵阳市维博电子有限责任公司 | Phase-sensitive rail signal detecting method and system |
| CN109061345A (en) * | 2018-07-13 | 2018-12-21 | 许继集团有限公司 | Effective value measurement method and system suitable for electric system |
| CN109061345B (en) * | 2018-07-13 | 2020-09-18 | 许继集团有限公司 | Effective value measuring method and system suitable for power system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103983852B (en) | 2017-04-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102508026B (en) | Harmonic wave analysis method for electric energy quality harmonic wave analyzer | |
| CN102435815B (en) | Operating method of resistive current on-line monitoring system of metal oxide arrester (MOA) | |
| CN102393488B (en) | Harmonic analysis method | |
| CN102435860B (en) | Working method of medium loss current tester | |
| CN103439566B (en) | Operating method of MOA resistive current tester with relatively high precision | |
| CN103869162A (en) | Dynamic signal phasor measurement method based on time domain quasi-synchronization | |
| Liu et al. | The evaluation of phasor measurement units and their dynamic behavior analysis | |
| CN106018956B (en) | A kind of power system frequency computational methods of adding window spectral line interpolation | |
| CN102879639A (en) | Real-time frequency measuring method in power system | |
| CN103257273A (en) | Method for measuring common-frequency periodic signal phase differences | |
| CN109061345B (en) | Effective value measuring method and system suitable for power system | |
| CN108896944B (en) | Laboratory calibrator of synchronous measuring device and synchronous phasor measuring method thereof | |
| CN106226590A (en) | A kind of synchronous phase measuring in power system method | |
| Jin et al. | A novel power harmonic analysis method based on Nuttall-Kaiser combination window double spectrum interpolated FFT algorithm | |
| CN109030957B (en) | Dielectric loss measuring method | |
| CN103543331B (en) | A kind of method calculating electric signal harmonic wave and m-Acetyl chlorophosphonazo | |
| CN103983852A (en) | Harmonic analysis method of power quality harmonic analyzer | |
| CN103592512A (en) | Harmonic analysis method of power quality harmonic analyzer | |
| CN103604989A (en) | Harmonic wave analysis method of power-quality harmonic wave analyzer | |
| CN102435816B (en) | Working method of MOA (Metal Oxide Arrester) resistive current tester | |
| CN103969507A (en) | Power quality harmonic analysis method | |
| CN104678170A (en) | Power harmonic analysis method based on harmonic analyzer and harmonic analyzer | |
| CN104849530A (en) | MOA resistive current cardinal wave measuring method | |
| CN109030942B (en) | Harmonic phase angle analysis method | |
| CN103267896B (en) | The measuring method of initial phase angle of periodic signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20230315 Address after: 201600 Floor 5, Building 5, No.168 Shuya Road, Songjiang District, Shanghai Patentee after: Shanghai Liudian Electric Co.,Ltd. Address before: 230000 Room 203, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province Patentee before: Hefei Jiuzhou Longteng scientific and technological achievement transformation Co.,Ltd. Effective date of registration: 20230315 Address after: 230000 Room 203, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province Patentee after: Hefei Jiuzhou Longteng scientific and technological achievement transformation Co.,Ltd. Address before: 213001, No. 1801, Wu Cheng Road, bell tower, Changzhou, Jiangsu Patentee before: JIANGSU University OF TECHNOLOGY |