CN1325580A - 尤其用于偏移调制模式的补偿方法 - Google Patents

尤其用于偏移调制模式的补偿方法 Download PDF

Info

Publication number
CN1325580A
CN1325580A CN99812847A CN99812847A CN1325580A CN 1325580 A CN1325580 A CN 1325580A CN 99812847 A CN99812847 A CN 99812847A CN 99812847 A CN99812847 A CN 99812847A CN 1325580 A CN1325580 A CN 1325580A
Authority
CN
China
Prior art keywords
sigma
dfe
receiver
received signal
coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN99812847A
Other languages
English (en)
Other versions
CN1144430C (zh
Inventor
J·奥尔迪斯
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.)
Ascom Intermediate AG
Original Assignee
Ascom Powerline Communications AG
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 Ascom Powerline Communications AG filed Critical Ascom Powerline Communications AG
Publication of CN1325580A publication Critical patent/CN1325580A/zh
Application granted granted Critical
Publication of CN1144430C publication Critical patent/CN1144430C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/01Equalisers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03012Arrangements for removing intersymbol interference operating in the time domain
    • H04L25/03114Arrangements for removing intersymbol interference operating in the time domain non-adaptive, i.e. not adjustable, manually adjustable, or adjustable only during the reception of special signals
    • H04L25/03146Arrangements for removing intersymbol interference operating in the time domain non-adaptive, i.e. not adjustable, manually adjustable, or adjustable only during the reception of special signals with a recursive structure

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Filters That Use Time-Delay Elements (AREA)

Abstract

本发明涉及在一台数字接收机中,用一种DFE结构对接收信号加以补偿的一种方法。接收信号根据一个一维信号区;或变为这样一种可转换的信号区(例如:BPSK、GMSK、OQPSK)。规定了DFE系数,使接收信号误差平方的实数部分希望值降至最小。而技术现状则相反,它不能用复数值的误差作为最佳化的基础。计算往往限于实数值。滤波系数同样可以用实数值来代替复数值。进步就在于根据这种原理上很简单的方法;能改善DFE结构的性能,与目前的技术相比,也可节省计算的成本。

Description

