CN100558096C - A kind of quadrature amplitude modulation demodulation method and device that is applied to communication system - Google Patents
A kind of quadrature amplitude modulation demodulation method and device that is applied to communication system Download PDFInfo
- Publication number
- CN100558096C CN100558096C CNB2007101195276A CN200710119527A CN100558096C CN 100558096 C CN100558096 C CN 100558096C CN B2007101195276 A CNB2007101195276 A CN B2007101195276A CN 200710119527 A CN200710119527 A CN 200710119527A CN 100558096 C CN100558096 C CN 100558096C
- Authority
- CN
- China
- Prior art keywords
- symbol
- bound
- mrow
- real
- demodulation
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 238000004891 communication Methods 0.000 title claims abstract description 29
- 238000007781 pre-processing Methods 0.000 claims abstract description 13
- 238000010586 diagram Methods 0.000 claims description 35
- 238000012937 correction Methods 0.000 claims description 8
- 238000013507 mapping Methods 0.000 claims description 6
- 238000009738 saturating Methods 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 abstract description 14
- 230000002950 deficient Effects 0.000 abstract 1
- 238000012545 processing Methods 0.000 description 7
- 238000010606 normalization Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 229920006395 saturated elastomer Polymers 0.000 description 5
- 238000010295 mobile communication Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 229910003460 diamond Inorganic materials 0.000 description 2
- 239000010432 diamond Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
The invention provides a kind of quadrature amplitude modulation demodulation method that is applied to communication system, may further comprise the steps: step 1, receiving terminal adopt planisphere that the symbol that receives is shone upon, and obtain the final symbol of described symbol, and to this final symbol operated in saturation; Step 2 is determined the predetermined thresholds at different levels of the soft bit demodulation of M unit quadrature amplitude modulation, goes out each bit of this final symbol according to the final sign computation after this predetermined threshold and the described operated in saturation.Symbol preprocessing process among the present invention only comprises phase place adjustment and symbol adjustment, and calculate in the process of each bit of described final symbol and adopt plus and minus calculation mostly, therefore demodulating process is simple, demodulation performance is good, having solved the receiving terminal Block Error Rate sharply increases and causes the data transmission efficiency of network to reduce, and the defective that reduces in the data transmission efficiency that the complicated demodulation computing that receiving terminal carries out is caused.
Description
Technical Field
The present invention relates to demodulation technologies in the field of communications, and in particular, to a quadrature amplitude modulation and demodulation method and apparatus applied to a communication system.
Background
With the rapid development of communication technology, the demand of people for transmitting more information on the limited communication network bandwidth is stimulated. In terms of Modulation schemes for communication, Quadrature Amplitude Modulation (QAM) can approach the limit of channel capacity more than channel resources occupied by other Modulation schemes when the bandwidth is limited, that is, QAM has higher transmission efficiency under the same bandwidth and transmission power.
For the above reasons, in high-speed data transmission, for example, in Wireless network communication such as Worldwide Interoperability for Microwave Access (WiMAX), satellite communication, optical cable communication, or in communication technologies under different standard communication protocols in the field of Mobile communication, for example, Global System for Mobile communication (GSM Global System for Mobile communication), Time Division synchronous Code Division Multiple Access (TD-S CDMA Time Division Multiple Access), Wireless Code Division Multiple Access (WCDMA Wireless CDMA), Code Division Multiple Access 2000(CDMA2000Code Division Multiple Access 2000), and various other technologies applied to high-speed data transmission, QAM methods such as 16QAM, 64QAM, even 256QAM, etc. are often used, and thus QAM is important, but the existing QAM methods or corresponding devices are due to the change of signal amplitude at the receiving end, therefore, the demodulation performance is reduced, the block error rate is increased rapidly, and the data transmission efficiency of the network is reduced; even though some prior art proposals for improving the above-mentioned drawbacks still require complex demodulation operations at the receiving end, these operations consume time, affect the demodulation capability of the receiving end, and still result in a reduction in the data transmission efficiency of the network.
Disclosure of Invention
The invention aims to provide a quadrature amplitude modulation and demodulation method and a quadrature amplitude modulation and demodulation device applied to a communication system, which solve the defects that the data transmission efficiency of a network is reduced due to the fact that the block error rate of a receiving end is increased sharply, and the data transmission efficiency is reduced due to the fact that too complex demodulation operation is carried out on the receiving end.
A quadrature amplitude modulation and demodulation method applied to a communication system, comprising the steps of: step one, the receiving end adopts a constellation diagram to map the received symbol, and adjusts the amplitude of a rotating symbol corresponding to the symbol according to a preset average power to obtain a final symbol
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mo>,</mo>
</mrow>
</math>
Wherein N is the magnitude of some or all of the selected rotation symbols; k is the value of the predetermined average power, and K > 0; saturating the final symbol; step two, determining each level of preset threshold of soft bit demodulation of M-element orthogonal amplitude modulation to obtain a preset average power value, wherein the value M of M in the M-element orthogonal amplitude modulation is 22*PP in the index is a positive integer; and determining the shortest distance a from the position of the rotation symbol to the coordinate axis, wherein the calculation formula of each level of threshold bound _ p is as follows: wherein, P is 1, 2., P-1; the bits of the final symbol obtained by calculating the ith symbol according to the threshold are respectively as follows:
|
|
|
bi(2P-1)=bound_(P-1)p-|bi(2P-3)|
bi(2P)=bound_(P-1)-|bi(2P-2)|。
the method, wherein the mapping the received symbols in the first step further includes: step A, judging if the phase adjustment needs to be carried out on the received symbol, firstly carrying out the phase adjustment on the symbol to obtain a rotating symbol, and turning to the step B; otherwise, directly turning to the step B; and B, adjusting the amplitude of the rotating symbol or the received symbol according to the preset average power to obtain a final symbol.
The method, when the method is applied to 16-ary qam demodulation in TD-SCDMA, the phase adjusting the received symbols further comprises: clockwise rotating the constellation diagram by 45 degrees to form a rectangular diagram under a rectangular plane coordinate system; adjusting the phase of the rotation sign to obtain a rotation sign; and the rotation symbol is obtained by adopting the following method: rotateReal part of symbol: imaginary part of the rotation sign: wherein,is the symbol that was received and is,is a sign of the rotation of the rotary shaft,is the real part of the received symbol,is the real part of the rotated symbol,is the imaginary part of the received symbol and,is the rotated sign imaginary part.
The method, wherein, in the step B, the final symbol is obtained by adjusting the amplitude of the rotation symbol according to the predetermined average powerAccording to the following method:
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
</mrow>
</math>
wherein N is the magnitude of some or all of the selected rotation symbols; said K is the number of said predetermined average powerAnd K > 0.
The method, wherein the obtaining of the final symbolThen, the final symbol is further processedA saturation operation comprising: if the final symbolIs greater than real or imaginary partThen the saturation isIf the final symbolIs less than the real or imaginary partThen the saturation is
The method, wherein if the value K of the predetermined average power is 1, the final symbol is determined to be the last symbolIf the real part or the imaginary part of (b) is greater than 1, the saturation is 1; if the final symbolIs less than-1, the saturation is-1.
The above-mentioned process, wherein,in the second step, the predetermined threshold is obtained according to the following method: obtaining a predetermined average power value, wherein M in M-ary quadrature amplitude modulation is 22*PP in the index is a positive integer; determining the shortest distance a from the position of the rotating symbol on the constellation diagram to the coordinate axis, wherein the calculation formula of each level of threshold bound _ p is as follows: wherein, P is 1, 2., P-1; the bits obtained by calculating the ith symbol according to the threshold are respectively:
|
|
|
bi(2P-1)=bound_(P-1)p-|bi(2P-3)|
bi(2P)=bound_(P-1)-|bi(2P-2)|。
in the above method, after calculating each bit of the final symbol in the second step, the method further includes: adopting hard bit demodulation to directly judge the obtained value, and carrying out judgment on the bi(1) To bi(2P) if the value is not less than 0, the value is judged to be 0, and if the value is less than 0, the value is judged to be 1; or b is the same asi(1) To biThe value of (2P) is 1 when it is 0 or more and 0 when it is less than 0.
The above method, wherein when M is 16, P is 2, then P is 1; according to the calculation formula To obtain Soft bit demodulation for 16-ary quadrature amplitude modulation is as follows:
when M is 64, P is 3, then P is 1 or 2; according to the calculation formula To obtain For 64-ary quadrature amplitude modulationThe soft bit demodulation of (a) is as follows:
when M is 256, P is 4, then P is 1, 2 or 3; according to the calculation formula To obtain The soft bit demodulation for 256-ary quadrature amplitude modulation is as follows:
the above method, wherein the calculation formula When the value K of the predetermined average power is 1: for soft bit demodulation of 16-ary quadrature amplitude modulation, then:
for soft bit demodulation of 64-ary quadrature amplitude modulation, then:
for soft bit demodulation of 256-ary quadrature amplitude modulation, then:
a quadrature amplitude modulation and demodulation device applied to a communication system comprises a symbol preprocessing module and a quadrature amplitude modulation and demodulation module; the symbol preprocessing module comprises a phase adjusting unit, an amplitude adjusting unit and a saturation correcting unit, and the three units are electrically connected in sequence; the phase adjusting unit is used for adjusting the constellation diagram into a rectangular diagram under a plane rectangular coordinate system; adjusting the received symbols to the corresponding positions of rectangular diagrams of which the constellation diagrams are positioned under a plane rectangular coordinate system; the amplitude adjusting unit adjusts the amplitude of the rotation symbol corresponding to the symbol according to the preset average power to obtain the final symbol
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mo>,</mo>
</mrow>
</math>
Wherein N is the magnitude of some or all of the selected rotation symbols; k is the value of the predetermined average power, and K > 0; the saturation correction unit performs saturation operation on the final symbol; the quadrature amplitude modulation and demodulation module is used for determining each level of preset threshold of the soft bit demodulation method of M-ary quadrature amplitude modulation and obtaining a preset average power value, wherein the value M of M in the M-ary quadrature amplitude modulation is 22*PP in the index is a positive integer; and determining the shortest distance a from the position of the rotation symbol to the coordinate axis, wherein the calculation formula of each level of threshold bound _ p is as follows: wherein, P is 1, 2., P-1; the bits of the final symbol obtained by calculating the ith symbol according to the threshold are respectively as follows:
|
|
|
bi(2P-1)=bound_(P-1)p-|bi(2P-3)|
bi(2P)=bound_(P-1)-|bi(2P-2)|。
the above device, wherein when the device is applied to 16-ary quadrature amplitude modulation and demodulation in TD-SCDMA, the phase adjusting unit rotates the constellation diagram clockwise by 45 degrees into a rectangular diagram under a rectangular plane coordinate system; adjusting the phase of the received symbol to a corresponding position in a rectangular graph under the rectangular coordinate system after the symbol rotates clockwise by 45 degrees; the phase adjustment unit of the symbol preprocessing module calculates the rotation symbol according to the following formulaReal part ofImaginary part
And rotating the symbolAnd real partImaginary partAnd sending the data to an amplitude adjusting unit.
The apparatus above, wherein the qam/qam module obtains a predetermined average power value K, and the value M of M in M-ary qam is 22*PCalculating a positive integer P; determining the shortest distance a from the position of the rotating symbol on the constellation diagram to the coordinate axis; according to the formula The values of the fixed thresholds bound _1, bound _2.. round _ (P-1) of each stage are determined, where P1, 2.., P-1.
The above apparatus, wherein the obtaining of each bit from the final symbol output by the saturation correction unit further includes: adopting a hard bit demodulation method to directly judge the acquired value, and carrying out the judgment on the bi(1) To bi(2P) if the value is not less than 0, the value is judged to be 0, and if the value is less than 0, the value is judged to be 1; or b is the same asi(1) To biThe value of (2P) is 1 when it is 0 or more and 0 when it is less than 0.
The apparatus above, wherein the qam/qam module demodulates qam with gray code characteristics according to a preset setting, and supports at least 16-ary qam/qam, 64-ary qam/qam, or 256-ary qam/qam.
The quadrature amplitude modulation and demodulation method and device applied to the communication system provided by the invention have the following steps that firstly, a receiving end adopts a constellation diagram to map a received symbol to obtain a final symbol of the symbol, and the final symbol is subjected to saturation operation; and step two, determining each level of preset threshold of soft bit demodulation of M-element quadrature amplitude modulation, and calculating each bit of the final symbol according to the preset threshold and the final symbol after saturation operation. The block error rate of the receiving end is greatly reduced, the data transmission efficiency of the network is improved, and meanwhile, the data transmission efficiency is improved because the demodulation operation is simpler.
Drawings
Fig. 1 is a 16-ary qam constellation according to an embodiment of the present invention;
fig. 2 is a process diagram of a 16-ary qam demodulation method according to an embodiment of the present invention;
fig. 3 is a schematic diagram of 16-ary qam-demod constellation conversion according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a 16-ary QAM/QAM apparatus according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of mapping 64-ary QAM constellation in the present invention;
fig. 6 is a schematic diagram of the qam demodulation method applied to a communication system according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and implementation effects of the present invention clearer, the technical solutions of the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
The invention provides a method and a device for quadrature amplitude modulation and demodulation in a communication system, which can be widely applied to different communication technologies in the communication field, and the specific technical scheme of the invention is described below by taking 16QAM adopted in TD-SCDMA standard in the mobile communication field as an embodiment. Table 1 is a 16QAM constellation mapping relationship defined in the 3GPP R525.223 protocol, where each 4 bits corresponds to a symbol, where a symbol is a complex symbol, and the real part of the symbol is also referred to as I-path and the imaginary part is also referred to as Q-path; by representing the bit-symbol mapping relationship reflected in table 1 in a rectangular plane coordinate system, as shown in fig. 1, the constellation diagram of 16QAM is in the shape of a diamond, and the internal angle of the diamond is 90 °. The value ranges of the I path and the Q path of the 16QAM symbol are all
Table 116 QAM constellation mapping table
Fig. 2 shows a 16QAM process in TD-SCDMA, which comprises the following specific steps:
And 203, forming data of the physical code channel at the transmitting end by the symbol, and transmitting the data to an antenna end through a radio frequency circuit to form a transmitting signal.
As can be seen from fig. 6, after the rf front end of the receiving end receives the signal of the transmitting end, the signal is processed by rf down-conversion and baseband, and then is demodulated by data to obtain a digital signal, the digital signal is demodulated by the QAM demodulator, all demodulated bits are connected in sequence to form a bit stream, and the bit stream is sent to the bit decoder by the quantizer, so that the required user information bit can be obtained finally. The demodulator demodulation process of QAM of the embodiment of the present invention comprises the following two main steps, wherein, the symbolThe symbol before the clockwise rotation is indicated,the symbol rotation after clockwise rotation is shown, and without loss of generality, 16QAM is taken as an example for description.
Step one, preprocessing data;
as shown in FIG. 3, (1) first, phase adjustment is performed to obtain a symbolIn which the constellation diagram is rotated by 45 degrees clockwise, i.e. for symbolsRotated 45 degrees clockwise, i.e. to the symbolMultiplication by e-jπ/4Due to said
<math>
<mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>jπ</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
</msup>
<mo>=</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>π</mi>
<mo>/</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>j</mi>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>π</mi>
<mo>/</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<mo>,</mo>
</mrow>
</math>
So symbolAnd e-jπ/4Product of (2)Is still a symbol of complex formTheHereinafter referred to as a rotation symbol, ofReal part ofImaginary partThe method can be obtained by the addition and subtraction operation in the following way:
(2) secondly, the rotating symbolCarrying out amplitude adjustment, namely carrying out power normalization on the symbol after phase adjustment to obtain a final symbolIf normalized to a certain real value K, the real value K describes a specific value of the power.
First, all or part of the rotation symbol is calculatedAverage power ofWherein i is the subscript or serial number corresponding to the ith symbol, and then the reciprocal is obtainedFinally multiplied by the rotation signTo obtainWhere N represents the number of said full or partial rotation symbols. In the embodiment of the invention, because the amplitude adjustment processing is carried out on the symbol after the phase rotation, namely the symbol is normalized to a real power adjustment value K, the obtained final symbolThe method comprises the following steps:
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mo>,</mo>
</mrow>
</math>
wherein,representing the final sign after the phase rotation and amplitude adjustment processing, wherein N used for calculating the average power of the rotation sign satisfies,
<math>
<mrow>
<mn>1</mn>
<mo>≤</mo>
<mi>N</mi>
<mo>≤</mo>
<mi>length</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</math>
wherein,representing all symbols to be processedAnd N adopts an integer power of 2 as much as possible, so that the receiving end reduces the complexity of calculation.
After the power normalization process, some final symbols may appearReal part ofOr imaginary partIs greater thanOr less thanAnd performing saturation processing on the real parts or the imaginary parts. The saturation process will be whatever the real part of the final symbolOr imaginary partIs greater thanIs uniformly saturated byAll is the real partOr imaginary partIs less thanIs uniformly saturated byFor example, when K is 1, and is The 1.23 saturation treatment is 1 and the-1.05 saturation treatment is-1.
Step two, for the final symbolPerforming 16QAM, or 64QAM, or 256 QAM; the methods of demodulation generally include two types: soft Bit (Bit) demodulation, hard Bit (Bit) demodulation.
Soft bit demodulationThe method is as follows, wherein bi(1) To bi(2P) then respectively represent the first bit to the 2P bit of the ith final symbol:
firstly, determining values of fixed thresholds bound _1 and bound _2.. round.. bound _ (P-1) of each level, wherein the values of the fixed thresholds of each level are only related to a specific QAM constellation diagram and a preset amplitude adjustment value K, and on the constellation diagram with the amplitude adjustment value K being 1, the shortest distance from a standard position of a symbol to a coordinate axis is a, then the fixed thresholds bound _ P of each level are calculated as follows: wherein P is 1, 2, 1, and P is a lower case letter and P is an upper case letter.
For M-16 quadrature amplitude modulation, according to the formula M-22*PTo obtain P ═ 2, then In the 16QAM method in TD-SCDMA according to the embodiment of the present invention, since the amplitude adjustment value K is 1, it is obvious that
When the embodiment of the invention is applied to a non-TD-SCDMA communication system, ifIf 64-ary quadrature amplitude modulation is required, then for M-64-ary quadrature amplitude modulation, the value is taken According to the formula M-22*PGet P to 3, as shown in fig. 5, there are two fixed thresholds, then The fixed threshold identified in the figure is the value when K is 1, and the value of b1b2 is consistent in each quadrant, so only the value of b3b4 b5b6 is identified in fig. 5.
When the embodiment of the invention is applied to a non-TD-SCDMA communication system, if 64-element quadrature amplitude modulation is required, for M-256-element quadrature amplitude modulation, the M-256-element quadrature amplitude modulation is selected According to the formula M-22*PTo get P to 4, there are three fixed thresholds:
then
After determining the respective fixed threshold, for soft bit demodulation using gray coded M-ary quadrature amplitude modulation, the respective bit of each final symbol is calculated using the following method:
the value of M in the M-element quadrature amplitude modulation is as follows: m22*PWherein P in the index is a positive integer, and in the formula, P1, 2
|
|
|
bi(2P-1)=bound_(P-1)p-|bi(2P-3)|
bi(2P)=bound_(P-1)-|bi(2P-2)|
The value obtained by the above formula is not usually an integer of 0 or 1, but is a decimal value, wherein | b in the formulai(2P-2)|、|bi(2P-3) | denotes the operation of taking the absolute value of a real number.
When M is 16, P is 2, then P is 1; soft bit demodulation for 16-ary quadrature amplitude modulation is as follows:
when M is 64, P is 3, then P is 1, 2; the soft bit demodulation for 64-ary quadrature amplitude modulation is as follows:
when M is 256, P is 4, then P is 1, 2, 3; the soft bit demodulation for 256-ary quadrature amplitude modulation is as follows:
the above describes the method using soft bit demodulation, and if the receiving end also uses the hard bit demodulation method, the value obtained by the soft bit demodulation method can be directly determined, taking fig. 5 as an example, for the bi(1) To bi(2P) taking an integer whenever the value falls within a predetermined certain range; for example, if the numerical value is 0 or more, it is determined as 0, and if it is less than 0, it is determined as 1; or if the numerical value is more than or equal to 0, the result is judged to be 1, and if the numerical value is less than 0, the result is judged to be 0; specifically, for example, the soft bit is 0.32, which can be judged as +1, and-0.04 can be judged as 0.
The hard bit demodulation method may be employed after the soft bit demodulation method, or only the hard bit demodulation method or only the soft bit demodulation method may be employed.
The invention also provides a device corresponding to the method, as shown in fig. 4, the device comprises: a symbol preprocessing module 401, a QAM demodulation module 405;
the preprocessing module is used for processing symbolsCarrying out pretreatment; in the TD-SCDMA communication system, the module rotates the constellation diagram of the received symbol clockwise by 45 degrees into a rectangular diagram under a rectangular plane coordinate system to obtain a rotating symbol; and adjusting the amplitude of the rotation symbol to a predetermined value according to the average power of the constellation to obtain a final symbolAnd then performing saturation treatment. The preprocessing module comprises a phase adjusting unit 402, an amplitude adjusting unit 403 and a saturation correcting unit 404, which are electrically connected in sequence. Wherein the phase adjustment unit 404 aligns the symbolsRotated 45 degrees clockwise, i.e. to the symbolMultiplication by e-jπ/4Due to said
<math>
<mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>jπ</mi>
<mo>/</mo>
<mn>4</mn>
</mrow>
</msup>
<mo>=</mo>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>π</mi>
<mo>/</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>j</mi>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>π</mi>
<mo>/</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<mo>,</mo>
</mrow>
</math>
So symbolAnd e-jπ/4Product of (2)Is still a symbol of complex formThe symbol may be referred to as a rotated symbol, the real part of whichImaginary partThe phase adjustment unit 404 obtains the following addition and subtraction operations:
the phase adjustment unit 404 obtains the rotation symbolAnd real partImaginary partSent to an amplitude adjustment unit 403, the amplitude adjustment unit 403 is used to adjust the symbolAmplitude adjustment, i.e. calculating the average power of the rotated sequenceThen, the reciprocal of the value is obtainedFinally multiplying by the rotated symbol sequence to obtainIn the invention, the data after phase rotation is normalized to a real value K representing power by carrying out power normalization processing, so that the obtained final symbolThe method comprises the following steps:
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mo>,</mo>
</mrow>
</math>
wherein,which represents the sign after phase rotation and normalization of the power, where N is satisfied for calculating the average power of the sequence,
<math>
<mrow>
<mn>1</mn>
<mo>≤</mo>
<mi>N</mi>
<mo>≤</mo>
<mi>length</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
</math>
here, N is an integer power of 2 as much as possible, so as to facilitate calculation and reduce implementation complexity. Symbol after power normalization processingThe phase adjustment unit 404 sends the phase adjustment unit 404 to the saturation correction unit 404.
The saturation correction unit 404 after normalization for saturation correctionFor those real partsOr imaginary partIs greater thanOr less thanIs saturated, i.e. the real partOr imaginary partIs greater thanIs uniformly saturated and corrected toWill partOr imaginary partIs less thanIs uniformly saturated and corrected to
The saturation correction unit 404 is connected to the QAM demodulation module 405.
The QAM demodulation module 405 is configured to demodulate and support at least 16QAM, 64QAM, 256 QAM; the methods of demodulation generally include two types: soft bit demodulation method, hard bit demodulation method.
In the soft bit demodulation method, the value of each fixed threshold bound _ p is determined first, the value of the threshold is only related to a specific QAM constellation diagram and a predetermined amplitude adjustment value K, as shown in fig. 5, the shortest distance from the constellation diagram with the amplitude adjustment value K of 1 to the coordinate axis is a, and then each fixed threshold is determined Wherein P is 1, 2. In the embodiment of the present invention, for the 16QAM method in TD-SCDMA, if K is 1, it is obvious that The specific implementation of M-16 quadrature amplitude modulation, M-64 quadrature amplitude modulation, and M-256 quadrature amplitude modulation is then the same as described in the method, and will not be described here in a repeated manner.
The QAM demodulation module 405 may select a module or a unit that supports M-16 quadrature amplitude modulation, M-64 quadrature amplitude modulation, and M-256 quadrature amplitude modulation. Meanwhile, the QAM demodulation module 405 provides a function of hard bit demodulation, and may directly determine the value obtained by the soft bit demodulation method, where b isi(1) To bi(2P) when the value falls within a certain range, e.g. for values greater than or equal toJudging the sample to be 0 if the sample is 0, and judging the sample to be 1 if the sample is less than 0; or if the numerical value is more than or equal to 0, the result is judged to be 1, and if the numerical value is less than 0, the result is judged to be 0; specifically, for example, the soft bit is 0.32, which can be judged as +1, and-0.04 can be judged as 0.
The hard bit demodulation method may be employed after the soft bit demodulation method, or only the hard bit demodulation method or the soft bit demodulation method may be employed.
The technical scheme described in the embodiment of the invention has simple demodulation process, wherein the symbol preprocessing process only comprises phase adjustment and amplitude adjustment, most of the symbol preprocessing process is addition operation, and the demodulation performance is good, thereby overcoming the defects of reduced network data transmission efficiency caused by the rapid increase of the block error rate of the receiving end and reduced data transmission efficiency caused by complex demodulation operation at the receiving end. And the technical problems solved by the technical scheme of the invention generally exist in the communication field, so the technical scheme of the invention can be applied in various communication fields.
Finally, it should be noted that the above embodiments are only used for illustrating the technical solutions of the present invention, and not for limiting, and all the values of the parameters can be adjusted according to the actual situation and are within the scope of the claims. It will be understood by those skilled in the art that various modifications and equivalent arrangements may be made without departing from the spirit and scope of the present invention and shall be covered by the appended claims.
Claims (12)
1. A qam demodulation method applied to a communication system, comprising the steps of:
step one, the receiving end adopts the constellation diagram to map the received symbols,
adjusting the amplitude of the rotation symbol corresponding to the symbol according to the preset average power to obtain the final symbol
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mo>,</mo>
</mrow>
</math>
Wherein N is the magnitude of some or all of the selected rotation symbols; k is the value of the predetermined average power, and K > 0; saturating the final symbol;
step two, determining each level of preset threshold of soft bit demodulation of M-element orthogonal amplitude modulation to obtain a preset average power value, wherein the value M of M in the M-element orthogonal amplitude modulation is 22*PP in the index is a positive integer; and determining the shortest distance a from the position of the rotation symbol to the coordinate axis, wherein the calculation formula of each level of threshold bound _ p is as follows: wherein, P is 1, 2., P-1;
the bits of the final symbol obtained by calculating the ith symbol according to the threshold are respectively as follows:
|
|
|
bi(2P-1)=bound_(P-1)p-|bi(2P-3)|
bi(2P)=bound_(P-1)-|bi(2P-2)|。
2. the method of claim 1, wherein said mapping the received symbols in step one further comprises:
if it is determined that phase adjustment is required for the received symbol, the symbol is first phase adjusted to obtain a rotated symbol.
3. The method of claim 2, wherein when the method is applied to 16-ary quadrature amplitude modulation and demodulation in TD-SCDMA, the phase adjusting the received symbols further comprises:
clockwise rotating the constellation diagram by 45 degrees to form a rectangular diagram under a rectangular plane coordinate system; adjusting the phase of the rotation sign to obtain a rotation sign; and the rotation symbol is obtained by adopting the following method:
real part of the rotation sign:
4. The method of claim 1, wherein the obtaining the final symbolThen, the final symbol is further processedA saturation operation comprising:
6. The method of claim 1, wherein the step two, after calculating the bits of the final symbol, further comprises:
adopting hard bit demodulation to directly judge the obtained value, and carrying out judgment on the bi(1) To bi(2P) if the value is not less than 0, the value is judged to be 0, and if the value is less than 0, the value is judged to be 1; or b is the same asi(1) To biThe value of (2P) is 1 when it is 0 or more and 0 when it is less than 0.
7. The method of claim 1, wherein when M is 16, P is 2, then P is 1; according to the calculation formula To obtain Soft bit demodulation for 16-ary quadrature amplitude modulation is as follows:
when M is 64, P is 3, then P is 1 or 2; according to the calculation formula To obtain The soft bit demodulation for 64-ary quadrature amplitude modulation is as follows:
when M is 256, P is 4, then P is 1, 2 or 3; according to the calculation formula To obtain The soft bit demodulation for 256-ary quadrature amplitude modulation is as follows:
8. the method of claim 7, wherein the computational formula When the value K of the predetermined average power is 1:
for soft bit demodulation of 16-ary quadrature amplitude modulation, then:
for soft bit demodulation of 64-ary quadrature amplitude modulation, then:
for soft bit demodulation of 256-ary quadrature amplitude modulation, then:
9. a quadrature amplitude modulation and demodulation device applied to a communication system is characterized by comprising a symbol preprocessing module and a quadrature amplitude modulation and demodulation module;
the symbol preprocessing module comprises a phase adjusting unit, an amplitude adjusting unit and a saturation correcting unit, and the three units are electrically connected in sequence;
the phase adjusting unit is used for adjusting the constellation diagram into a rectangular diagram under a plane rectangular coordinate system; adjusting the received symbols to the corresponding positions of rectangular diagrams of which the constellation diagrams are positioned under a plane rectangular coordinate system;
the amplitude adjusting unit adjusts the amplitude of the rotation symbol corresponding to the symbol according to the preset average power to obtain the final symbol
<math>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>norm</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mi>K</mi>
</msqrt>
<mo>/</mo>
<msqrt>
<mrow>
<mo>(</mo>
<munderover>
<mi>Σ</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>N</mi>
</msqrt>
<mo>)</mo>
</mrow>
<mo>×</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mi>rot</mi>
</msub>
<mo>,</mo>
</mrow>
</math>
Wherein N is the magnitude of some or all of the selected rotation symbols; k is the value of the predetermined average power, and K > 0;
the saturation correction unit performs saturation operation on the final symbol;
the quadrature amplitude modulation and demodulation module is used for determining each level of preset threshold of the soft bit demodulation method of M-ary quadrature amplitude modulation and obtaining a preset average power value, wherein the value M of M in the M-ary quadrature amplitude modulation is 22*PP in the index is a positive integer; and determining the shortest distance a from the position of the rotation symbol to the coordinate axis, wherein the calculation formula of each level of threshold bound _ p is as follows: wherein, P is 1, 2., P-1;
the bits of the final symbol obtained by calculating the ith symbol according to the threshold are respectively as follows:
|
|
|
bi(2P-1)=bound_(P-1)p-|bi(2P-3)|
bi(2P)=bound_(P-1)-|bi(2P-2)|。
10. the apparatus of claim 9, wherein when the apparatus is applied to 16-ary quadrature amplitude modulation and demodulation in TD-SCDMA, the phase adjustment unit rotates the constellation diagram clockwise by 45 degrees into a rectangular diagram under a rectangular plane coordinate system; adjusting the phase of the received symbol to a corresponding position in a rectangular graph under the rectangular coordinate system after the symbol rotates clockwise by 45 degrees;
the phase adjustment unit of the symbol preprocessing module calculates the rotation symbol according to the following formulaReal part ofImaginary part
11. The apparatus as claimed in claim 9, wherein the final symbol output by the saturation correction unit further comprises after obtaining each bit: adopting a hard bit demodulation method to directly judge the acquired value, and carrying out the judgment on the bi(1) To bi(2P) if the value is not less than 0, the value is judged to be 0, and if the value is less than 0, the value is judged to be 1; or b is the same asi(1) To biThe value of (2P) is 1 when it is 0 or more and 0 when it is less than 0.
12. The apparatus of claim 9, wherein the qam demodulation block demodulates qam having gray code characteristics according to a preset setting, and supports at least 16-ary qam demodulation, 64-ary qam demodulation, or 256-ary qam demodulation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101195276A CN100558096C (en) | 2007-07-25 | 2007-07-25 | A kind of quadrature amplitude modulation demodulation method and device that is applied to communication system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101195276A CN100558096C (en) | 2007-07-25 | 2007-07-25 | A kind of quadrature amplitude modulation demodulation method and device that is applied to communication system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101102301A CN101102301A (en) | 2008-01-09 |
CN100558096C true CN100558096C (en) | 2009-11-04 |
Family
ID=39036392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007101195276A Expired - Fee Related CN100558096C (en) | 2007-07-25 | 2007-07-25 | A kind of quadrature amplitude modulation demodulation method and device that is applied to communication system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100558096C (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102546518B (en) * | 2010-12-21 | 2014-12-17 | 中兴通讯股份有限公司 | Method and device used for demodulating 64 quadrature amplitude modulation (QAM) and used in high speed downlink packet access (HSDPA) business |
CN106713213B (en) * | 2015-07-25 | 2018-09-28 | 华为技术有限公司 | A kind of Ditital modulation method, demodulation method, relevant apparatus and system |
EP3490207B1 (en) | 2016-08-15 | 2021-03-31 | Huawei Technologies Co., Ltd. | Target constellation diagram determination method and apparatus |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5838728A (en) * | 1994-11-23 | 1998-11-17 | At&T Wireless Services, Inc. | High rate reed-solomon concatenated trellis coded 16 star QAM system for transmission of data over cellular mobile radio |
CN1443416A (en) * | 2001-05-17 | 2003-09-17 | 连宇通信有限公司 | A Method of Applying Quadrature Amplitude Modulation to Wireless Mobile Communication System |
CN1481634A (en) * | 2001-09-18 | 2004-03-10 | 三星电子株式会社 | Appts and method for calculating channel decoder input soft decision value in data communications system |
CN1697442A (en) * | 2005-07-08 | 2005-11-16 | 中兴通讯股份有限公司 | Demodulation method of 16QAM in time division synchronized CDMA system |
-
2007
- 2007-07-25 CN CNB2007101195276A patent/CN100558096C/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5838728A (en) * | 1994-11-23 | 1998-11-17 | At&T Wireless Services, Inc. | High rate reed-solomon concatenated trellis coded 16 star QAM system for transmission of data over cellular mobile radio |
CN1443416A (en) * | 2001-05-17 | 2003-09-17 | 连宇通信有限公司 | A Method of Applying Quadrature Amplitude Modulation to Wireless Mobile Communication System |
CN1481634A (en) * | 2001-09-18 | 2004-03-10 | 三星电子株式会社 | Appts and method for calculating channel decoder input soft decision value in data communications system |
CN1697442A (en) * | 2005-07-08 | 2005-11-16 | 中兴通讯股份有限公司 | Demodulation method of 16QAM in time division synchronized CDMA system |
Also Published As
Publication number | Publication date |
---|---|
CN101102301A (en) | 2008-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6434200B1 (en) | TCM revisiting system and method | |
US8340202B2 (en) | Method and apparatus for efficient soft modulation for gray-mapped QAM symbols | |
JP4470377B2 (en) | Propagation path estimation method in mobile communication system | |
AU2002321932A1 (en) | Demodulation apparatus and method in a communication system employing 16-ary QAM | |
CN112019298B (en) | Coding modulation method, demodulation decoding method, device and equipment | |
JP2008512036A (en) | Method and apparatus for calculating log approximation rate for decoding in receiver of mobile communication system | |
WO2003017610A1 (en) | Demodulation apparatus and method in a communication system employing 16-ary qam | |
WO1993007701A1 (en) | System for determining absolute phase of differentially-encoded phase-modulated signal | |
WO2006076599A1 (en) | Multiresolution modulation using variable euclidean distance ratio and blind receiver | |
AU2002329092A1 (en) | Apparatus and method for calculating soft decision value input to channel decoder in a data communication system | |
WO2003026243A1 (en) | Apparatus and method for calculating soft decision value input to channel decoder in a data communication system | |
CN104283835A (en) | Quadrature amplitude modulation qubit demodulation method and device | |
CN112350970A (en) | Multiphase frequency shift keying modulation and demodulation method and equipment | |
CN104956636B (en) | The method and apparatus that frequency orthogonal amplitude modulation is supported in wireless communication system | |
CN100558096C (en) | A kind of quadrature amplitude modulation demodulation method and device that is applied to communication system | |
EP3565208B1 (en) | Method and apparatus for demodulating signal | |
CN113364717B (en) | 32APSK soft information calculation method suitable for DVBS2 modulation | |
EP2073477B1 (en) | Combined 16QAM De-Mapping, 16QAM Constellation Re-Arrangement and TD-SCDMA De-Scrambling | |
EP0987863B1 (en) | Soft decision method and apparatus for 8PSK demodulation | |
CN103905370B (en) | Soft demodulation method and device for quadrature amplitude modulation (QAM) signal and mobile terminal | |
JP3851143B2 (en) | MODULATION SYSTEM IDENTIFICATION CIRCUIT, RECEPTION DEVICE EQUIPPED WITH SAME, WIRELESS STATION, AND MODULATION SYSTEM IDENTIFICATION METHOD | |
CN105530215B (en) | 8PSK soft differential decoding method based on CORDIC | |
Mueller et al. | A DSP implemented dual 9600/7200 bps TCM modem for mobile communications over FM voice radios | |
JP4139706B2 (en) | Receiver | |
CN113411279A (en) | Reusable Q power demapping method and system based on DVB-S2 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091104 |