Channel Estimation for High Speed Wireless Systems using Gaussian Particle Filter and Auxiliary Particle Filter
论文来源
International Conference on Communication and Signal Processing, April 4-6, 2019, India,Gargi Rajam, P. Sandeeptha and Sudheesh P
摘要
无线通信系统是从一个设备到另一个设备的传输介质。由于多径和多普勒频移,无误差传播变得不可能被实现。信道估计作为一种提高无线系统传输可靠性的技术变得越来越受欢迎。本篇论文通过使用GPF(Gaussian Particle Filter)和APF(Auxiliary Particle Filter)来进行信道估计。
基本概念
1.时变信道
在移动通信系统中由于收发双方的信道随着周围环境的变化而改变,但是由于系统的高移动性,使得信道发生快速变化.在接收端就会接收到由于多径产生的不同的频率响应.而使用PF(Particle Filter)或者SMC(Sequential Monte Carlo)可以解决这个问题.状态方程如下:
\]
\]
其中, A是协方差矩阵如下所示:
\left\{
\begin{matrix}
a_{11} & \cdots & a_{1n}\\
\vdots & \ddots & \vdots \\
a_{n1} & \cdots & a_{nn}
\end{matrix}
\right\}
\]
\(x_t= \left[ \begin{matrix} x_t(0) & \cdots & x_t(m-1) \end{matrix} \right]\) 是第t个时刻的输入向量.
\(h_t= \left[ \begin{matrix} h_t(0) & \cdots & h_t(n-1) \end{matrix} \right]\) 是第t个时刻的信道系数.
\(y_t= \left[ \begin{matrix} y_t(0) & \cdots & y_t(m+n-1) \end{matrix} \right]\) 是第t个时刻的输出向量.
\(w_t= \left[ \begin{matrix} w_t(0) & \cdots & w_t(n-1) \end{matrix} \right]\) 是第t个时刻的系统噪声.
2.粒子滤波
粒子滤波器广泛应用于通信、信号处理、图像处理、视觉和目标跟踪等领域.具体算法如下:
1)初始化N个权重相等的随机粒子.
2)根据贝叶斯公式:\(zw_t \propto w^i_{t-1}* \frac{p(y_t|x^i_t)p(y_t|x^i_{t-1})}{q(x^i_t|x^i_{t-1},y_t)}\) .权重近似为:\(w_t = p(y_t|h^i_t)\)
3) 重采样
具体如下图所示:

3.高斯粒子滤波
GPF(Gaussian Particle Filter)比EKF(Extend Kalman Filter)和UKF(Unscented Kalman Filter)有更快的收敛性.
4.辅助粒子滤波
APF(Auxiliary Particle Filter)均匀的改变权重,防止问题过快退化.但是更好的适用于少量数据的估计.算法具体如下

比较
PF,GPF和APF估计的信道均方误差如下所示:

均方误差公式为$$ RMSE = \sqrt{ \frac{ \sum_{i=1}^T(h_t-\hat{h_t})^2}{T}} $$
GPF克服了这一问题,因为收敛速度更快,但需要更多的粒子数.
APF显示出更好的结果,不会在粒子数和通道长度之间做任何妥协.
借鉴之处
1.将同一问题的不同解决方法放在一块比较.(可提出新的优化方法)
2.使用均方误差来对不同方法估计出的信道进行评判好坏.
Channel Estimation for High Speed Wireless Systems using Gaussian Particle Filter and Auxiliary Particle Filter的更多相关文章
- 信道估计(channel estimation)图解——从SISO到MIMO原理介绍
1. 引言 在所有通信中,信号都会通过一个介质(称为信道),并且信号会失真,或者在信号通过信道时会向信号中添加各种噪声.正确解码接收到的信号而没有太多错误的方法是从接收到的信号中消除信道施加的失真和噪 ...
- Introducing Microsoft Sync Framework: Sync Services for File Systems
https://msdn.microsoft.com/en-us/sync/bb887623 Introduction to Microsoft Sync Framework File Synchro ...
- Cesium中级教程8 - Introduction to Particle Systems 粒子系统入门
Cesium中文网:http://cesiumcn.org/ | 国内快速访问:http://cesium.coinidea.com/ What is a particle system? 什么是粒子 ...
- usb.ids
# # List of USB ID's # # Maintained by Vojtech Pavlik <vojtech@suse.cz> # If you have any new ...
- ACM会议列表与介绍(2014/05/06)
Conferences ACM SEACM Southeast Regional Conference ACM Southeast Regional Conference the oldest, co ...
- 简明Python3教程 1.翻译
这里有许许多多本书不同语言的译本,感谢那些不知疲倦的志愿者们! 如果你想帮助这些译本,请参看下面的关于志愿者和语言的列表,然后决定是要开始一个新的翻译项目还是帮助改进现有的翻译项目. 如果你计划开始一 ...
- 斯坦福CS课程列表
http://exploredegrees.stanford.edu/coursedescriptions/cs/ CS 101. Introduction to Computing Principl ...
- LTE Air interface Channels-----http://www.rfwireless-world.com/Tutorials/LTE-logical-transport-physical-channels.html
LTE technology works based on three channel types viz. logical channel,transport channel and physica ...
- linux运维中的命令梳理(四)
----------管理命令---------- ps命令:查看进程 要对系统中进程进行监测控制,查看状态,内存,CPU的使用情况,使用命令:/bin/ps (1) ps :是显示瞬间进程的状态,并不 ...
随机推荐
- 一张图看懂三维GIS
- 五、linux-mysql 下mysql的管理(二)
1.插入语句: create table test( id int(4) not null auto_increment, name char(20) not null, primary ...
- scala编程(七)——内建控制结构
几乎所有的 Scala 的控制结构都会产生某个值.这是函数式语言所采用的方式,程序被看成是计算值的活动,因此程序的控件也应当这么做.另外,指令式语言经常具有三元操作符(如 C,C++和 Java 的? ...
- 用windows下的Anaconda搭建Django虚拟环境
Django 是一个Python定制框架,可用于简便.快速的开发数据库驱动的web站点. 要使用Django,首先要建立一个虚拟工作环境. 那么,为什么要搭建虚拟环境呢?我们来看以下的开发场景: 假设 ...
- winform显示word和ppt文档
最近所做的项目中需要在Winform窗体中显示Office文档.刚开始就使用webBrowser控件实现的,但是后来发现这个控件在显示Office文档的时候有个限制:只支持Office2003之前的版 ...
- CentOS-DHCP服务搭建
title date tags layout CentOS6.5 DHCP服务器搭建 2018-08-26 Centos6.5服务器搭建 post 1.安装dhcp软件包 yum install -y ...
- SMTP错误码/建议解决方法
SMTP错误码/建议解决方法 错误总表 420 1. Timeout Communication Problem Encountered During Transmission. Thie Is a ...
- Leetcode_开篇碎碎念
已经计划写leetcode3天了,但是到目前为止一道题还没有写,执行力我太差了. 写leetcode的原因有两个,一是锻炼自己代码水平,提高实力:二是考研据了解到的都会有上机考试,基本是c++,jav ...
- HTML常用数据类型
.数学函数: Math.ceil():天花板数 //大于当前小数的最小整数 Math.floor():地板数 //小于当前小数的最大整数 Math.round():四舍五入取整数 Math.rando ...
- [2015普及组-D]推销员 思维que
题:https://www.cometoj.com/problem/0221 #include<iostream> #include<cstring> #include< ...