尤其用于偏移调制模式的补偿方法
技术领域
本发明涉及在一台数字接收机中,用一种DFE结构(决策反馈均衡器)对接收信号加以补偿的一种方法,这时接收信号系根据一维信号区或一种可转换的信号区。
技术现状
典型的情况和GSM(全球移动通信系统)、HIPERLAN高性能无线电局域网)、DECT(无线电话数据交换)等一样,出现的传输信道之特点是对多路传送有干扰作用。
众所周知,使用一个决策反馈均衡器(DFE,Decision Feedback Equalizer)是为了补偿在一个数字通信系统中受到一个线性选频过程阳如在一个无线电通道中多路传送)干扰的信号。
一个DFE的性能取决于用于计算或确定前馈和反馈部分中滤波系数的质量。对于未知的通道,一般是通过适合性训练来确定系数的。如果相反,一个频道的脉冲响应是已知的,那么则可以从中推导出DEF的最佳系数。
一个DEF的结构本身是非常简单的,所以使用也很随意。但总是不能达到所希望的性能。
对本发明的表述
本发明的任务是说明以上所述方式的一种方法,根据已知或事先估计的通道脉冲响应,以最小的计算成本来确定最佳的系数,同时;按照技术现状来实现高于常见DFE的性能。
由专利权限1的特点确定了本发明任务的实现。根据本发明,规定了DFE的系数,使误差二次实数部分的期望值最小。
与技术现状不同;也可以不用复杂的误差作为最佳化的基础。计算更多地限于实数值。反馈滤波器的滤波系数不太复杂,前馈滤波器的滤波系数也是一般的。跳越点仅在于以原理上很简单的方式应能改善DFE结构的性能,甚至还能降低对技术现状的计算成本。
对于一个H进制的BPSK信号,最好根据以下作进一步解释的公式(Ⅰ)和(11)来计算系数。
本发明不仅适用于BPSK信号(BPSI一二进位相键穆控法),还可适用于GMSK和OQPSK调制法(GMSK一高斯最小移动键控法,OQPSK一偏移正交相位移动键控法)。就是虽然具有一个二维信号区;但(用一个适当的变换)可以以一个(至少是近似的)等量一维表达式表示的也可视作是一维的调制法。
本发明方法在电路上实现转换毫无困难。一般可以在一个处理器或ASIC中进行编程计算。
本发明例如适用于一个HIPERLAN系统。(例如由Ascom Tech公司的EP0795976A2专利中得到一种特别不错的系统结构)。一个欧洲通迅标准(ETS)规定了一种高性能无线局域网(HIPERLAN)的技术特点。HIPERLAN是一种具有大数据流通量的短程通讯子系统(参见ETSI 1995,ETS 300652,UDC:621396)。ETS-HIPERLAN标准规定的频带为5.15~5.30GHZ。
从下列详细说明和整个专利权限中可以看出本发明的更为优越的设计型式和性能的组合。
对图纸的简述
图1:一个DFE的示意图;
图2:一种结构示例的示意图;
图3:本发明方法的性能与技术现状的对比图;
图4a~c:技术现状与本发明的误差情况对照表;
图5:一台BPSK接收机的示意图;
图6:一台GMSK接收机的示意图。
实现本发明的方法
以下通过与技术现状的比较来介绍本发明的原理。
图1为一种常见的DFE分程序结构。由载体向下调制的接收信号I被送入DFE的一个前馈滤波器FF中。然后再同通过反馈滤波器FB送回DD(决策设备)的估计信号结合(添加)。于是将信号I加在决策器DD输入端上。根据技术现状,对系数f和g(在以上情况下是指有多个系数分量的矢量)计算如下: min E f . g { | I ~ - I ^ | 2 } - - - - ( A ) 与比相反,本发明采用下式进行优化: min E f . g { ( Re ( I ~ - I ^ ) ) 2 } - - - - ( B )
与技术现状的差别在于对滤波系数的计算方式。而DEF的其他结构保持不变。然后通过结构示例作详细介绍。
图2表明了一个DFE的具体示例。经DFE处理的信号,同现代相干数字接收机一样;都是通过复数来表示。这时实数代表同相成分;而虚数则代表正交分量。图2所示的DEF,根据一般理解,是一个复数系数和复数数据。
如果对误差的实数值按MMSE准则进行最优化(MMSE=最小均方差),则应通过下列方程求出前馈滤波系教:(1) h M + 1 - i R = σ 2 2 f i R + Σ m = 1 M f m R Σ n = 1 M h n + 1 - i R h n + 1 - m R - Σ m = 1 M f m 1 Σ n = 1 M h n + 1 - i R h n + 1 - m 1 - h M + 1 - i 1 = σ 2 2 f i 1 - Σ m = 1 M f m R Σ n = 1 M h n + 1 - i 1 h n + 1 - m R + Σ m = 1 M f m 1 Σ n = 1 M h n + 1 - i 1 h n + 1 - m 1
这是wh 2M实数方程(1≤i≤M)。指数太大或大小的系数在此取0值。对长度L的矢量,指数从1到L。滤波系数值可用常用的线性方程来解。以下再对这个标准方法作进一步阐述。
通过下式求出反馈滤波系数: ( II ) g i - M R = - Σ m = 1 M f m R h i + 1 - m R - f m 1 h i + 1 - m 1 这是个N-1的方式;因为M+1≤i≤M+N-1。式(Ⅰ)和(Ⅱ)根据下列惯用法:N通道脉冲响应;M前馈滤波器的长度;hi R通道脉冲响应的实数部分,1≤i≤M;hi 1通道脉冲响应的虚数部分,1≤i≤N;fi R DFE前馈部件滤波系数的实数部分,1≤i≤M;fi 1 DFE前馈部件滤波系数的虚数部分,1≤i≤M;gi R DFE反馈部件滤波系数的实数部分,1≤i≤N-1;
σ2 DFE输入的噪声功率(噪声功率的实数部分和虚数部分的组合)。如果这
个值是未知的,则可将它设为一个常数;才不会大大降低性能。
在大多数情况下,M=N。N<M并无益处。如果N>M,这个复数则会降低到性能负载。但尽管如此;本发明的计算依据平均方差仍能得出最佳滤波系数。
反馈滤波器的长度等长或比通道响应的长度要短(即N-1)。如果选择较大的长度,附力。抽头的系数都会变为o。而在决策器的输入上,较短的长度则会引起码间干扰。由于对反馈滤波器增加了抽头;并不会大大提高总体的复杂程度;所以一般都采用全长。
反馈滤波器的系数没有虚数部分。因为反馈滤波器的输入为实数,输出也是实数。(不考虑决策器输入的虚数部分)。
本发明对滤波器系数的计算适用于不同的用途。以下则表明如可能改进一个HPERLAN接收机的性能。这时可将一个已知复数值的MMSE方法与本发明实值的*MSE方法进行对比。此外还需假定接收器要进行一次三根天线的选择分集。对相应接收器的模拟可以估计出程序包的误差率。
假定在接收器中参数σ2为10dB,并力。在接收到的信号功率上。此外,以45ns或75 ns的延迟扩展从无线电通道中输出。DFE便有8个前情抽头和7个反馈抽头。
图3所示的结果表明,使用本发明的二种计算方法有了显著地改善。长时间延迟扩展(75ns)的误差率较高。在20dB信噪比和45ns延迟扩展时;可确定在可测阈值以下的误差率。
本发明方法的效应可从图4a-4c上看清。如用QPSK作为调制方法,决策器就会根据哪个是决策器下一个输入值来输入四个复数值1+j,1-j,-1+j;-1-j中的-个复数值。由于噪声和未消除的简便码间干扰会使输入值引起失真。在图4a~c中将其以云状分布的形式表示。
降低复数方差可使得在每个区点四周产出环状分布,如图4a和4b所示。与止匕相反;本发明减少方差实数部分则是形成准椭圆形分布(图4C)。从复平面上看,(复数值)方差的平均值要大于技术现状(图4a,b)。但误差却在虚数轴上。在实数轴中;它要小于技术现状。但当决策器的输出只能是实数值时;虚数轴上较大的误差则不起作用。
图5表明,在一个BPSK接收器中如何来采用本发明。在一台发射器中;数据1受进行载体振荡的一个BPSK调制器2调制。在一台接收器中,解调器3用于将接收信号转换成频率基带并进行相应滤波。然后以代码速度测试信号(扫描器4)。扫描器的输出一方面要经通道评估器5;另一方面又要经DFE 7的处理。本发明对系数的计算在系数计算器6中进行。DEF 7的输出端加上传递的数据8。接收器的结构是常见的那种。而新的则是上述在系数计算器6中如何测定系数的方式方法。
原则上本发明也可用于一种QPSK方法(这时必须进行相应的进行调制解调)。DFE与BPSK接收器不同,在任何情况下都必须采用复数。
图6为一个GMSK传递方法的一般示意图。数据9在发射器一侧经一个预编码器10,以常用方式进行预编码,再用一个载体振荡的GMSK调制器11进行调制。在一台接收器中,解调器12用于将接收信号转换成频率基带以及进行相应的过滤。然后以码速扫描信号(扫描器13)。将扫描器的输出乘以一个相位系数ji(移相器14,乘法器15);然后一方面经通道评估器16,另一方面又经DFEIS的处理。本发明对系数的计算是在系数计算器门中进行的。DEFIS的输出端加上传递的数据19。接收器的结构也是常见的那种。而新的则是上述在系数计算器6中如何测定系数的方式方法。
以下就本发明如何用于GMSK和OQPSK调制方法,其中哪一个在第一时刻出现二维信号区加以说明。
对于具有bk∈[-1,+1],k=…-1,0,1,2…代码的一个二进制数据流,在复数基带表示法中可按下式来表示GMSK调制信号: ( III ) S 0 ( t ) = Aexp [ jπ 2 Σ k b k ∫ - ∝ 1 - kT g ( τ ) dτ + φ 0 ]
A和φ0表示振幅或载波相位;g(τ)是(高斯分响应)脉冲,它决定了相位的调制,T代表代码或位的时间。
根据脉冲g(τ),通过下列线性部分响应QAM信号能很好地逼近调制的信号。 ( IV ) S o ~ ( t ) = Aexp ( jφ o ) Σ k α k g ~ ( t - kT )
式中,αk项是复数值数据代号,它只与bk有关;数值范围为+1,-1,+j;-j。g(t)为一个部分响应的脉冲生成函数,即: ( V ) α k = exp ( jπ 2 Σ m = - ∝ k b N )
众所周知,(Baier,A等所著“Bit Synchronization and Timing SenitivltyIn Adaptivev Viterbi Equallzers for Narrowband-TDMA Digital Mobile RadioSystems”,IEEE 1988,ch2622-9/8/0000-0377);上述逼近式能很好地用于采用GSM和HIPERLAN,以0.3时间-宽带积进行的GMSK调制。这个逼近精确地等于从数值范围[+1,-1,+j,-jl以数据代码进行的线性QAM调制。它的和 交替地为偶数和奇数,从而使传递的代码ak也交替地为实数和虚数。这种调制称作为OQPSK偏移正交相键控)。ak和bk之间的过渡极为简单。必须指出,从ak向bk的过渡不易出现误差,而反相传递则不是如此。bk序列的唯一误差只可能是由推导ak。序列中很多(可能是无穷多的)误差所致。
在接收器中必须反馈传递的ak。以下假定发射机与接收器中具有同一种框架同步。从第一个ak中我们知道;它为实数值(也就是说非+1,即-1)。如果第一个为虚数,则要求对下列公式略作修正。传递的信号为 。接收信号为r(t),需用通道脉冲响应和接收器的模拟滤波顺的卷积积分来表示: ( VI ) r ( t ) = A Σ k α k h ( t - kT ) 式中,h(t)为传递信号。
Figure A9981284700084
、初始相移、通道脉冲响应和接收器一侧滤波器总体脉冲响应的的卷积积分。
在接收器中,根据通道的码率扫描复数值的基带信号,产生一个时间离散的信号。它可以写成下列: ( VII ) r ^ i = A Σ k α k h ( iT + λ - kT )
设定扫描相位为λ。对总体不加限制时,可设λ=0,因为时间延迟总是包括在通道脉冲响应中的。
向DFE送入信号之间,应将它乘以j-i相位:
Figure A9981284700086
Figure A9981284700087
Figure A9981284700088
ck为从ak中推导出的数序。必须注意;j-1相只能取[+1,-1,+j,-j]值。于是只要将接收信号乘以这个相即可(参见图6中乘法器14)。
Figure A9981284700091
如果具有一个框架同步,则可避免其中的一种情况。因此忽略第二种可能。可以看出,从接收器一侧测得的信号值是具有某一函数…)实数数序C内卷积积分,它包括:
一调制脉冲的生成;
一通道脉冲的响应;
一载波信号的初始相;
一扫描的时间偏移;
一接收器中以1湘的旋转。
例如在接收器中,可用一个训练序列和一个关系计算式来求出这个函数。在接收器中就是用这个函数来计算DFE的滤波系数。DFE只能产生一个实数输出,因为依据的数据都为实数(ck)。最后(知道指数k)求出初始数据ak
如上所述,通过OQPSK调制可以很好地逼近GMSK调洲前提是时间一带宽积是已知的并在ak和bk之间进行数据转换)。用这种方法,本发明的DFE也可用于GMSK和OPQSK。它只需要对数据进行一次附加的、但极为简单的转换即可。如果在发射器中,预编码用于GMSK调制前,那应更为简单了。
当时间一带宽的积处于最不利的情况时,本发明对GMSK的补偿所产生的、使性能变差要比OQPSK来得轻,尽管在数据转换后GMSK并不是严格的线性。当然,如果时间带宽的积处于一般的数量级时,这种劣化则可忽略不计。
总而言之;实验中采用本发明可以成为用途极广的一维调制方法,采用这种杰出的DFB结构可以改善补偿的效果。可用实数值(英语:real values)代替复数值(英语:complex values)在反馈滤波器中进行计算。即使前馈滤波器的输出才只需要用实数值。在这种滤波器中只需要相应地对输出的实数值进行计算。本发明的接收机例如可以是GSM电话或无绳DECT电话机,或影民据HIPERLAN,用于计算机之间61数据通讯。

