CN101394688B - Non-thermal resistance prediction method for power compression of coil driven loudspeaker - Google Patents

Non-thermal resistance prediction method for power compression of coil driven loudspeaker Download PDF

Info

Publication number
CN101394688B
CN101394688B CN2008101964795A CN200810196479A CN101394688B CN 101394688 B CN101394688 B CN 101394688B CN 2008101964795 A CN2008101964795 A CN 2008101964795A CN 200810196479 A CN200810196479 A CN 200810196479A CN 101394688 B CN101394688 B CN 101394688B
Authority
CN
China
Prior art keywords
voice coil
power
power compression
thermal resistance
temperature
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.)
Expired - Fee Related
Application number
CN2008101964795A
Other languages
Chinese (zh)
Other versions
CN101394688A (en
Inventor
沈勇
吴志成
王旭
王晓楠
孙飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University filed Critical Nanjing University
Priority to CN2008101964795A priority Critical patent/CN101394688B/en
Publication of CN101394688A publication Critical patent/CN101394688A/en
Application granted granted Critical
Publication of CN101394688B publication Critical patent/CN101394688B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Audible-Bandwidth Dynamoelectric Transducers Other Than Pickups (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

The invention relates to a prediction method for power compression of a coil driven loudspeaker without thermal resistance. For arbitrary coil driven loudspeaker, an impedance curve and the temperature variation of a voice coil are measured, the power compression value of a characteristic frequency point is calculated through the computing formula of the power compression without thermal resistance, then a modified coefficient is calculated through the modified formula, and finally power compression curves under different input power conditions are obtained. Based on the measured impedance curve and the temperature variation of the voice coil, the computing formula of the power compression without thermal resistance is adopted for calculating, and the computing formula is as follows: PC=10*log10[(ZTa + kRe DeltaT) (1+k DeltaT)/ZTa], wherein the power compression indicates the relative reducing of the frequency response amplitude of the coil driven loudspeaker under different input power conditions.

Description

Moving-coil speaker power compression without thermal resistance Forecasting Methodology
Technical field
The present invention relates to the Forecasting Methodology of moving-coil speaker power compression, is a kind of non-thermal resistance prediction method based on moving-coil speaker cell impedance curved measurement and temperature rise curve measurement of novelty.
Background technology
Loud speaker is the electroacoustic transducer that electrical signal conversion is become acoustical signal and medium radiation towards periphery.Loud speaker can be divided into moving-coil type, electromagnetic type, electrostatic, piezoelectric type, ionic formula etc. by the transducing mode, and wherein moving-coil speaker is widely used, and plays a significant role in life, practical applications field.
Moving-coil speaker has different temperature rise situation when under the different capacity initial conditions, working, the reason that temperature raises is the voice coil loudspeaker voice coil heating.Based on the temperature characterisitic of voice coil loudspeaker voice coil material, its direct current resistance meeting raises along with voice coil temperature and increases, and direct current resistance increase can cause reducing of input electric power, finally causes the reduction of sound output.The relative decrease of the frequency response of moving-coil speaker under the different capacity initial conditions is exactly the power compression value.
In order to measure the power compression value; Someone has proposed to calculate entire thermal resistance based on the quadratic equation that the hot equivalent model of linearity is found the solution about Δ T under the constant voltage initial conditions, and utilization and entire thermal resistance, temperature rise then, the power compression computing formula that input power is relevant are calculated the power compression value.But calculate so several drawbacks are arranged, the one, linear hot equivalent model lost efficacy in low-frequency range, and just there is very mistake in the thermal resistance of calculating with this method itself, and the power compression that calculates must be forbidden; The 2nd, this method thinks that voice coil loudspeaker voice coil is in same steady temperature at full frequency band, thereby the value that calculates and actual result deviation are very big.
The relevant method of having published at present:
1, based on the power compression computational methods of the hot equivalent model thermal resistance of linearity.
(1992) have proposed the computing formula of moving-coil speaker power compression,
PC = 10 Log 10 [ P e R 0 ( 1 + KΔ T ) Δ T ] , R 0Be entire thermal resistance, P eBe input electric power.
Calculate R 0Then need to obtain through separating following two formulas:
k 2ΔT 2R e+k 1R eΔT-V 2=0
V 2 = ΔT k 1 R e R 0 + k 2 Δ T 2 R e R 0
K wherein 1Be the ratio of loud speaker minimum impedance and the resistance of voice coil loudspeaker voice coil direct current, k 2Be temperature coefficient, R 0Be entire thermal resistance.
This method based on hot equivalent model lost efficacy at low frequency it be restricted.
Referring to DOUGLAS J.BUTTON.Heat dissipation and power compression in loudspeakers [J] .J.AudioEng.Soc., 1992,40 (1/2): 32-41.
2, calculate that based on direct current resistance variation temperature rise comes the rated output compression
(2002) have proposed to be similar to the temperature rise of reckoning voice coil loudspeaker voice coil based on the linear hot equivalent model in three rank, calculate the power compression value again.
ΔT ≈ Z min ( T a ) 2 δ R e ( T a ) ( 1 + 4 u 2 R t δ R e ( T a ) Z min 2 ( T a ) - 1 )
PC = 10 log [ R t ( 1 + δΔT ) ΔT ]
Wherein Δ T is the front and back temperature difference, and δ is a temperature coefficient, R tBe entire thermal resistance.
This method is based on the power compression calculating of thermal resistance equally, because itself is not based on the temperature data of actual measurement, last resultant error is very big, has method 1 restricted condition equally simultaneously.Referring to WOLFGANG KLIPPEL.Nonlinear modeling ofthe heat transfer in loudspeakers [J] .J.Audio Eng.Soc., 2004,52 (1/2): 3-21.
Situation to sum up, existing technical method is measured the calculating power compression based on thermal resistance, and is not only loaded down with trivial details but also can not get accurate result.The power compression computational methods that propose no thermal resistance are not also arranged so far.More and more universal at current high-power high-fidelity loudspeaker, moving-coil speaker power compression situation is more and more by the epoch that people paid close attention to, and a simple and easy quick power compression Forecasting Methodology all benefits for research and development, production.
Summary of the invention
The objective of the invention is: a kind of moving-coil speaker power compression without thermal resistance Forecasting Methodology is provided.For any moving-coil speaker, as long as change, utilize the criterion of no thermal resistance power compression directly to extrapolate the power compression value through impedance curve and the voice coil temperature of measuring moving-coil speaker, dope the power compression curve.
The objective of the invention is to realize like this: the impedance curve of at first measuring moving-coil speaker; The variation of respective frequencies point measurement voice coil temperature then; Utilize the power compression without thermal resistance computing formula; Obtain the power compression value of characteristic frequency point, revise through correction formula again, dope the power compression curve under the different capacity condition at last.
Moving-coil speaker power compression without thermal resistance Forecasting Methodology: for any moving-coil speaker; The at first impedance curve of testing of loudspeakers and temperature curve; Extract the modulus of impedance value of needed Frequency point then; Utilize the power compression without thermal resistance computing formula directly to calculate the power compression value again, frequency range is predicted, directly make the power compression curve with the frequency range algorithm like needs.The curve of making reduces curve for the relative amplitude of the frequency response before and after the rising of moving-coil speaker voice coil temperature.
The present invention chooses two characteristic frequency points and measures: valley point behind resonant frequency point and the impedance curve peak.Measuring system is simple, and the measuring equipment cost is low, can express delivery dope power compression, and the inventive method also can guarantee certain relative accuracy.The present invention also aims to carry out the single-point prediction for the characteristic frequency point, the power compression value of characteristic frequency point can characterize the quality of the heat dispersion of moving-coil speaker to a certain extent.Do not adopt the thermal resistance calculation power compression, need not linear or nonlinear hot equivalent model, do not receive the material and the effect on structure of moving-coil speaker basin frame, vibrating diaphragm etc.
The inventive method also is included on the basis of the variations in temperature of measuring impedance curve and voice coil loudspeaker voice coil, utilizes the power compression without thermal resistance computing formula to calculate, and this calculating formula is: PC = 10 Log 10 [ ( Z Ta + k R e Δ T ) ( 1 + KΔ T ) Z Ta ] .
Improvement of the present invention also comprises: moving-coil speaker power compression without thermal resistance Forecasting Methodology; Correction formula is
Figure G2008101964795D00022
correction factor a; B, c are all based on the measurement of impedance curve and variations in temperature.
Z wherein TaBe the loudspeaker impedance mould value before the temperature rising, Z TcBe the loudspeaker impedance mould value after the temperature rising, f SOBe valley point frequency behind the impedance curve of loudspeaker peak, f is the Frequency point of required correction, R eBe the resistance of voice coil loudspeaker voice coil direct current, k is a voice coil loudspeaker voice coil material temperature coefficient, and Δ T is the temperature rise of voice coil loudspeaker voice coil, and QT is the total quality factor of loud speaker, P rBe the rated power of loudspeaker, P is an input power.
Characteristics of the present invention are: the moving-coil speaker power compression computational methods that proposed not rely on thermal resistance calculation first.As long as the impedance curve and the variations in temperature of known loudspeaker just can calculate the power compression value and dope the compression curve under the different capacity initial conditions.
Description of drawings
Fig. 1 is the flow process of the inventive method,
Fig. 2 is the inventive method block diagram
Fig. 3 is corresponding for the impedance curve of loudspeaker of the present invention and compression curve
Fig. 4 is for measuring the circuit of voice coil temperature.
Fig. 5 has provided and has utilized correction back power compression not have theoretical sound pressure level curve and measured curve contrast that the thermal resistance calculation formula calculates.
Embodiment
Fig. 1,2 is a basic skills of the present invention.
Flow process of the present invention is as shown in Figure 1.The resistance of voice coil loudspeaker voice coil direct current increases along with the increase of voice coil temperature, and relational expression is:
R ETc=R ETa(1+k Δ T) R ETcBe the voice coil loudspeaker voice coil direct current resistance after heating up, R ETaBe the resistance of the voice coil loudspeaker voice coil direct current before heating up, Δ T is the front and back temperature difference, and k is a voice coil loudspeaker voice coil material temperature coefficient, is the moving-coil speaker of copper for the voice coil loudspeaker voice coil material, is 0.00393K -1
The power compression effect is meant that the voice coil temperature rising makes the resistance of voice coil loudspeaker voice coil direct current increase when speaker operation, under loud speaker constant voltage input condition, according to P=U 2/ R, the power of the loud speaker of feeding can descend, and so just can finally cause the decline of sensitivity and radiation efficiency.It is generally acknowledged that the power compression effect is in particular in the variation of direct current resistance, electrical quality factor, acoustic radiation efficiency and BL, direct result is exactly to cause the amplitude of the frequency response of loud speaker to descend relatively.Loud speaker electric power sound system (H 1And the hot (H of system (s)) 2(s)) all regard nonlinear and time-varying system continuous time (Fig. 2) as, after signal of telecommunication X (s) feeds loud speaker, the thermal power of system output Y 2(s) through the H of system 2(s), output temperature difference Y 3(s), negative feedback is to the H of system 1(s), cause last acoustics output Y 1(s) minimizing.
Before the calculating of power compression all had been based on the method for thermal resistance; Owing to count thermal resistance; Will certainly cause the last result who in different frequency range, calculates difference again; Because the air forced convection that moving-coil speaker produces when low frequency, it is non-linear to make that thermal resistance occurs, and the power compression result who obtains is inaccurate.
The present invention obtains the more convenient quick very computational methods of pinpoint accuracy that have again from impedance curve.
Under the constant voltage initial conditions, the output sound pressure level before voice coil temperature raises:
SPL 1=SPL 0+10log 10η 1P 1=SPL 0+10log 10η 1+10log 10P 1
Output sound pressure level after temperature raises:
SPL 2=SPL 0+10log 10η 2P 2=SPL 0+10log 10η 2+10log 10P 2
η 1, η 2Be divided into and be the loudspeaker enclosure of temperature after raising with temperature before raising.P 1And P 2Be respectively before the temperature rising and the loudspeaker enclosure after the temperature rising.SPL 0Be the output sound pressure level under the 1W input power condition.
At ka<<under 1 condition, &eta; = B 2 l 2 R e &rho; 0 2 &pi; c 0 S D 2 ( M MS + M MR ) 2 w 2 / w s 2 1 Q T 2 + ( w w s - w s w ) 2
After voice coil temperature raise, (the 1+k Δ T) before the direct current resistance becomes greatly times, then η became 1/ original (1+k Δ T), therefore
&eta; 2 = 1 1 + k&Delta;T &eta; 1
PC = SPL 1 - SPL 2 = 10 log 10 ( P 1 P 2 ) + 10 log 10 ( &eta; 1 &eta; 2 ) = 10 log 10 ( P 1 P 2 ) + 10 log 10 ( 1 + k&Delta;T )
= 10 log 10 ( 1 + k&Delta;T ) P 1 P 2
P 1 = V 2 Z Ta
P 2 = V 2 Z Tc
Z Tc=Z Ta+kR eΔT
PC = 10 log 10 [ ( Z Ta + k R e &Delta;T ) ( 1 + k&Delta;T ) Z Ta ]
Z wherein TaBe the loudspeaker impedance mould value before the temperature rising, Z TcBe the loudspeaker impedance mould value after the temperature rising, R eBe the resistance of voice coil loudspeaker voice coil direct current, k is a voice coil loudspeaker voice coil material temperature coefficient, and Δ T is the temperature rise of voice coil loudspeaker voice coil.
Can see that from following formula the voice coil loudspeaker voice coil temperature rise reaches under the situation of stable state, the PC value is at Z TaReach maximum when reaching minimum.So the selected characteristic Frequency point, promptly there is good representativeness the valley point behind the peak of impedance curve, simultaneously since loud speaker at the resonance frequency place frequency of occurrences peak value of response, choosing resonance frequency is another feature Frequency point.From Fig. 3, can obviously find out the corresponding relation of power compression peak value and impedance curve characteristic frequency point.
Fig. 4 is for measuring the circuit of voice coil temperature.The basic principle of circuit is to utilize precision resister and the moving-coil speaker dividing potential drop of 1K.Input is that signal is exported through power amplifier.C6, C7 are capacitance (2.2mF/35V), and the butt joint of bipolarity electric capacity can realize jumbo nonpolar electric capacity.R1 is the cement resistor of 0.22ohm/5W, is used to monitor the voice coil loudspeaker voice coil electric current, and R2 is that temperature coefficient is 25ppm/ ℃ a precision resistance, and C1, C2 are used to protect direct-flow voltage regulation source, measuring DC voltage behind the 3 rank LPFs.Since added cement resistor, so calculating formula is:
V measure 5 = R eTc + 0.22 1000 + R eTc + 0.22
R eTc = 1000.22 V measure 5 - V measure - 0.2
&Delta;T = R eTc - R eTa R eTa &delta;
According to the temperature-measuring results of voice coil loudspeaker voice coil temperature measurement circuit with from the modulus of impedance value that impedance curve obtains, directly calculate the power compression value with formula.
Because under the actual measurement condition, the measured deviation of steady temperature can cause the prediction margin error.Need make correction to the result; Correction term is
Figure G2008101964795D00044
correction factor a; B, c obtains through the power compression data computation at characteristic frequency point place.Q wherein TBe the total quality factor of loud speaker, f S0Be valley point frequency behind the impedance curve of loudspeaker peak, P rBe the rated power of loudspeaker, P is an input power.
The concrete calculated examples of measuring
Be example with 1 inch moving-coil speaker below, the process of utilizing non-thermal resistance prediction method to calculate the moving-coil speaker power compression is described.
According to measuring flow process; At first measure the impedance curve of moving-coil speaker; Characteristic of correspondence Frequency point place is with the variations in temperature of circuit measuring voice coil loudspeaker voice coil shown in Figure 3 again; The power compression formula of the no thermal resistance of utilization calculates the power compression value of single-point then, revises, and calculates the prediction curve that then can obtain power compression through frequency range again.The actual measurement parameter of this moving-coil speaker:
Re=4.31Ω,f 0=212Hz,f s0=530Hz,Q m=9.30,Q e=4.35,Qt=2.96,Z max=13.53Ω。Under 212HZ and the 530HZ initial conditions; When the different input power situation; Its temperature raises when being respectively 50 ℃ and 80 ℃; Not having the thermal resistance calculation formula through power compression, to calculate the power compression theoretical value be 1.03dB/1.65dB (70%Pr/100%Pr), and the compression measured value is 0.76dB/1.32dB (70%Pr/100%Pr), and the valley point is 530Hz behind the impedance peak; Its corresponding compression theory value is 1.08dB/2.26dB (70%Pr/100%Pr), and the compression measured value is 0.88dB/1.92dB (70%Pr/100%Pr).
Through the theoretical power (horse-power) compressed value that the characteristic frequency point is calculated and the comparison of actual measured value, the correction term parameter is carried out at 2 ask average computation, obtain a=316*10 -4, b=0.09, so the c=0.51 correction term is:
3.16 &times; 10 - 4 &times; P Pr + 0.09 &times; f - f s 0 f s 0 + 0.51 &times; 1 Q T .
Fig. 5 has provided and has utilized correction back power compression not have theoretical sound pressure level curve and measured curve contrast that the thermal resistance calculation formula calculates.Can see that promptly below the 750HZ, the power compression curve and the theoretical curve that utilize Forecasting Methodology provided by the invention to obtain have the good goodness of fit under low frequency condition (ka < < 1).

Claims (2)

1. moving-coil speaker power compression without thermal resistance Forecasting Methodology; It is characterized in that for any moving-coil speaker; Measure the variations in temperature of impedance curve and voice coil loudspeaker voice coil; Calculate the power compression value of characteristic frequency point through the power compression without thermal resistance computing formula, calculate the correction formula correction factor then, obtain the power compression curve under the different input power condition at last; On the basis of the variations in temperature of measuring impedance curve and voice coil loudspeaker voice coil; Utilize the power compression without thermal resistance computing formula to calculate, this calculating formula is:
Correction formula does
Figure FSB00000712808300012
Correction factor a, b, c are all based on the measurement of impedance curve and variations in temperature; Z TaBe the loudspeaker impedance mould value before the temperature rising, f S0Be valley point frequency behind the impedance curve of loudspeaker peak, f is the Frequency point of required correction, R eBe the resistance of voice coil loudspeaker voice coil direct current, k is a voice coil loudspeaker voice coil material temperature coefficient, and Δ T is the temperature rise of voice coil loudspeaker voice coil, Q TBe the total quality factor of loud speaker, P rBe the rated power of loudspeaker, P is an input power; Through the theoretical power (horse-power) compressed value that the characteristic frequency point is calculated and the comparison of actual measured value, the correction term parameter is carried out at 2 ask average computation, obtain correction factor a, b, c; Correction term is:
3.16 &times; 10 - 4 &times; P Pr + 0.09 &times; f - f s 0 f s 0 + 0.51 &times; 1 Q T .
2. moving-coil speaker power compression without thermal resistance Forecasting Methodology as claimed in claim 1, it is characterized in that getting earlier resonant frequency point and impedance curve peak after two Frequency points in valley point measure.
CN2008101964795A 2008-09-09 2008-09-09 Non-thermal resistance prediction method for power compression of coil driven loudspeaker Expired - Fee Related CN101394688B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101964795A CN101394688B (en) 2008-09-09 2008-09-09 Non-thermal resistance prediction method for power compression of coil driven loudspeaker

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101964795A CN101394688B (en) 2008-09-09 2008-09-09 Non-thermal resistance prediction method for power compression of coil driven loudspeaker

Publications (2)

Publication Number Publication Date
CN101394688A CN101394688A (en) 2009-03-25
CN101394688B true CN101394688B (en) 2012-05-09

Family

ID=40494702

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101964795A Expired - Fee Related CN101394688B (en) 2008-09-09 2008-09-09 Non-thermal resistance prediction method for power compression of coil driven loudspeaker

Country Status (1)

Country Link
CN (1) CN101394688B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101924054B (en) * 2010-05-12 2014-10-22 上海华虹宏力半导体制造有限公司 Method for measuring change of resistance of high-resistance semiconductor substrate with change of thermal budgets
CN107135458B (en) * 2017-06-13 2022-11-04 广州市迪士普音响科技有限公司 Method and device for measuring power borne by loudspeaker
KR102524292B1 (en) * 2018-08-06 2023-04-21 현대자동차주식회사 Speaker temperature controlling device for vehicle and method for controlling the same and vehicle including the same
US11333622B2 (en) * 2019-05-01 2022-05-17 Cirrus Logic, Inc. Thermal model of transducer for thermal protection and resistance estimation
CN111343560B (en) * 2020-03-17 2021-06-18 厦门傅里叶电子有限公司 Method for testing and tracking resonant frequency fo of mobile phone loudspeaker

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1589065A (en) * 2004-08-13 2005-03-02 黄新民 System capable of eliminating heat effect of speaker

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1589065A (en) * 2004-08-13 2005-03-02 黄新民 System capable of eliminating heat effect of speaker

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
李海锋.扬声器热仿真分析.《电声技术》.2008,第32卷(第5期),32-34. *
沈勇等.扬声器中涡流的热效应.《应用声学》.2006,第25卷(第2期),96-101.
沈勇等.扬声器的热效应与功率压缩.《电声技术》.2005,(第2005年第7期),18-21.
邬宁
邬宁;沈勇等.扬声器中涡流的热效应.《应用声学》.2006,第25卷(第2期),96-101. *
邬宁;沈勇等.扬声器的热效应与功率压缩.《电声技术》.2005,(第2005年第7期),18-21. *

Also Published As

Publication number Publication date
CN101394688A (en) 2009-03-25

Similar Documents

Publication Publication Date Title
CN101394688B (en) Non-thermal resistance prediction method for power compression of coil driven loudspeaker
US10080085B2 (en) Loudspeaker nonlinear compensation method and apparatus
CN106068007A (en) Nonlinear loudspeaker system identifying method
CN105120414B (en) Intelligent audio management system
CN103780997B (en) audio power management system and method thereof
CN101044793B (en) System and method for automatically adjusting hearing aid based on acoustic reflectance
CN102118678B (en) Method and system for measuring speaker parameters by using current sensor
US11057718B2 (en) Load change diagnostics for acoustic devices and methods
CN102158793A (en) Method utilizing laser sensor to measure speaker parameters and system
JP4917650B2 (en) Audio system with synthesized positive impedance
KR20130001162A (en) Active audio transducer protection
CN106797520A (en) The method and hearing aid device system of operating hearing aid system
CN102742300A (en) Control of a loudspeaker output
CN103500574B (en) A kind of intelligent electronic horn and its implementation
Klippel Adaptive stabilization of electro-dynamical transducers
US20190238994A1 (en) Armature-based acoustic receiver having improved output and method
CN104344886A (en) Method for evaluating high sound intensity characteristics of sound intensity measuring instrument
US20190260333A1 (en) Controlling parameters of an amplifier system based on a measured physical quantity
JPWO2015012401A1 (en) Speaker control device
CN101296533A (en) Loudspeaker diaphragm equivalent mass and equivalent force measuring method
US20160261954A1 (en) High-fidelity sound reproduction equipment
CN214310734U (en) Audio signal detection pen
Chen et al. Investigation of the Use of Loudspeaker as a Liner Alternator for Thermoacoustic Application
CN205883501U (en) Multi -mode microphone
CN217738342U (en) Acoustic control radar current surveying system

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120509

Termination date: 20150909

EXPY Termination of patent right or utility model