CN103124190B - A kind of relaying mimo system method for precoding based on nonlinear organization - Google Patents

A kind of relaying mimo system method for precoding based on nonlinear organization Download PDF

Info

Publication number
CN103124190B
CN103124190B CN201310090457.1A CN201310090457A CN103124190B CN 103124190 B CN103124190 B CN 103124190B CN 201310090457 A CN201310090457 A CN 201310090457A CN 103124190 B CN103124190 B CN 103124190B
Authority
CN
China
Prior art keywords
relaying
matrix
signal
base station
user
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
CN201310090457.1A
Other languages
Chinese (zh)
Other versions
CN103124190A (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.)
Shandong University
Original Assignee
Shandong 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 Shandong University filed Critical Shandong University
Priority to CN201310090457.1A priority Critical patent/CN103124190B/en
Publication of CN103124190A publication Critical patent/CN103124190A/en
Application granted granted Critical
Publication of CN103124190B publication Critical patent/CN103124190B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)

Abstract

The invention discloses a kind of relaying mimo system method for precoding based on nonlinear organization, comprise the steps: the first step: signal a, after precoding processing, sends to relaying by base station; The signal a sent from base station adds noise n 1after, form signal x; Second step: relaying obtains channel H between base station and relaying according to the signal x received 1channel condition information, simultaneously relaying asks for the optimum pre-coding matrix W in base station based on Zero Forcing; 3rd step: relaying by the signal x received through delivery and feedback, carry out non-linear preliminary treatment after send to user, signal x adds noise n 2after, form signal y; 4th step: user obtains the channel condition information of the equivalent channel H between relaying and user according to the signal y received, calculates the optimum pre-coding matrix F of relaying, relaying feedback matrix B, user side weighting diagonal matrix G; The advantage of nonlinear precoding can be given full play to, the bit error rate performance of further improvement system in relay.

Description

