CN108594256B - 一种基于脉冲编码技术的相干激光雷达 - Google Patents
一种基于脉冲编码技术的相干激光雷达 Download PDFInfo
- Publication number
- CN108594256B CN108594256B CN201810335326.8A CN201810335326A CN108594256B CN 108594256 B CN108594256 B CN 108594256B CN 201810335326 A CN201810335326 A CN 201810335326A CN 108594256 B CN108594256 B CN 108594256B
- Authority
- CN
- China
- Prior art keywords
- pulse
- signal
- optical
- signals
- atmosphere
- 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
Links
- 230000001427 coherent effect Effects 0.000 title claims abstract description 68
- 238000005516 engineering process Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 claims abstract description 23
- 239000013307 optical fiber Substances 0.000 claims abstract description 21
- 108091026890 Coding region Proteins 0.000 claims abstract description 7
- 230000003287 optical effect Effects 0.000 claims description 51
- 238000001228 spectrum Methods 0.000 claims description 38
- 238000012545 processing Methods 0.000 claims description 17
- 239000000835 fiber Substances 0.000 claims description 9
- 238000000691 measurement method Methods 0.000 claims description 4
- 230000010355 oscillation Effects 0.000 claims description 4
- 238000012163 sequencing technique Methods 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims 1
- 230000000295 complement effect Effects 0.000 abstract description 31
- 238000001514 detection method Methods 0.000 abstract description 18
- 238000009825 accumulation Methods 0.000 abstract description 10
- 230000006835 compression Effects 0.000 abstract description 4
- 238000007906 compression Methods 0.000 abstract description 4
- 238000000605 extraction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 21
- 239000011159 matrix material Substances 0.000 description 12
- 238000005311 autocorrelation function Methods 0.000 description 9
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000012935 Averaging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/50—Systems of measurement based on relative movement of target
- G01S17/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/26—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/483—Details of pulse systems
- G01S7/484—Transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Multimedia (AREA)
- Aviation & Aerospace Engineering (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开了一种基于脉冲编码技术的相干激光雷达。本发明针对相干激光雷达回波信号信噪比低,提取困难,提出一种采用脉冲编码技术改善系统信噪比从而增大系统动态范围的方法。本发明通过电光调制器输出脉冲编码序列编码,降低了各发射波形之间的时域互相关,并利用Golay互补序列的互补性质减小自相关旁瓣,降低了波形脉冲压缩后的旁瓣。另外,本发明按照解码规则进行解码,解决了光纤激光器脉冲峰值功率受限的问题,显著地提高了相干激光雷达信噪比,有效提高了探测距离,且不需要脉冲积累时间,不牺牲探测时间和距离分辨率。
Description
技术领域
本发明涉及激光雷达领域,尤其涉及一种基于脉冲编码技术的相干激光雷达。
背景技术
得益于高空分辨率和高测量精度,相干测风激光雷达是激光遥感的热点研究领域。特别值得注意的是1.5μm-2μm的全光纤相干测风激光雷达具有体积小,眼安全、易于集成等优点,是世界各国争相发展的领域。国际上从事脉冲相干激光雷达研究的单位包括NASA、NOAA、美国洛克希德马丁公司(Lockheed Martin)、美国雷神公司(Raytheon)、美国相干公司(CTI)、日本三菱公司(Mitsubishi Electric)、法国Leosphere公司、法国宇宙航空航天研究中心(ONERA)、英国SgurrEnergy公司、英国QinetiQ公司等。国内相干激光雷达的起步较晚,主要单位包括哈尔滨工业大学、电子科技大学、中国电子科技集团公司第二十七研究所、中科院上海光学精密机械研究所、中国科大、中国海洋大学等单位。
本发明的发明人经研究发现:与传统的测风手段相比,由于激光在大气传输过程中,受大气衰减及各种复杂天气状况影响,相干激光雷达接收的回波信号具有很低信噪比和和很大随机性。多普勒激光雷达信号可以描述为一个零均值复高斯随机过程的信号光附加一个统计独立的白噪声。目前,最常用的多普勒频移估计算法是周期图最大值法,即直接提取每个距离门激光雷达时域信号功率谱的最大值对应的频率值。当信噪比(SNR)较低时,频率估计是有偏的,从而增大了测量误差。提高雷达回波信号的信噪比是雷达接收机从噪声中提取回波信号的关键,也是进行后续数据分析处理的前提。
本发明的发明人进一步经研究发现:为提高信噪比,目前采用方法有增加光纤激光器的发射脉冲功率,增大望远镜口径等。采用增加发射脉冲功率的方法,在一定程度上能够提高信噪比,但目前的光纤激光器功率受限,且增大发射激光脉宽,会牺牲探测距离分辨率。采取增大望远镜口径的方法,会大大提高成本。
发明内容
为了解决上述技术问题,本发明采用通信中某些码型对探测脉冲光进行调制,并依据一定解码规则进行处理,可以提高探测距离,获得相同的距离分辨率的同时,有效提高系统的信噪比。
本发明的目的之一是提供一种基于脉冲编码技术的相干激光雷达,在不牺牲探测距离分辨率的情况下,提高系统的信噪比。
本发明的目的是通过以下技术方案实现的:
一种基于脉冲编码技术的相干激光雷达,包括:
激光器、第一分束器、电光调制器、声光调制器、信号发生器、放大器、光学收发模块、耦合器、探测器、采集卡、数据处理装置;其中:
激光器输出的激光信号进入经第一分束器后分为本振光和信号光,信号光经电光调制器后进入声光调制器,信号发生器用于分别向电光调制器和声光调制器输入预设的控制信号;所述电光调制器用于输出脉冲编码序列;从声光调制器输出的信号经放大器放大后进入光学收发模块,光学收发模块用于将输入的信号输出到大气中,并用于接收大气回波信号;从光学收发模块输出的大气回波信号和本振光经耦合器耦合后依次输出到探测器、采集卡和数据处理装置。
进一步地,所述第一分束器为光纤分束器;所述放大器为光纤放大器,所述耦合器为光纤耦合器。
进一步地,所述光学收发模块包括发射望远镜和接收望远镜,所述发射望远镜用于将放大器放大的信号输出到大气中;接收望远镜用于接收大气回波信号;
所述激光器的输出端与分束器的输入端连接,分束器的第一输出端与电光调制器的输入端连接,所述的分束器的第二输出端与耦合器的第一输入端连接;
所述信号发生器分别与电光调制器和声光调制器连接;声光调制器的输出端与放大器的输入端连接,放大器的输出端口与发射望远镜的输入端连接,接收望远镜的输出端与耦合器的第二输入端连接,耦合器的输出端与探测器、采集卡、和数据处理装置依次连接。
进一步地,所述光学收发模块包括环行器和收发望远镜,
所述激光器的输出端与分束器的输入端连接,分束器的第一输出端与电光调制器的输入端连接,所述的分束器的第二输出端与耦合器的第一输入端连接;
所述信号发生器分别与电光调制器和声光调制器连接;声光调制器的输出端与放大器的输入端连接,放大器的输出端口与环行器的输入端连接,所述环行器的收发端与收发望远镜的输入端连接,所述环行器的输出端与耦合器的第二输入端连接,耦合器的输出端与探测器、采集卡、和数据处理装置依次连接。
进一步地,所述采集卡为双通道采集卡,其中,所述双通道采集卡的两个通道中,一个通道用于记录信号发生器调制产生的信号,另一个通道用于记录大气后向散射信号。
本发明还提供了一种基于脉冲编码技术的相干激光雷达风速测量方法,基于脉冲编码技术的相干激光雷达,方法包括:
S1、调制电光调制器,得到M组光脉冲编码序列,分别与大气相互作用,得到M组大气回波信号;M为正整数;
S2、对得到的M组大气回波信号进行距离门的划分,在每个距离门内进行快速傅里叶变换,得到M组大气回波信号的频域信号;
S3、将得到的频域信号按照对应的光脉冲编码序列进入到光纤中的顺序进行排序,并与采集卡采集到的信号发生器调制信号按照时间先后顺序对齐,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,并根据所述频率信息得到多普勒频移;
S4、根据预设的多普勒频移和风速及发射脉冲波长之间的关系,得到大气风速信息。
进一步地,步骤S3之前,还包括:
多次重复步骤S1、S2,得到多组的频域信号;
对得到的频域信号进行累加,得到平均后的多组激光雷达的频域信号。
进一步地,M为4;
所述步骤S3中,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,包括:
Ak、Bk为Golay脉冲编码信号;
其中,N表示码长,k表示对应编码序列中的序号,a、b表示由“+1”和“-1”构成N阶序列;
本发明提供的一种基于脉冲编码技术的相干激光雷达,激光器输出的激光信号进入经第一分束器后分为本振光和信号光,信号光经电光调制器后进入声光调制器,信号发生器用于分别向电光调制器和声光调制器输入预设的控制信号;电光调制器输出脉冲编码序列;从声光调制器输出的信号经放大器放大后进入光学收发模块,光学收发模块用于将输入的信号输出到大气中,并用于接收大气回波信号;从光学收发模块输出的大气回波信号和本振光经耦合器耦合后依次输出到探测器、采集卡和数据处理装置。本发明通过电光调制器输出脉冲编码序列编码,降低了各发射波形之间的时域互相关,并利用Go l ay互补序列的互补性质减小自相关旁瓣,降低了波形脉冲压缩后的旁瓣。另外,本发明按照预定的解码规则进行解码,解决了光纤激光器脉冲峰值功率受限的问题,显著地提高了相干激光雷达信噪比,有效提高了探测距离,且不需要脉冲积累时间,不牺牲探测时间和距离分辨率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案和优点,下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它附图。
图1为本发明公开的一种基于脉冲编码技术的相干激光雷达的一个示意图;
图2为本发明公开的一种基于脉冲编码技术的相干激光雷达的另一示意图;
图3为本发明公开的编码前后相干激光雷达的脉冲发射示意图;
图4(a)为序列Ak的自相关函数的示意图;
图4(b)为序列Bk的自相关函数的示意图;
图4(c)为序列Ak和Bk自相关函数之和的示意图;
图5为本发明中脉冲编码相干激光雷达信号解码流程图;
图6(a)为模拟的相干激光雷达信号距离与风速对应关系的频谱图;
图6(b)为模拟的相干激光雷达信号距离与电流对应关系的频谱图;
图7(a)为相干激光雷达编码前的反演得到的频谱图;
图7(b)为相干激光雷达编码后的反演得到的频谱图;
图8(a)为传统的相干激光雷达系统经过脉冲累加反演后得到的风速结果;
图8(b)为脉冲编码系统累加后的反演的到的风速结果。
具体实施方式
下面介绍的是本发明的多个可能实施例中的一些,旨在提供对本发明的基本了解,并不旨在确认本发明的关键或决定性的要素或限定所要保护的范围。容易理解,根据本发明的技术方案,在不变更本发明的实质精神下,本领域的一般技术人员可以提出可相互替换的其他实现方式。因此,以下具体实施方式以及附图仅是对本发明的技术方案的示例性说明,而不应当视为本发明的全部或者视为对本发明技术方案的限定或限制。
下面的描述中,为描述的清楚和简明,并没有对图中所示的所有多个部件进行描述。附图中示出了多个部件为本领域普通技术人员提供本发明的完全能够实现的公开内容。对于本领域技术人员来说,许多部件的操作都是熟悉而且明显的。
实施例1:
如图1所示,本发明提供了一种基于脉冲编码技术的相干激光雷达,包括:激光器1、第一分束器、电光调制器3、声光调制器4、信号发生器5、放大器、光学收发模块、耦合器、探测器10、采集卡11、数据处理装置;其中:
激光器1输出的激光信号进入经第一分束器后分为本振光和信号光,信号光经电光调制器3后进入声光调制器4,信号发生器5用于分别向电光调制器3和声光调制器4输入预设的控制信号;所述电光调制器3用于输出脉冲编码序列;从声光调制器4输出的信号经放大器放大后进入光学收发模块,光学收发模块用于将输入的信号输出到大气中,并用于接收大气回波信号;从光学收发模块输出的大气回波信号和本振光经耦合器耦合后依次输出到探测器10、采集卡11和数据处理装置。
电光调制器也称为EOM,声光调制器也称为AOM。本发明的电光调制器EOM用来产生编码脉冲,是核心器件。在本发明实施例中,AOM用来产生80MHz的频移。
所述激光器1为连续光激光器;进一步地,所述第一分束器为光纤分束器2,当然,第一分束器还可以是其他分束器,例如棱镜;所述放大器为光纤放大器EDFA6,所述耦合器为光纤耦合器9。
进一步地,所述光学收发模块包括发射望远镜7和接收望远镜8,所述发射望远镜7用于将放大器放大的信号输出到大气中;接收望远镜8用于接收大气回波信号;
所述激光器1的输出端与分束器2的输入端连接,分束器2的第一输出端与电光调制器3的输入端连接,所述的分束器2的第二输出端与耦合器的第一输入端连接;
所述信号发生器5分别与电光调制器3和声光调制器4连接;声光调制器4的输出端与放大器的输入端连接,放大器的输出端口与发射望远镜7的输入端连接,接收望远镜8的输出端与耦合器的第二输入端连接,耦合器的输出端与探测器10、采集卡11、和数据处理装置依次连接。数据处理装置可以是电脑9或计算机。
进一步地,所述采集卡11为双通道采集卡,其中,所述双通道采集卡的两个通道中,一个通道用于记录信号发生器调制产生的信号,另一个通道用于记录大气后向散射信号。
本发明还提供了一种基于脉冲编码技术的相干激光雷达风速测量方法,基于脉冲编码技术的相干激光雷达,方法包括:
S1、调制电光调制器,得到M组光脉冲编码序列,分别与大气相互作用,得到M组大气回波信号;M为正整数;
S2、对得到的M组大气回波信号进行距离门的划分,在每个距离门内进行快速傅里叶变换,得到M组大气回波信号的频域信号;
S3、将得到的频域信号按照对应的光脉冲编码序列进入到光纤中的顺序进行排序,并与采集卡采集到的信号发生器调制信号按照时间先后顺序对齐,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,并根据所述频率信息得到多普勒频移;
S4、根据预设的多普勒频移和风速及发射脉冲波长之间的关系,得到大气风速信息。
进一步地,步骤S3之前,还包括:
多次重复步骤S1、S2,得到多组的频域信号;
对得到的频域信号进行累加,得到平均后的多组激光雷达的频域信号。
进一步地,M为4;
所述步骤S3中,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,包括:
Ak、Bk为Golay脉冲编码信号;
其中,N表示码长,k表示对应编码序列中的序号,a、b表示由“+1”和“-1”构成N阶序列;
下面以一个具体实施场景对本发明的方案进一步阐释。
如图1所示,激光器1的输出端与光纤分束器2的输入端连接,光纤分束器2输出端a与EOM3的输入端连接,所述的光纤分束器2的输出端与光纤耦合器9的第一输入端连接。
所述的EOM3的输出端和信号发生器5的输出端分别和电光调制器3和声光调制器4连接;光开关4的输出端与EDFA6的输入端连接,EDFA6的输出端口与发射望远镜连接,激光经由大气,返回到接收望远镜8,并与光纤耦合器9的第二输入端口连接,本振光和探测光经由探测器10,探测器10的输出端和采集卡11的输入端连接,采集卡11的输出端和电脑12连接。
图3为本发明公开的编码前后相干激光雷达的脉冲发射示意图。图3的上面一行为编码前的相干激光雷达的脉冲发射示意图,图3的下面一行为编码后的相干激光雷达的脉冲发射示意图。本发明实施例中采用Golay互补脉冲序列。Golay互补脉冲序列的结构方式通常由低阶的互补序列构成高阶的互补序列,即将一序列b接在一序列a的后面形成一高阶序列,再将序列b的反码接在a的后面形成一个高阶序列,从而构成一对互补序列。即:
Golay互补码包含一对双极性码,然而,由于在光纤中只能传输正的光脉冲,因此不能像在电域或声域中使用双极性的互补序列。本发明为了解决上述问题,通过偏置的方法,通过传送4组单极性正脉冲,从而达到与Golay双极性脉冲相同的效果。即:
图4(a)、图4(b)和图4(c)展示了采用的Golay脉冲序列相关函数及其和的函数示意。图4(a)为序列Ak的自相关函数的示意图;图4(b)为序列Bk的自相关函数的示意图;图4(c)为序列Ak和Bk自相关函数之和的示意图。从图中可以看出Golay互补码具有很好的自相关性,降低了波形脉冲压缩后的旁瓣。
图5为本发明中脉冲编码相干激光雷达信号解码流程图。
首先对得到的四组激光雷达信号进行距离门的划分,在每个距离门内进行快速傅里叶变换,得到4组激光雷达的频域信号。重述上述步骤N次,对得到的频域信号进行累加,得到平均后的4组激光雷达的频域信号。
进一步,将第一和第二组单极性脉冲产生的后向散射信号频谱相减,并将第三和第四组单极性脉冲产生的后向散射信号频谱相减。
进一步的,将Golay互补序列的第一码Ak与上一步中第一组和第二组单极性码相减结果取相关。Golay互补序列的第二码Bk与上一步中第三组和第四组单极性码相减结果取相关。
最后将上一步得到的两组结果相加,即:
式中表示编码脉冲Uk、Wk、产生的频谱。corr表示取相关。在得到后向散射频谱后,采用周期图最大值法,对每个距离门的功率谱的最大值所对应的频率值进行提取,该频率值即为多普勒频移值,得到每个距离门处的多普勒频移值后,根据多普勒频移和大气风速之间的关系Δν=2V/λ(Δν表示多普勒频移,V表示风速),反演出最后的风速信息,λ为发射激光波长。
需要说明的是,大气的后向散射谱是时域信号做FFT变换后的得到的,假设时域信号有10000个点,每100个点划分为一个距离门,共100个距离门。每个距离门进行FFT,得到的就是大气的后向散射谱,每个频谱中都有一个峰值,提取这个峰值所对应的值,这个值就是对应该距离门的多普勒信号,即大气中该位置的频率信息。大气中频率信息是根据大气后向散射谱提取得到的。
图7(a)为相干激光雷达编码前的反演得到的频谱图;图7(b)为相干激光雷达编码后的反演得到的频谱图。可见,采用脉冲编码技术后,信噪比得到了提高。为了评估提出的脉冲编码技术对信噪比的影响,对传统的单脉冲相干激光雷达信号进行累加平均。对编码系统的频域信号取2500次累加平均,由于解码一组编码信号需要四组时域信号,为了使传统的单脉冲相干激光雷达和编码的相干激光雷达测量次数相同,对传统相干激光雷达的频域信号进行10000次叠加,得到了如图6(a)和图6(b)所示的频谱图,图6(a)为模拟的相干激光雷达信号距离与风速对应关系的频谱图;图6(b)为模拟的相干激光雷达信号距离与电流对应关系的频谱图,可见,使用编码技术后,提高了探测距离,系统信噪比提高。
图8(a)为传统的相干激光雷达系统经过脉冲累加反演后得到的风速结果;图8(b)为脉冲编码系统累加后的反演的到的风速结果。从图8(a)和图8(b)可见,采用脉冲编码可提高风速反演的信噪比。
本发明的保护范围不限于Golay脉冲编码,其他编码形式,包括但不限于Simplex码、CCPOMS码等也在本发明的保护范围内。其中Simplex码是从归一化的Hadamard矩阵中得到的,其中Hadamard矩阵满足HHT=nIn,其中In为单位矩阵。去掉Hadamard矩阵的第一行与第一列,并将其中的元素“0”替换为“1”,“1”替换成“-1”,则矩阵中每一行码都为Simplex码,Simplex码为单极性码,可以用于激光雷达系统中。CCPOMS源于Golay和Shapiro编码思想,其构造方式如下,
本发明实施例提供的一种基于脉冲编码技术的相干激光雷达具有如下有益效果:
1)Golay互补序列构成的脉冲串进行空时编码,降低各发射波形之间的时域互相关,并利用Golay互补序列的互补性质减小自相关旁瓣,降低了波形脉冲压缩后的旁瓣。
2)使用Golay编码脉冲作为相干激光雷达的探测脉冲,解决了光纤激光器脉冲峰值功率受限的问题,显著地提高了相干激光雷达信噪比,有效提高了探测距离,且不需要脉冲积累时间,不牺牲探测时间和距离分辨率。
实施例2:
如图2所示,本发明提供了一种基于脉冲编码技术的相干激光雷达,包括:激光器1、第一分束器、电光调制器3、声光调制器4、信号发生器5、放大器、光学收发模块、耦合器、探测器10、采集卡11、数据处理装置;其中:
激光器1输出的激光信号进入经第一分束器后分为本振光和信号光,信号光经电光调制器3后进入声光调制器4,信号发生器5用于分别向电光调制器3和声光调制器4输入预设的控制信号;所述电光调制器3用于输出脉冲编码序列;从声光调制器4输出的信号经放大器放大后进入光学收发模块,光学收发模块用于将输入的信号输出到大气中,并用于接收大气回波信号;从光学收发模块输出的大气回波信号和本振光经耦合器耦合后依次输出到探测器10、采集卡11和数据处理装置。
电光调制器也称为EOM,声光调制器也称为AOM。
所述激光器1为连续光激光器;进一步地,所述第一分束器为光纤分束器2,当然,第一分束器还可以是其他分束器,例如棱镜;所述放大器为光纤放大器EDFA6,所述耦合器为光纤耦合器9。
进一步地,所述光学收发模块包括环行器71和收发望远镜81,所述激光器1的输出端与分束器2的输入端连接,分束器2的第一输出端与电光调制器3的输入端连接,所述的分束器2的第二输出端与耦合器的第一输入端连接;
所述信号发生器5分别与电光调制器3和声光调制器4连接;声光调制器4的输出端与放大器的输入端连接,放大器的输出端口与环行器71的输入端连接,所述环行器71的收发端与收发望远镜81的输入端连接,所述环行器71的输出端与耦合器的第二输入端连接,耦合器的输出端与探测器10、采集卡11、和数据处理装置依次连接。
进一步地,所述采集卡11为双通道采集卡,其中,所述双通道采集卡的两个通道中,一个通道用于记录信号发生器调制产生的信号,另一个通道用于记录大气后向散射信号。
本发明还提供了一种基于脉冲编码技术的相干激光雷达风速测量方法,基于脉冲编码技术的相干激光雷达,方法包括:
S1、调制电光调制器,得到M组光脉冲编码序列,分别与大气相互作用,得到M组大气回波信号;M为正整数;
S2、对得到的M组大气回波信号进行距离门的划分,在每个距离门内进行快速傅里叶变换,得到M组大气回波信号的频域信号;
S3、将得到的频域信号按照对应的光脉冲编码序列进入到光纤中的顺序进行排序,并与采集卡采集到的信号发生器调制信号按照时间先后顺序对齐,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,并根据所述频率信息得到多普勒频移;
S4、根据预设的多普勒频移和风速及发射脉冲波长之间的关系,得到大气风速信息。
进一步地,步骤S3之前,还包括:
多次重复步骤S1、S2,得到多组的频域信号;
对得到的频域信号进行累加,得到平均后的多组激光雷达的频域信号。
进一步地,M为4;
所述步骤S3中,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,包括:
Ak、Bk为Golay脉冲编码信号;
其中,N表示码长,k表示对应编码序列中的序号,a、b表示由“+1”和“-1”构成N阶序列;
下面以一个具体实施场景对本发明的方案进一步阐释。
图3为本发明公开的编码前后相干激光雷达的脉冲发射示意图。图3的上面一行为编码前的相干激光雷达的脉冲发射示意图,图3的下面一行为编码后的相干激光雷达的脉冲发射示意图。本发明实施例中采用Golay互补脉冲序列。Golay互补脉冲序列的结构方式通常由低阶的互补序列构成高阶的互补序列,即将一序列b接在一序列a的后面形成一高阶序列,再将序列b的反码接在a的后面形成一个高阶序列,从而构成一对互补序列。即:
Golay互补码包含一对双极性码,然而,由于在光纤中只能传输正的光脉冲,因此不能像在电域或声域中使用双极性的互补序列。本发明为了解决上述问题,通过偏置的方法,通过传送4组单极性正脉冲,从而达到与Golay双极性脉冲相同的效果。即:
图4为本发明中采用的Golay脉冲序列相关函数及其和的函数示意图。图4(a)为序列Ak的自相关函数;图4(b)为序列Bk的自相关函数;图4(c)为序列Ak和Bk自相关函数之和。从图中可以看出Golay互补码具有很好的自相关性,降低了波形脉冲压缩后的旁瓣。
图5为本发明中脉冲编码相干激光雷达信号解码流程图。
首先对得到的四组激光雷达信号进行距离门的划分,在每个距离门内进行快速傅里叶变换,得到4组激光雷达的频域信号。重述上述步骤N次,对得到的频域信号进行累加,得到平均后的4组激光雷达的频域信号。
进一步,将第一和第二组单极性脉冲产生的后向散射信号频谱相减,并将第三和第四组单极性脉冲产生的后向散射信号频谱相减。
进一步的,将Golay互补序列的第一码Ak与上一步中第一组和第二组单极性码相减结果取相关。Golay互补序列的第二码Bk与上一步中第三组和第四组单极性码相减结果取相关。
最后将上一步得到的两组结果相加,即:
式中表示编码脉冲Uk、Wk、产生的频谱。corr表示取相关。在得到后向散射频谱后,采用周期图最大值法,对每个距离门的功率谱的最大值所对应的频率值进行提取,该频率值即为多普勒频移值,得到每个距离门处的多普勒频移值后,根据多普勒频移和大气风速之间的关系Δν=2V/λ(Δν表示多普勒频移,V表示风速),反演出最后的风速信息,λ为发射激光波长。
需要说明的是,大气的后向散射谱是时域信号做FFT变换后的得到的,假设时域信号有10000个点,每100个点划分为一个距离门,共100个距离门。每个距离门进行FFT,得到的就是大气的后向散射谱,每个频谱中都有一个峰值,提取这个峰值所对应的值,这个值就是对应该距离门的多普勒信号,即大气中该位置的频率信息。大气中频率信息是根据大气后向散射谱提取得到的。
图7为相干激光雷达编码前后的反演得到的频谱图。图7a为相干激光雷达编码前的反演得到的频谱图;图7b为相干激光雷达编码后的反演得到的频谱图。可见,采用脉冲编码技术后,信噪比得到了提高。为了评估提出的脉冲编码技术对信噪比的影响,对传统的单脉冲相干激光雷达信号进行累加平均。对编码系统的频域信号取2500次累加平均,由于解码一组编码信号需要四组时域信号,为了使传统的单脉冲相干激光雷达和编码的相干激光雷达测量次数相同,对传统相干激光雷达的频域信号进行10000次叠加,得到了如图6所示的频谱图,图6a为距离和风速对应关系的频谱图;图6b从图6发现,使用编码技术后,提高了探测距离,系统信噪比提高。
图8为相干激光雷达编码前后反演得到的风速图。图8(a)为传统的相干激光雷达系统经过脉冲累加反演后得到的风速结果;图8(b)为脉冲编码系统累加后的反演的到的风速结果。从图8可见,采用脉冲编码可提高风速反演的信噪比。
本发明的保护范围不限于Golay脉冲编码,其他编码形式,包括但不限于Simplex码、CCPOMS码等也在本发明的保护范围内。其中Simplex码是从归一化的Hadamard矩阵中得到的,其中Hadamard矩阵满足HHT=nIn,其中In为单位矩阵。去掉Hadamard矩阵的第一行与第一列,并将其中的元素“0”替换为“1”,“1”替换成“-1”,则矩阵中每一行码都为Simplex码,Simplex码为单极性码,可以用于激光雷达系统中。CCPOMS源于Golay和Shapiro编码思想,其构造方式如下:
本发明实施例提供的一种基于脉冲编码技术的相干激光雷达具有如下有益效果:
1)Golay互补序列构成的脉冲串进行空时编码,降低各发射波形之间的时域互相关,并利用Golay互补序列的互补性质减小自相关旁瓣,降低了波形脉冲压缩后的旁瓣。
2)使用Golay编码脉冲作为相干激光雷达的探测脉冲,解决了光纤激光器脉冲峰值功率受限的问题,显著地提高了相干激光雷达信噪比,有效提高了探测距离,且不需要脉冲积累时间,不牺牲探测时间和距离分辨率。
以上例子主要说明了本发明提供的系统和方法。尽管只对其中一些本发明的实施方式进行了描述,但是本领域普通技术人员应当了解,本发明可以在不偏离其主旨与范围内以许多其他的形式实施。因此,所展示的例子与实施方式被视为示意性的而非限制性的,在不脱离如所附各权利要求所定义的本发明精神及范围的情况下,本发明可能涵盖各种的修改与替换。
Claims (6)
1.一种基于脉冲编码技术的相干激光雷达,其特征在于,包括:
激光器(1)、第一分束器(2)、电光调制器(3)、声光调制器(4)、信号发生器(5)、放大器(6)、光学收发模块、耦合器(9)、探测器(10)、采集卡(11)、数据处理装置(12);其中:
所述激光器(1)输出的激光信号进入经所述第一分束器(2)后分为本振光和信号光,所述信号光经电光调制器(3)后进入声光调制器(4),所述信号发生器(5)用于分别向电光调制器(3)和声光调制器(4)输入预设的控制信号;所述电光调制器(3)用于输出脉冲编码序列;从声光调制器(4)输出的信号经放大器(6)放大后进入光学收发模块,光学收发模块用于将输入的信号输出到大气中,并用于接收大气回波信号;从光学收发模块输出的大气回波信号和本振光经耦合器(9)耦合后依次输出到探测器(10)、采集卡(11)和数据处理装置(12);
所述脉冲编码序列包括:
Ak、Bk为Golay脉冲编码信号;
其中,N表示码长,k表示对应编码序列中的序号,a、b表示由“+1”和“-1”构成N阶序列;
一种基于脉冲编码技术的相干激光雷达风速测量方法,基于上述一种基于脉冲编码技术的相干激光雷达,包括:
S1、调制电光调制器,得到M组光脉冲编码序列,分别与大气相互作用,得到M组大气回波信号;M为正整数;
S2、对得到的M组大气回波信号进行距离门的划分,在每个距离门内进行快速傅里叶变换,得到M组大气回波信号的频域信号;
S3、将得到的频域信号按照对应的光脉冲编码序列进入到光纤中的顺序进行排序,并与采集卡采集到的信号发生器调制信号按照时间先后顺序对齐,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,并根据所述频率信息得到多普勒频移;
S4、根据预设的多普勒频移和风速及发射脉冲波长之间的关系,得到大气风速信息;
M为4;
所述步骤S3中,根据预设的解码规则,对得到的M组频域信号进行解码,得到大气中不同位置处的频率信息,包括:
Ak、Bk为Golay脉冲编码信号;
其中,N表示码长,k表示对应编码序列中的序号,a、b表示由“+1”和“-1”构成N阶序列;
2.根据权利要求1所述的基于脉冲编码技术的相干激光雷达,其特征在于,所述第一分束器(2)为光纤分束器;所述放大器(6)为光纤放大器,所述耦合器(9)为光纤耦合器。
3.根据权利要求2所述的基于脉冲编码技术的相干激光雷达,其特征在于,所述光学收发模块包括发射望远镜(7)和接收望远镜(8),所述发射望远镜(7)用于将放大器(6)放大的信号输出到大气中;接收望远镜(8)用于接收大气回波信号;
所述激光器(1)的输出端与分束器(2)的输入端连接,分束器(2)的第一输出端与电光调制器(3)的输入端连接,所述的分束器(2)的第二输出端与耦合器(9)的第一输入端连接;
所述信号发生器(5)分别与电光调制器(3)和声光调制器(4)连接;声光调制器(4)的输出端与放大器(6)的输入端连接,放大器(6)的输出端口与发射望远镜(7)的输入端连接,接收望远镜(8)的输出端与耦合器(9)的第二输入端连接,耦合器(9)的输出端与探测器(10)、采集卡(11)、和数据处理装置(12)依次连接。
4.根据权利要求2所述的基于脉冲编码技术的相干激光雷达,其特征在于,所述光学收发模块包括环行器(71)和收发望远镜(81),
所述激光器(1)的输出端与分束器(2)的输入端连接,分束器(2)的第一输出端与电光调制器(3)的输入端连接,所述的分束器(2)的第二输出端与耦合器(9)的第一输入端连接;
所述信号发生器(5)分别与电光调制器(3)和声光调制器(4)连接;声光调制器(4)的输出端与放大器(6)的输入端连接,放大器(6)的输出端口与环行器(71)的输入端连接,所述环行器(71)的收发端与收发望远镜(81)的输入端连接,所述环行器(71)的输出端与耦合器(9)的第二输入端连接,耦合器(9)的输出端与探测器(10)、采集卡(11)、和数据处理装置(12)依次连接。
5.根据权利要求1所述的基于脉冲编码技术的相干激光雷达,其特征在于,所述采集卡(11)为双通道采集卡,其中,所述双通道采集卡的两个通道中,一个通道用于记录信号发生器调制产生的信号,另一个通道用于记录大气后向散射信号。
6.根据权利要求1所述的一种基于脉冲编码技术的相干激光雷达,其特征在于,步骤S3之前,还包括:
多次重复步骤S1、S2,得到多组的频域信号;
对得到的频域信号进行累加,得到平均后的多组激光雷达的频域信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810335326.8A CN108594256B (zh) | 2018-04-16 | 2018-04-16 | 一种基于脉冲编码技术的相干激光雷达 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810335326.8A CN108594256B (zh) | 2018-04-16 | 2018-04-16 | 一种基于脉冲编码技术的相干激光雷达 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108594256A CN108594256A (zh) | 2018-09-28 |
CN108594256B true CN108594256B (zh) | 2021-10-12 |
Family
ID=63622673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810335326.8A Active CN108594256B (zh) | 2018-04-16 | 2018-04-16 | 一种基于脉冲编码技术的相干激光雷达 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108594256B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11614543B2 (en) * | 2018-10-09 | 2023-03-28 | GM Global Technology Operations LLC | Transimpedance amplifier for Lidar system |
GB2578788B (en) * | 2018-11-09 | 2022-10-05 | Toshiba Kk | An investigation system and method |
CN110082778B (zh) * | 2019-05-27 | 2023-04-21 | 中国科学技术大学 | 基于单光子探测的相干测风激光雷达 |
CN111190189B (zh) * | 2020-01-13 | 2022-11-08 | 中国科学院上海光学精密机械研究所 | 多功能双调频相干激光雷达 |
CN112327272B (zh) * | 2020-11-06 | 2021-06-29 | 深圳煜炜光学科技有限公司 | 一种脉冲可编码的激光雷达装置及其检测方法 |
CN112965086B (zh) * | 2021-02-08 | 2022-07-15 | 中国科学技术大学 | 基于脉冲编码的1.5μm波长气溶胶探测激光雷达及信号解码方法 |
CN113406603B (zh) * | 2021-08-19 | 2021-12-03 | 武汉镭晟科技有限公司 | 一种用于相干激光雷达的激光器模块 |
CN115963479A (zh) * | 2021-10-13 | 2023-04-14 | 华为技术有限公司 | 一种信号发射、处理方法及相关装置 |
CN113985447B (zh) * | 2021-12-30 | 2022-04-08 | 青岛镭测创芯科技有限公司 | 一种相干测风激光雷达及测量方法 |
CN114428229B (zh) * | 2022-01-02 | 2024-09-06 | 西安电子科技大学 | 任意时间间隔脉冲编码雷达的信噪比提升方法 |
CN114924255A (zh) * | 2022-04-20 | 2022-08-19 | 厦门大学 | 基于编码脉冲的水体激光雷达及编码脉冲激光生成方法 |
CN115236697B (zh) * | 2022-09-23 | 2022-12-13 | 中国人民解放军国防科技大学 | 分时多频脉冲测风激光雷达系统和风速测量方法 |
CN115453574B (zh) * | 2022-10-27 | 2023-04-25 | 南京信息工程大学 | 一种大气多参数探测多功能激光雷达 |
CN115951332B (zh) * | 2023-03-15 | 2023-05-23 | 中国人民解放军国防科技大学 | 二元脉冲幅度编码测风激光雷达风速测量方法及系统 |
CN118011416B (zh) * | 2024-04-09 | 2024-06-04 | 北京中科国光量子科技有限公司 | 一种基于稳定相位编解码的激光雷达及相位补偿方法 |
CN118244285B (zh) * | 2024-05-28 | 2024-08-30 | 北醒(北京)光子科技有限公司 | 一种激光雷达及其测距方法与控制系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4893923A (en) * | 1986-10-06 | 1990-01-16 | Laser Science, Inc. | Doppler laser radar system |
US5815250A (en) * | 1997-05-27 | 1998-09-29 | Coherent Technologies, Inc. | Doublet pulse coherent laser radar for precision range and velocity measurements |
CN104865564A (zh) * | 2015-05-19 | 2015-08-26 | 王明军 | 一种大尺寸目标激光雷达成像技术 |
CN106226778A (zh) * | 2016-08-23 | 2016-12-14 | 成都信息工程大学 | 一种高分辨率测量远程目标的相干激光雷达系统 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2001280758A1 (en) * | 2000-07-24 | 2002-02-05 | University Of Kansas Center For Research, Inc. | Laser radar system and method |
US7463368B2 (en) * | 2003-09-10 | 2008-12-09 | Metris Canada Inc | Laser projection system, intelligent data correction system and method |
WO2011016892A2 (en) * | 2009-05-15 | 2011-02-10 | Michigan Aerospace Corporation | Range imaging lidar |
CN101750608B (zh) * | 2009-12-25 | 2012-07-11 | 哈尔滨工业大学 | 扫描相干激光成像雷达程控装置及方法 |
US9599714B2 (en) * | 2012-06-25 | 2017-03-21 | Mitsubishi Electric Corporation | Wind measurement coherent lidar |
WO2015189915A1 (ja) * | 2014-06-10 | 2015-12-17 | 三菱電機株式会社 | レーザレーダ装置 |
CN104820223B (zh) * | 2015-04-10 | 2017-07-14 | 中国科学院上海光学精密机械研究所 | 基于m序列相位编码的光学域匹配滤波测距装置 |
CN106556881A (zh) * | 2015-09-28 | 2017-04-05 | 东莞前沿技术研究院 | 风场监测系统 |
US10401499B2 (en) * | 2016-03-16 | 2019-09-03 | Raytheon Company | Laser beam projection system with dynamic phase compensation |
CN105866795B (zh) * | 2016-05-19 | 2018-08-03 | 中国科学技术大学 | 基于全光纤扫描f-p干涉仪的大动态风场探测激光雷达 |
CN106019313B (zh) * | 2016-08-05 | 2018-09-07 | 中国科学技术大学 | 基于偏振双边缘的单像素探测测风激光雷达 |
CN206114903U (zh) * | 2016-08-23 | 2017-04-19 | 成都信息工程大学 | 一种高分辨率测量远程目标的相干激光雷达系统 |
CN106154289B (zh) * | 2016-09-23 | 2018-09-07 | 中国科学技术大学 | 基于差分受激布里渊增益效应的直接测风激光雷达 |
CN106940444B (zh) * | 2017-02-27 | 2019-09-13 | 南京红露麟激光雷达科技有限公司 | 基于微波差分增益的相干多普勒测风激光雷达 |
CN107144829B (zh) * | 2017-06-29 | 2019-11-19 | 南京信息工程大学 | 一种高效的激光雷达回波信号去噪方法 |
CN107783144A (zh) * | 2017-10-30 | 2018-03-09 | 南京牧镭激光科技有限公司 | 测风激光雷达装置 |
-
2018
- 2018-04-16 CN CN201810335326.8A patent/CN108594256B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4893923A (en) * | 1986-10-06 | 1990-01-16 | Laser Science, Inc. | Doppler laser radar system |
US5815250A (en) * | 1997-05-27 | 1998-09-29 | Coherent Technologies, Inc. | Doublet pulse coherent laser radar for precision range and velocity measurements |
EP1042688A2 (en) * | 1997-05-27 | 2000-10-11 | Coherent Technologies Inc. | Doublet pulse coherent laser radar for precision range and velocity measurements |
CN104865564A (zh) * | 2015-05-19 | 2015-08-26 | 王明军 | 一种大尺寸目标激光雷达成像技术 |
CN106226778A (zh) * | 2016-08-23 | 2016-12-14 | 成都信息工程大学 | 一种高分辨率测量远程目标的相干激光雷达系统 |
Non-Patent Citations (2)
Title |
---|
Fiber-based coherent pulsed Doppler lidar for atmospheric monitoring;Pearson;《PROCEEDINGS OF SPIE》;20021231;50-57 * |
基于概率统计的光子激光雷达性能提高的理论与实验研究;张子静;《中国博士学位论文全文数据库 信息科技辑》;20160315;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108594256A (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108594256B (zh) | 一种基于脉冲编码技术的相干激光雷达 | |
CN106092305B (zh) | 分布式光纤传感系统及其振动检测定位方法 | |
CN104931973B (zh) | 用于激光雷达系统的非对称相位编码测距方法 | |
CN102636251B (zh) | 一种相位敏感光时域反射系统及提高其信噪比的方法 | |
CN106054204A (zh) | 一种面向长距离高精度的复合式激光测距方法及系统 | |
CN110455324B (zh) | 一种基于cp-φotdr的高重复率准分布式传感系统及其实现方法 | |
CN104181509B (zh) | 基于跳频和多相交替码的非相干散射雷达信号处理方法 | |
CN112965086B (zh) | 基于脉冲编码的1.5μm波长气溶胶探测激光雷达及信号解码方法 | |
CN111998933B (zh) | 一种基于脉冲编码的光纤光栅振动测量装置及方法 | |
CN110161521B (zh) | 一种基于真随机编码的光子计数激光雷达 | |
CN106603149A (zh) | 高速率激光通信方法和高精度激光测距一体化方法 | |
CN112697257B (zh) | 无衰落多波长分布式声波传感系统和差分旋转矢量叠加法 | |
CN110794421A (zh) | 伪随机码时延自差分干涉三维成像激光雷达方法及装置 | |
CN104158587B (zh) | 基于周期开关键控混沌信号的光时域反射方法 | |
CN113848533B (zh) | 一种多载频多相位编码探测干扰共享信号生成方法 | |
CN110646779A (zh) | 相位编码不饱和调制方法、装置、激光雷达测距测速方法及激光雷达系统 | |
CN104931975A (zh) | 基于微波光子信号处理的相位编码激光成像雷达 | |
Shiloh et al. | Highly-sensitive distributed dynamic strain sensing via perfect periodic coherent codes | |
CN211577419U (zh) | 相位编码不饱和调制装置和激光雷达系统 | |
CN109507684B (zh) | 基于类噪声脉冲的高空间分辨率探测系统及探测方法 | |
EP4024085A1 (en) | Phase encoding unsaturated modulation method and device, laser radar ranging and speed measurement method and laser radar system | |
CN111412935A (zh) | 基于时分复用的高重复率准分布式传感系统 | |
Zhu et al. | High anti-interference 3D imaging LIDAR system based on digital chaotic pulse position modulation | |
CN107171716B (zh) | 一种基于相关编码的在线链路监测系统及方法 | |
CN114111860B (zh) | 基于多频脉冲编码的分布式φ-otdr传感方法及系统 |
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 | ||
TA01 | Transfer of patent application right | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20191031 Address after: No. 164, Changguo South Road, Yanling Town, Danyang City, Zhenjiang City, Jiangsu Province Applicant after: Xia Hedi Address before: 231241 Anhui Hefei Shushan Economic Development Zone, Wen Shui Road, electricity supplier park three phase 2 building A 8 floor 8428 room. Applicant before: Hefei Finel Photoelectric Technology Co., Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |