CN115314092B - Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method - Google Patents

Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method Download PDF

Info

Publication number
CN115314092B
CN115314092B CN202210965546.5A CN202210965546A CN115314092B CN 115314092 B CN115314092 B CN 115314092B CN 202210965546 A CN202210965546 A CN 202210965546A CN 115314092 B CN115314092 B CN 115314092B
Authority
CN
China
Prior art keywords
led
ogsm
photodetector
illuminance
positioning area
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.)
Active
Application number
CN202210965546.5A
Other languages
Chinese (zh)
Other versions
CN115314092A (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.)
Xian Technological University
Original Assignee
Xian Technological 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 Xian Technological University filed Critical Xian Technological University
Priority to CN202210965546.5A priority Critical patent/CN115314092B/en
Publication of CN115314092A publication Critical patent/CN115314092A/en
Application granted granted Critical
Publication of CN115314092B publication Critical patent/CN115314092B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0602Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using antenna switching
    • H04B7/0608Antenna selection according to transmission parameters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B10/00Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/11Arrangements specific to free-space transmission, i.e. transmission through air or vacuum
    • H04B10/114Indoor or close-range type systems
    • H04B10/116Visible light communication
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Communication System (AREA)

Abstract

The invention discloses a Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method, which specifically comprises the following steps: building a visible light indoor communication link system model; establishing a mathematical model of a channel comprising a direct line-of-sight link LOS and a first order reflection link NLOS; determining the mounting interval of LEDs and dividing positioning areas; collecting the received illuminance RSS value of each LED and the actual coordinates of the reference point at the positioning point, and storing the RSS value and the actual coordinates into a fingerprint library; determining a fingerprint library of the positioning areas, and calculating a Pearson coefficient between any LED combination and the photoelectric detector in any one of the positioning areas according to the fingerprint library of the positioning areas; and taking the corresponding LED combination with high correlation as an actually activated antenna combination. The invention does not depend on the system channel characteristics, realizes the multiplexing of the time domain and the space domain, improves the error code performance, ensures the reliability and simultaneously improves the application range of OGSM-MIMO systems.

Description

一种基于Pearson系数的多阵列可见光OGSM-MIMO天线选择 方法A multi-array visible light OGSM-MIMO antenna selection method based on Pearson coefficient

技术领域Technical Field

本发明属于可见光室内通信技术领域,涉及一种基于Pearson系数的多阵列可见光OGSM-MIMO天线选择方法。The invention belongs to the technical field of visible light indoor communication and relates to a multi-array visible light OGSM-MIMO antenna selection method based on Pearson coefficient.

背景技术Background technique

在可见光室内照明场景下,为了同时兼顾照明与通信,通常采用多阵列光源布局模型,尤其在大场景下,LED阵列数量会急剧增加,从而使传统的OGSM-MIMO系统可靠性得不到保障,基于此,保证可靠性的同时进一步提升OGSM-MIMO系统的适用范围,研究者们引入了天线选择方法。In visible light indoor lighting scenarios, in order to take into account both lighting and communication, a multi-array light source layout model is usually adopted. Especially in large scenes, the number of LED arrays will increase sharply, making the reliability of the traditional OGSM-MIMO system unable to be guaranteed. Based on this, in order to ensure reliability while further improving the applicability of the OGSM-MIMO system, the researchers introduced an antenna selection method.

目前传统的天线选择方法主要包括随机天线选择方法、基于范数的天线选择方法。随机选择方法为评价天线选择方法系统性能的下限;基于范数选择的天线选择方法是基于已知信道,然而实际由于外界环境干扰,信道是时变信道。At present, traditional antenna selection methods mainly include random antenna selection method and norm-based antenna selection method. The random selection method is the lower limit of evaluating the system performance of the antenna selection method; the norm-based antenna selection method is based on a known channel, but in reality, due to external environmental interference, the channel is a time-varying channel.

发明内容Summary of the invention

为了解决上述问题,本发明提供一种基于Pearson系数的多阵列可见光OGSM-MIMO天线选择方法,不依赖系统信道特性,实现时间域与空间域的复用,改善了误码性能,解决了现有技术中存在的问题。In order to solve the above problems, the present invention provides a multi-array visible light OGSM-MIMO antenna selection method based on Pearson coefficient, which does not rely on system channel characteristics, realizes time domain and space domain multiplexing, improves bit error performance, and solves the problems existing in the prior art.

本发明所采用的技术方案是,一种基于Pearson系数的多阵列可见光OGSM-MIMO天线选择方法,具体按照以下步骤进行:The technical solution adopted by the present invention is a multi-array visible light OGSM-MIMO antenna selection method based on the Pearson coefficient, which is specifically performed in the following steps:

步骤1、确定室内空间参数,搭建可见光室内通信链路系统模型,在室内空间建立三维坐标系;Step 1: Determine the indoor space parameters, build a visible light indoor communication link system model, and establish a three-dimensional coordinate system in the indoor space;

步骤2、建立包含直射视距链路LOS与一阶反射链路NLOS的信道数学模型;Step 2, establish a channel mathematical model including the direct line-of-sight link LOS and the first-order reflection link NLOS;

步骤3、基于包含LOS与NLOS的信道数学模型,通过接收平面的光照度标准差以及接收平面光照度RSS最小值两个条件确定LED的安装间距,进而进行定位区域的划分;Step 3: Based on the channel mathematical model including LOS and NLOS, the installation spacing of the LEDs is determined by the two conditions of the standard deviation of the illuminance on the receiving plane and the minimum RSS value of the illuminance on the receiving plane, and then the positioning area is divided;