A kind of relaying mimo system method for precoding based on nonlinear organization
Technical field
The present invention relates to a kind of relaying MIMO based on nonlinear organization (Multiple-Input Multiple-Output, multiple-input, multiple-output) system precoding method, belong to digital mobile communication system and technical field.
Background technology
The advantages such as it is large that MIMO communication system has capacity, and antijamming capability is strong, high speed propagation, all have a wide range of applications in future wireless system.Introduce relaying technique in mimo systems, effectively can resist wireless channel decline, improve power system capacity, expand system and cover, become the important candidate technologies in future broadband wireless communication systems.
In practical communication situation, transmitting terminal (base station) is fixed to relay relative position, and can be placed in higher position, therefore can build the line-of-sight propagation channel between base station and relay station, channel condition information changes less relatively, even can be considered as constant-parameter channel.And the communication channel between relay and user side, due to the mobility of user side and the complexity of reception environment, channel condition information change more complicated.Now there are some researches show, nonlinear precoding can realize the efficient channel power larger than linear predictive coding, so adopt nonlinear organization to have very great help to the communication performance improved between lifting relay to user side at relaying.The current research to relaying mimo system precoding technique is mostly based on linear link structure, the i.e. relay structure pattern of base station (linear/non-linear precoding) → relay (linear predictive coding) → user side, and adopt the research of nonlinear precoding to be also one piece of blank substantially to relay.
Summary of the invention
Object of the present invention is exactly the problems referred to above existed to solve prior art, a kind of relaying mimo system method for precoding based on nonlinear organization is provided, namely the precoding algorithms of THP (Tomlinson-Harashima precoding) is adopted based on relay, contrast with traditional relay structure, method of the present invention can give full play to the advantage of nonlinear precoding in relay, the bit error rate performance of further improvement system.
To achieve these goals, the present invention adopts following technical scheme:
Based on a relaying mimo system method for precoding for nonlinear organization, concrete steps are:
The first step: signal a, after precoding processing, sends to relaying by base station.What send from base station is the signal of a after the precoding processing of base station, adds noise n 1after, form signal x;
Second step: relaying obtains the channel condition information of channel between base station and relaying according to the signal x received, relaying is based on ZF (Zero Forcing simultaneously, ZF) criterion asks for the optimum pre-coding matrix W in base station, wherein, the channel condition information matrix H of channel between base station and relaying 1represent; ;
3rd step: relaying sends to user after the signal x received is carried out non-linear preliminary treatment (delivery and feedback).What therefrom secondary went out is x, through the non-linear pretreated signal of relaying, adds noise n 2after, form signal y;
4th step: user obtains the channel condition information of the equivalent channel between relaying and user according to the signal y received, calculate the optimum pre-coding matrix F of relaying, relaying feedback matrix B, user side weighting diagonal matrix G, wherein, the channel condition information of the equivalent channel between relaying and user represents by matrix H.
Wherein, the signal x in the described first step adopts following formula to obtain:
x=H 1Wa+n 1(1)
Wherein, a is initialize signal data flow, and W is the optimum pre-coding matrix in base station, H 1, n 1base station and relay well channel condition information and base station and relay well interchannel noise respectively.A ∈ C m × 1, W ∈ C m × M, H 1∈ C l × M, n 1∈ C l × 1, x ∈ C l × 1, M, L are the antenna number of base station and relaying respectively.
The processing method asking for the optimum pre-coding matrix W in base station in described second step is: relaying, according to ZF (ZF) criterion, is tried to achieve:
W = K T r ( H 1 - 1 H 1 - T ) H 1 - 1
Wherein, K=min (M, L, N), N are user side antenna number, H 1 -1to matrix H 1get inverse operation, H 1 -Tto matrix H 1get inverse after transpose operation again, Tr (H 1 -1h 1 -T) be to matrix (H 1 -1h 1 -T) track taking computing.
Wherein, the non-linear preliminary treatment of relaying to signal in the 3rd step obtains according to following formula:
x-(B-I)v=v (2)
Can be obtained by formula (2)
v=B -1x (3)
Can be obtained by formula (3), (1)
y=H 2Fv+n 2=H 2FB -1x+n 2
=H 2FB -1(H 1Wa+n 1)+n 2(4)
=H 2FB -1H 1Wa+H 2FB -1n 1+n 2
Wherein, v is not yet through the signal of relay precoding, and y is the signal that user side receives, and B is relaying feedback matrix, and I is L rank unit matrixs, and F is the optimum pre-coding matrix of relaying, H 2, n 2channel condition information and interchannel noise between relaying and user between relaying and user respectively, v ∈ C l × 1, y ∈ C n × 1, B ∈ C l × L, F ∈ C l × L, H 2∈ C n × L, n 2∈ C n × 1.
Ask for the optimum pre-coding matrix F of relaying, relaying feedback matrix B in described 4th step, user side weighting diagonal matrix G is obtained by following steps:
1. ask for the equivalent channel matrix H between relaying and user, equivalent system is non-linear single-hop structure;
y=H 2FB -1H 1Wa+H 2FB -1n 1+n 2=HFB -1a+n (5)
Wherein, equivalent channel matrix is H=H 2h 1w, the equivalent noise between relaying and user is n=H 2fB -1n 1+ n 2
2. the equivalent channel matrix H couple between relaying and user carries out ORTHOGONAL TRIANGULAR DECOMPOSITION (Orthogonal-Triangular Decomposition, QR) computing, matrix H is resolved into an orthonomal matrix Q and upper triangular matrix R, ask for the optimum pre-coding matrix F of relaying, relaying feedback matrix B, user side weighting diagonal matrix G.
H'=QR (6)
Wherein, H' is the conjugate transpose operation to matrix H, and Q is an orthogonal matrix, and R is a upper triangular matrix.
G = diag [ r 11 - 1 , r 22 - 1 . . . , r kk - 1 ] - - - ( 7 )
B=GR' (8)
F=Q' (9)
R 11, r 22,, r kKbe the 1st row the 1st row of upper triangular matrix R, the 2nd row the 2nd arranges ..., the element that line K K arranges, diag [] builds diagonal matrix computing, and R' is the conjugate transpose operation to matrix R, and Q' is the conjugate transpose operation to matrix Q.
Beneficial effect of the present invention:
1. the present invention changes base station constant in comparatively slow or a period of time to relay well channel H relatively at channel condition information 1adopt linear predictive coding, and be relayed to channel H between user channel condition information is fast-changing 2adopt nonlinear precoding, because the movement of user can cause the changeable of communication channel state information, nonlinear precoding can realize the launching beam and the own channel vector matching that make each user, thus strengthens signal power, realizes larger efficient channel power;
2. under identical signal to noise ratio condition, than traditional relay structure, adopt relaying pre-coding scheme of the present invention, the systematic bits error rate is lower, namely improves communication system performance further.
Accompanying drawing explanation
Fig. 1 is non-linear relaying pre-coding scheme-LT relaying (the Linear Source and Tomlinson – Harashima Relay) structure principle chart that the present invention adopts;
Fig. 2 is linear link structure principle chart;
Fig. 3 is traditional non-linear relaying (THP relaying) structure principle chart;
Fig. 4 is the system equivalent structure figure of the present invention the 4th step;
Fig. 5 illustrates simulated conditions of the present invention;
Fig. 6-Fig. 7 is fixing H 1under channel SNRs condition, linear link, THP relaying, the comparison diagram of LT relaying three kinds of precoding algorithms bit error rate performances.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described with enforcement.
Shown in composition graphs 1, Fig. 2, Fig. 3, the structure principle chart of the LT relaying pre-coding scheme that Fig. 1 carries for the present invention.Wherein, H 1, n 1base station and relay well channel and this interchannel noise respectively, H 2, n 2be the noise of channel and this channel between relaying and user respectively, W is base station pre-coding matrix, and F is relaying pre-coding matrix, and B is feedback matrix, and I is unit matrix, and G is weighting diagonal matrix.A is system input signal, and x is the signal that hop receives, and v is not yet through the signal of relay precoding, and y is the signal that user side receives.
In the present embodiment, base station, relaying have M, L root antenna respectively, there is N number of single-antenna subscriber, M=L=N=2, symbol waiting for transmission is the QPSK modulation symbol of stochastic generation, and channel is smooth Ruili (Rayleigh) decline, and channel condition information is completely known, the reception noise of relaying and user's receiving terminal is the white complex gaussian noise of zero mean unit variance
The present embodiment comprises the following steps:
The first step: signal a, after precoding processing, sends to relaying by base station.
Signal precoding process in the described first step is:
x=H 1Wa+n 1
Wherein, a is initialize signal data flow, and W is base station pre-coding matrix, H 1, n 1base station and relay well channel and this interchannel noise respectively, a ∈ C m × 1, W ∈ C m × M, H 1∈ C l × M, n 1∈ C l × 1, M, L are the antenna number of base station and relaying respectively.
A ∈ C in this example 2 × 1, H 1the signal to noise ratio of channel is SNR 1={ 20,30} (unit: decibel dB), M=L=2.
Second step: relaying obtains channel matrix H between base station and relaying according to the signal x received 1, ask for the optimum pre-coding matrix W in base station based on ZF (Zero Forcing, ZF) criterion.
The method asking for the optimum pre-coding matrix W in base station in described second step is: relaying, according to ZF (ZF) criterion, is tried to achieve:
W = K T r ( H 1 - 1 H 1 - T ) H 1 - 1
Wherein, K=min (M, L, N), N are user side antenna number, H 1 -1to matrix H 1get inverse operation, H 1 -Tto matrix H 1get inverse after transpose operation again, Tr (H 1 -1h 1 -T) be to matrix (H 1 -1h 1 -T) track taking computing.
In this example, K=min (M, L, N)=min (2,2,2)=2
3rd step: relaying sends to user after the signal x received is carried out non-linear preliminary treatment (delivery and feedback).
The non-linear preliminary treatment of relaying to signal in described 3rd step obtains according to following formula:
x-(B-I)v=v
v=B -1x
y=H 2Fv+n 2=H 2FB -1x+n 2
=H 2FB -1(H 1Wa+n 1)+n 2
=H 2FB -1H 1Wa+H 2FB -1n 1+n 2
Wherein v is not yet through the signal of relay precoding, and y is the signal that user side receives, and B is feedback matrix, and I is L rank unit matrixs, and F is the optimum pre-coding matrix of relaying, H 2, n 2channel condition information and channel interchannel noise between relaying and user between relaying and user respectively, v ∈ C l × 1, y ∈ C n × 1, B ∈ C l × L, F ∈ C l × L, H 2∈ C n × L, n 2∈ C n × 1, N is the antenna number of user side.
N=2, H in the present embodiment 2the signal to noise ratio of channel is SNR 2=[0:5:30] (unit: decibel dB).
4th step: user obtains the equivalent channel matrix H between relaying and user according to the signal y received, calculates relaying optimum code matrix F, relaying feedback matrix B, user side weighting diagonal matrix G.
The method asking for the optimum pre-coding matrix F of relaying in described 4th step is:
1. ask for the equivalent channel matrix H between relaying and user, equivalent system is non-linear single-hop structure, as shown in Figure 4;
y=H 2FB -1H 1Wa+H 2FB -1n 1+n 2=HFB -1a+n,
Wherein, H=H 2h 1w, n be between relaying and user equivalent noise, n=H 2fB -1n 1+ n 2;
2. the equivalent channel matrix H couple between relaying and user carries out ORTHOGONAL TRIANGULAR DECOMPOSITION (Orthogonal-Triangular Decomposition, QR) computing, matrix H is resolved into an orthonomal matrix Q and upper triangular matrix R, ask for the optimum pre-coding matrix F of relaying, relaying feedback matrix B, user side weighting diagonal matrix G.
H'=QR
Wherein, H' is the conjugate transpose operation to matrix H, and Q is an orthogonal matrix, and R is a upper triangular matrix.
G = diag [ r 11 - 1 , r 22 - 1 . . . , r kk - 1 ]
B=GR'
F=Q'
R 11, r 22,, r kKbe the 1st row the 1st row of upper triangular matrix R, the 2nd row the 2nd arranges ..., the element that line K K arranges, diag [] builds diagonal matrix computing, and R' is the conjugate transpose operation to matrix R, and Q' is the conjugate transpose operation to matrix Q.
It is the simulated conditions that the present embodiment adopts shown in Fig. 5.In the present embodiment, have stochastic generation 10000 secondary channels altogether and realize, every secondary channel all transmits 2000 QPSK modulation symbols in realizing.The effect that the present embodiment obtains can be further illustrated by the concrete data obtained in Fig. 6, Fig. 7 emulation experiment.In order to illustrate that the LT relay structure adopting the present embodiment is relative to the superiority of traditional linear link and the precoding of THP relay structure, Fig. 6, Fig. 7 sets forth the comparison of performance under three kinds of system configurations, namely fix H 1under channel SNRs condition, error rate of system is with H 2the Performance comparision of channel SNRs change.As shown in the figure, along with H 1the increase of channel SNRs, adopts the precoding algorithms of LT relay structure to manifest advantage (ber curve is minimum) gradually.In addition, at H 2channel low signal-to-noise ratio region (SNR 2≤ 20dB), adopt the systematic function of LT relay structure precoding algorithms to be all better than other two schemes.Consider in practical communication environment, H 2channel signal transmission process is disturbed comparatively large, not easily obtains comparatively high s/n ratio, so the method tool of this example has significant practical applications.

