CN109768816A - 一种非高斯噪声3d-mimo系统数据检测方法 - Google Patents
一种非高斯噪声3d-mimo系统数据检测方法 Download PDFInfo
- Publication number
- CN109768816A CN109768816A CN201811552253.4A CN201811552253A CN109768816A CN 109768816 A CN109768816 A CN 109768816A CN 201811552253 A CN201811552253 A CN 201811552253A CN 109768816 A CN109768816 A CN 109768816A
- Authority
- CN
- China
- Prior art keywords
- function
- gaussian
- order
- mixed
- variance
- 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.)
- Granted
Links
Landscapes
- Radio Transmission System (AREA)
Abstract
一种非高斯噪声3D‑MIMO系统数据检测方法,首先根据接收信号的特征得到各高斯函数方差及其混合系数,再进行混合高斯模型的阶次选择计算;然后利用EM算法计算出该阶次下噪声的方差与混合系数,从而近似得到接收信号的概率密度函数,利用变分推断算法和接收信号的概率密度函数找到一个与待检测数据后验概率密度函数最相似的狄拉克函数,其相似程度用Kullback‑Leibler散度函数表征;最后利用梯度下降法得到该函数上升最快方向,求出狄拉克函数的参数,该参数即为待检测信号的估计值;本发明解决了传统数据检测算法针对非高斯噪声恢复误差过高的问题,有效地提高了3D‑MIMO通信系统的可靠性。
Description
技术领域
本发明涉及通信技术领域,特别涉及一种非高斯噪声3D-MIMO系统数据检测方法。
背景技术
随着无线通信技术的发展,无线网络的需求带动了无线数据业务的迅速增长。尤其是在以大数据为基础的万物互联时代,无线通信系统的传输速率需要进一步提升。在这种背景下,多输入多输出(MIMO)技术成为业界关注的热点。MIMO技术通过在基站上装置大规模天线整列以显著增加频谱效率,尤其是在容量需求较大或覆盖范围较广时,其可以显著改善数据传输质量。3D-MIMO技术以多天线为基础,通过采用二维天线阵列和先进的信号处理算法,动态地调整发送波束的下倾角,增加了信号在垂直方向上的自由度,可以实现精确的三维波束成形,使得信号能量更集中、信号方向更精确。由于3D-MIMO可以实现更好的干扰抑制和空间多用户复用,是提高系统容量和传输效率的有效手段,因此该技术成为了5G通信系统的核心技术之一。3D-MIMO系统数据检测技术已有一套完备的体系,然而这些方法都是以高斯噪声信道为基础来实现精确的估值。这是由于高斯噪声便于数学表达和计算,且这一假设符合中心极限定理。然而已有文献表明,实际信道噪声并不都是服从高斯分布。例如,电磁干扰、突发脉冲、大气噪声、闪电等环境下,信道噪声表现出非高斯特性。
发明内容
为了克服上述现有技术的缺陷,解决3D-MIMO系统在非高斯噪声下的信号检测问题,本发明的目的在于提供一种非高斯噪声3D-MIMO系统数据检测方法,将变分推断与梯度下降相结合,从而使接收端准确地检测出非高斯噪声干扰下的发送数据。
为了达到上述目的,本发明是通过下述技术方案来实现的。
一种非高斯噪声3D-MIMO系统数据检测方法,包括下述步骤:
(1)、根据接收信号y的特征,初始化混合高斯函数的阶次,利用Kmeans聚类算法计算该阶次下混合高斯函数中各高斯函数的均值,由均值及聚类结果计算出各高斯函数方差及其混合系数;
(2)、由步骤(1)中得到的各高斯函数方差及其混合系数进行混合高斯模型的阶次选择计算,若不满足判定条件,则更新混合高斯函数的阶次,并返回执行步骤(1)中的Kmeans算法继续计算阶次更新后的混合高斯函数参数,直到满足判定条件为止;
(3)、根据步骤(2)中最终得到的阶次,利用EM算法计算该阶次下各高斯函数的方差和混合高斯系数,在计算时,每次迭代过程中假设各高斯函数的均值为零;
(4)、结合步骤(3)得出的方差和混合高斯系数,用狄拉克函数近似待检测数据的后验概率,计算两个函数的Kullback-Leibler散度以衡量它们的相似程度;
(5)计算步骤(4)中Kullback-Leibler散度函数关于待检测数据的偏导,得到该函数上升最快方向;
(6)利用梯度下降算法迭代地计算出待检测数据。
所述步骤(1)中,混合高斯函数的各参数计算方法为:
给定阶次K正整数初始值,利用Kmeans聚类算法计算接收信号y在该阶次下各高斯函数的均值μi,i=1,2,…,K,利用均值及聚类结果得到各高斯函数的方差及混合系数:
上式中,Si为接收信号y中聚类结果为第i类的元素下标集合,Ni为Si中元素的个数,L表示接收信号y的维度,yk表示y中的第k个元素,λi分别表示第i个高斯函数的方差和混合系数。
所述步骤(2)中,阶次计算过程为:
对若
不满足此条件,则K=K-1。上式中,η1和η2分别表示在均值和方差方面使两个高斯分量合并为一个的阈值。
所述步骤(4)中,待检测信号x的后验概率密度与参数为x′的狄拉克函数δ(x-x′)的Kullback-Leibler散度为:
上式中,Q表示狄拉克函数δ(x-x′),p表示x的后验概率,H表示信道矩阵,ρ2I为发送信号x的协方差矩阵,c为某一常数。
所述步骤(5)中,偏导计算结果为:
所述步骤(6)中,迭代过程为:
上式中,t表示迭代次数,x′t为第t次迭代结果,μ(t)为第t步的迭代步长,为F(Q,p)的梯度,当满足迭代次数达到规定程度或者迭代结果趋于稳定时,该迭代终止,迭代结果便是该算法下待检测数据的值。
本发明提出了一种3D-MIMO信道非高斯噪声下的数据检测算法,该算法将信道噪声建模为混合高斯函数,根据接收信号的特征求出该函数的参数,得到了接收信号的概率密度函数,利用变分推断算法找到一种与该密度函数最相似的狄拉克函数。由于很难找到该函数参数的闭式解,因此本算法利用梯度下降法逐次迭代计算这一数据,该数据即为待检测信号。本发明解决了传统数据检测算法针对非高斯噪声恢复误差过高的问题,有效地提高了3D-MIMO通信系统的可靠性。
附图说明
图1为本算法的实施流程图。
图2为本发明所述算法与其它算法的对比。
具体实施方式
为了更加清晰说明本发明的目的、技术方案及优点,以下结合附图及实施例对本发明作进一步详细说明。应当理解,此处所述的具体实施例仅用以解释本发明,并不用于限定本发明。
考虑一个单蜂窝下行OFDM通信系统,其基站阵列上装置N个天线单元,假设信道是频率非选择性,则3D-MIMO系统接收信号模型可表示为:
y=Hx+n
上式中,H表示信道矩阵,表示用户接收信号,L表示多径数,其元素yi为第i条路径上的值,表示基站发送信号,为信道噪声,其分布未知。本发明考虑用均值为零的混合高斯函来近似代替其密度函数,即假设n中的任一元素n的概率密度函数为:
为了在接收端准确检测出发送数据x,本发明提出了一种变分推断与最陡下降法相结合的算法,实施步骤为:
(1)根据接收信号y的特征,初始化混合高斯函数的阶次,利用Kmeans聚类算法计算该阶次下混合高斯函数中各高斯函数的均值,由均值及聚类结果计算出各高斯函数方差及其混合系数;
具体而言,所述步骤(1)中,混合高斯函数的各参数计算方法为:
给定阶次K正整数初始值,利用Kmeans聚类算法计算接收信号y在该阶次下各高斯函数的均值μi,i=1,2,…,K,利用均值及聚类结果得到各高斯函数的方差及混合系数:
上式中,Si为接收信号y中聚类结果为第i类的元素下标集合,Ni为Si中元素的个数,yk表示y中的第k个元素,λi分别表示第i个高斯函数的方差和混合系数。
(2)由步骤(1)中得到的参数进行混合高斯模型的阶次选择计算,若不满足判定条件,则更新混合高斯函数的阶次,并返回执行步骤(1)中的Kmeans算法继续计算阶次更新后的混合高斯函数参数,直到满足判定条件为止;
阶次计算过程为:对若
不满足此条件,则K=K-1。上式中,η1和η2分别表示在均值和方差方面使两个高斯分量合并为一个的阈值。
(3)根据步骤(2)中最终得到的阶次,利用EM算法计算该阶次下各高斯函数的方差和混合系数,从而得到接收信号的概率密度函数,在计算时,每次迭代过程中假设各高斯函数的均值为零。
(4)结合步骤(3)得出的方差和混合系数,用狄拉克函数近似待检测数据的后验概率,计算两个函数的Kullback-Leibler散度以衡量它们的相似程度;
待检测信号x的后验概率密度函数与参数为x′的狄拉克函数δ(x-x′)的Kullback-Leibler散度为:
上式中,Q表示狄拉克函数δ(x-x′),p表示x的后验概率,H表示信道矩阵,ρ2I为发送信号x的协方差矩阵,c为某一常数。
(5)计算步骤(4)中Kullback-Leibler散度函数关于待检测数据的偏导,得到该函数上升最快方向;
偏导计算过程为:
(6)利用梯度下降算法迭代地计算出待检测数据。
所述步骤(6)中,迭代过程为:
上式中,t表示迭代次数,x′t为第t次迭代结果,μ(t)为第t步的迭代步长,为F(Q,p)的梯度。当满足迭代次数达到一定程度或者迭代结果趋于稳定时,该迭代终止,迭代结果便是该算法下待检测数据的值。
图2给出了按照本发明所述算法得到的检测信号归一化均方误差,图中其余两条曲线为最小二乘法(LS)与最小均方误差(MMSE)算法得的结果。需要说明的是,MMSE算法适用于高斯噪声信道,应用时需要提前得到噪声的方差。由于本发明所述模型为非高斯噪声,准确方差未知,为了便于应用该算法,仿真过程中将噪声按高斯分布处理,用接收信号方差代替了噪声方差。本次仿真设置参数为:路径数L=50,基站天线个数N=128,阶次初始值K=10,下降步长μ(t)=-0.5。为由图可见,本发明所述算法的归一化均方误差明显优于其它算法,且在低信噪比的情况下依然有着较为理想的检测性能。
Claims (6)
1.一种非高斯噪声3D-MIMO系统数据检测方法,其特征在于,包括下述步骤:
(1)、根据接收信号y的特征,初始化混合高斯函数的阶次,利用Kmeans聚类算法计算该阶次下混合高斯函数中各高斯函数的均值,由均值及聚类结果计算出各高斯函数方差及其混合系数;
(2)、由步骤(1)中得到的各高斯函数方差及其混合系数进行混合高斯模型的阶次选择计算,若不满足判定条件,则更新混合高斯函数的阶次,并返回执行步骤(1)中的Kmeans算法继续计算阶次更新后的混合高斯函数参数,直到满足判定条件为止;
(3)、根据步骤(2)中最终得到的阶次,利用EM算法计算该阶次下各高斯函数的方差和混合系数,在计算时,每次迭代过程中假设各高斯函数的均值为零;
(4)、结合步骤(3)中得出的方差和混合系数,用狄拉克函数近似待检测数据的后验概率,计算两个函数的Kullback-Leibler散度以衡量它们的相似程度;
(5)计算步骤(4)中Kullback-Leibler散度函数关于待检测数据的偏导,得到该函数上升最快方向;
(6)利用梯度下降算法迭代地计算出待检测数据。
2.根据权利要求1所述的一种非高斯噪声3D-MIMO系统数据检测方法,其特征在于,所述步骤(1)中,混合高斯函数的各参数计算方法为:
给定阶次K正整数初始值,利用Kmeans聚类算法计算接收信号y在该阶次下各高斯函数的均值μi,i=1,2,…,K,利用均值及聚类结果得到各高斯函数的方差及混合系数:
上式中,Si为接收信号y中聚类结果为第i类的元素下标集合,Ni为Si中元素的个数,L表示接收信号y的维度,yk表示y中的第k个元素,λi分别表示第i个高斯函数的方差和混合系数。
3.根据权利要求2所述的一种非高斯噪声3D-MIMO系统数据检测方法,其特征在于,所述步骤(2)中,阶次计算过程为:
对若
不满足此条件,则K=K-1,上式中,η1和η2分别表示在均值和方差方面使两个高斯分量合并为一个的阈值。
4.根据权利要求3所述的一种非高斯噪声3D-MIMO系统数据检测方法,其特征在于,所述步骤(4)中,待检测信号x的后验概率密度与参数为x′的狄拉克函数δ(x-x′)的Kullback-Leibler散度为:
上式中,Q表示狄拉克函数δ(x-x′),p表示x的后验概率,H表示信道矩阵,ρ2I表示发送信号x的协方差矩阵,c为某一常数。
5.根据权利要求4所述的一种非高斯噪声3D-MIMO系统数据检测方法,其特征在于,所述步骤(5)中,偏导计算过程为:
6.根据权利要求5所述的一种非高斯噪声3D-MIMO系统数据检测方法,其特征在于,所述步骤(6)中,迭代过程为:
上式中,t表示迭代次数,x′t为第t次迭代结果,μ(t)为第t步的迭代步长,为F(Q,p)的梯度,当满足迭代次数达到规定程度或者迭代结果趋于稳定时,该迭代终止,迭代结果便是该算法下待检测数据的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811552253.4A CN109768816B (zh) | 2018-12-19 | 2018-12-19 | 一种非高斯噪声3d-mimo系统数据检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811552253.4A CN109768816B (zh) | 2018-12-19 | 2018-12-19 | 一种非高斯噪声3d-mimo系统数据检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109768816A true CN109768816A (zh) | 2019-05-17 |
CN109768816B CN109768816B (zh) | 2020-11-17 |
Family
ID=66451541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811552253.4A Active CN109768816B (zh) | 2018-12-19 | 2018-12-19 | 一种非高斯噪声3d-mimo系统数据检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109768816B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110233808A (zh) * | 2019-05-31 | 2019-09-13 | 西安交通大学 | 一种ftn系统信号检测方法 |
CN110568415A (zh) * | 2019-07-22 | 2019-12-13 | 广东工业大学 | 混合高斯模型下基于Arctan函数的信号检测方法 |
CN112953607A (zh) * | 2021-02-22 | 2021-06-11 | 西安交通大学 | 一种mimo-ofdm系统量化噪声消除方法、介质及设备 |
CN113641957A (zh) * | 2021-06-29 | 2021-11-12 | 南京邮电大学 | 一种基于归一化流的信号检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103227761A (zh) * | 2013-05-24 | 2013-07-31 | 北京邮电大学 | 一种基于经验似然方法的多径非高斯噪声信道的估计方法 |
US20170264391A1 (en) * | 2016-03-09 | 2017-09-14 | Samsung Electronics Co., Ltd. | System and method for multiple input multiple output (mimo) detection with soft slicer |
CN107276646A (zh) * | 2017-07-03 | 2017-10-20 | 东南大学 | 大规模mimo‑ofdm上行链路迭代检测方法 |
CN108512787A (zh) * | 2018-04-13 | 2018-09-07 | 电子科技大学 | 大规模mimo系统的超精细信道估计方法 |
CN108736934A (zh) * | 2018-05-18 | 2018-11-02 | 东南大学 | 一种高效的大规模mimo系统信号检测方法 |
-
2018
- 2018-12-19 CN CN201811552253.4A patent/CN109768816B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103227761A (zh) * | 2013-05-24 | 2013-07-31 | 北京邮电大学 | 一种基于经验似然方法的多径非高斯噪声信道的估计方法 |
US20170264391A1 (en) * | 2016-03-09 | 2017-09-14 | Samsung Electronics Co., Ltd. | System and method for multiple input multiple output (mimo) detection with soft slicer |
CN107276646A (zh) * | 2017-07-03 | 2017-10-20 | 东南大学 | 大规模mimo‑ofdm上行链路迭代检测方法 |
CN108512787A (zh) * | 2018-04-13 | 2018-09-07 | 电子科技大学 | 大规模mimo系统的超精细信道估计方法 |
CN108736934A (zh) * | 2018-05-18 | 2018-11-02 | 东南大学 | 一种高效的大规模mimo系统信号检测方法 |
Non-Patent Citations (2)
Title |
---|
冯讯: "一种MIMO雷达非高斯杂波中的信号分离方法", 《信号处理》 * |
汪霜玲,何茜,何子述: "非高斯杂波下分置天线MIMO 雷达的动目标检测", 《计算机应用研究》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110233808A (zh) * | 2019-05-31 | 2019-09-13 | 西安交通大学 | 一种ftn系统信号检测方法 |
CN110233808B (zh) * | 2019-05-31 | 2020-07-28 | 西安交通大学 | 一种ftn系统信号检测方法 |
CN110568415A (zh) * | 2019-07-22 | 2019-12-13 | 广东工业大学 | 混合高斯模型下基于Arctan函数的信号检测方法 |
CN110568415B (zh) * | 2019-07-22 | 2022-04-08 | 广东工业大学 | 混合高斯模型下基于Arctan函数的信号检测方法 |
CN112953607A (zh) * | 2021-02-22 | 2021-06-11 | 西安交通大学 | 一种mimo-ofdm系统量化噪声消除方法、介质及设备 |
CN113641957A (zh) * | 2021-06-29 | 2021-11-12 | 南京邮电大学 | 一种基于归一化流的信号检测方法 |
CN113641957B (zh) * | 2021-06-29 | 2024-08-13 | 南京邮电大学 | 一种基于归一化流的信号检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109768816B (zh) | 2020-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210377079A1 (en) | Time-frequency block-sparse channel estimation method based on compressed sensing | |
Paul et al. | A novel beamspace channel estimation technique for millimeter wave massive MIMO systems | |
CN109768816A (zh) | 一种非高斯噪声3d-mimo系统数据检测方法 | |
CN108832976B (zh) | 一种大规模mimo系统的上行链路信道估计方法 | |
US8032184B2 (en) | Method for generating downlink beamforming weighting vectors | |
Benzaghta et al. | Massive MIMO systems for 5G: A systematic mapping study on antenna design challenges and channel estimation open issues | |
Awan et al. | Detection for 5G-NOMA: An online adaptive machine learning approach | |
CN110971547B (zh) | 一种基于毫米波/太赫兹的宽带大规模终端检测和信道估计方法 | |
CN106716866B (zh) | 乒乓波束成形 | |
CN113315560B (zh) | 一种紧凑型平面阵列Massive MIMO系统的波束赋形方法 | |
CN108199987A (zh) | 一种fdd大规模mimo系统中的下行链路信道估计方法 | |
WO2017126511A1 (en) | Method and computer implemented method for decoding symbol transmitted over mmwave channel, receiver for receiving and decoding symbols transmitted over mmwave channel | |
CN108199990B (zh) | 一种非高斯噪声3d-mimo信道估计方法 | |
Kim et al. | Spatial wideband channel estimation for mmWave massive MIMO systems with hybrid architectures and low-resolution ADCs | |
Cheng et al. | Channel path identification in mmWave systems with large-scale antenna arrays | |
CN105553526B (zh) | 大规模mimo系统导频长度及功率联合分配方法 | |
CN106788631A (zh) | 一种基于局部校准的大规模mimo互易性校准方法 | |
US20150222334A1 (en) | Monotonic Optimization Method for Achieving the Maximum Weighted Sum-Rate in Multicell Downlink MISO Systems | |
CN106357309A (zh) | 基于非理想信道下大规模mimo线性迭代检测方法 | |
CN105610479B (zh) | 一种大规模mu-mimo系统信道估计方法 | |
CN105812032B (zh) | 基于波束块结构压缩感知的信道估计方法 | |
CN105610484B (zh) | 大规模mimo低复杂度迭代接收方法 | |
CN108572347A (zh) | 基于通信信号信道状态响应信息的面阵二维测角方法及系统 | |
CN117614781A (zh) | 基于ris的两阶段超分辨率参数信道估计方法及装置 | |
CN101083498A (zh) | 优化阵列的接收天线间的间距的方法 |
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 |