步骤4、在定位区域中选取若干定位点,在定位点采集接收到的各LED的光照度RSS值和该参考点的实际坐标,并存入指纹库;Step 4: Select several positioning points in the positioning area, collect the received RSS value of the illuminance of each LED and the actual coordinates of the reference point at the positioning point, and store them in the fingerprint library;

步骤5、确定定位区域的指纹库;Step 5: Determine the fingerprint library of the positioning area;

步骤6、根据定位区域的指纹库计算任意一个定位区域内任意LED组合与光电检测器间的Pearson系数;Step 6: Calculate the Pearson coefficient between any LED combination and the photodetector in any positioning area according to the fingerprint library of the positioning area;

步骤7、基于任意LED组合与光电检测器间的Pearson系数,将相关性高的对应LED组合作为实际激活的天线组合。Step 7: Based on the Pearson coefficient between any LED combination and the photodetector, the corresponding LED combination with high correlation is used as the actually activated antenna combination.

进一步的,所述步骤2中,任意一个光电检测器接收到的任意M个LED光源的光照度E为:Furthermore, in step 2, the illuminance E of any M LED light sources received by any photodetector is:

式中,D表示任意一个光电检测器接收端与LED光源之间的距离,D′表示任意一个光电检测器接收端与LED光源产生的次级光源之间的距离;Ej为多光源下第j个光电检测器的直射水平照度,E'j为多次级光源下第j个光电检测器的反射光照度,I(φi)为单光源到任意点的光强,φi表示第i个LED的发射角,表示第j个光电检测器的入射角,Ii'为单光源到任意点的反射光强,β表示次级光源的发射角。Where D represents the distance between the receiving end of any photodetector and the LED light source, D′ represents the distance between the receiving end of any photodetector and the secondary light source generated by the LED light source; Ej is the direct horizontal illuminance of the jth photodetector under multiple light sources, E′j is the reflected light illuminance of the jth photodetector under multiple secondary light sources, I( φi ) is the light intensity from a single light source to any point, φi represents the emission angle of the i-th LED, represents the incident angle of the jth photodetector, I i ' is the reflected light intensity from a single light source to any point, and β represents the emission angle of the secondary light source.

进一步的,所述步骤3中,人眼接收光照度舒适范围光照度RSS是300lx~1500lx。Furthermore, in step 3, the comfort range illuminance RSS of the illuminance received by the human eye is 300lx to 1500lx.

进一步的,所述步骤4中,指纹库为:Furthermore, in step 4, the fingerprint library is:

fn={Sn,(xn,yn,zn),tn|n=1,2,…,N}f n ={S n ,(x n , yn ,z n ),t n |n=1,2,…,N}

式中,表示第n个光电检测器PDn处来自各LEDm的接收光照度;(xn,yn,zn)表示第n个光电检测器的坐标位置,单位:m;tn表示第n个光电检测器接收的光照次数;In the formula, represents the received light illuminance from each LED m at the nth photodetector PD n ; (x n , y n , z n ) represents the coordinate position of the nth photodetector, in m; t n represents the number of light exposures received by the nth photodetector;

所述LEDm表示第m个LED,m=1,2,…,M;LEDm表示为:LED m represents the mth LED, m=1, 2, ..., M; LED m is represented by:

LEDm={Rm,(xm,ym,zm)|m=1,2,…,M}LED m = {R m ,(x m ,y m ,z m )|m = 1,2,…,M}

其中,Rm表示第m个LED的中心发光强度,单位:lx;(xm,ym,zm)为第m个LED的坐标位置,单位:m。Wherein, R m represents the central luminous intensity of the mth LED, unit: lx; (x m , y m , z m ) is the coordinate position of the mth LED, unit: m.

进一步的,所述步骤5中,定位区域的指纹库fn,k根据步骤4中的指纹库fn确定,即任意一个定位区域内第n个光电检测器接收到LED组合的平均光照度的集合,见下式:Furthermore, in step 5, the fingerprint library fn ,k of the positioning area is determined according to the fingerprint library fn in step 4, that is, the set of average illuminances of the LED combination received by the nth photodetector in any positioning area, as shown in the following formula:

fn,k={Sn,k,(xn,k,yn,k,zn,k),tn,k|n=1,2,…,Nk}f n,k = {S n,k , (x n,k , yn ,k , z n,k ), t n,k | n = 1, 2, …, N k }

式中,Sn,k表示第k个定位区域内第n个光电检测器接收到来自第k个定位区域内第m个LED的接收光照度,(xn,k,yn,k,zn,k)表示第k个区域第n个光电检测器PD的坐标位置,单位:m;Where, Sn ,k represents the received light illumination received by the nth photodetector in the kth positioning area from the mth LED in the kth positioning area, (xn ,k , yn,k , zn,k ) represents the coordinate position of the nth photodetector PD in the kth area, unit: m;

第k个定位区域内第m个LED的位置LEDm,k表示为:The position LED m,k of the mth LED in the kth positioning area is expressed as:

LEDm,k={Rm,k,(xm,k,ym,k,zm,k)|m=1,2,…,Mk}LED m,k = {R m,k , (x m,k , y m,k , z m,k ) | m = 1, 2, …, M k }

始终,Rm,k表示第k个区域的第m个LED的中心发光强度,单位:lx;(xm,k,ym,k,zm,k)表示第k个区域内第m个LED的坐标位置,单位:m。Throughout, R m,k represents the central luminous intensity of the mth LED in the kth region, unit: lx; (x m,k , y m,k , z m,k ) represents the coordinate position of the mth LED in the kth region, unit: m.