Claims (1)

1., based on a relaying mimo system method for precoding for nonlinear organization, it is characterized in that, comprise the steps:
The first step: signal a, after precoding processing, sends to relaying by base station; The signal a sent from base station adds noise n 1after, form signal x;
Second step: relaying obtains the channel condition information of channel between base station and relaying according to the signal x received, relaying asks for the optimum pre-coding matrix W in base station based on Zero Forcing simultaneously, wherein, the channel condition information matrix H of channel between base station and relaying 1represent;
3rd step: relaying by the signal x received through delivery and feedback, carry out non-linear preliminary treatment after send to user, signal x adds noise n 2after, form signal y;
4th step: user obtains the channel condition information of the equivalent channel between relaying and user according to the signal y received, calculates the optimum pre-coding matrix F of relaying, relaying feedback matrix B, user side weighting diagonal matrix G; Wherein, the channel condition information of the equivalent channel between relaying and user represents by matrix H;
Wherein, the signal x in the first step adopts following formula to obtain:
x=H 1Wa+n 1
Wherein, a is initialize signal data flow, and W is the optimum pre-coding matrix in base station, H 1, n 1base station and relay well channel condition information and base station and relay well interchannel noise respectively; A ∈ C m × 1, W ∈ C m × M, H 1∈ C l × M, n 1∈ C l × 1, x ∈ C l × 1, M, L are the antenna number of base station and relaying respectively;
Wherein, the method asking for the optimum pre-coding matrix W in base station in second step is, relaying is according to Zero Forcing, try to achieve according to following formula:
W = K Tr ( H 1 - 1 H 1 - T ) H 1 - 1
Wherein, K=min (M, L, N), N are user side antenna number, H 1 -1to matrix H 1get inverse operation, H 1 -Tto matrix H 1get inverse after transpose operation again, Tr (H 1 -1h 1 -T) be to matrix (H 1 -1h 1 -T) track taking computing;
Wherein, the non-linear preliminary treatment of relaying to signal in the 3rd step obtains according to following formula:
x-(B-I)v=v
v=B -1x
y=H 2Fv+n 2=H 2FB -1x+n 2
=H 2FB -1(H 1Wa+n 1)+n 2
=H 2FB -1H 1Wa+H 2FB -1n 1+n 2
Wherein, v is not yet through the signal of relay precoding, and y is the signal that user side receives, and B is relaying feedback matrix, and I is L rank unit matrixs, and F is the optimum pre-coding matrix of relaying, H 2, n 2channel condition information and interchannel noise between relaying and user between relaying and user respectively, v ∈ C l × 1, y ∈ C n × 1, B ∈ C l × L, F ∈ C l × L, H 2∈ C n × L, n 2∈ C n × 1;
Ask for the optimum pre-coding matrix F of relaying, relaying feedback matrix B in 4th step, user side weighting diagonal matrix G is obtained by following steps:
5-1) ask for the equivalent channel matrix H between relaying and user, equivalent system is non-linear single-hop structure:
y=H 2FB -1H 1Wa+H 2FB -1n 1+n 2=HFB -1a+n
Wherein, H=H 2h 1w, n are the equivalent noise between relaying and user, n=H 2fB -1n 1+ n 2;
5-2) ORTHOGONAL TRIANGULAR DECOMPOSITION computing is carried out to the equivalent channel matrix H between relaying and user, matrix H is resolved into an orthonomal matrix Q and upper triangular matrix R, asks for the optimum pre-coding matrix F of relaying, relaying feedback matrix B, user side weighting diagonal matrix G:
H'=QR
G = diag [ r 11 - 1 , r 22 - 1 . . . , r KK - 1 ]
B=GR'
F=Q';
Wherein, H' is the conjugate transpose operation to matrix H, and Q is an orthogonal matrix, and R is a upper triangular matrix; r 11, r 22..., r kKbe the 1st row the 1st row of upper triangular matrix R, the 2nd row the 2nd arranges ..., the element that line K K arranges, diag [] builds diagonal matrix computing, and R' is the conjugate transpose operation to matrix R, and Q' is the conjugate transpose operation to matrix Q.
CN201310090457.1A 2013-03-20 2013-03-20 A kind of relaying mimo system method for precoding based on nonlinear organization Expired - Fee Related CN103124190B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310090457.1A CN103124190B (en) 2013-03-20 2013-03-20 A kind of relaying mimo system method for precoding based on nonlinear organization

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310090457.1A CN103124190B (en) 2013-03-20 2013-03-20 A kind of relaying mimo system method for precoding based on nonlinear organization