Claims (4)

1.在一台数字接收机中用一种DFE即决策反馈均衡器的结构来补偿接收信号的方法,这时接收信号系根据一个一维信号区;或变为这样一种可转换的信号区,其特征在于,规定了DFE系数;使接收信号误差平方的实数部分希望值降至最小。
2.如权利要求1所述的方法,其特征在于,信号区相当于一次BPSK调制;可按下式确定DFE的系数: ( I ) h M + 1 - i R = σ 2 2 f i R + Σ m = 1 M f m R Σ n = 1 M h m + 1 - i R h n + 1 - m R - Σ m = 1 M f m 1 Σ n = 1 M h n + 1 - i R h n + 1 - m 1 - h M + 1 - i 1 = σ 2 2 f i 1 - Σ m = 1 M f m R Σ n = 1 M h n + 1 - i 1 h n + 1 - m R + Σ m = 1 M f m 1 Σ n = 1 M h m + 1 - i 1 h n + 1 - m 1 ( II ) g i - M R = - Σ m = 1 M f m R h i + 1 - m R - f m 1 h i + 1 - m 1
3.如权利要求1所述的方法,其特征在于,信号区相当于一次GMSK调制或一次OQPSK调制,在接收器中以j-1相旋转测试值,在此i表示测试值的顺序代号。
4.用一个DFE来补偿一台数字接收器中接收信号的电路,这时接收信号根据一个一维信号区,或变为这样一种可转换的信号区,其特征在于,具有一个用于计算DFE系数的电路,使接收信号误差平方的实数部分希望值降至最小。
CNB998128473A 1998-10-30 1999-10-29 尤其用于偏移调制模式的补偿方法 Expired - Fee Related CN1144430C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP98811090A EP0998083A1 (de) 1998-10-30 1998-10-30 Verfahren zum Entzerren, insbesondere für Offset-Modulationsarten
EP98811090.4 1998-10-30