进一步的,所述步骤6中,第k个定位区域内任意LED组合与光电检测器间的Pearson系数ρ(LEDm,k,fn,k)表示为:Furthermore, in step 6, the Pearson coefficient ρ(LED m,k ,f n,k ) between any LED combination and the photodetector in the kth positioning area is expressed as:

式中,是第k个定位区域指纹库的第n个光电检测器接收实际LED组合的平均光照度;σ表示样本标准差,/>表示第k个定位区域第n个光电检测器接收到第m个LED的光照度。In the formula, is the average illuminance of the actual LED combination received by the nth photodetector of the kth localization area fingerprint library; σ represents the sample standard deviation, /> It means that the nth photodetector in the kth positioning area receives the light intensity of the mth LED.

本发明的有益效果是:The beneficial effects of the present invention are:

本发明实施例结合了可见光通信链路模型、LED信道漫反射模型、可见光OGSM-MIMO系统与Pearson系数的相关性数学模型,拟合室内信道环境,并采用基于Pearson系数的天线选择方法对OGSM-MIMO系统进行天线选择,利用不同位置光电检测器(Photoelectric Detector,PD)端与激活的发射天线之间的Pearson系数的相关性来选择最优天线组合,不依赖系统信道特性,从而实现时间域与空间域的复用,提升了系统的误码性能。The embodiment of the present invention combines the visible light communication link model, the LED channel diffuse reflection model, and the mathematical model of the correlation between the visible light OGSM-MIMO system and the Pearson coefficient to fit the indoor channel environment, and adopts an antenna selection method based on the Pearson coefficient to select antennas for the OGSM-MIMO system. The correlation of the Pearson coefficients between the photoelectric detector (PD) ends at different positions and the activated transmitting antennas is used to select the optimal antenna combination, which does not rely on the system channel characteristics, thereby realizing multiplexing of the time domain and the space domain, and improving the system's bit error performance.

本发明实施例还可以应用于大规模OGSM-MIMO系统中,减少了计算量,不局限于室内光通信,激光等光领域也可以应用。The embodiments of the present invention can also be applied to large-scale OGSM-MIMO systems, reducing the amount of calculation, and are not limited to indoor optical communications, but can also be applied to optical fields such as lasers.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

图1是本发明实施例中可见光通信链路模型。FIG. 1 is a visible light communication link model in an embodiment of the present invention.

图2是本发明实施例中LED信道漫反射模型。FIG. 2 is a diffuse reflection model of an LED channel in an embodiment of the present invention.

图3是本发明实施例基于天线选择的OGSM-MIMO系统模型。FIG. 3 is an OGSM-MIMO system model based on antenna selection according to an embodiment of the present invention.

图4是不同天线选择方法的OGSM-MIMO系统的误码率仿真图。FIG4 is a simulation diagram of the bit error rate of the OGSM-MIMO system with different antenna selection methods.

图5是本发明实施例基于Pearson系数的OGSM-MIMO系统的理论误码性能仿真图。FIG5 is a simulation diagram of theoretical bit error performance of an OGSM-MIMO system based on the Pearson coefficient according to an embodiment of the present invention.

图6是本发明实施例不同条件下基于Pearson系数的OGSM-MIMO系统误码率分析。FIG6 is an analysis of the bit error rate of the OGSM-MIMO system based on the Pearson coefficient under different conditions in an embodiment of the present invention.

具体实施方式Detailed ways

下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will be combined with the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

实施例Example

一种基于Pearson系数的多阵列可见光OGSM-MIMO天线选择方法,包括如下步骤:A multi-array visible light OGSM-MIMO antenna selection method based on Pearson coefficient comprises the following steps:

步骤1、确定室内空间参数,如图1所示,以4m×4m×3m室内房间为实验仿真模型,搭建可见光室内通信链路系统模型,以室内空间的一个墙角作为坐标原点建立三维坐标系。Step 1, determine the indoor space parameters, as shown in Figure 1, take a 4m×4m×3m indoor room as the experimental simulation model, build a visible light indoor communication link system model, and establish a three-dimensional coordinate system with a corner of the indoor space as the coordinate origin.

步骤2、建立包含直射视距链路LOS与一阶反射链路NLOS的信道数学模型,如图2所示;其中,单光源到任意点的光强I(φi)定义为:Step 2: Establish a channel mathematical model including the direct line-of-sight link LOS and the first-order reflection link NLOS, as shown in FIG2 ; wherein the light intensity I(φ i ) from a single light source to any point is defined as:

I(φi)=I0cosrφi (1)I(φ i )=I 0 cos r φ i (1)

式中,I0表示LED的中心发光强度,φi表示第i个LED的发射角,r为光源的辐射模式,r的取值范围是0~2,数值表示反射次数,0表示不考虑反射,2表示反射2次。Where I0 represents the central luminous intensity of the LED, φi represents the emission angle of the i-th LED, r is the radiation pattern of the light source, and the value range of r is 0 to 2. The value represents the number of reflections, 0 means no reflection is considered, and 2 means 2 reflections.

则多光源下第j个PD的直射水平照度Ej为:Then the direct horizontal illumination Ej of the jth PD under multiple light sources is:

式中,表示第j个PD的入射角,Eij表示第j个PD端接收到第i个LED的光照度,(xj,yj,0)表示第j个PD的坐标,(Xi,Yi,Zi)表示第i个LED的坐标。图2中Dj表示第j个PD的坐标,W表示图1中任意一个区域空间宽度,H=3表示室内高度为3m。In the formula, represents the incident angle of the jth PD, Eij represents the illuminance of the i-th LED received by the j-th PD, ( xj , yj , 0) represents the coordinates of the j-th PD, and ( Xi , Yi , Zi ) represents the coordinates of the i-th LED. In Figure 2, Dj represents the coordinates of the j-th PD, W represents the width of any area in Figure 1, and H=3 represents the indoor height of 3m.

