CN105871440A - 混合载波多天线分量传输的信号接收方法 - Google Patents

混合载波多天线分量传输的信号接收方法 Download PDF

Info

Publication number
CN105871440A
CN105871440A CN201610428749.5A CN201610428749A CN105871440A CN 105871440 A CN105871440 A CN 105871440A CN 201610428749 A CN201610428749 A CN 201610428749A CN 105871440 A CN105871440 A CN 105871440A
Authority
CN
China
Prior art keywords
signal
mixed carrier
alpha
multiple antennas
components
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
Application number
CN201610428749.5A
Other languages
English (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201610428749.5A priority Critical patent/CN105871440A/zh
Publication of CN105871440A publication Critical patent/CN105871440A/zh
Pending legal-status Critical Current

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/0413MIMO systems
    • H04B7/0417Feedback systems
    • H04B7/0421Feedback systems utilizing implicit feedback, e.g. steered pilot signals
    • 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/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/0848Joint weighting
    • H04B7/0854Joint weighting using error minimizing algorithms, e.g. minimum mean squared error [MMSE], "cross-correlation" or matrix inversion
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00Arrangements affording multiple use of the transmission path
    • H04L5/003Arrangements for allocating sub-channels of the transmission path
    • H04L5/0048Allocation of pilot signals, i.e. of signals known to the receiver

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Radio Transmission System (AREA)

Abstract

混合载波多天线分量传输的信号接收方法,涉及无线通信领域。本发明是为了实现混合载波信号在接收端难以进行正确恢复的问题。其接收方法:采用接收天线接收发射的信号,并利用不同的导频序列分别估计出不同子信道的信道状态信息;根据变换阶数α和信号的组合方式,在接收端生成发射信号变换矩阵;根据获得的信道状态信息和得到的发射信号变换矩阵,计算接收信号生成矩阵;利用接收到的信号和收信号生成矩阵,采用最小均方误差(MMSE)方法恢复原始信号。本发明适用于无线通信过程中。

Description

混合载波多天线分量传输的信号接收方法
技术领域
本发明涉及无线通信领域。
背景技术
随着科技的进步,移动通信的信道环境变得越来越复杂,在许多新的通信场景下,信道所呈现出的衰落特性已经不止是单一的时域选择性衰落或频域选择性衰落,而是一种时域与频域同时衰落的时频双选特性。在这种信道下,无论是单载波体制所采用的时域信号,还是多载波体制采用的频域信号都难以获得较好的性能,而混合载波体制所采用的时频协同信号却能够获得较好的效果。当时频协同信号应用于多天线场景时,由于不同分量对不同衰落特性的信道具有不同的适应性,因而,可以将混合载波信号的不同分量映射到与之匹配的子信道上,从而获得更好的通信效果。但由于不同子信道之间的衰落是随机的,导致混合载波信号不同分量之间的数学约束关系遭到破坏,如果在接收端直接进行反变换将难以获得正确的信号。
发明内容
本发明是为了实现混合载波信号在接收端难以进行正确恢复的问题,从而提供一种混合载波多天线分量传输系统接收方法。
混合载波多天线分量传输的信号接收方法,
混合载波多天线分量传输的信号发射方法,它包括以下步骤:
步骤一、将待发送信号进行调制,并分成长度为N的数据块;N为正数;
步骤二、将步骤一所述的长度为N的数据块进行α阶加权分数傅里叶变换;α为正数;
步骤三、将步骤二变换后不同的信号分量或分量的组合形式映射到不同的发射天线上;
步骤四、将步骤三中映射到不同天线上的信号插入不同的导频序列,并送入信道进行传输;
其特征是:混合载波多天线分量传输的信号接收方法,它包括以下步骤:
步骤五、采用接收天线接收步骤四发射的信号,并利用不同的导频序列分别估计出不同子信道的信道状态信息;
步骤六、根据变换阶数α和信号的组合方式,在接收端生成发射信号变换矩阵;
步骤七、根据步骤五中获得的信道状态信息和步骤六中得到的发射信号变换矩阵,计算接收信号生成矩阵;
步骤八、利用接收到的信号和收信号生成矩阵,采用最小均方误差(MMSE)方法恢复原始信号。
本发明通过结合信道状态信息与加权分数傅里叶变换的联合逆变换方法为混合载波多天线分量传输系统提供了一种可行的接收方案,使得接收端能够正确恢复发送端所发送的信号。在各个子信道衰落特性不同的情况下,获得更好的系统性能。
附图说明
图1是混合载波多天线分量传输系统示意图;
图2是混合载波多天线分量传输系统接收机结构示意图;
具体实施方式
具体实施方式一、结合图1和图2说明本具体实施方式,混合载波多天线分量传输的信号接收方法,该接收方法主要通过以下步骤实现:
步骤一、在发送端,待发送信号经过调制后,分成长度为N的数据块;
步骤二、将长度为N的数据块送入加权分数傅里叶变换模块,根据选择的变换阶数α进行变换。如果需要进行变换的原始数据为X0,则加权分数傅里叶变换的矩阵表示形式如下:
FαX0=(w0(α)I+w1(α)F+w2(α)PI+w3(α)PF)X0
其中:{w0(α),w1(α),w2(α),w3(α)}为加权系数,其具体表达式为:
w l ( α ) = c o s [ ( α - l ) π 4 ] c o s [ ( α - l ) π 2 ] exp [ 3 j ( α - l ) π 4 ] , ( l = 0 , 1 , 2 , 3 )
I为单位阵,而F为离散傅里叶变换矩阵,其具体的归一化表达式为:
其中q=e-j2π/N
P为置换矩阵,其具体形式为:
P = 1 0 0 ... 0 0 0 0 0 0 ... 0 0 1 0 0 0 ... 0 1 0 · · · · · · · · · · · · · · · · · · 0 0 1 ... 0 0 0 0 1 0 ... 0 0 0 .
步骤三、将不同的信号分量或分量的组合形式映射到不同的发射天线上。
以两发一收的多天线系统为例,可将混合载波信号中的时域分量与时域反转分量合并,映射到多普勒较大的子信道上,同时将频域分量与频域反转分量合并,映射到多径较大的信道上,即:
X 1 = ( w 0 ( α ) I + w 2 ( α ) P I ) X 0 X 2 = ( w 1 ( α ) F + w 3 ( α ) P F ) X 0
步骤四、将映射到不同天线上的信号插入不同的导频序列,并送入信道进行传输。
步骤五、接收天线接收到信号后,将信号分为两路,一路用于最终的逆变换,另一路则送入信道估计器,利用不同的导频序列分别估计出不同子信道的信道状态信息。
以两发一收的多天线系统为例,经过信道估计后,可以得到两条子信道的信道增益矩阵:
H1=[h1(1) h1(2) … h1(n)]
H2=[h2(1) h2(2) … h2(n)]
其中,h1(n)和h2(n)分别为子信道1和子信道2在第n时刻的信道增益系数,n∈[1,N]。
步骤六、根据变换阶数α和信号的组合方式,在接收端生成发射信号变换矩阵。假设发送端信号按照步骤三中的映射规则进行映射则其发射信号变换矩阵为:w0(α)I+w2(α)PI和w1(α)F+w3(α)PF。
事实上由于发射信号的变换矩阵只与变换阶数α和信号的组合方式有关,因而当收发两端事先约定好变换阶数α和信号的组合方式时,发射信号变换矩阵可提前计算好后存储在接收端备用。
步骤七、根据步骤五中获得的信道状态信息和步骤六中得到的发射信号变换矩阵计算接收信号生成矩阵,其表达式为:
G=H1(w0(α)I+w2(α)PI)+H2(w1(α)F+w3(α)PF)
步骤八、利用接收到的信号Y和收信号生成矩阵G,MMSE联合逆变换器依据步骤七计算出的接收信号生成矩阵,利用公式:
X ~ M M S E = ( G H G + σ 2 I ) - 1 G H Y
计算出原始信号的估计值,而后通过解星座映射和信号判决从而得到最终结果;式中σ2为噪声功率。
本发明通过一种结合信道状态信息与加权分数傅里叶变换的联合逆变换方法为混合载波多天线分量传输系统提供了一种可行的接收方案,使得接收端能够正确恢复发送端所发送的信号。在各个子信道衰落特性不同的情况下,获得更好的系统性能。
具体实施方式二、本具体实施方式与具体实施方式一所述的混合载波多天线分量传输的信号接收方法的区别在于,还包括4发1收的情况。其特征在于,步骤三中产生的四个分量分别映射到4个发射天线上,即:
X 1 = w 0 ( α ) I X 0 X 2 = w 1 ( α ) FX 0 X 3 = w 2 ( α ) PIX 0 X 4 = w 3 ( α ) PFX 0
步骤五中需要估计出的子信道的信道状态信息由两个变为四个,可分别表示为H1=[h1(1) h1(2) … h1(n)]、H2=[h2(1) h2(2) … h2(n)]、H3=[h3(1) h3(2) … h3(n)]和H4=[h4(1) h4(2) … h4(n)]。
步骤六中的发射信号变换矩阵也由两个变为4个,分别是w0(α)I、w1(α)F、w2(α)PI和w3(α)PF。步骤七中的接收信号生成矩阵表达式变为:G=H1w0(α)I+H2w1(α)F+H3w2(α)PI+H4w3(α)PF。其他步骤不变。

Claims (10)

1.混合载波多天线分量传输的信号接收方法,
混合载波多天线分量传输的信号发射方法,它包括以下步骤:
步骤一、将待发送信号进行调制,并分成长度为N的数据块;N为正数;
步骤二、将步骤一所述的长度为N的数据块进行α阶加权分数傅里叶变换;α为正数;
步骤三、将步骤二变换后不同的信号分量或分量的组合形式映射到不同的发射天线上;
步骤四、将步骤三中映射到不同天线上的信号插入不同的导频序列,并送入信道进行传输;
其特征是:混合载波多天线分量传输的信号接收方法,它包括以下步骤:
步骤五、采用接收天线接收步骤四发射的信号,并利用不同的导频序列分别估计出不同子信道的信道状态信息;
步骤六、根据变换阶数α和信号的组合方式,在接收端生成发射信号变换矩阵;
步骤七、根据步骤五中获得的信道状态信息和步骤六中得到的发射信号变换矩阵,计算接收信号生成矩阵;
步骤八、利用接收到的信号和收信号生成矩阵,采用最小均方误差(MMSE)方法恢复出原始信号。
2.根据权利要求1所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤二所述将步骤一所述的长度为N的数据块进行α阶加权分数傅里叶变换,具体为:
设长度为N的数据块为X0,则加权分数傅里叶变换的矩阵表示形式如下:
FαX0=(w0(α)I+w1(α)F+w2(α)PI+w3(α)PF)X0
式中:{w0(α),w1(α),w2(α),w3(α)}为加权系数,具体表达式为:
w l ( α ) = c o s [ ( α - l ) π 4 ] c o s [ ( α - l ) π 2 ] exp [ 3 j ( α - l ) π 4 ] , ( l = 0 , 1 , 2 , 3 )
I为单位阵,而F为离散傅里叶变换矩阵,具体的归一化表达式为:
F = 1 N 1 1 ... 1 1 q ... q N - 1 . . . . . . ... . . . 1 q N - 1 ... q ( N - 1 ) ( N - 1 )
其中:q=e-j2π/N
P为置换矩阵,具体形式为:
P = 1 0 0 ... 0 0 0 0 0 0 ... 0 0 1 0 0 0 ... 0 1 0 . . . . . . . . . . . . . . . . . . 0 0 1 ... 0 0 0 0 1 0 ... 0 0 0 .
3.根据权利要求1所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤三中,将步骤二变换后不同的信号分量或分量的组合形式映射到不同的发射天线上,具体方法为:
对于两发一收的多天线系统,将混合载波信号中的时域分量与时域反转分量合并,映射到多普勒大于预设阈值的子信道上,同时将频域分量与频域反转分量合并,映射到多径大于预设值的信道上,即:
X 1 = ( w 0 ( α ) I + w 2 ( α ) P I ) X 0 X 2 = ( w 1 ( α ) F + w 3 ( α ) P F ) X 0 .
4.根据权利要求3所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤五中,采用接收天线接收步骤四发射的信号,并利用不同的导频序列分别估计出不同子信道的信道状态信息,具体为:
经过信道估计后,得到两条子信道的信道增益矩阵:
H1=[h1(1) h1(2) … h1(n)]
H2=[h2(1) h2(2) … h2(n)]
式中,h1(n)和h2(n)分别为子信道1和子信道2在第n时刻的信道增益系数,n∈[1,N]。
5.根据权利要求4所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤六中,在接收端生成发射信号变换矩阵,分别为:w0(α)I+w2(α)PI和w1(α)F+w3(α)PF。
6.根据权利要求5所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤七中,计算接收信号生成矩阵,结果为:
G=H1(w0(α)I+w2(α)PI)+H2(w1(α)F+w3(α)PF)。
7.根据权利要求5所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤八中,采用最小均方误差(MMSE)方法恢复原始信号,具体表达式为:
X ~ M M S E = ( G H G + σ 2 I ) - 1 G H Y
式中:Y为接收到的信号;G为接收信号生成矩阵;σ2为噪声功率。
8.根据权利要求1所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤三中,将步骤二变换后不同的信号分量或分量的组合形式映射到不同的发射天线上,具体方法为:
对于四发一收的多天线系统,将产生的四个分量分别映射到4个发射天线上,即:
X 1 = w 0 ( α ) IX 0 X 2 = w 1 ( α ) FX 0 X 3 = w 2 ( α ) PIX 0 X 4 = w 3 ( α ) PFX 0 .
9.根据权利要求8所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤五中,采用接收天线接收步骤四发射的信号,并利用不同的导频序列分别估计出不同子信道的信道状态信息,具体为:
经过信道估计后,得到四条子信道的信道增益矩阵:
分别表示为:
H1=[h1(1) h1(2) … h1(n)]
H2=[h2(1) h2(2) … h2(n)]
H3=[h3(1) h3(2) … h3(n)]
H4=[h4(1) h4(2) … h4(n)]。
10.根据权利要求9所述的混合载波多天线分量传输的信号接收方法,其特征在于步骤六中,在接收端生成发射信号变换矩阵,分别为:
w0(α)I、w1(α)F、w2(α)PI和w3(α)PF;
步骤七中,计算接收信号生成矩阵,结果为:
G=H1w0(α)I+H2w1(α)F+H3w2(α)PI+H4w3(α)PF。
CN201610428749.5A 2016-06-15 2016-06-15 混合载波多天线分量传输的信号接收方法 Pending CN105871440A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610428749.5A CN105871440A (zh) 2016-06-15 2016-06-15 混合载波多天线分量传输的信号接收方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610428749.5A CN105871440A (zh) 2016-06-15 2016-06-15 混合载波多天线分量传输的信号接收方法

Publications (1)

Publication Number Publication Date
CN105871440A true CN105871440A (zh) 2016-08-17

Family

ID=56649765

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610428749.5A Pending CN105871440A (zh) 2016-06-15 2016-06-15 混合载波多天线分量传输的信号接收方法

Country Status (1)

Country Link
CN (1) CN105871440A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017054632A1 (zh) * 2015-09-29 2017-04-06 华为技术有限公司 信号发送方法、信号接收方法、发射端及接收端
CN106656427A (zh) * 2017-02-24 2017-05-10 哈尔滨工业大学 基于奇偶对称序列的双天线信号传输方法
CN106953825A (zh) * 2017-03-21 2017-07-14 哈尔滨工业大学 基于加权分数傅立叶变换的多用户多址通信发射和接收方法
CN108924077A (zh) * 2018-06-25 2018-11-30 哈尔滨工业大学 广义混合载波系统时间选择性衰落信道下的传输方法
CN112332930A (zh) * 2020-11-06 2021-02-05 西北工业大学 一种移动条件下的时间反转方法
CN116015372A (zh) * 2022-12-29 2023-04-25 国家工业信息安全发展研究中心 大规模mimo数字信息传输方法、系统、设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101651649A (zh) * 2009-09-14 2010-02-17 哈尔滨工业大学 基于四项加权分数傅立叶变换与分布式发送天线的载波协同通信方法
WO2015184875A1 (zh) * 2014-06-06 2015-12-10 华为技术有限公司 一种双选信道的补偿方法、系统及相关装置
CN105262528A (zh) * 2015-09-18 2016-01-20 哈尔滨工业大学 基于加权分数傅立叶变换域的4天线发射分集方法
CN105322991A (zh) * 2015-10-15 2016-02-10 哈尔滨工业大学 基于wfrft预编码的多输入多输出传输系统及传输方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101651649A (zh) * 2009-09-14 2010-02-17 哈尔滨工业大学 基于四项加权分数傅立叶变换与分布式发送天线的载波协同通信方法
WO2015184875A1 (zh) * 2014-06-06 2015-12-10 华为技术有限公司 一种双选信道的补偿方法、系统及相关装置
CN105262528A (zh) * 2015-09-18 2016-01-20 哈尔滨工业大学 基于加权分数傅立叶变换域的4天线发射分集方法
CN105322991A (zh) * 2015-10-15 2016-02-10 哈尔滨工业大学 基于wfrft预编码的多输入多输出传输系统及传输方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KUN WANG等: ""Performance Analysis of Hybrid Carrier System with MMSE Equalization over Doubly-Dispersive Channels"", 《IEEE COMMUNICATION LETTERS》 *
吴玉龙: ""基于加权分数傅立叶变换的双弥散信道内信号波形设计"", 《硕士学位论文》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017054632A1 (zh) * 2015-09-29 2017-04-06 华为技术有限公司 信号发送方法、信号接收方法、发射端及接收端
CN106656427A (zh) * 2017-02-24 2017-05-10 哈尔滨工业大学 基于奇偶对称序列的双天线信号传输方法
CN106953825A (zh) * 2017-03-21 2017-07-14 哈尔滨工业大学 基于加权分数傅立叶变换的多用户多址通信发射和接收方法
CN108924077A (zh) * 2018-06-25 2018-11-30 哈尔滨工业大学 广义混合载波系统时间选择性衰落信道下的传输方法
CN108924077B (zh) * 2018-06-25 2020-10-30 哈尔滨工业大学 广义混合载波系统时间选择性衰落信道下的传输方法
CN112332930A (zh) * 2020-11-06 2021-02-05 西北工业大学 一种移动条件下的时间反转方法
CN112332930B (zh) * 2020-11-06 2021-09-07 西北工业大学 一种移动条件下的时间反转方法
CN116015372A (zh) * 2022-12-29 2023-04-25 国家工业信息安全发展研究中心 大规模mimo数字信息传输方法、系统、设备及介质
CN116015372B (zh) * 2022-12-29 2023-09-12 国家工业信息安全发展研究中心 大规模mimo数字信息传输方法、系统、设备及介质

Similar Documents

Publication Publication Date Title
CN105871440A (zh) 混合载波多天线分量传输的信号接收方法
CN102144377B (zh) 用于硬限信号的信道估计与均衡
CN1756248B (zh) 多入多出正交频分复用移动通信系统及信道估计方法
CN103701730B (zh) 信道时域相关性低复杂度压缩感知的信道估计方法及装置
CN105978674A (zh) 基于压缩感知的fdd下大规模mimo信道估计的导频优化方法
CN103634038B (zh) 基于多天线的联合doa估计和波束形成的多径信号接收方法
CN102624666B (zh) 稀疏信道模型下多路收发的正交多载波水声通信循环译码方法
CN101222458B (zh) Mimo-ofdm信道的低阶递推最小均方误差估计
CN104410590A (zh) 一种基于压缩感知的短波ofdm抑制干扰联合信道估计方法
CN101494528A (zh) 发射分集块传输系统的训练序列设计及其信道估计方法
CN103607262B (zh) 空时分组编码mimo系统中的两阶段预编码方法
CN106302298A (zh) 一种消除正交频分复用水声通信系统限幅噪声的方法
CN103297111A (zh) Mimo上行多用户信号检测方法、检测装置及接收系统
CN106161328A (zh) 基于载波索引调制的mimo‑ofdm系统的检测方法
CN102647259A (zh) 基于奇异值分解的多用户mimo系统上行链路传输方法
CN108259397A (zh) 基于自适应正则化子空间追踪压缩感知算法的大规模mimo系统信道估计
CN102227098A (zh) 一种多模mimo-scfde自适应传输系统频域承载点选取方法
CN103701571B (zh) 用于td‑lte‑a中继系统的八天线双码本设计方法
CN104065462A (zh) 中继干扰信道下具有分集增益的信号传输处理方法
CN100499611C (zh) 无线通信系统空域最大后验概率检测方法
CN104639473B (zh) 一种基于完全互补序列和压缩感知的mimo信道估计方法
CN106130613A (zh) 一种基于酉空时码获得灵活发射分集的空间调制方法
CN101662441B (zh) 信号估计方法、装置及检测系统
CN102075220B (zh) 一种基于时域降噪的信道估计装置及方法
CN102832986B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20160817