Publications (2)

Publication Number Publication Date
CN1325580A true CN1325580A (zh) 2001-12-05
CN1144430C CN1144430C (zh) 2004-03-31

Family

ID=8236418

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB998128473A Expired - Fee Related CN1144430C (zh) 1998-10-30 1999-10-29 尤其用于偏移调制模式的补偿方法

Country Status (12)

Country Link
EP (2) EP0998083A1 (zh)
JP (1) JP2002529964A (zh)
KR (1) KR20010089314A (zh)
CN (1) CN1144430C (zh)
AU (1) AU6322599A (zh)
BR (1) BR9914950A (zh)
CA (1) CA2347747A1 (zh)
HK (1) HK1040025A1 (zh)
ID (1) ID29202A (zh)
IL (1) IL142795A0 (zh)
NO (1) NO20012104D0 (zh)
WO (1) WO2000027083A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101836407B (zh) * 2007-12-10 2013-03-27 中兴通讯股份有限公司 终端射频性能的优化方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0614627B2 (ja) * 1985-06-04 1994-02-23 富士通株式会社 モデムのトレーニング方法
CA2268272A1 (en) * 1996-10-10 1998-04-16 Statistical Signal Processing, Inc. Signal processing apparatus employing the spectral property of the signal

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101836407B (zh) * 2007-12-10 2013-03-27 中兴通讯股份有限公司 终端射频性能的优化方法