假设任意LED灯源产生的次级光源为P点,P点坐标可以表示为P(X1,0,Z1)。同理,单光源到任意点的反射光强I′i定义为:Assume that the secondary light source generated by any LED light source is point P, and the coordinates of point P can be expressed as P(X 1 ,0,Z 1 ). Similarly, the reflected light intensity I′ i from a single light source to any point is defined as:

I′i=kIi cosm-1β (3)I′ i = kI i cos m-1 β (3)

式中,Ii表示第i个LED直射到P点的光照强度,k为墙面的反射系数,β表示次级光源的发射角,m为光源的辐射模式。则多次级光源下第j个PD的反射光照度E'j为:In the formula, Ii represents the light intensity of the i-th LED directly irradiating to point P, k is the reflection coefficient of the wall, β represents the emission angle of the secondary light source, and m is the radiation mode of the light source. Then the reflected light illuminance E'j of the j-th PD under multiple secondary light sources is:

E′ij表示第j个PD接收到第i个LED的光照度。 E′ij represents the illuminance of the i-th LED received by the j-th PD.

则任意一个PD端接收到的任意M个LED光源的光照度E为:Then the illuminance E of any M LED light sources received by any PD end is:

式中,D表示任意一个PD端与LED光源之间的距离,D′表示任意一个PD端与LED光源产生的次级光源之间的距离。Wherein, D represents the distance between any PD end and the LED light source, and D′ represents the distance between any PD end and the secondary light source generated by the LED light source.

实施例中,即任意一个PD端接收到的任意M个LED光源的光照度E为:In the embodiment, That is, the illuminance E of any M LED light sources received by any PD end is:

通过公式(5)使得PD接收到不同LED的光照度差异更加明显,Pearson系数相关性区别更大,原因是现实生活墙面也会反射光,可以更好的模拟现实环境。Formula (5) makes the difference in the illuminance received by the PD from different LEDs more obvious, and the Pearson coefficient correlation is more different. This is because walls in real life also reflect light, which can better simulate the real environment.

步骤3、结合步骤2建立的LED信道漫反射模型(包括LOS与NLOS的信道数学模型)确定各LED的安装间距以及根据接收底面RSS值(光照度)进行定位区域的划分。Step 3: Determine the installation spacing of each LED based on the LED channel diffuse reflection model (including the channel mathematical model of LOS and NLOS) established in step 2 and divide the positioning area according to the RSS value (illuminance) of the receiving bottom surface.

LED的安装间距不同导致RSS值不同,最后通过接收平面的光照度标准差以及接收平面RSS最小值两个条件确定了LED的安装间距,进而确定定位区域的划分。国际照明标准规定人眼接收光照度舒适范围RSS是300lx~1500lx,根据不同区域的LED产生的光照度,选择最优的光照度区域,模拟的环境更加贴近于现实生活,考虑因素更为全面。Different installation spacing of LEDs leads to different RSS values. Finally, the installation spacing of LEDs is determined by two conditions: the standard deviation of the illumination on the receiving plane and the minimum RSS value on the receiving plane, and then the division of the positioning area is determined. The international lighting standard stipulates that the comfortable range of illumination received by the human eye is 300lx to 1500lx. According to the illumination generated by LEDs in different areas, the optimal illumination area is selected. The simulated environment is closer to real life and the factors considered are more comprehensive.

步骤4、在定位区域中选取若干定位点,在定位点采集接收到的各LED的RSS值和该参考点的实际坐标,并存入指纹库;其中,包括M个LED,第m个LED记为LEDm,m=1,2,…,M;可表示为:Step 4: Select several positioning points in the positioning area, collect the received RSS values of each LED and the actual coordinates of the reference point at the positioning points, and store them in the fingerprint library; wherein, there are M LEDs, the mth LED is recorded as LED m , m = 1, 2, ..., M; which can be expressed as:

LEDm={Rm,(xm,ym,zm)|m=1,2,…,M} (6)LED m = {R m ,(x m ,y m ,z m )|m = 1,2,…,M} (6)

其中,Rm表示第m个LED的中心发光强度,单位:lx;(xm,ym,zm)为第m个LED的坐标位置,单位:m。定位区域有N个PD,第n个PD记为PDn,n=1,2…,N,则在PDn处接收到来自LEDm信号的光照强度可记为因此,指纹库可定义为:Where R m represents the central luminous intensity of the mth LED, unit: lx; (x m ,y m ,z m ) is the coordinate position of the mth LED, unit: m. There are N PDs in the positioning area, and the nth PD is recorded as PD n , n = 1, 2…, N, then the light intensity received from the LED m signal at PD n can be recorded as Therefore, the fingerprint library can be defined as:

fn={Sn,(xn,yn,zn),tn|n=1,2,…,N} (7)f n = {S n , (x n , yn , z n ), t n | n = 1, 2, …, N} (7)

式中,fn表示不同PD接收端接收到不同LED光照度的集合,表示PDn处来自各LED的接收光照度;(xn,yn,zn)表示该PD的坐标位置,单位:m;tn表示PD接收的光照次数。In the formula, fn represents the set of different LED light intensities received by different PD receiving ends. represents the received light illuminance from each LED at PD n ; ( xn , yn , zn ) represents the coordinate position of the PD, unit: m; tn represents the number of times the PD receives light.

