CN101626354A - 多径信道的噪声估计方法和装置 - Google Patents
多径信道的噪声估计方法和装置 Download PDFInfo
- Publication number
- CN101626354A CN101626354A CN200910084280A CN200910084280A CN101626354A CN 101626354 A CN101626354 A CN 101626354A CN 200910084280 A CN200910084280 A CN 200910084280A CN 200910084280 A CN200910084280 A CN 200910084280A CN 101626354 A CN101626354 A CN 101626354A
- Authority
- CN
- China
- Prior art keywords
- channel
- footpath
- noise estimation
- noise
- multipath channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Mobile Radio Communication Systems (AREA)
- Noise Elimination (AREA)
Abstract
本发明提供了一种多径信道的噪声估计方法,包括以下步骤:对多径信道进行信道估计,得到多径信道的各径信道的时域响应估计值;设置噪声门限;选取各径信道的所有径中其功率低于噪声门限的径,并根据选取的径的时域响应估计值计算多径信道的噪声估计值。本发明克服了相关技术中利用信号子空间的噪声估计方法需要对接收信号进行高阶统计,复杂度较高的问题,进而达到了降低信道噪声估计复杂度的效果。
Description
技术领域
本发明涉及无线通信领域,具体而言,涉及一种多径信道的噪声估计方法和装置。
背景技术
无线通信系统中,信噪比为接收端输出的信号平均功率与噪声平均功率之比,是衡量接收性能的一个重要指标,可用于去除噪声、提高软解调精度、功率控制、自适应编码调制、信道译码等。其中,噪声平均功率可通过对信道进行噪声估计来得到,故信道的噪声估计是无线通信中一项至关重要的技术。
针对多径时变信道的无线通信系统,相关技术提供了一种基于OFDM(Orthogonal Frequency Division Multiplexing,正交频分复用技术)的利用信号子空间的噪声估计方法。
在实现本发明过程中,发明人发现相关技术中利用信号子空间的噪声估计方法需要对接收信号进行高阶统计,复杂度较高。
发明内容
本发明旨在提供一种多径信道的噪声估计方法,能够解决相关技术中利用信号子空间的噪声估计方法需要对接收信号进行高阶统计,复杂度较高的问题。
在本发明的实施例中,提供了一种多径信道的噪声估计方法,包括以下步骤:对多径信道进行信道估计,得到多径信道的各径信道的时域响应估计值;设置噪声门限;选取各径信道的所有径中其功率低于噪声门限的径,并根据选取的径的时域响应估计值计算多径信道的噪声估计值。
另一方面,在本发明的实施例中,还提供了一种多径信道的噪声估计装置,包括:信道估计模块,用于对多径信道进行信道估计,得到多径信道的各径信道的时域响应估计值;设置模块,用于设置噪声门限;计算模块,用于选取各径信道的所有径中其功率低于噪声门限的径,并根据选取的径的时域响应估计值计算多径信道的噪声估计值。
因为上述实施例利用多径信道的各径信道的时域响应估计值,并通过设置噪声门限来完成对信道的噪声估计,所以克服了相关技术中利用信号子空间的噪声估计方法需要对接收信号进行高阶统计,复杂度较高的问题,进而达到了降低信道噪声估计复杂度的效果。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1示出了根据本发明第一实施例的噪声估计方法的流程图;
图2示出了根据本发明第二实施例的噪声估计方法的流程图;
图3示出了根据本发明第三实施例的噪声估计方法的流程图;
图4示出了根据本发明第四实施例的噪声估计方法的流程图;
图5示出了根据本发明第五实施例的噪声估计装置的结构图。
具体实施方式
下面将参考附图并结合实施例,来详细说明本发明。
图1示出了根据本发明第一实施例的噪声估计方法的流程图,包括以下步骤:
步骤S10,对多径信道进行信道估计,得到多径信道的各径信道的时域响应估计值h(l),其中l=0,1,……,L-1;
步骤S20,设置噪声门限PT;
步骤S30,选取各径信道的所有径中其功率|h(l)|2低于噪声门限PT的径,并根据选取的径的时域响应估计值计算多径信道的噪声估计值PN。
上述实施例首先对多径信道进行信道估计来得到多径信道的各径信道的时域响应估计值h(l),即使信道估计结果为信道的频域响应H(k),也可通过IFFT(Inverse Fast Fourier Transform,快速傅立叶逆变换)转化为时域响应,然后设置噪声门限PT,并与各径信道的功率|h(l)|2进行比较,其功率|h(l)|2高于噪声门限PT的径被认为是有效多径,其功率|h(l)|2低于噪声门限PT的径被认为是噪声多径,最后,根据噪声多径的时域响应估计值h(l)计算多径信道的噪声估计值PN。其中,时域响应估计值h(l)的模值|h(l)|可视为径的幅度,多径信道时域响应估计值h(l)的模值平方|h(l)|2可视为径的功率,故通过比较径的功率|h(l)|2与噪声门限PT的大小关系,可判断出该径是否属于噪声。本实施例基于二阶统计特性,根据多径信道的时域响应估计值h(l)得到多径信道的噪声估计值PN,克服了相关技术中利用信号子空间的噪声估计方法需要对接收信号进行高阶统计,复杂度较高的问题,进而达到了降低信道噪声估计复杂度的效果。
优选地,如图2所示,步骤S20具体包括:步骤S201,对时域响应估计值h(l)进行滑动窗处理,计算滑动窗内的各径信道的平均功率;步骤S202,从多个位置的滑动窗的平均功率中选择最小值P1;步骤S203,将最小值P1与第一系数α的乘积αP1设置为噪声门限PT。
滑动窗用于在多径信道的各径信道的时域响应估计值h(l)中滑动,对于其滑动经过的多个位置,均计算此时滑动窗内的各径信道的平均功率,从多个位置的滑动窗的平均功率中选择最小值P1,即
其中,l=k和l=k+Nw-1分别表示平均功率达到最小值P1时滑动窗的起始和结束位置;Nw为滑动窗的窗长,其选取原则是:在保证平均功率达最小时对应的滑动窗不包含有效多径的前提下,Nw的值应尽可能地大;|h(l)|表示对复数h(l)取模。此时得到的最小值P1可视为多径信道噪声估计的初步结果。再将最小值P1与第一系数α的乘积αP1设置为噪声门限PT,即PT=αP1,其中,第一系数α的值,若取得过大会导致将有效多径误判为噪声多径,若取得过小则会导致将噪声多径误判为有效多径。通常情况下,α的取值范围为3~8,其具体取值根据具体情况可作适当调整,并不局限于上述给定的取值范围,比如当窗长Nw的值较小时α可取得略大一些。本实施例利用对多径信道的时域响应估计值进行滑动窗处理来设置噪声门限PT,简单易行,同时通过利用统计特性又兼顾了噪声判定的准确性。
优选地,步骤S20具体包括:将上一次的噪声估计值与第二系数α′的乘积设置为噪声门限PT。
利用上一次的噪声估计值,也可确定噪声门限的设置值,乘积因子调整为第二系数α′,通常情况下,第二系数α′的取值范围为2~5,其具体取值根据具体情况可作适当调整,并不局限于上述给定的取值范围。本实施例直接利用上一次的噪声估计值来设置噪声门限,进一步降低了噪声估计过程的复杂度。
优选地,步骤S30具体包括:选取各径信道的所有径中其功率|h(l)|2低于噪声门限PT的径,计算选取的径的平均功率P2,即
其中,A2为所有满足|h(l)|2<PT的l,N2为A2中元素的个数,并设置P2为多径信道的噪声估计值PN,即PN=P2。
本实施例通过计算被判为噪声多径的各径信道的平均功率P2,来确定多径信道的噪声估计值PN的方案,简单易行的同时也保证了较高的噪声估计精度。
优选地,如图3所示,步骤S30具体包括:步骤S301,选取各径信道的所有径中其功率|h(l)|2低于噪声门限PT的径,并计算选取的径的平均功率P2;步骤S302,对平均功率P2进行低通滤波;步骤S303,设置滤波后的结果P3为噪声估计值PN,即PN=P3。
本实施例中,对被判为噪声多径的各径信道的平均功率P2进行低通滤波,并设置滤波后的结果为噪声估计值PN,由于屏蔽了噪声估计过程中产生的高频干扰,从而提高了噪声估计精度。
优选地,步骤S30具体包括:选取各径信道的所有径中其功率|h(l)|2低于噪声门限PT的径,并计算选取的径的平均功率P2;对平均功率P2与之前一次或多次的噪声估计值进行加权平均;设置加权平均后的结果P′3为噪声估计值PN,即PN=P′3。
本实施例中,计算得到噪声多径的各径信道的平均功率P2后,将其与之前一次或多次的噪声估计值进行加权平均,并设置加权平均后的结果为噪声估计值PN,在噪声估计的其他条件一定时,对多次估计结果进行加权平均,可减小估计过程中产生的偶然误差,从而提高噪声估计精度。
优选地,如图4所示,步骤S10具体包括:步骤S101,对多径信道进行信道估计;步骤S102,对信道估计结果进行低通滤波,并将低通滤波后的结果设置为多径信道的各径信道的时域响应估计值h(l);步骤S30具体包括:步骤S304,选取各径信道的所有径中其功率|h(l)|2低于噪声门限PT的径,并计算选取的径的平均功率P′2;步骤S305,设置平均功率P′2与第三系数β的乘积为噪声估计值PN。
本实施例中,对信道估计结果进行低通滤波,并将低通滤波后的结果设置为多径信道的各径信道的时域响应估计值h(l),同时为了补偿对信道估计结果进行低通滤波而造成的增益损失,需要对判为噪声多径的各径信道的平均功率P′2作相应的调整,即将P′2与第三系数β的乘积设置为噪声估计值PN,即PN=βP′2,其中,β>1,β的具体取值需根据所采用的低通滤波器的特性来确定。本实施例通过低通滤波屏蔽了来自信道估计过程中产生的高频干扰,从而实现了提高噪声估计精度的效果。
优选地,步骤S10具体包括:对多径信道进行信道估计;将当前的信道估计结果与之前一次或多次的信道估计结果进行加权平均,并将加权平均后的结果设置为多径信道的各径信道的时域响应估计值h(l);步骤S30具体包括:选取各径信道的所有径中其功率|h(l)|2低于噪声门限PT的径,并计算选取的径的平均功率P″2;设置平均功率P2与第四系数β′的乘积为噪声估计值PN。
本实施例中,将当前的信道估计结果与之前一次或多次的信道估计结果进行加权平均,并将加权平均后的结果设置为多径信道的各径信道的时域响应估计值h(l),同时为了补偿对信道估计结果进行加权平均造成的增益损失,需要对判为噪声多径的各径信道的平均功率P″2作相应的调整,即将P′2与第四系数β′的乘积设置为噪声估计值PN,即PN=β′P″2,其中,β′>1,β′的具体取值需参考加权系数来确定,本实施例通过将当前的信道估计结果与之前一次或多次的信道估计结果进行加权平均,减小了信道估计过程中的偶然误差,从而提高噪声估计精度。
另一方面,图5示出了根据本发明第五实施例的噪声估计装置50的结构图,该装置包括:
信道估计模块501,用于对多径信道进行信道估计,得到多径信道的各径信道的时域响应估计值;
设置模块502,用于设置噪声门限;
计算模块503,用于选取各径信道的所有径中其功率低于噪声门限的径,并根据选取的径的时域响应估计值计算多径信道的噪声估计值。
上述实施例首先通过信道估计模块501对多径信道进行信道估计来得到多径信道的各径信道的时域响应估计值h(l),即使信道估计结果为信道的频域响应H(k),也可通过IFFT(Inverse Fast FourierTransform,快速傅立叶逆变换)转化为时域响应,然后通过设置模块502设置噪声门限PT,并与各径信道的功率|h(l)|2进行比较,其功率|h(l)|2高于噪声门限PT的径被认为是有效多径,其功率|h(l)|2低于噪声门限PT的径被认为是噪声多径,最后,由计算模块503根据噪声多径的时域响应估计值h(l)计算多径信道的噪声估计值PN。其中,时域响应估计值h(l)的模值|h(l)|可视为径的幅度,多径信道时域响应估计值h(l)的模值平方|h(l)|2可视为径的功率,故通过比较径的功率|h(l)|2与噪声门限PT的大小关系,可判断出该径是否属于噪声。本实施例基于统计特性,根据多径信道的时域响应估计值h(l)得到多径信道的噪声估计值PN,克服了相关技术中利用信号子空间的噪声估计方法需要对接收信号进行高阶统计,复杂度较高的问题,进而达到了降低信道噪声估计复杂度的效果。
从以上的描述中,可以看出,本发明上述的实施例以较低的复杂度实现了较高的噪声估计精度。
显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种多径信道的噪声估计方法,其特征在于,包括以下步骤:
对所述多径信道进行信道估计,得到所述多径信道的各径信道的时域响应估计值;
设置噪声门限;
选取所述各径信道的所有径中其功率低于所述噪声门限的径,并根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值。
2.根据权利要求1所述的噪声估计方法,其特征在于,设置噪声门限,具体包括:
对所述时域响应估计值进行滑动窗处理,计算所述滑动窗内的各径信道的平均功率;
从多个位置的所述滑动窗的所述平均功率中选择最小值;
将所述最小值与第一系数的乘积设置为所述噪声门限。
3.根据权利要求1所述的噪声估计方法,其特征在于,设置噪声门限,具体包括:
将上一次的所述噪声估计值与第二系数的乘积设置为所述噪声门限。
4.根据权利要求1至3中任一项所述的噪声估计方法,其特征在于,根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值,具体包括:
设置选取的径的平均功率为所述多径信道的噪声估计值。
5.根据权利要求1至3中任一项所述的噪声估计方法,其特征在于,根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值,具体包括:
计算选取的径的平均功率;
对选取的径的所述平均功率进行低通滤波;
设置滤波后的结果为所述噪声估计值。
6.根据权利要求1至3中任一项所述的噪声估计方法,其特征在于,根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值,具体包括:
计算选取的径的平均功率;
对选取的径的所述平均功率与之前一次或多次的所述噪声估计值进行加权平均;
设置加权平均后的结果为所述噪声估计值。
7.根据权利要求2所述的噪声估计方法,其特征在于,
对所述多径信道进行信道估计,得到所述多径信道的各径信道的时域响应估计值,具体包括:
对所述多径信道进行信道估计;
对信道估计结果进行低通滤波,并将低通滤波后的结果设置为所述多径信道的各径信道的时域响应估计值;
根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值,具体包括:
计算选取的径的平均功率;
设置选取的径的所述平均功率与第三系数的乘积为所述噪声估计值。
8.根据权利要求2所述的噪声估计方法,其特征在于,
对所述多径信道进行信道估计,得到所述多径信道的各径信道的时域响应估计值,具体包括:
对所述多径信道进行信道估计;
将当前的信道估计结果与之前一次或多次的所述信道估计结果进行加权平均,并将加权平均后的结果设置为所述多径信道的各径信道的时域响应估计值;
根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值,具体包括:
计算选取的径的平均功率;
设置选取的径的所述平均功率与第四系数的乘积为所述噪声估计值。
9.一种多径信道的噪声估计装置,其特征在于,包括:
信道估计模块,用于对所述多径信道进行信道估计,得到所述多径信道的各径信道的时域响应估计值;
设置模块,用于设置噪声门限;
计算模块,用于选取所述各径信道的所有径中其功率低于所述噪声门限的径,并根据选取的径的所述时域响应估计值计算所述多径信道的噪声估计值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910084280 CN101626354B (zh) | 2009-05-20 | 2009-05-20 | 多径信道的噪声估计方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 200910084280 CN101626354B (zh) | 2009-05-20 | 2009-05-20 | 多径信道的噪声估计方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101626354A true CN101626354A (zh) | 2010-01-13 |
CN101626354B CN101626354B (zh) | 2013-06-05 |
Family
ID=41522045
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 200910084280 Active CN101626354B (zh) | 2009-05-20 | 2009-05-20 | 多径信道的噪声估计方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101626354B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102075465A (zh) * | 2011-02-14 | 2011-05-25 | 电信科学技术研究院 | 一种噪声功率的估计方法及设备 |
CN102208930A (zh) * | 2010-03-30 | 2011-10-05 | 清华大学 | 快速信号子空间估计方法 |
CN102255839A (zh) * | 2010-05-18 | 2011-11-23 | 中兴通讯股份有限公司 | 一种适用于cmmb系统的信道估计方法和装置 |
CN103338165A (zh) * | 2013-05-24 | 2013-10-02 | 华为技术有限公司 | 一种信道估计方法和装置 |
CN103428119A (zh) * | 2012-05-24 | 2013-12-04 | 中兴通讯股份有限公司 | 一种消除信道估计奇异值影响的方法和基站 |
CN103916340A (zh) * | 2014-04-04 | 2014-07-09 | 电信科学技术研究院 | 一种噪声功率估计方法及网络侧设备 |
CN112910803A (zh) * | 2019-12-03 | 2021-06-04 | 深圳市中兴微电子技术有限公司 | 信道估计方法、装置、设备及存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7532864B2 (en) * | 2003-02-17 | 2009-05-12 | Panasonic Corporation | Noise power estimation method and noise power estimation device |
CN1288854C (zh) * | 2003-12-01 | 2006-12-06 | 北京天碁科技有限公司 | 一种用于无线通信系统的去噪方法 |
CN100349383C (zh) * | 2004-04-14 | 2007-11-14 | 华为技术有限公司 | 一种信道估计的方法和装置 |
CN101127745B (zh) * | 2006-08-16 | 2011-09-14 | 大唐移动通信设备有限公司 | 一种信道估计方法及装置 |
-
2009
- 2009-05-20 CN CN 200910084280 patent/CN101626354B/zh active Active
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102208930A (zh) * | 2010-03-30 | 2011-10-05 | 清华大学 | 快速信号子空间估计方法 |
CN102208930B (zh) * | 2010-03-30 | 2014-01-22 | 清华大学 | 快速信号子空间估计方法 |
CN102255839A (zh) * | 2010-05-18 | 2011-11-23 | 中兴通讯股份有限公司 | 一种适用于cmmb系统的信道估计方法和装置 |
CN102255839B (zh) * | 2010-05-18 | 2014-02-19 | 中兴通讯股份有限公司南京分公司 | 一种适用于cmmb系统的信道估计方法和装置 |
CN102075465A (zh) * | 2011-02-14 | 2011-05-25 | 电信科学技术研究院 | 一种噪声功率的估计方法及设备 |
CN103428119B (zh) * | 2012-05-24 | 2018-02-02 | 中兴通讯股份有限公司 | 一种消除信道估计奇异值影响的方法和基站 |
CN103428119A (zh) * | 2012-05-24 | 2013-12-04 | 中兴通讯股份有限公司 | 一种消除信道估计奇异值影响的方法和基站 |
CN103338165A (zh) * | 2013-05-24 | 2013-10-02 | 华为技术有限公司 | 一种信道估计方法和装置 |
CN103338165B (zh) * | 2013-05-24 | 2017-04-26 | 华为技术有限公司 | 一种信道估计方法和装置 |
CN103916340B (zh) * | 2014-04-04 | 2017-03-15 | 电信科学技术研究院 | 一种噪声功率估计方法及网络侧设备 |
CN103916340A (zh) * | 2014-04-04 | 2014-07-09 | 电信科学技术研究院 | 一种噪声功率估计方法及网络侧设备 |
CN112910803A (zh) * | 2019-12-03 | 2021-06-04 | 深圳市中兴微电子技术有限公司 | 信道估计方法、装置、设备及存储介质 |
WO2021109639A1 (zh) * | 2019-12-03 | 2021-06-10 | 深圳市中兴微电子技术有限公司 | 信道估计方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN101626354B (zh) | 2013-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101626354B (zh) | 多径信道的噪声估计方法和装置 | |
US8107517B2 (en) | Average-tap energy based thresholding for channel estimation in multi antenna systems | |
EP1414164A1 (en) | Method for cross-talk cancellation | |
US10153921B2 (en) | Self-interference channel estimation method and device | |
US9071480B2 (en) | Unified interference power estimation | |
JP2008507218A (ja) | Ofdmaのための専用パイロットトーンを用いた反復的なチャネルおよび干渉推定 | |
CN102647387B (zh) | 同频干扰的消除方法及装置 | |
US20160142228A1 (en) | Channel measurement method for large-scale antennas, and user terminal | |
CN101268667B (zh) | 估计发送数据向量的方法、码元确定设备和系统 | |
CN102111363B (zh) | 一种信道估计的方法和装置 | |
CN103346983A (zh) | 一种基于梳状导频的ofdm自适应复插值信道估计方法 | |
Sethuraman et al. | Low-SNR capacity of noncoherent fading channels | |
CN102594737A (zh) | 一种邻区干扰检测方法及系统 | |
CN102035768A (zh) | 一种载波干扰噪声比的测量方法及装置 | |
CN102238114A (zh) | 有效信噪比确定方法和装置 | |
CN103428126B (zh) | 一种基于梳状导频的ofdm自适应信道估计方法 | |
US8780956B2 (en) | Method and device for determining smooth window length in channel estimation | |
TW202402016A (zh) | 無線通訊的窄頻干擾管理方法及裝置 | |
CN103138820A (zh) | Vamos模式下的信号检测方法及装置 | |
CN102163990A (zh) | 多用户干扰消除方法及基站 | |
CN101356758B (zh) | Ofdm接收装置和使用其的ofdm接收设备 | |
Koch et al. | On noncoherent fading relay channels at high signal-to-noise ratio | |
CN104468426A (zh) | Lte上行信道估计方法及系统 | |
US20080298480A1 (en) | Asynchronous MIMO-OFDM spatial covariance estimation | |
CN101777885A (zh) | 滤波器系数的确定方法和装置 |
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 | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: 100089 No. A1, A2 and D, 8th floor, building 2, yard 1, Wangzhuang Road, Haidian District, Beijing Patentee after: Gaotuoxunda (Beijing) Microelectronics Co.,Ltd. Address before: 100084 science building, Tsinghua Science Park, Beijing, Haidian District, B702 Patentee before: ALTO BEAM (CHINA) Inc. |