Publications (2)

Publication Number Publication Date
CN103124190A CN103124190A (en) 2013-05-29
CN103124190B true CN103124190B (en) 2015-09-09

Family

ID=48455061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310090457.1A Expired - Fee Related CN103124190B (en) 2013-03-20 2013-03-20 A kind of relaying mimo system method for precoding based on nonlinear organization

Country Status (1)

Country Link
CN (1) CN103124190B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102594524A (en) * 2012-04-07 2012-07-18 西安电子科技大学 Orthogonal space-time block code transmission method based on an optimal relay linear weighting matrix

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE0403218D0 (en) * 2004-12-30 2004-12-30 Ericsson Telefon Ab L M Method and apparatus related to communication

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102594524A (en) * 2012-04-07 2012-07-18 西安电子科技大学 Orthogonal space-time block code transmission method based on an optimal relay linear weighting matrix

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
On the Design of ZF and MMSE Tomlinson-Harashima Precoding in Multiuser MIMO Amplify-and-Forward Relay System;Hyungjong Kim等;《IEEE 20th International Symposium on Personal,Indoor and Mobile Radio Communications》;20091231;第2509-2513页 *

Also Published As

Publication number Publication date
CN103124190A (en) 2013-05-29

Similar Documents

Publication Publication Date Title
KR102508858B1 (en) Method and Apparatus for Transmitting Diversity
CN101409576A (en) Pre-encoding method and system for managing multiuser radio communication system
US8594215B2 (en) MIMO system having a plurality of service antennas for data transmission thereof
CN102714506B (en) System and method for transceivers in a wireless network
US8218668B2 (en) Multi-user MIMO method, apparatus and communication system for transmitting signal
Dampahalage et al. Weighted-sum-rate maximization for an reconfigurable intelligent surface aided vehicular network
El Hassan et al. Performance assessment of linear precoding for multi-user massive MIMO systems on a realistic 5G mmWave channel
Maimaiti et al. A low-complexity algorithm for the joint antenna selection and user scheduling in multi-cell multi-user downlink massive MIMO systems
CN104852878A (en) Downlink multiuser MIMO (Multiple-Input Multiple-Output) system pre-encoding method capable of lowering complexity based on sum mean square error minimum principle
US9072058B2 (en) Method and apparatus for power optimization in wireless systems with large antenna arrays
CN105262528A (en) Four-antenna transmit diversity method based on weighted fractional Fourier transformation domain
CN103124190B (en) A kind of relaying mimo system method for precoding based on nonlinear organization
Upadhya et al. Tracking abruptly changing channels in mmWave systems using overlaid data and training
CN108683440A (en) The improved modular algebra method for precoding of greedy algorithm is based under a kind of high-speed mobile
CN101521531B (en) Method of enhancing down-link performance of TD-LTE by antenna selection
EP4264849A1 (en) Transceiver method between receiver (rx) and transmitter (tx) in an overloaded communication channel
KR20230013542A (en) MIMO system including intelligent reflective surface and optimal phase transformation matrix acquisition method for intelligent reflective surface to improve sum capacity of channels applied to it
Devi et al. A coherent hybrid precoding for homogenize millimeter-wave multiple-in multiple-out systems for 5G communication
Zemen et al. Combating massive MIMO channel aging by orthogonal precoding
Hawisa et al. Modeling and Analyzing the BER and SER Performances of MIMO System using Multi Level-DPSK Modulation Technique Under Rayleigh Fading Channel
CN112910524A (en) Millimeter wave beam tracking method based on probability model
Wang et al. To cooperate or not: A capacity perspective
CN109067444A (en) Improved modular algebra precoding algorithms are dispatched based on user under high-speed mobile scene
Sun Spectral efficiency and energy efficiency in massive MIMO systems
EP1985051A1 (en) Method of multiple-antenna communication having improved utilization of channel correlations

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: 20150909

Termination date: 20160320

CF01 Termination of patent right due to non-payment of annual fee