步骤5、确定定位区域的指纹库。根据步骤3划分的LED光源布局模型,选择划分的任意一个区域;假设室内底面划分为k个定位区域,记为fn,k,k∈{1,2,…,K},则定位区域的指纹库定义为:Step 5: Determine the fingerprint library of the positioning area. According to the LED light source layout model divided in step 3, select any area divided; assuming that the indoor bottom surface is divided into k positioning areas, denoted as f n,k , k∈{1,2,…,K}, then the fingerprint library of the positioning area is defined as:

fn,k={Sn,k,(xn,k,yn,k,zn,k),tn,k|nk=1,2,…,Nk} (8)f n,k = {S n,k ,(x n,k ,yn ,k ,z n,k ),t n,k | nk = 1,2,…, Nk } (8)

fn,k表示第k个定位区域中第n个PD接收到LED组合的平均光照度的集合;Sn,k表示第k个定位区域内PDn接收到来自LEDm,k的接收光照度,(xn,k,yn,k,zn,k)表示第k个区域第n个光电检测器PD的坐标位置,单位:m;tn,k表示第k个定位区域中第n个光电检测器PD接收的光照次数。f n,k represents the set of average light illumination received by the nth PD in the kth positioning area from the LED combination; Sn ,k represents the received light illumination received by PD n in the kth positioning area from LED m,k , ( xn,k , yn,k , zn ,k ) represents the coordinate position of the nth photodetector PD in the kth area, unit: m; tn,k represents the number of light illuminations received by the nth photodetector PD in the kth positioning area.

第k个定位区域内第m个LED的位置LEDm,k表示为:The position LED m,k of the mth LED in the kth positioning area is expressed as:

LEDm,k={Rm,k,(xm,k,ym,k,zm,k)|m=1,2,…,Mk} (9)LED m,k = {R m,k ,(x m,k ,y m,k ,z m,k )|m = 1,2,…,M k } (9)

式中,Rm,k表示第k个区域的第m个LED的中心发光强度,单位:lx;(xm,k,ym,k,zm,k)表示第k个区域PD的坐标位置,单位:m。Wherein, R m,k represents the central luminous intensity of the mth LED in the kth region, unit: lx; (x m,k ,y m,k ,z m,k ) represents the coordinate position of the kth region PD, unit: m.

步骤6、选择步骤5划分的任意一个区域,根据步骤4的指纹库计算区域内任意LED组合与光电检测器间的Pearson系数。Step 6: Select any area divided in step 5, and calculate the Pearson coefficient between any LED combination and the photodetector in the area according to the fingerprint library in step 4.

Pearson相关系数是一种线性相关系数,定义为两点的协方差与标准差乘积的比值,Pearson系数的值介于-1到1之间,值越大,相关度越高,因此可以作为天线组合的选择标准,对于第k个分区内任意LED组合与指纹库PD间的相关系数ρ(LEDm,k,fn,k)表示为:The Pearson correlation coefficient is a linear correlation coefficient, which is defined as the ratio of the covariance of two points to the product of the standard deviation. The value of the Pearson coefficient is between -1 and 1. The larger the value, the higher the correlation. Therefore, it can be used as a selection criterion for antenna combinations. The correlation coefficient ρ(LED m,k ,f n,k ) between any LED combination in the kth partition and the fingerprint library PD is expressed as:

式中,是第k个分区指纹库的第n个PD的接收实际LED组合的平均光照度;将相似度高的m个LED作为发射天线组合;σ表示样本标准差,是Pearson系数的原始公式的参数;表示第k个区域第n个PD接收到第m个LED的光照度。In the formula, is the average illuminance of the actual LED combination received by the nth PD in the kth partition fingerprint library; m LEDs with high similarity are used as the transmitting antenna combination; σ represents the sample standard deviation, which is the parameter of the original formula of the Pearson coefficient; It means that the nth PD in the kth area receives the illuminance of the mth LED.

根据VLC通信链路模型假设第k个分区指纹库中有6个LED,其中心发光强度为21(lx),则LED信息可以根据式(9)表示为:According to the VLC communication link model, it is assumed that there are 6 LEDs in the kth partition fingerprint library, and the central luminous intensity is 21 (lx). Then the LED information can be expressed according to formula (9):

通过可见光通信模型可得各PD端接收光照度,假设PD端有2个,坐标信息表示为(0.8,0.8,0),(0.8,1.6,0),则接收端分别接收6个LED可得接收光照度有12种,分区指纹库根据公式(8)表示为:Through the visible light communication model, the received light intensity of each PD can be obtained. Assuming there are two PDs, the coordinate information is expressed as (0.8, 0.8, 0) and (0.8, 1.6, 0). Then, the receiving end receives 6 LEDs respectively, and there are 12 types of received light intensity. The partition fingerprint library is expressed according to formula (8):

式中,考虑到现实生活中噪声干扰,因此测量3次取平均,tn=3时平均光照值,例如由于已知6个LED,2个PD接收端,在OGSM-MIMO系统中,即发射天线Nt=6,接收天线Nr=2,假设激活的发射天线(LED)Na=2,则OGSM-MIMO系统中总共有15种天线组合,然而OGSM-MIMO实际天线组合必须满足2的整幂次方,因此实际天线组合便是从15种任选8种,根据公式(10)、(11)、(12)可得Pearson系数集合:In the formula, considering the noise interference in real life, the average of three measurements is taken, and the average illumination value when t n = 3 is, for example Since there are 6 LEDs and 2 PD receiving ends, in the OGSM-MIMO system, that is, the transmitting antenna Nt = 6, the receiving antenna Nr = 2, and assuming that the activated transmitting antenna (LED) Na = 2, there are a total of 15 antenna combinations in the OGSM-MIMO system. However, the actual antenna combination of OGSM-MIMO must satisfy the integral power of 2, so the actual antenna combination is to select 8 of the 15. According to formulas (10), (11), and (12), the Pearson coefficient set can be obtained:

