CN107911182A - 一种无线信道环境特征参数突变检测方法 - Google Patents

一种无线信道环境特征参数突变检测方法 Download PDF

Info

Publication number
CN107911182A
CN107911182A CN201711164439.8A CN201711164439A CN107911182A CN 107911182 A CN107911182 A CN 107911182A CN 201711164439 A CN201711164439 A CN 201711164439A CN 107911182 A CN107911182 A CN 107911182A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
interval
munderover
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
Application number
CN201711164439.8A
Other languages
English (en)
Other versions
CN107911182B (zh
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.)
Hefei University of Technology
Original Assignee
Hefei University 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 Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201711164439.8A priority Critical patent/CN107911182B/zh
Publication of CN107911182A publication Critical patent/CN107911182A/zh
Application granted granted Critical
Publication of CN107911182B publication Critical patent/CN107911182B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/309—Measuring or estimating channel quality parameters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3911—Fading models or fading generators
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H04B17/30—Monitoring; Testing of propagation channels
    • H04B17/391—Modelling the propagation channel
    • H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明提供一种无线信道环境特征参数突变检测方法，包括如下步骤：利用无线信道大尺度模型计算阴影衰落序列；计算阴影衰落序列的增量序列；对增量序列做赋值、取值划分处理；根据增量序列的概率分布，计算无线信道的先验状态转移概率矩阵；计算无线信道的先验状态转移概率矩阵的相似度，统计正常通信相似度波动范围；统计监测样本信号下一采样时刻点实际所处的状态区间的个数，计算监测数据状态转移概率矩阵；计算无线信道的先验状态转移概率矩阵与监测数据状态转移概率矩阵的相似度，判断所监测信道环境特征参数是否发生突变。本发明的有益效果在于：运算简单快速，便于实现实时监测；对异常的发生较为敏感，能够准确确定异常的发生。

Description