Also Published As

Publication number Publication date
HK1040025A1 (zh) 2002-05-17
JP2002529964A (ja) 2002-09-10
IL142795A0 (en) 2002-03-10
EP1125404A1 (de) 2001-08-22
WO2000027083A1 (de) 2000-05-11
EP0998083A1 (de) 2000-05-03
CN1144430C (zh) 2004-03-31
BR9914950A (pt) 2001-11-06
NO20012104L (no) 2001-04-27
CA2347747A1 (en) 2000-05-11
KR20010089314A (ko) 2001-09-29
NO20012104D0 (no) 2001-04-27
AU6322599A (en) 2000-05-22
ID29202A (id) 2001-08-09

Similar Documents

Publication Publication Date Title
CN112737985A (zh) 基于深度学习的大规模mimo信道联合估计和反馈方法
US8290082B2 (en) Method for signal transmission-reception in a multi-user multiple transmit and multiple receive antenna radio communication system
CN1256828A (zh) 通信系统中确定信号对干扰加噪声功率比的方法和系统
US7342872B1 (en) Differential OFDM using multiple receiver antennas
CN113866752B (zh) 一种svd预编码的mimo-otfs雷达目标检测方法
Phan-Huy et al. Receive antenna shift keying for time reversal wireless communications
KR20040078151A (ko) 멀티캐리어 확산-스펙트럼 신호들의 수신
US20070206697A1 (en) Signal receiving method and signal receiving equipment for multiple input multiple output wireless communication system
CN110278167B (zh) 一种对iq不平衡进行连续估计与补偿的无线通信方法
US11996899B2 (en) Method of discrete digital signal recovery in noisy overloaded wireless communication systems in the presence of hardware impairments
Yeh et al. New parallel algorithm for mitigating the frequency offset of OFDM systems
Golden et al. Transmitter design for data transmission in the presence of a data-like interferer
US20050073946A1 (en) Transmitter and receiver for use with an orthogonal frequency division multiplexing system
CN1144430C (zh) 尤其用于偏移调制模式的补偿方法
Sundberg Block error probability for noncoherent FSK with diversity for very slow Rayleigh fading in Gaussian noise
US20070183536A1 (en) Method of differential-phase/absolute-amplitude polar qam
WO2021176532A1 (ja) 無線通信システム、無線通信方法及び送信装置
US7440527B2 (en) Method of differential-phase/absolute-amplitude QAM
Sangirov et al. An enhanced semi-blind channel estimation for MIMO-OFDM systems
Babulkar A Comprehensive Review on Channel Estimation Techniques in MIMO-OFDM
CN111917441B (zh) 一种基于大规模多输入多输出毫米波系统的信道估计方法
Muthukumaran et al. ANN based Quick and Accurate Estimate of Channel Availability by using Hybrid M1-MMSE Equalizer
Kafedziski Joint source channel coding of images over frequency selective fading channels with feedback using dct and multicarrier block pulse amplitude modulation
Fines et al. Synchronization techniques for all digital 16-ary QAM receivers operating over land mobile satellite links
AU668166B2 (en) Optimum estimation of data and parameters for communication systems

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
C19 Lapse of patent right due to non-payment of the annual fee
CF01 Termination of patent right due to non-payment of annual fee
REG Reference to a national code

Ref country code: HK

Ref legal event code: WD

Ref document number: 1040025

Country of ref document: HK