步骤7、通过步骤6计算的Pearson系数相关性选择其对应的天线组合作为实际激活的天线组合,从而获得最优通信效果;天线组合的数量受到2的整幂次约束,原因是经过约束条件限制,每种天线组合都会对应相同数目的调制信号组合;本发明实施例由激活所有天线改变为激活部分天线,只需要考虑发射天线端最优的几根天线,不仅降低多径干扰对通信质量的影响,也减少了天线数量对误码性能的影响,减少多径干扰、提升误码性能,克服了天线数与信道间干扰的关系制约对通信质量的影响,同时计算量小,对于天线数更多的大规模OGSM-MIMO系统计算量仍然很小,适用性更强。Step 7, select the corresponding antenna combination as the actually activated antenna combination through the Pearson coefficient correlation calculated in step 6, so as to obtain the optimal communication effect; the number of antenna combinations is constrained by an integer power of 2, because after the constraint condition, each antenna combination will correspond to the same number of modulation signal combinations; the embodiment of the present invention changes from activating all antennas to activating some antennas, and only needs to consider the best antennas at the transmitting antenna end, which not only reduces the impact of multipath interference on communication quality, but also reduces the impact of the number of antennas on bit error performance, reduces multipath interference, improves bit error performance, overcomes the impact of the relationship between the number of antennas and inter-channel interference on communication quality, and at the same time, the amount of calculation is small, and the amount of calculation for large-scale OGSM-MIMO systems with more antennas is still very small, and the applicability is stronger.

本发明实施例选择的天线组合应用于OGSM-MIMO系统中,如图3所示,已知实际天线组合,同时采用L-PAM调制方式,调制阶数L=2,则OGSM-MIMO系统映射表如表1所示。The antenna combination selected in the embodiment of the present invention is applied to the OGSM-MIMO system, as shown in FIG3 . The actual antenna combination is known, and the L-PAM modulation mode is adopted, and the modulation order L=2, then the OGSM-MIMO system mapping table is shown in Table 1.

表1基于天线选择的OGSM-MIMO映射Table 1 OGSM-MIMO mapping based on antenna selection

图3中,xu表示传输的信号,如表1中的输入比特组成的数据流,x1,u表示天线组合映射比特数,x2,u星座映射比特数,x1,u、x2,u都可以通过OGSM-MIMO系统理论公式得到,yu表示接收端接收到的信号。Nt表示发射端有Nt个发射天线,Nr表示接收端有Nr个接收天线。ML表示最大似然检测算法,s/p、p/s分别表示串并转换和并串转换,OGSM是光广义空间调制(Optical Generalized Spatial Modulation),Symbol modulation表示符号调制,Symboldemodulation表示符号解调,Antenna modulation表示天线调制,Antenna demodulation表示天线解调。In Figure 3, xu represents the transmitted signal, such as the data stream composed of input bits in Table 1, x1,u represents the number of bits mapped by the antenna combination, x2 ,u represents the number of bits mapped by the constellation, and x1 ,u and x2 ,u can be obtained through the theoretical formula of the OGSM-MIMO system. yu represents the signal received by the receiving end. Nt represents that the transmitting end has Nt transmitting antennas, and Nr represents that the receiving end has Nr receiving antennas. ML represents the maximum likelihood detection algorithm, s/p and p/s represent serial-to-parallel conversion and parallel-to-serial conversion respectively, OGSM is optical generalized spatial modulation, Symbol modulation represents symbol modulation, Symboldemodulation represents symbol demodulation, Antenna modulation represents antenna modulation, and Antenna demodulation represents antenna demodulation.

分析基于Pearson系数的多阵列可见光OGSM-MIMO天线选择方法相比其他传统天线选择方法的性能。The performance of the multi-array visible light OGSM-MIMO antenna selection method based on the Pearson coefficient is analyzed compared with other traditional antenna selection methods.

如附图4所示,发射端激活天线数Na=2,接收PD端Nr=2,调制方式为2PAM,从仿真结果可以看出,相比基于范数的天线选择方法(Norm selection)以及随机选择方法(Random selection),基于本发明Pearson系数天线选择方法的OGSM-MIMO系统误码率性能更优;当误码率(BER)达到10-3时,采用基于Pearson系数的天线选择方法,OGSM-MIMO系统误码性能(37.8dB)分别比Norm selection(40dB)、Random selection(41.8dB)改善了2.2dB、4dB。As shown in FIG4 , the number of activated antennas at the transmitting end is Na = 2, the number of receiving PD end is Nr = 2, and the modulation mode is 2PAM. It can be seen from the simulation results that, compared with the antenna selection method based on the norm (Norm selection) and the random selection method (Random selection), the OGSM-MIMO system bit error rate performance based on the Pearson coefficient antenna selection method of the present invention is better; when the bit error rate (BER) reaches 10 -3 , the OGSM-MIMO system bit error performance (37.8dB) is improved by 2.2dB and 4dB compared with the Norm selection (40dB) and the Random selection (41.8dB), respectively.

如附图5所示,调制方式采用2PAM,发射端激活天线数Na分别为2根、3根、4根,接收天线Nr=2。从仿真结果可以看出:1)当信噪比比较低时,OGSM-MIMO系统的理论误码率高于实际误码率,当信噪比比较高时,理论误码率与实际误码率基本重合。2)相同的调制方式下,增加发射天线数,会降低误码性能;说明本发明实施例在增加激活发射天线数保证可靠性的同时进一步提升OGSM-MIMO系统的适用范围。As shown in Figure 5, the modulation mode adopts 2PAM, the number of activated antennas Na at the transmitting end is 2, 3, and 4 respectively, and the receiving antenna Nr = 2. From the simulation results, it can be seen that: 1) when the signal-to-noise ratio is relatively low, the theoretical bit error rate of the OGSM-MIMO system is higher than the actual bit error rate, and when the signal-to-noise ratio is relatively high, the theoretical bit error rate and the actual bit error rate basically coincide. 2) Under the same modulation mode, increasing the number of transmitting antennas will reduce the bit error performance; this shows that the embodiment of the present invention further improves the scope of application of the OGSM-MIMO system while increasing the number of activated transmitting antennas to ensure reliability.