一种无线信道环境特征参数突变检测方法
技术领域
本发明属于无线信道特征提取技术领域，具体涉及一种无线信道环境特征参数突变检测方法。
背景技术
在移动通信中，发射端和接收端之间通过电磁波来传输信号，我们可以想象两者之间有一些看不见的电磁通路，并把这些电磁通路称为无线信道。无线信道与周围的环境密切相关，不同环境下的无线信道具有一些差异化的特征，如何发现并提取这些特征并将其应用于优化无线网络，是当前的一个研究热点。
无线信道的大尺度衰减特性是无线通信网络规划和覆盖范围预测的重要研究内容。长期以来，移动信道和地空信道的传播特性在相关领域受到重视并得到广泛研究，近地面信道由于缺乏应用需求的推动而研究较少。近年来，随着无线传感器网络技术的不断发展，其在地面或低空环境下的广泛应用前景开始引起人们对近地面无线信道研究的重视，但当前的研究重点仍集中于上层路由、组网以及接入协议，涉及信道传播特性的研究工作很少。在极低天线高度下，无线信号传播对信道的考察较以往的研究要求更细微化。地面的微小起伏都会阻挡收发信机间的传播路径而形成非视距(non-line-of-sight,NLOS)通信，不同地形地貌引起的阴影衰落更是千差万别，该种环境即使能进行有效的信号强度测量，也难以形成可借鉴的信道模型，当前的信号检测技术大部分已无法满足高精度检测场景的需求。
发明内容
为了解决现有信号检测技术无法满足高精度检测场景需求的技术缺陷，本发明提供一种利用无线信道大尺度衰减采集阴影衰落序列，建立马尔可夫转移概率矩阵，准确快速的检测信道环境是否发生突变的方法。
本发明是通过以下技术方案实现的：
一种无线信道环境特征参数突变检测方法，包括如下步骤：
步骤S1，在正常通信情况下，采集无线信道内的功率信号序列，并利用无线信道大尺度模型计算阴影衰落序列；
步骤S2，计算所述阴影衰落序列的增量序列；对所述增量序列做赋值、取值划分处理；根据所述取值划分后所述增量序列的概率分布，计算无线信道的先验状态转移概率矩阵；
步骤S3，在正常通信情况下，重复所述步骤S1和步骤S2若干次；计算所述无线信道的先验状态转移概率矩阵的相似度，称为第一相似度，统计正常通信相似度波动范围；
步骤S4，统计监测样本信号下一采样时刻点实际所处的状态区间的个数，计算监测数据状态转移概率矩阵；计算所述无线信道的先验状态转移概率矩阵与所述监测数据状态转移概率矩阵的相似度，称为第二相似度，并根据所述正常通信相似度波动范围和所述第二相似度判断所监测信道环境特征参数是否发生突变。
本发明相对于现有技术的有益效果在于：
第一，运算简单快速，便于实现实时监测；
第二，状态转移概率矩阵对异常的发生较为敏感，能够准确确定异常的发生。
附图说明
图1为实施例1的无线信道环境特征参数突变检测方法的总流程图；
图2为阴影衰落序列采集的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例1：
首先，对本发明提供的无线信道环境特征参数突变检测方法所需的硬件环境做简要说明：本发明提供的检测方法针对的是无线信道(即无线通信中发送端和接收端之间通路)的环境特征参数突变检测，一般情况下，该方法基于的硬件环境包括：一台信号发射机、一台信号接收机、一组发射天线和一组接收天线，在检测过程中，利用信号发射机通过发射天线向外发射信号，信号接收机通过接收天线接收上述信号。
如图1所示，无线信道环境特征参数突变检测方法包括如下步骤：
步骤S1，在正常通信情况下，采集无线信道内的功率信号序列，并利用无线信道大尺度模型计算阴影衰落序列。具体包括如下步骤：
步骤S101，在正常通信情况下，在信号接收机端以远小于原信号频率fo的采样频率fc 采集功率信号序列；
步骤S102，利用无线信道大尺度模型，见下式：
其中，dn为信号接收机到信号发射机的距离；PL(dn)[dB]为路径功率损耗对数形式，dB为单位；为参考点路径功率损耗，d0为已知参考距离；为路径损耗因子；
计算阴影衰落序列
在本实施例中，步骤S102中的无线信道大尺度模型中计算阴影衰落序列的方法进一步包括：
步骤S1021，计算路径功率损耗PL(dn)[dB]，如下式所示：
其中Pt,Gt,分别为发射信号功率、发射增益、dn位置处的接收增益和接收信号功率；
步骤S1022，计算参考位置d0处的路径功率损耗如下式所示：
其中Pt,Gt,分别为发射信号功率、发射增益、参考点d0位置处的接收增益和接收信号功率；
步骤S1023，计算值，如下式所示：
步骤S1024，计算阴影衰落序列，如下式所示：
其中，PL(dn,m)表示为路径功率损耗序列PL(dn)的第m个值，m＝1,2,...,M。
步骤S2，计算所述阴影衰落序列的增量序列；对所述增量序列做赋值、取值划分处理；根据所述取值划分后所述增量序列的概率分布，计算无线信道的先验状态转移概率矩阵。
步骤S201，计算所述阴影衰落序列的增量序列：m＝1,2,...,M-1；其中，为阴影衰落序列；
步骤S202，对所述增量序列做赋值处理，如下式：
其中，“←”为赋值运算，为的均值；
再对所述增量序列做取值划分处理，即，将的取值划分为5个区间，计算每个区间的均值 表示第i个区间的均值，i＝1,2,...,5；若表示落在第i,i＝1,2,...,5个区间内的第m'个值，M'为落在第i区间的个数，则计算其中，表示在原序列中序列数为mi所对应的值，为的方差，按以下规则(1)-(5)判定所属区间：
规则(1)，当s∈[0,1)时，与属于同一区间；
规则(2)，当s∈[1,2)时，若且i≤4，则位于右侧相邻区间，若i＝5，则也位于第5区间；若且i≥2，位于左侧相邻区间，若i＝1，则也位于第1区间；规则(3)，当s∈[2,3)时，若且i≤3，则位于右侧第2个区间，若i≥4，则位于第5区间；若且i≥3，位于左侧第2个区间，若i≤2，则位于第1区间；规则(4)，当s∈[3,4)时，若且i≤2，则位于右侧第3个区间，若i≥3，则位于第5区间；若且i≥4，位于左侧第3个区间，若i≤3，则位于第1区间；规则(5)，当s∈[4,+∞)时，若则位于第5区间；若则位于第1区间。
其中，所述增量序列的5个状态区间的划分方法具体为：
以为第1区间，其中，为序列的标准差；以为第2区间；以为第3区间；以为第4区间；以为第5区间。
步骤S203，统计落入各区间的个数Nij，i＝1,2,...,5，j＝1,2,...,5；
计算前一时刻在第i区间，后一时刻落在第j区间的概率pij，如下式所示：
则先验状态转移概率矩阵Pth为：
步骤S3，在正常通信情况下，重复所述步骤S1和步骤S2若干次；计算所述无线信道的先验状态转移概率矩阵的相似度，称为第一相似度，统计正常通信相似度波动范围。
步骤S301，在正常通信情况下，重复所述步骤S1和步骤S2L次；其中，L为大于100的整数；
步骤S302，计算所述无线信道的先验状态转移概率矩阵Pth的相似度el,l＝1,2,...,L，称为第一相似度，如下式：
统计正常通信相似度波动范围：
其中，为第l次概率矩阵，其元素为如果满足则波动范围e≥elmin，其中elmin为el,l＝1,2,...,L的最小值。
步骤S4，统计监测样本信号下一采样时刻点实际所处的状态区间的个数，计算监测数据状态转移概率矩阵；计算所述无线信道的先验状态转移概率矩阵与所述监测数据状态转移概率矩阵的相似度，称为第二相似度，并根据所述正常通信相似度波动范围和所述第二相似度判断所监测信道环境特征参数是否发生突变。
步骤S401，统计监测样本信号下一采样时刻点实际所处的状态区间为j,j＝1,2,...,5的个数
步骤S402，记Pre元素为 计算监测数据状态转移概率矩阵：
步骤S403，计算所述无线信道的先验状态转移概率矩阵与所述监测数据状态转移概率矩阵Pre的相似度ere，称为第二相似度，如下式所示：
其中，若ere≥elmin，则监测信道无异常，否则所监测信道环境特征参数发生突变。
本发明的有益效果如下：
本发明提供的无线信道环境特征参数突变检测方法基于马尔可夫性原理实现。马尔可夫过程(Markov process)是一类随机过程，具有马尔可夫性特性，即在已知“现在”的条件下，“将来”与“过去”无关，因此也称作无后效性。本发明引入无线信道大尺度衰减功率的马尔可夫转移概率，可以简化监测模型，准确快速的确定监测信道环境是否发生突变。
现有的信号检测方法一般基于频谱检测、压缩感知检测法，而本发明引入马尔可夫概率矩阵检测后，计算复杂度降低。
并且作为一种能量检测法，本发明具有通用性、盲检测和高性能。
相比传统的检测参数，本发明采用的马尔可夫概率矩阵检测参数所需的先验知识较少。本发明对突变的反应主要是概率计算，比数值计算的可靠性更高。
本发明对大尺度衰减的阴影衰落序列使用马尔可夫概率矩阵检测主要是从非视距因素出发，因此对于不同的环境具有良好的适应性。
本发明引入的马尔可夫概率矩阵检测研究的是数值的波动变化程度，因此其检测应用范围较广，可用于检测原信号的切换、环境的变换及异常信号介入等等，例如，所处检测环境不变，正在发送的信号变成另一个不同的信号，或发送信号不变，检测范围内物理环境受到破坏；再或者，检测范围内有异常信号介入影响正常通信时，所检测特征参数都将发生突变。本发明提供的检测方法能够快速有效准确检测上述无线信道环境突变情况。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种无线信道环境特征参数突变检测方法，其特征在于，包括如下步骤：
步骤S1，在正常通信情况下，采集无线信道内的功率信号序列，并利用无线信道大尺度模型计算阴影衰落序列；
步骤S2，计算所述阴影衰落序列的增量序列；对所述增量序列做赋值、取值划分处理；根据所述取值划分后所述增量序列的概率分布，计算无线信道的先验状态转移概率矩阵；
步骤S3，在正常通信情况下，重复所述步骤S1和步骤S2若干次；计算所述无线信道的先验状态转移概率矩阵的相似度，称为第一相似度，统计正常通信相似度波动范围；
步骤S4，统计监测样本信号下一采样时刻点实际所处的状态区间的个数，计算监测数据状态转移概率矩阵；计算所述无线信道的先验状态转移概率矩阵与所述监测数据状态转移概率矩阵的相似度，称为第二相似度，并根据所述正常通信相似度波动范围和所述第二相似度判断所监测信道环境特征参数是否发生突变。
2.根据权利要求1所述的一种无线信道环境特征参数突变检测方法，其特征在于，所述步骤S2进一步包括：
步骤S201，计算所述阴影衰落序列的增量序列：m＝1,2,...,M-1；其中，为阴影衰落序列；
步骤S202，对所述增量序列做赋值处理，如下式：
<mrow> <msub> <mi>&amp;Delta;&amp;epsiv;</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>&amp;LeftArrow;</mo> <mo>|</mo> <msub> <mi>&amp;Delta;&amp;epsiv;</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;Delta;&amp;mu;</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mo>|</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，“←”为赋值运算，为的均值；
再对所述增量序列做取值划分处理，即，将的取值划分为5个区间，计算每个区间的均值 表示第i个区间的均值，i＝1,2,...,5；若表示落在第i,i＝1,2,...,5个区间内的第m'个值，M'为落在第i区间的个数，则计算其中，表示在原序列中序列数为mi所对应的值，为的方差，按以下规则(1)-(5)判定所属区间：
规则(1)，当s∈[0,1)时，与属于同一区间；
规则(2)，当s∈[1,2)时，若且i≤4，则位于右侧相邻区间，若i＝5，则也位于第5区间；若且i≥2，位于左侧相邻区间，若i＝1，则也位于第1区间；规则(3)，当s∈[2,3)时，若且i≤3，则位于右侧第2个区间，若i≥4，则位于第5区间；若且i≥3，位于左侧第2个区间，若i≤2，则位于第1区间；规则(4)，当s∈[3,4)时，若且i≤2，则位于右侧第3个区间，若i≥3，则位于第5区间；若且i≥4，位于左侧第3个区间，若i≤3，则位于第1区间；规则(5)，当s∈[4,+∞)时，若则位于第5区间；若则位于第1区间；
步骤S203，统计落入各区间的个数Nij，i＝1,2,...,5，j＝1,2,...,5；
计算前一时刻在第i区间，后一时刻落在第j区间的概率pij，如下式所示：
<mrow> <msub> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>N</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mrow> <msub> <mi>N</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> </msub> <mo>+</mo> <msub> <mi>N</mi> <mrow> <mi>i</mi> <mn>2</mn> </mrow> </msub> <mo>+</mo> <mo>...</mo> <mo>+</mo> <msub> <mi>N</mi> <mrow> <mi>i</mi> <mn>5</mn> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
则先验状态转移概率矩阵Pth为：
<mrow> <msub> <mi>P</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>p</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mn>12</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>p</mi> <mn>15</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mn>22</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>p</mi> <mn>25</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mrow></mrow> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>p</mi> <mn>51</mn> </msub> </mtd> <mtd> <msub> <mi>p</mi> <mn>12</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>p</mi> <mn>55</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mrow> <mn>5</mn> <mo>&amp;times;</mo> <mn>5</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
3.根据权利要求1所述的一种无线信道环境特征参数突变检测方法，其特征在于，所述步骤S1进一步包括：
步骤S101，在正常通信情况下，在信号接收机端以远小于原信号频率fo的采样频率fc 采集无线信道内的功率信号序列；
步骤S102，利用无线信道大尺度模型，见下式：
<mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;lsqb;</mo> <mi>d</mi> <mi>B</mi> <mo>&amp;rsqb;</mo> <mo>=</mo> <mover> <mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>+</mo> <mn>10</mn> <msub> <mi>k</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>d</mi> <mi>n</mi> </msub> <msub> <mi>d</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;epsiv;</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，dn为信号接收机到信号发射机的距离；PL(dn)[dB]为路径功率损耗对数形式，dB为单位；为参考点路径功率损耗，d0为已知参考距离；为路径损耗因子；
计算阴影衰落序列
4.根据权利要求1所述的一种无线信道环境特征参数突变检测方法，其特征在于，所述步骤S3进一步包括：
步骤S301，在正常通信情况下，重复所述步骤S1和步骤S2L次；其中，L为大于100的整数；
步骤S302，计算所述无线信道的先验状态转移概率矩阵Pth的相似度el,l＝1,2,...,L，称为第一相似度，如下式：
<mrow> <msub> <mi>e</mi> <mi>l</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>1</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> <mn>1</mn> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>l</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> <mi>l</mi> </msubsup> <mo>)</mo> </mrow> </mrow> <mrow> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>1</mn> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>l</mi> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mfrac> <mo>,</mo> <mi>l</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>L</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
统计正常通信相似度波动范围：
其中，为第l次概率矩阵，其元素为如果满足则波动范围e≥elmin，其中elmin为el,l＝1,2,...,L的最小值。
5.根据权利要求1所述的一种无线信道环境特征参数突变检测方法，其特征在于，所述步骤S4进一步包括：
步骤S401，统计监测样本信号下一采样时刻点实际所处的状态区间为j,j＝1,2,...,5的个数
步骤S402，记Pre元素为 计算监测数据状态转移概率矩阵：
步骤S403，计算所述无线信道的先验状态转移概率矩阵与所述监测数据状态转移概率矩阵Pre的相似度ere，称为第二相似度，如下式所示：
<mrow> <msub> <mi>e</mi> <mrow> <mi>r</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>1</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>t</mi> <mi>h</mi> </mrow> <mn>1</mn> </msubsup> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mover> <msub> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <msub> <mi>&amp;mu;</mi> <mrow> <mi>r</mi> <mi>e</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <msup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mn>1</mn> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>5</mn> </munderover> <msup> <mover> <msub> <mi>p</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;OverBar;</mo> </mover> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，若ere≥elmin，则监测信道无异常，否则所监测信道环境特征参数发生突变。
6.根据权利要求3所述的一种无线信道环境特征参数突变检测方法，其特征在于，在所述步骤S102中，无线信道大尺度模型中计算阴影衰落序列的方法进一步包括：
步骤S1021，计算路径功率损耗PL(dn)[dB]，如下式所示：
<mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;lsqb;</mo> <mi>d</mi> <mi>B</mi> <mo>&amp;rsqb;</mo> <mo>=</mo> <mi>log</mi> <mi> </mi> <msub> <mi>P</mi> <mrow> <mi>r</mi> <mo>,</mo> <msub> <mi>d</mi> <mi>n</mi> </msub> </mrow> </msub> <mo>-</mo> <mi>log</mi> <mi> </mi> <msub> <mi>G</mi> <mi>t</mi> </msub> <mo>-</mo> <mi>log</mi> <mi> </mi> <msub> <mi>G</mi> <mrow> <mi>r</mi> <mo>,</mo> <msub> <mi>d</mi> <mi>n</mi> </msub> </mrow> </msub> <mo>-</mo> <mi>log</mi> <mi> </mi> <msub> <mi>P</mi> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中分别为发射信号功率、发射增益、dn位置处的接收增益和接收信号功率；
步骤S1022，计算参考位置d0处的路径功率损耗如下式所示：
<mrow> <mover> <mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;lsqb;</mo> <mi>d</mi> <mi>B</mi> <mo>&amp;rsqb;</mo> <mo>=</mo> <mi>log</mi> <mi> </mi> <msub> <mi>P</mi> <mrow> <mi>r</mi> <mo>,</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> </mrow> </msub> <mo>-</mo> <mi>log</mi> <mi> </mi> <msub> <mi>G</mi> <mi>t</mi> </msub> <mo>-</mo> <mi>log</mi> <mi> </mi> <msub> <mi>G</mi> <mrow> <mi>r</mi> <mo>,</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> </mrow> </msub> <mo>-</mo> <mi>log</mi> <mi> </mi> <msub> <mi>P</mi> <mi>t</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中分别为发射信号功率、发射增益、参考点d0位置处的接收增益和接收信号功率；
步骤S1023，计算值，如下式所示：
<mrow> <msub> <mi>k</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mo>=</mo> <mfrac> <mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mover> <mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> </mrow> <mrow> <mn>10</mn> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>d</mi> <mi>n</mi> </msub> <msub> <mi>d</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤S1024，计算阴影衰落序列，如下式所示：
<mrow> <msub> <mi>&amp;epsiv;</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mrow> <mo>(</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>&amp;lsqb;</mo> <mi>d</mi> <mi>B</mi> <mo>&amp;rsqb;</mo> <mo>=</mo> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mi>n</mi> </msub> <mo>,</mo> <mi>m</mi> <mo>)</mo> </mrow> <mo>-</mo> <mover> <mrow> <mi>P</mi> <mi>L</mi> <mrow> <mo>(</mo> <msub> <mi>d</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> <mo>&amp;OverBar;</mo> </mover> <mo>-</mo> <mn>10</mn> <msub> <mi>k</mi> <msub> <mi>d</mi> <mi>n</mi> </msub> </msub> <mi>l</mi> <mi>o</mi> <mi>g</mi> <mrow> <mo>(</mo> <mfrac> <msub> <mi>d</mi> <mi>n</mi> </msub> <msub> <mi>d</mi> <mn>0</mn> </msub> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，PL(dn,m)表示为路径功率损耗序列PL(dn)的第m个值，m＝1,2,...,M。
7.根据权利要求2所述的一种无线信道环境特征参数突变检测方法，其特征在于，在所述步骤S202中，所述增量序列的5个状态区间的划分方法具体为：
以为第1区间，其中，为序列的标准差；以为第2区间；以为第3区间；以为第4区间；以为第5区间。
CN201711164439.8A 2017-11-21 2017-11-21 一种无线信道环境特征参数突变检测方法 Active CN107911182B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711164439.8A CN107911182B (zh) 2017-11-21 2017-11-21 一种无线信道环境特征参数突变检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711164439.8A CN107911182B (zh) 2017-11-21 2017-11-21 一种无线信道环境特征参数突变检测方法

Publications (2)

Publication Number Publication Date
CN107911182A true CN107911182A (zh) 2018-04-13
CN107911182B CN107911182B (zh) 2020-07-14

Family

ID=61846905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711164439.8A Active CN107911182B (zh) 2017-11-21 2017-11-21 一种无线信道环境特征参数突变检测方法

Country Status (1)

Country Link
CN (1) CN107911182B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006005791A (ja) * 2004-06-18 2006-01-05 Tokyo Univ Of Science 通信路推定及びデータ検出方法
CN1852543A (zh) * 2006-01-09 2006-10-25 华为技术有限公司 一种物理随机接入信道异常的检测方法
US8354925B1 (en) * 2008-07-14 2013-01-15 Lockheed Martin Corporation Monitoring using RF communication technology

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006005791A (ja) * 2004-06-18 2006-01-05 Tokyo Univ Of Science 通信路推定及びデータ検出方法
CN1852543A (zh) * 2006-01-09 2006-10-25 华为技术有限公司 一种物理随机接入信道异常的检测方法
US8354925B1 (en) * 2008-07-14 2013-01-15 Lockheed Martin Corporation Monitoring using RF communication technology

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
周作成, 张延华: "基于动态规划( DP) 的信道模型变化检测", 《电子学报》 *
陈典铖: "基于突变检测的频谱感知算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Also Published As

Publication number Publication date
CN107911182B (zh) 2020-07-14

Similar Documents

Publication Publication Date Title
Zou et al. Freedetector: Device-free occupancy detection with commodity wifi
CN102404058B (zh) 电磁波区分设备、电磁波区分方法以及电磁波区分程序
CN100518012C (zh) 认知无线电系统的授权用户信号检测方法
Alsindi et al. Analysis of time of arrival estimation using wideband measurements of indoor radio propagations
Wei et al. Analysis of scaled largest eigenvalue based detection for spectrum sensing
Liu et al. A research on CSI-based human motion detection in complex scenarios
CN109698724A (zh) 入侵检测方法、装置、设备及存储介质
CN103532647A (zh) 基于WiFi物理层时域特征的视距传播路径判断方法和系统
CN109143157B (zh) 基于信号强度值混合滤波的动态更新参数的测距方法
CN104683050B (zh) 一种能有效对抗噪声不确定性的多天线全盲频谱感知方法
Nguyen et al. Improved localization accuracy using machine learning: Predicting and refining RSS measurements
CN103346847A (zh) 基于迭代攻击检测的联合压缩频谱感知方法
CN103905129B (zh) 基于谱型分析的信号检测及信号信息判读方法
Xin-Di et al. The improvement of RSS-based location fingerprint technology for cellular networks
CN103346849A (zh) 一种抵抗模仿主用户信号的恶意攻击的频谱感知方法
CN107911182B (zh) 一种无线信道环境特征参数突变检测方法
Choi et al. Simultaneous crowd estimation in counting and localization using WiFi CSI
Kotian et al. Assessment of proactive transmission power control for wireless sensor networks
GB2476379A (en) Measuring a radio wave reception environment using RSSI values
CN102546116A (zh) 基于压缩感知的跳频信号盲检测方法
CN106686719B (zh) 一种rssi拟合方法
Korolkov et al. Received-signal-strength-based approach for detection and 2D indoor localization of evil twin rogue access point in 802.11
Deng et al. Data pre-processing in compressive sensing based indoor fingerprinting positioning
CN118590989B (zh) 一种基于随机森林的Wi-Fi指纹室内定位方法及系统
Gao et al. Threshold selection for energy detection in Cognitive Radio with a non-line-Of-sight environment

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