CN103913723A - 一种处理非相干散射雷达散射信号的波形设计方法 - Google Patents
一种处理非相干散射雷达散射信号的波形设计方法 Download PDFInfo
- Publication number
- CN103913723A CN103913723A CN201410098372.2A CN201410098372A CN103913723A CN 103913723 A CN103913723 A CN 103913723A CN 201410098372 A CN201410098372 A CN 201410098372A CN 103913723 A CN103913723 A CN 103913723A
- Authority
- CN
- China
- Prior art keywords
- codes
- code
- design method
- long pulse
- alternate
- 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.)
- Pending
Links
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种处理非相干散射雷达散射信号的波形设计方法,在非相干散射雷达信号处理的基础上,结合巴克码,长脉冲编码以及交替码对回波信号进行处理,根据三者自相关函数表达式,利用MATlAB软件求得三种编码的自相关函数以及其功率谱。本发明的优点是:可探测到电离层多种参数,结合巴克码、长脉冲编码以及交替码的信号处理过程,利用MATLAB软件求出其自相关函数和功率谱。
Description
技术领域
本发明涉及一种波形设计方法,尤其涉及一种处理非相干散射雷达散射信号的波形设计方法。
背景技术
电离层对无线电通信、导航、测量以及人类空间活动等都有着非常重要的影响,从而对电离层的研究具有十分重要的意义,非相干散射雷达是探测电离层的强有力工具,是唯一能够测量整个电离层高度的设备,可覆盖从六十公里至两千公里左右的范围,从而探测到电离层多种参数,包括电子密度、电子温度、离子密度和离子温度等。由于非相干散射雷达的建造和运行费用都比较昂贵,所以并没有广泛运用,目前全世界大概只有十部非相干散射雷达,集中于美洲和北欧。我国于2008年开工建设的空间科学重大项目——子午工程,其中涉及了领头设备非相干散射雷达技术,本发明是一种处理非相干散射雷达回波信号的方法,具体给出了巴克码、长脉冲编码和交替码的信号处理算法。
探测电离层是指用直接或间接的探测方法,获得电离层物理参量及其变化规律的工作,之前就有用连续波法来进行探测,后来发明电离层垂直探测装置,直到戈登提出的电离层非相干散射探测理论,从此拉开了非相干散射雷达发展的序幕。之后人们开始研究一种新的电离层探测手段——非相干散射雷达,是目前测量空间大气环境参数最先进的手段。
研究非相干散射雷达的过程中,有两个问题需要解决,一是散射信号功率谱的计算,二是对等离子体散射谱进行建模分析并且反演得到电离层的参量。第一个就是涉及到信号处理的部分,在没有实测数据的情况下,可以在非相干散射雷达信号处理的基础上,结合巴克码,长脉冲编码以及交替码对回波信号进行处理,从而有了非相干散射雷达散射信号的波形设计方法。
在深刻理解非相干散射基本原理的基础上,熟悉回波信号处理的一般过程,研究基本雷达调制方式:巴克码,长脉冲编码以及交替码,根据其自相关函数表达式,利用MATlAB软件可求得三种编码的自相关函数以及其功率谱。
发明内容
本发明的目的在于提供一种处理非相干散射雷达散射信号的波形设计方法,当我们探测电离层时,不同的编码方式和回波处理过程将会影响其探测谱的形状,为了消除这两方面所带来的不平衡性,所以需要对求得的功率谱进行修正,由此我们引入了模糊函数这一概念,在整个波形设计过程中,其核心就是模糊函数,包括幅度模糊函数、二维模糊函数、距离幅度模糊函数以及距离模糊函数。非相干散射雷达散射信号处理主要有巴克码、长脉冲编码以及交替码三种编码方式。
本发明是这样来实现的,其特征在于:在非相干散射雷达信号处理的基础上,结合巴克码,长脉冲编码以及交替码对回波信号进行处理,根据三者自相关函数表达式,利用MATlAB软件求得三种编码的自相关函数以及其功率谱。
本发明所述巴克码是一种相位调制,是二相编码信号的一种,在相位编码中,发射脉冲被划分为一系列单元,每个单元具有一个确定的相位,一般在实际处理中,相位有 和0,这种相位编码对应的两个包络值是-1和1或者0和1。巴克码是一种脉冲压缩编码波形,可获得较高的距离分辨率和最大可能探测距离;其特征在于:若一个n 位的巴克码{},每个码元只可能取值+1或-1,则它必然满足条件:
在处理巴克码的过程中采用的是匹配滤波,使其输出信噪比最大,其滤波系数是巴克码包络的镜像函数,然后对经过滤波器卷积后的回波信号进行采样,得到一系列的数值,由上式可知,信号的自相关可以求出。但是在实际应用中,直接进行卷积会造成很大的计算量,其中卷积了很多不必要的点,在计算自相关函数很多点是没有用的。所以为了简化以及更加快速的计算,直接将滤波器系数和回波信号相乘,根据lag时延间隔的要求,对乘积结果进行相关求和,就可以得到自相关函数。对于巴克码编码信号的自相关函数就是其距离幅度模糊函数,并且在lag0处的中心有一个较大的峰值,较小的旁瓣,可以得到较高的距离分辨率,但是随着lag的增大,会导致旁瓣的增大以及距离分辨率的下降,所以一般是求得lag0处的自相关,即巴克码只能求某一高度电离层等离子体零时延处的自相关。巴克码的处理方式比较简单,匹配滤波后每个高度得到一个滤波结果,将其平方得到功率,有比较高的距离分辨率,但是只能求得每个高度上的功率剖面,含有较少的信息量,所以,一般在实际应用中,将其作为一种辅助手段结合其它编码进行探测实验。
本发明所述长脉冲编码调制是最简单的调制方式,尽管调制方式简单,但是数据采样和距离门的确定比较困难,脉冲宽度较长,其高度分辨率在10km~30km之间,因此通常适用于探测电离层F层或低仰角和低空间分辩率探测,而在E层等低电离层,因为各电离层参数的垂向梯度较大,需要较小的距离分辨率,这种调制方式不能满足探测要求。长脉冲的自相关函数可以表示为:
其中,i 表示第 i 个高度的距离门,n 是脉冲宽度与 lag 时延间隔。
本发明所述交替码是将单个长脉冲划分为一系列小单元,所有的单元之间可以进行相关操作,且单元之间不影响,可以分为强交替码和弱交替码,因为强交替码和弱交替码的工作原理是相同的,所以其工作原理可简单的由四位弱交替码给出。
本发明的优点是:可探测到电离层多种参数,结合巴克码、长脉冲编码以及交替码的信号处理过程,利用MATLAB软件求出其自相关函数和功率谱。
附图说明
图1为本发明巴克码编码回波信号功率谱图。
图2为本发明长脉冲编码回波信号功率谱图。
图3为本发明交替码编码回波信号功率谱图。
具体实施方式
以下结合附图说明对本发明的实施例作进一步详细描述,但本实施例并不用于限制本发明,凡是采用本发明的相似结构及其相似变化,均应列入本发明的保护范围。
如图1所示,巴克码处采用的是匹配滤波,对于巴克码编码信号的自相关函数就是其距离幅度模糊函数,在处理巴克码调制的过程中运用了脉冲压缩和非相干累积,首先将发射信号经过13位巴克码编码调制,然后将其进行脉冲压缩得到回波信号的自相关函数,再做傅里叶变换得到功率谱。
如图2所示,长脉冲编码调制过程中,假设脉冲宽度是lag时延间隔的20倍,即n为20,阴影部分表示一个距离门,反映同一个高度的信息,每个单元格表示的是滤波点的乘积,将这些乘积沿对角线相加便可以得到相应时延的自相关。为了消除传统方法下的不平衡性,v的值是可以变动的,主要是使得所有 lag 求解涉及到的滤波点数尽量接近,可以有效的提高信噪比,这里假设v为3,根据自相关函数的表达式即可以得到其自相关函数,进而求得功率谱。
如图3所示,交替码调制过程中,求出其分延时的自相关函数,将其合并,可以得到回波信号的自相关函数,并求得功率谱。
Claims (4)
1.一种处理非相干散射雷达散射信号的波形设计方法,其特征在于:在非相干散射雷达信号处理的基础上,结合巴克码,长脉冲编码以及交替码对回波信号进行处理,根据三者自相关函数表达式,利用MATlAB软件求得三种编码的自相关函数以及其功率谱。
2.根据权利要求1所述的一种处理非相干散射雷达散射信号的波形设计方法,其特征在于:所述巴克码是一种相位调制,是二相编码信号的一种,若一个n 位的巴克码{ },每个码元只可能取值+1或-1,则它满足条件:
。
3.根据权利要求1所述的一种处理非相干散射雷达散射信号的波形设计方法,其特征在于:所述长脉冲编码调制通常适用于探测电离层F层或低仰角和低空间分辩率探测,长脉冲的自相关函数可以表示为:
其中,i 表示第 i 个高度的距离门,n 是脉冲宽度与 lag 时延间隔。
4.根据权利要求1所述的一种处理非相干散射雷达散射信号的波形设计方法,其特征在于:所述交替码是将单个长脉冲划分为一系列小单元,所有的单元之间可以进行相关操作,且单元之间不影响,可分为强交替码和弱交替码。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410098372.2A CN103913723A (zh) | 2014-03-18 | 2014-03-18 | 一种处理非相干散射雷达散射信号的波形设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410098372.2A CN103913723A (zh) | 2014-03-18 | 2014-03-18 | 一种处理非相干散射雷达散射信号的波形设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103913723A true CN103913723A (zh) | 2014-07-09 |
Family
ID=51039549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410098372.2A Pending CN103913723A (zh) | 2014-03-18 | 2014-03-18 | 一种处理非相干散射雷达散射信号的波形设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103913723A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898101A (zh) * | 2015-05-11 | 2015-09-09 | 南昌大学 | 电离层非相干散射雷达探测威力仿真方法 |
CN105137400A (zh) * | 2015-09-06 | 2015-12-09 | 哈尔滨工业大学 | 一种瞬态极化雷达波形获得方法及基于该方法的雷达信号传输方法 |
CN105182305A (zh) * | 2015-06-04 | 2015-12-23 | 南昌大学 | 电离层非相干散射雷达仿真系统 |
CN108254730A (zh) * | 2018-04-11 | 2018-07-06 | 西安电子科技大学 | 基于阻尼拟合的雷达零时延处自相关函数处理方法 |
CN109782234A (zh) * | 2019-01-28 | 2019-05-21 | 南昌大学 | 用于雷达探测的复合编码方法及设备 |
CN110412536A (zh) * | 2019-08-15 | 2019-11-05 | 南京信息工程大学 | 基于变尺度中值比率的目标检测方法 |
CN115356720A (zh) * | 2022-07-20 | 2022-11-18 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种非合作电离层探测信号重塑方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010013839A1 (en) * | 1999-12-13 | 2001-08-16 | University Corporation For Atmospheric Research | Bistatic radar system for centralized, near-real-time synchronized, processing of data to identify scatterers |
CN103592646A (zh) * | 2013-11-01 | 2014-02-19 | 南昌大学 | 相控阵非相干散射雷达系统 |
-
2014
- 2014-03-18 CN CN201410098372.2A patent/CN103913723A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20010013839A1 (en) * | 1999-12-13 | 2001-08-16 | University Corporation For Atmospheric Research | Bistatic radar system for centralized, near-real-time synchronized, processing of data to identify scatterers |
CN103592646A (zh) * | 2013-11-01 | 2014-02-19 | 南昌大学 | 相控阵非相干散射雷达系统 |
Non-Patent Citations (3)
Title |
---|
张梁民: "非相干散射雷达谱分析与信号处理方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
杜鑫: "非相干散射雷达信号处理软件设计与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
王玮: "非相干散射雷达信号处理算法研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898101A (zh) * | 2015-05-11 | 2015-09-09 | 南昌大学 | 电离层非相干散射雷达探测威力仿真方法 |
CN105182305A (zh) * | 2015-06-04 | 2015-12-23 | 南昌大学 | 电离层非相干散射雷达仿真系统 |
CN105137400A (zh) * | 2015-09-06 | 2015-12-09 | 哈尔滨工业大学 | 一种瞬态极化雷达波形获得方法及基于该方法的雷达信号传输方法 |
CN108254730A (zh) * | 2018-04-11 | 2018-07-06 | 西安电子科技大学 | 基于阻尼拟合的雷达零时延处自相关函数处理方法 |
CN108254730B (zh) * | 2018-04-11 | 2021-07-20 | 西安电子科技大学 | 基于阻尼拟合的雷达零时延处自相关函数处理方法 |
CN109782234A (zh) * | 2019-01-28 | 2019-05-21 | 南昌大学 | 用于雷达探测的复合编码方法及设备 |
CN110412536A (zh) * | 2019-08-15 | 2019-11-05 | 南京信息工程大学 | 基于变尺度中值比率的目标检测方法 |
CN110412536B (zh) * | 2019-08-15 | 2022-06-24 | 南京信息工程大学 | 基于变尺度中值比率的目标检测方法 |
CN115356720A (zh) * | 2022-07-20 | 2022-11-18 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种非合作电离层探测信号重塑方法 |
CN115356720B (zh) * | 2022-07-20 | 2024-05-10 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种非合作电离层探测信号重塑方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103913723A (zh) | 一种处理非相干散射雷达散射信号的波形设计方法 | |
JP6177467B1 (ja) | ドップラー能力を高めたレーダー動作 | |
CN107683422B (zh) | 车辆雷达系统 | |
CN110531336A (zh) | 一种物体检测识别方法及系统 | |
Liu et al. | Observation of compact intracloud discharges using VHF broadband interferometers | |
CN103217673B (zh) | 一种非均匀Weibull杂波背景下的CFAR检测方法 | |
Shin et al. | Ultrasonic distance measurement method with crosstalk rejection at high measurement rate | |
JP2018205174A (ja) | レーダ装置及びそのレーダ信号処理方法 | |
RU2559165C1 (ru) | Устройство для определения направления и дальности до источника сигнала | |
US20190162812A1 (en) | System and method for determining the source location of a firearm discharge | |
CN111427018A (zh) | 一种雷达干扰装备干扰效果评估方法 | |
Wu et al. | A highly accurate ultrasonic ranging method based on onset extraction and phase shift detection | |
CN111474524A (zh) | 一种雷达干扰装备干扰效果监测与决策支持系统 | |
CN104833970A (zh) | 一种适用于存在遮挡效应的旋转目标稀疏多普勒成像方法 | |
CN104316160A (zh) | 基于小波脊的水下声信号瞬时频率解调方法 | |
CN105629224A (zh) | 调频连续波雷达高精度测距方法 | |
CN107346022B (zh) | 基于微波干涉仪的高精度测船雷达及测速方法 | |
CN104751675A (zh) | 基于有限新息率理论脉冲信号采样的车位检测方法 | |
Hyun et al. | Method to improve range and velocity error using de-interleaving and frequency interpolation for automotive FMCW radars | |
CN106066472B (zh) | 一种二维振速梯度水听器的被动目标相关检测方法 | |
Chen et al. | Properties of “site error” of lightning direction-finder (DF) and its modeling | |
RU2545068C1 (ru) | Способ измерения изменения курсового угла движения источника зондирующих сигналов | |
RU2688921C2 (ru) | Способ измерения дальности и радиальной скорости в РЛС с зондирующим составным псевдослучайным ЛЧМ импульсом | |
CN109581350A (zh) | 基于时频积分插值的雷达测距测速方法及装置 | |
CN110832341A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140709 |
|
WD01 | Invention patent application deemed withdrawn after publication |