如附图6所示,发射端激活天线数Na=2,接收天线Nr=6,OGSM-MIMO系统调制方式采用2PAM,4PAM、8PAM。从仿真结果以及OGSM-MIMO系统可以看出:随着调制阶数的增加,传输速率会得到提升,但误码性能会下降;当误码率达到10-3时,2PAM调制方式相比4PAM、8PAM调制方式,传输速率分别降低了2bpcu,4bpcu,但误码性能改善了约2dB、4dB。As shown in Figure 6, the number of active antennas at the transmitting end is Na = 2, the number of receiving antennas is Nr = 6, and the OGSM-MIMO system modulation uses 2PAM, 4PAM, and 8PAM. From the simulation results and the OGSM-MIMO system, it can be seen that: as the modulation order increases, the transmission rate will be improved, but the bit error performance will decrease; when the bit error rate reaches 10-3 , the 2PAM modulation mode reduces the transmission rate by 2bpcu and 4bpcu respectively compared with the 4PAM and 8PAM modulation modes, but the bit error performance is improved by about 2dB and 4dB.

以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。The above description is only a preferred embodiment of the present invention and is not intended to limit the protection scope of the present invention. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention are included in the protection scope of the present invention.

Claims (5)

1. The multi-array visible light OGSM-MIMO antenna selection method based on the Pearson coefficient is characterized by comprising the following steps of:
step 1, determining indoor space parameters, building a visible light indoor communication link system model, and building a three-dimensional coordinate system in an indoor space;
Step 2, establishing a channel mathematical model comprising a direct line-of-sight link LOS and a first-order reflection link NLOS;
Step 3, determining the mounting interval of the LEDs based on a channel mathematical model comprising LOS and NLOS through two conditions of the standard deviation of illuminance of a receiving plane and the minimum value of illuminance RSS of the receiving plane, and dividing a positioning area;
step 4, selecting a plurality of positioning points in the positioning area, collecting the received illuminance RSS value of each LED and the actual coordinates of the positioning points at the positioning points, and storing the RSS value and the actual coordinates of the positioning points into a fingerprint library;
Step 5, determining a fingerprint library of the positioning area;
Step 6, calculating Pearson coefficients between any LED combination and the photoelectric detector in any one positioning area according to the fingerprint library of the positioning area;
Step 7, based on the Pearson coefficient between any LED combination and the photoelectric detector, the corresponding LED combination with high relativity is used as an actually activated antenna combination;
In the step 6, the Pearson coefficient ρ (LED m,k,fn,k) between any LED combination and the photodetector in the kth positioning region is expressed as:
Where R m,k represents the center emission intensity of the mth LED of the kth region, The nth photodetector which is the kth positioning area fingerprint database receives the average illuminance of the actual LED combination; sigma represents the sample standard deviation,/>Indicating that the kth positioning area nth photodetector received the illumination of the mth LED.
2. The method for selecting the multi-array visible light OGSM-MIMO antenna based on the Pearson coefficients according to claim 1, wherein the illuminance E of any M LED light sources received by any one of the photodetectors in step 2 is:
wherein D represents the distance between the receiving end of any one photoelectric detector and the LED light source, and D' represents the distance between the receiving end of any one photoelectric detector and the secondary light source generated by the LED light source; e j is the direct horizontal illumination of the jth photodetector under multiple light sources, E j' is the reflected illumination of the jth photodetector under multiple light sources, I (phi i) is the light intensity from a single light source to any point, phi i represents the emission angle of the ith LED, Representing the angle of incidence of the jth photodetector, I i' is the reflected intensity of the single light source to any point, and β represents the emission angle of the secondary light source.
3. The method for selecting the multi-array visible light OGSM-MIMO antenna based on the Pearson coefficients according to claim 1, wherein in the step 3, the illuminance RSS of the human eye in the comfortable range of illuminance is 300 lx-1500 lx.
4. The method for selecting the multi-array visible light OGSM-MIMO antenna based on the Pearson coefficients according to claim 1, wherein in the step 4, the fingerprint library is:
fn={Sn,(xn,yn,zn),tn|n=1,2,L,N}
In the method, in the process of the invention, Representing the received illuminance from each LED m at the nth photodetector PD n; (x n,yn,zn) represents the coordinate position of the nth photodetector, in units: m; t n represents the number of times of illumination received by the nth photodetector;
The LED m represents an mth LED, m=1, 2, l, m; LED m is represented as:
LEDm={Rm,(xm,ym,zm)|m=1,2,…,M}
Wherein, R m represents the central luminous intensity of the mth LED, unit: lx; (x m,ym,zm) is the coordinate position of the mth LED in units: m.
5. The method for selecting a multi-array visible light OGSM-MIMO antenna based on Pearson coefficients according to claim 1, wherein in step 5, the fingerprint library f n,k of the positioning area is determined according to the fingerprint library f n in step 4, that is, the average illuminance set received by the nth photodetector in any positioning area from the LED combinations is expressed by the following formula:
fn,k={Sn,k,(xn,k,yn,k,zn,k),tn,k|n=1,2,L,Nk}
Where S n,k denotes the received illuminance from the nth LED in the kth positioning area received by the nth photodetector in the kth positioning area, (x n,k,yn,k,zn,k) denotes the coordinate position of the nth photodetector PD in the kth area in units of: m;
The position of the mth LED in the kth positioning area LED m,k is denoted as:
LEDm,k={Rm,k,(xm,k,ym,k,zm,k)|m=1,2,…,Mk}
Wherein R m,k represents the center emission intensity of the m-th LED of the k-th region, in units of: lx; (x m,k,ym,k,zm,k) represents the coordinate position of the mth LED in the kth region, in units of: m.
CN202210965546.5A 2022-08-12 2022-08-12 Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method Active CN115314092B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210965546.5A CN115314092B (en) 2022-08-12 2022-08-12 Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210965546.5A CN115314092B (en) 2022-08-12 2022-08-12 Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method

