CN104155648B - 基于阵列数据重排的高频地波雷达单次快拍music测向方法 - Google Patents
基于阵列数据重排的高频地波雷达单次快拍music测向方法 Download PDFInfo
- Publication number
- CN104155648B CN104155648B CN201410424616.1A CN201410424616A CN104155648B CN 104155648 B CN104155648 B CN 104155648B CN 201410424616 A CN201410424616 A CN 201410424616A CN 104155648 B CN104155648 B CN 104155648B
- Authority
- CN
- China
- Prior art keywords
- wave radar
- ground wave
- data
- array
- frequency
- 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.)
- Expired - Fee Related
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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
-
- 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
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/74—Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
-
- 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
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/0218—Very long range radars, e.g. surface wave radar, over-the-horizon or ionospheric propagation systems
- G01S2013/0227—OTH, Over-The-Horizon radar
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法,包括对高频地波雷达的回波信号进行二维FFT处理后,利用CFAR检测技术,检测出上述待侧向目标所在的距离单元、多普勒单元,得到阵列频域一次快拍数据,将阵列接收的一次快拍数据重新排列成Toeplitz矩阵,进而估计协方差矩阵,对协方差矩阵进行特征值分解,得到噪声子空间UN,然后利用MUSIC算法构造空间谱方式,进行目标方位估计,克服了降维MUSIC算法构造协方差矩阵引起的阵列天线孔径损失问题。实测数据表明,基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法具有较高的方位估计精度和分辨力。
Description
技术领域
本发明涉及一种高频地波雷达单次快拍MUSIC测向方法,具体涉及一种基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法。
背景技术
高频地波雷达,又称为高频表面波雷达,工作于高频(3~30MHz)频段,利用垂直极化高频电磁波沿海面(地表面)绕射传播特性,实现对海面舰船目标和低空飞行目标的超视距探测。高频地波雷达由于自身体制的限制,发射信号波长较长,会导致其对目标方位信息不敏感,方位向估计的精度较低,而目标的方位信息对于后期的数据处理又有十分重要的作用,因而对于目标方位的估计具有非常现实的意义。为保证***得到较高的信噪比与良好的角度分辨力,一般需要使用大型接收阵列天线,而在实际工程实践中会带来一系列的问题和麻烦,因此小型化阵列天线备受亲睐。然而,小型化阵列孔径宽度有限,常规波束形成方法经常会造成目标不能正常分辨,测角误差较大。因此,为了改善测角型能,Schmidt(1986)提出的多重信号分类(MUSIC)算法、Paulraj等(1987)提出旋转不变子空间(ESPRIT)算法、Barabell(1989)提出Root-MUSIC算法等超分辨方法应用到高频地波雷达中。
利用MUSIC算法进行DOA估计时,其性能与快拍数有关,快拍数越多,DOA估计精度越高,但大量的快拍积累给实际的工程应用带来了不便:一是它无法对快速时变信号进行测向;二是较多的快拍积累不利于实时计算。对于一次快拍而言,如果用XXH来估计协方差矩阵,估计的效果会非常差。协方差矩阵估计的不准,特征分解时得到的噪声矢量与阵列导向矢量的正交性就会降低,导致空间谱估计不准。
针对这个问题,位寅生等提出了一种降维单次快拍MUSIC算法,这种方法存在阵元损失,致使天线自由度降低的问题,尤其在阵列数比较小的情况下,阵元损失会使得角度估计性能降低。相关的参考文献如下:
[1]R.O.Schmidt.Multiple Emitter Location and Signal ParameterEstimation[J].IEEE Transaction on Antenna and Propagation.1986,34(3):276-280.
[2]B.D.Rao and K.V.S.Hari.Performance analysis of root-music[J].IEEETrans.Aacoust,Speech,Signal Process,1989,37(12):1939-1949.
[3]WEI Yin-sheng,TONG Peng,GUO Xiao-jiang.Single snapshotsuperresolution algorithm for HFSWR based on noise eigenvector reconstruction
[4]Q.S.Ren,A.J.Wlllis.Extending MUSIC to signal snapshot and on linedirection finding applications[J].Proc of IEEE Radar,Edinburgh,UK,1997:783-787.
发明内容
本发明的目的是提供一种基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法,该方法将阵列接收的一次快拍数据重新排列成Toeplitz矩阵,进而估计协方差矩阵,提高了方位估计精度高,而且有更高的分辨力。
为了实现上述目的,本发明的技术方案如下:
利用高频地波雷达朝海面发射无线电磁波对海面舰船目标的超视距探测,在中波和短波段海水表面的地波传播衰减很小,而且地波在一定程度上会沿着弯曲的地球表面传播,到达地平线以下很远的地方,即实现超视距传播,通过目标对高频无线电波的反射实现对数百公里内的目标进行探测,其特征在于还包括以下步骤:
步骤1:对高频地波雷达的回波信号进行二维FFT处理后,利用CFAR检测技术,检测出上述待侧向目标所在的距离单元、多普勒单元,得到阵列频域一次快拍数据,其结构为xi,i为1至M的自然数,M为地波雷达天线阵元数;
步骤2:利用步骤1得到的一次快拍数据进行阵列数据重排,排列成Toeplize矩阵X;
其中,*为共轭,M为地波雷达天线阵元数(行数);
步骤3:利用步骤2得到的矩阵X,估计协方差矩阵Rx;
RX=X*XH
其中H表示共轭转置;
步骤4:利用步骤3得到的协方差矩阵Rx进行特征值分解,得到噪声子空间UN,然后利用MUSIC算法构造空间谱方式,进行待侧向目标方位估计;其空间谱函数表达式为:
其中a(θ)为方向导向矢量,H表示共轭转置,UN为噪声子空间;
步骤5:利用步骤4构造的空间谱函数进行谱峰搜索,得到目标的方位估计。
上述步骤1:首先对高频地波雷达的回波信号进行二维FFT处理后,将时域信号转换到频域信号,得到距离-多普勒谱图,即纵坐标为距离、横坐标为多普勒的二维数据;利用CFAR检测技术,检测出待测向目标所在的距离单元、多普勒单元,将待测向目标所在距离单元与多普勒单元的数据提取出来,即可得到阵列频域一次快拍数据;
上述步骤4:利用步骤3得到的协方差矩阵Rx进行特征值分解,可得到N(信号个数,一个距离单元及其多普勒单元内一般为一个信号)个大特征值和M-N个小特征值,它们对应的特征向量分别为u1,…uN,uN+1,uM,则UN=[uN+1,…,uM]的各列张成噪声子空间,然后利用MUSIC算法构造空间谱方式,进行目标方位估计,其空间谱表达式为:
其中a(θ)为方向导向矢量,H表示共轭转置,UN为噪声子空间。
与现有技术相比,本发明的创新之处体现在以下方面:
MUSIC算法是高频地波雷达目标超分辨测向方法中经典的方法,其方位估计性能与雷达阵列接收数据的快拍数有关,快拍数越多,方位估计精度越高。但大量的快拍积累使得算法计算量大,不利于实时估计目标方位,给实际的工程应用带来不便。在实际应用中,高频地波雷达通常采用频域一次快拍数据进行目标方位估计,利用降维MUSIC算法估计协方差矩阵,但这种方法相当于减小了阵列天线的有效孔径,降低了阵列天线的自由度,导致目标方位估计精度下降。
针对这一问题,本发明提出了一种基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法,该方法将阵列接收的一次快拍数据重新排列成Toeplitz矩阵,进而估计协方差矩阵,克服了降维MUSIC算法构造协方差矩阵引起的阵列天线孔径损失问题。实测数据表明,基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法具有较高的方位估计精度和分辨力。
附图说明
图1为本发明的基本流程示意图。
图2为距离多普勒(R-D谱)图。
图3为本发明的单目标方位估计图。
图4为本发明的目标方位估计误差统计图。
图5为本发明的多目标方位估计图。
具体实施方式
下面结合附图,对本发明的方法做进一步说明:
利用高频地波雷达朝海面发射无线电磁波对海面舰船目标的超视距探测,在中波和短波段海水表面的地波传播衰减很小,而且地波在一定程度上会沿着弯曲的地球表面传播,到达地平线以下很远的地方,即实现超视距传播,通过目标对高频无线电波的反射实现对数百公里内的目标进行探测。
以雷达天线阵元数为8元等间距直线阵,频率为4.7MHz,积累时间为291s,目标方位以阵列法向为0°,顺时针为正的高频地波雷达为例。如图1所示,其具体步骤如下:
步骤1:对高频地波雷达的回波信号进行二维FFT处理后,利用CFAR检测技术,检测出上述待测向目标所在的距离单元、多普勒单元,如图2所示;得到阵列频域一次快拍数据,其结构为xi,i为1至M的自然数,M为高频地波雷达天线阵元数;
步骤2:利用步骤1得到的一次快拍数据进行阵列数据重排,排列成Toeplize矩阵X;
其中,*为共轭,M为阵元数(行数),为高频地波雷达天线阵元数;
步骤3:利用步骤2得到的矩阵X,估计协方差矩阵Rx;
RX=X*XH
其中H表示共轭转置;
步骤4:利用步骤3得到的协方差矩阵Rx进行特征值分解,得到噪声子空间UN,然后利用MUSIC算法构造空间谱方式,进行待测向目标方位估计;其空间谱函数表达式为:
其中a(θ)为方向导向矢量,H表示共轭转置,UN为噪声子空间;
步骤5:利用步骤4构造的空间谱函数进行谱峰搜索,得到目标的方位估计,如图3所示。然后统计多个单目标的方位估计偏差结果,如图4所示,可见本方法的偏差较小,因此具有较高的测向精度。
对于待测向多目标的情况,以两个待侧向目标为例,两个待侧向目标信号具有大致相同的距离单元和多普勒单元,即具有大致相同的距离和速度,但方位与原目标方位有一定的间隔,重复步骤2-步骤5,得到多目标方位估计结果。如图5所示可见,本方法对具有基本相同的距离和速度的目标也具有很好的分辨效果。
实施例
上述步骤1:对高频地波雷达的回波信号进行二维FFT处理后,利用CFAR检测技术,检测出上述待测向目标所在的距离单元为51、多普勒单元为97,得到阵列频域一次快拍数据为
上述步骤2:利用步骤1得到的一次快拍数据进行阵列数据重排,排列成Toeplize矩阵X;
上述步骤3:利用步骤2得到的矩阵X,估计协方差矩阵Rx;
上述步骤4:利用步骤3得到的协方差矩阵Rx进行特征值分解,可得到1个大特征值和7个小特征值,噪声子空间为UN,然后利用MUSIC算法构造空间谱方式,进行目标方位估计。
步骤5:利用步骤4构造的空间谱函数进行谱峰搜索,
得到max(Pmusic)对应的角度为-7.7°。
Claims (3)
1.基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法,包括:
利用高频地波雷达朝海面发射无线电磁波对海面舰船目标的超视距探测,通过目标对高频无线电波的反射实现对数百公里内的目标进行探测,其特征在于还包括以下步骤:
步骤1:对高频地波雷达的回波信号进行二维FFT处理后,利用CFAR检测技术,检测出上述待测向目标所在的距离单元、多普勒单元,得到阵列频域一次快拍数据,其结构为xi,i为1至M的自然数,M为地波雷达天线阵元数;
步骤2:利用步骤1得到的一次快拍数据进行阵列数据重排,排列成Toeplize矩阵X;
其中,*为共轭,M为地波雷达天线阵元数;
步骤3:利用步骤2得到的矩阵X,估计协方差矩阵Rx;
RX=X*XH
其中H表示共轭转置;
步骤4:利用步骤3得到的协方差矩阵Rx进行特征值分解,得到噪声子空间UN,然后利用MUSIC算法构造空间谱方式,进行待测向目标方位估计;其空间谱函数表达式为:
其中a(θ)为方向导向矢量,H表示共轭转置,UN为噪声子空间;
步骤5:利用步骤4构造的空间谱函数进行谱峰搜索,得到目标的方位估计。
2.如权利要求1所述的基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法,其特征在于上述步骤1:首先对高频地波雷达的回波信号进行二维FFT处理后,将时域信号转换到频域信号,得到距离-多普勒谱图,是纵坐标为距离、横坐标为多普勒频率的二维数据;利用CFAR检测技术,检测出待测向目标所在的距离单元、多普勒单元,将待测向目标所在距离单元与多普勒单元的数据提取,即可得到阵列频域一次快拍数据。
3.如权利要求1所述的基于阵列数据重排的高频地波雷达单次快拍MUSIC测向方法,其特征在于上述步骤4:利用步骤3得到的协方差矩阵Rx进行特征值分解,可得到N个大特征值和M-N个小特征值,它们对应的特征向量分别为u1,…uN,uN+1,uM,则UN=[uN+1,…,uM]的各列张成噪声子空间,然后利用MUSIC算法构造空间谱方式,进行目标方位估计,其空间谱表达式为:
其中a(θ)为方向导向矢量,H表示共轭转置,UN为噪声子空间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410424616.1A CN104155648B (zh) | 2014-08-26 | 2014-08-26 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410424616.1A CN104155648B (zh) | 2014-08-26 | 2014-08-26 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104155648A CN104155648A (zh) | 2014-11-19 |
CN104155648B true CN104155648B (zh) | 2016-08-17 |
Family
ID=51881198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410424616.1A Expired - Fee Related CN104155648B (zh) | 2014-08-26 | 2014-08-26 | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104155648B (zh) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104698433B (zh) * | 2015-03-16 | 2017-08-29 | 电子科技大学 | 基于单快拍数据的相干信号doa估计方法 |
CN105182283B (zh) * | 2015-08-17 | 2017-08-25 | 周口师范学院 | 一种无源雷达固定目标时域测向方法 |
CN105182285A (zh) * | 2015-10-14 | 2015-12-23 | 中国电子科技集团公司第二十八研究所 | 一种基于声矢量二维嵌套阵列的目标测向方法 |
CN105372635A (zh) * | 2015-11-19 | 2016-03-02 | 哈尔滨工业大学 | 基于改进的降维空时自适应处理的船载高频地波雷达海杂波抑制方法 |
CN109283536B (zh) * | 2018-09-01 | 2023-02-14 | 哈尔滨工程大学 | 一种多波束测深声呐水体成像波束形成方法 |
CN109212526B (zh) * | 2018-10-17 | 2022-07-01 | 哈尔滨工业大学 | 用于高频地波雷达的分布式阵列目标角度测量方法 |
CN109856627A (zh) * | 2019-01-11 | 2019-06-07 | 中国船舶重工集团公司第七二四研究所 | 基于空间补零插值的三角阵相控阵2d-fft测向方法 |
CN110334638B (zh) * | 2019-06-28 | 2021-09-10 | 西安理工大学 | 基于快速music算法的道路双黄线检测方法 |
CN110310667A (zh) * | 2019-07-12 | 2019-10-08 | 莫毓昌 | 一种基于声学特征的设备故障检测方法 |
CN113391257B (zh) * | 2020-03-13 | 2023-04-18 | 光宝科技新加坡私人有限公司 | 用于物件的角度估测的运算装置及物件的角度估测方法 |
CN111521968B (zh) * | 2020-05-22 | 2022-05-20 | 南京理工大学 | 基于目标空间分集的欠定doa估计方法 |
WO2021243529A1 (zh) * | 2020-06-01 | 2021-12-09 | 华为技术有限公司 | 一种波达角aoa估计方法和装置 |
CN112346005B (zh) * | 2020-10-30 | 2022-07-12 | 哈尔滨工程大学 | 一种应用于均匀圆水听器阵的空域旋转方位估计方法 |
CN113203997B (zh) * | 2021-04-15 | 2023-11-03 | 杭州加速科技有限公司 | 基于fpga的雷达超分辨测向方法、***及应用 |
CN114442084B (zh) * | 2022-01-21 | 2024-06-14 | 哈尔滨工业大学(威海) | 便携式高频地波雷达阵列、雷达***性能分析方法及雷达 |
CN115061109A (zh) * | 2022-06-09 | 2022-09-16 | 山东大学 | 一种基于被动雷达的低慢小多目标检测方法及*** |
CN116256716B (zh) * | 2023-01-18 | 2023-11-24 | 珠海微度芯创科技有限责任公司 | 基于加权线性拟合的毫米波雷达超分辨方法和*** |
CN116819429B (zh) * | 2023-06-30 | 2024-02-27 | 中国人民解放军军事科学院***工程研究院 | 一种联合时域快拍检测的弱信号测向方法 |
-
2014
- 2014-08-26 CN CN201410424616.1A patent/CN104155648B/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
MUSIC空间谱估计并行运算算法;位寅生等;《***工程与电子技术》;20120301;第34卷(第1期);12-16 * |
一种不考虑信源数目的修正MUSIC算法;张志军等;《电讯技术》;20071030;第47卷(第5期);58-62 * |
三角形Toeplize矩阵的三角本原指数;张景晓等;《数学的实践与认识》;20040131;第34卷(第1期);139-144 * |
基于四阶累积量的被动测向声纳高分辨力算法研究;王良等;《中国海洋大学学报》;20070331;第37卷(第2期);211-214 * |
Also Published As
Publication number | Publication date |
---|---|
CN104155648A (zh) | 2014-11-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104155648B (zh) | 基于阵列数据重排的高频地波雷达单次快拍music测向方法 | |
US12032091B2 (en) | Signal detection and denoising systems | |
Wang et al. | Conjugate ESPRIT for DOA estimation in monostatic MIMO radar | |
CN103616661B (zh) | 一种稳健的远场窄带信号源个数估计方法 | |
CN103353596A (zh) | 基于压缩感知的波束空间域米波雷达测高方法 | |
CN103926573B (zh) | 基于四阶累积量的单基地mimo雷达分布式目标角度估计方法 | |
CN104155632A (zh) | 一种基于局部相关性的改进子空间海杂波抑制方法 | |
CN109633558A (zh) | 一种基于极化时频分布的波达方向估计算法 | |
CN102135617A (zh) | 双基地多输入多输出雷达多目标定位方法 | |
CN103323827B (zh) | 基于快速傅里叶变换的mimo雷达***角度估计方法 | |
CN102393525A (zh) | 子空间投影的导航干扰抑制与信号增强方法 | |
CN104749554A (zh) | 一种基于递归秩损的幅相误差校准和波达方向估计方法 | |
CN109254272B (zh) | 一种共点式极化mimo雷达的两维角度估计方法 | |
CN101644760B (zh) | 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法 | |
CN109375154A (zh) | 一种冲击噪声环境下基于均匀圆阵的相干信号参数估计方法 | |
CN103885049B (zh) | 基于最小冗余线性稀疏子阵的米波雷达低仰角估计方法 | |
CN113032721B (zh) | 一种低计算复杂度的远场和近场混合信号源参数估计方法 | |
CN105403871A (zh) | 一种双基地mimo雷达阵列目标角度估计和互耦误差校准方法 | |
CN109765521A (zh) | 一种基于子阵划分的波束域成像方法 | |
CN105572630B (zh) | 基于多波位联合处理的单脉冲目标doa估计方法 | |
CN107290732A (zh) | 一种量子大***的单基地mimo雷达测向方法 | |
CN104459713A (zh) | 卫星导航接收机对欺骗式干扰来波方向的估计方法 | |
CN114779236A (zh) | 一种基于空间平滑music改进的米波雷达低仰角测高方法 | |
CN105572642B (zh) | 一种基于两级架构的空时自适应处理方法 | |
Dakulagi | A new Nystrom approximation based efficient coherent DOA estimator for radar applications |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160817 Termination date: 20170826 |