Publications (2)

Publication Number Publication Date
CN115314092A CN115314092A (en) 2022-11-08
CN115314092B true CN115314092B (en) 2024-05-10

Family

ID=83862729

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210965546.5A Active CN115314092B (en) 2022-08-12 2022-08-12 Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method

Country Status (1)

Country Link
CN (1) CN115314092B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116131948B (en) * 2023-04-17 2023-07-25 广东工业大学 Gray-GSM constellation construction method, device, terminal and medium based on visible light communication

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080114132A (en) * 2007-06-27 2008-12-31 재단법인서울대학교산학협력재단 Method and apparatus for selecting a transmitting antenna in uplink of a multiple input multiple output wireless communication system
CN103728608A (en) * 2013-12-26 2014-04-16 电子科技大学 Antenna arrangement method for improving MIMO-OTH radar detecting performance in ionized layer double-Gaussian model
CN109831253A (en) * 2019-03-19 2019-05-31 东莞信大融合创新研究院 Visible light MIMO antenna selection method based on multi-color CSK modulation
CN110289908A (en) * 2019-06-18 2019-09-27 西安工业大学 An adaptive visible light MIMO communication method considering lighting

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080114132A (en) * 2007-06-27 2008-12-31 재단법인서울대학교산학협력재단 Method and apparatus for selecting a transmitting antenna in uplink of a multiple input multiple output wireless communication system
CN103728608A (en) * 2013-12-26 2014-04-16 电子科技大学 Antenna arrangement method for improving MIMO-OTH radar detecting performance in ionized layer double-Gaussian model
CN109831253A (en) * 2019-03-19 2019-05-31 东莞信大融合创新研究院 Visible light MIMO antenna selection method based on multi-color CSK modulation
CN110289908A (en) * 2019-06-18 2019-09-27 西安工业大学 An adaptive visible light MIMO communication method considering lighting

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
MIMO系统的发射天线选择及自适应链路实现;王文杰;张莹;;系统工程与电子技术;20070815(08);全文 *
Through-the-Debris Imagery of Vital Signs Using UWB MIMO Bioradar;Zhao Li;《2019 Photonics & Electromagnetics Research Symposium - Fall (PIERS - Fall)》;20200305;全文 *
兼顾照明的可见光MIMO通信系统模型;朱彤;赵黎;刘智港;张峰;;红外与激光工程;20180825(第08期);全文 *
基于光广义空间调制的VLC-MIMO系统研究;赵黎;中国激光;20220718;全文 *

Also Published As

Publication number Publication date
CN115314092A (en) 2022-11-08

Similar Documents

Publication Publication Date Title
Zeng et al. Improvement of date rate by using equalization in an indoor visible light communication system
US10289094B2 (en) System and method for the optimization of radiance modelling and controls in predictive daylight harvesting
CN108732537B (en) An indoor visible light localization method based on neural network and received signal strength
Zhu et al. A novel 3D non-stationary channel model for 6G indoor visible light communication systems
CN105634593B (en) A kind of indoor visible light communication LED array layout optimization method based on genetic algorithm
CN115314092B (en) Pearson coefficient-based multi-array visible light OGSM-MIMO antenna selection method
CN1292955A (en) Apparatus and method for improved connectivity in wireless optical communication systems
CN111664853B (en) Anti-NLOS interference visible light localization method and system based on linear regression model
CN113949447B (en) Geometric random channel modeling method for indoor visible light communication
CN112468954A (en) Visible light indoor stereo positioning method based on neural network
CN115657620A (en) Multi-objective coordination optimization control method and system based on indoor visible light communication
Gfeller et al. Dynamic cell planning for wireless infrared in-house data transmission
CN113219408B (en) Improved RBF neural network indoor visible light positioning method and system
Ding et al. Impact of LED array simplification on indoor visible light communication channel modeling
CN116087877B (en) Point classification indoor visible light positioning method based on artificial intelligence algorithm
Sayed et al. A neural network-based vlc indoor positioning system for moving users
CN117440554B (en) Scene modeling method and system for realizing LED light source based on digitization
CN106911387B (en) Visible light communication multiple cell interference management method
Al-Ghamdi et al. Analysis of optical wireless links employing a beam clustering method and diversity receivers
Ala’F et al. On the effect of light emitting diodes positions on the performance of an indoor visible light communication system
CN119210531B (en) Microwave signal transmission method and system based on non-line-of-sight microwave equipment
US20240260162A1 (en) System and method to optimize a light emitting diode power allocation framework
Serpi et al. Machine learning assisted indoor visible light communications radio environment maps
Al-Ghamdi et al. Performance analysis of mobile optical wireless systems employing a novel beam clustering method and diversity detection
CN114786267A (en) Resource allocation optimization method, controller and communication system for multi-power terminal communication

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant