CN108566256B - 一种频谱地图的构建方法 - Google Patents
一种频谱地图的构建方法 Download PDFInfo
- Publication number
- CN108566256B CN108566256B CN201810245329.2A CN201810245329A CN108566256B CN 108566256 B CN108566256 B CN 108566256B CN 201810245329 A CN201810245329 A CN 201810245329A CN 108566256 B CN108566256 B CN 108566256B
- Authority
- CN
- China
- Prior art keywords
- expression
- center point
- specifically
- point position
- value
- 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.)
- Active
Links
Images
Classifications
-
- 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
-
- 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
- H04B17/318—Received signal strength
- H04B17/327—Received signal code power [RSCP]
-
- 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
- H04B17/345—Interference values
-
- 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
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W16/00—Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
- H04W16/18—Network planning tools
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W84/00—Network topologies
- H04W84/18—Self-organising networks, e.g. ad-hoc networks or sensor networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Quality & Reliability (AREA)
- Radar Systems Or Details Thereof (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供一种频谱地图的构建方法,包括以下步骤:第一步:利用区域内随机部署的多个频谱传感器采集获得关于功率的观测值;第二步:预先设置中心点数量和衰减指数初始值;采用一种改进的算法来选取中心点位置初始值;第三步:采用最小二乘拟合方法结合迭代交替最小化方法得到权重、中心点位置和衰减指数三者的估计值;第四步:通过特定表达式获得区域内任意位置处的功率估计值,从而得到区域内关于功率的频谱地图。应用本发明方法,效果是:本发明方法无需利用关于发射机和传播环境的相关信息,从而应用场景方面更具备普适性;所需的观测点数量更少,从而对传感器数量和部署要求更低,更经济;本发明方法无需进行参数调节,更利于实现。
Description
技术领域
本发明涉及通信技术领域,特别地,涉及一种频谱地图的构建方法。
背景技术
频谱地图就是将频谱信息叠加在地图上,通常表示某一无线电参数(如接收信号功率、信道增益、干扰功率等等)在感兴趣区域内的分布状况。由于考虑到了空闲频谱的空间复用,频谱地图信息更有利于提高频率利用率。
最常用的频谱地图构建方法是各种确定性的和地理统计的空间插值方法,如距离反比法、自然邻点插值法、改进Shepard法和克里金插值法等等。在上述插值方法中,虽然克里金插值法最为复杂,且需要较多的观测值,但因为它是地理统计学上最佳的线性无偏估计,能得到较为准确的结果,故而更为常用。
此外,可用于频谱地图构建的方法还有以下方法:(1)基于压缩感知的方法:考虑到发射机在空间上的稀疏性,将频谱地图构建表述为一个压缩感知问题,然后使用正交匹配追踪方法进行求解;(2)基于字典学习的方法:使用字典学习和压缩感知框架来获得射频功率在空间和时间上的分布状况;(3)基于矩阵补全的方法:将频谱数据表述为一个多维的频谱张量,然后考虑到频谱数据的低秩特性,利用矩阵补全的方法进行求解;(4)基于径向基函数的方法:提出了一种基于自适应高斯径向基函数的方法,将频谱地图构建表述为一个包含最小二乘代价函数和二阶正则化项的优化问题,并使用梯度下降方法进行求解;(5)基于发射机位置估计的方法:提出了一种基于发射机位置估计的方法,首先估计出发射机位置和发射功率,然后选择合适的传播模型来构建频谱地图。
以上方法存在以下缺陷:(1)基于压缩感知的方法、基于字典学习的方法以及基于矩阵补全的方法需要大量的观测值,从而要求感兴趣区域内密集部署大量传感器,这种部署方式耗费巨大,且在实际应用中有时是无法实现的;(2)基于矩阵补全的方法和基于径向基函数的方法仅仅从数学的角度来考虑频谱地图构建问题,而没有考虑其物理含义,没有利用到电波传播的知识,从而导致构建精度不佳;(3)基于字典学习的方法和基于径向基函数的方法中的参数需要精心选择,而这种调参过程会带来需要大量额外的工作量;(4)基于发射机位置估计的方法需要用到关于发射机和传播环境的相关信息,而这些信息在实际应用中是无法获得的。
因此,设计一种新的频谱地图的构建方法具有重要意义。
发明内容
本发明的目的在于提供一种无需参数调整、观测点数量少且具有普遍适用性的频谱地图的构建方法,具体技术方案如下:
一种频谱地图的构建方法,包括以下步骤:
其中:观测值估计值权重Ξ=[1M×1|Ξ′],Ξ′∈RM×I,且矩阵Ξ′的第(m,i)个元素Ξ′mi=(||xm-ci||κ+τ)-1;d取值为2;τ的取值为10-6-10-3;中心点数量I选取1-100个。
其中:xm表示第m(=1,…,M)个频谱传感器的位置,ym表示第m(=1,…,M)个频谱传感器采集获得的关于功率的观测值。
其中:测量值y=[y1,y2,…,yM]T,‖·‖表示向量的2范数;
αj=(ΞTΞ)-1ΞTy 5);
其中:δ表示步长;
应用本发明的技术方案,具有以下效果:(1)本发明方法无需利用关于发射机和传播环境的相关信息,从而应用场景方面更具备普适性;(2)本发明方法在建模过程中考虑到了电波传播的知识,在相同传感器数量和部署情况下能得到更高精度的频谱地图构建结果,或者在相同的频谱地图构建精度条件下,所需的观测点数量更少,从而对传感器数量和部署要求更低,更经济;(3)本发明方法无需进行参数调节,是一种非参数的频谱地图构建方法,更利于实现;(4)本发明使用幂函数的线性组合进行建模(即采用表达式2)对观测值进行拟合),结合使用迭代交替最小化方法对模型中权重、中心点位置和指数的优化求解方法,得到更为精准的频谱地图。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照附图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是真实的频谱地图的示意图;
图2是实施例1所得频谱地图的示意图;
图3是本发明方法与基于自适应高斯径向基函数(RBF)方法的效果比较图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以根据权利要求限定和覆盖的多种不同方式实施。
实施例1:
一种频谱地图的构建方法,具体包括以下步骤:
其中:观测值估计值权重α=[α0,α1,…,αI]T;Ξ=[1M×1|Ξ′],Ξ′∈RM×I,且矩阵Ξ′的第(m,i)个元素Ξ′mi=(||xm-ci||κ+τ)-1;d取值为2;τ的取值为10-6-10-3;中心点数量I选取1-100个(实际应用时,一般选择大于预估的辐射源个数)。
其中:xm表示第m(=1,…,M)个频谱传感器的位置,ym表示第m(=1,…,M)个频谱传感器采集获得的关于功率的观测值。
其中:测量值y=[y1,y2,…,yM]T,||·||表示向量的2范数;
αj=(ΞTΞ)-1ΞTy 5);
其中:δ表示步长;
优选的,所述迭代交替最小化方法的迭代终止条件为:迭代次数j大于预先设置的最大迭代次数Jmax(最大迭代次数可根据实际应用而定),或者足够小或者待优化变量收敛(此处为权重处于收敛、中心点位置处于收敛和衰减指数处于收敛)。
考虑一个50m×50m的区域R,区域内有3个辐射源和M个随机部署的传感器。区域R内传播环境实际的衰减指数为3,真实的频谱地图详见图1,其中:米字符表示传感器位置,正方形符号表示由改进的K-means++算法选取的中心点初始位置,圆圈表示本发明实施例1所提方法求得的中心点最优位置。
采用实施例的方法使用50个传感器(图1中米字符位置)和5个中心点构建的频谱地图详见图2,取值d取2,τ的取值为10-4,中心点数量I取5。对比图1和图2可以看出,所提算法能够有效地优化中心点位置和衰减指数,从而使得构建的频谱地图在视觉上与真实的频谱地图基本相同。
为了量化本发明方法的性能,将其与基于自适应高斯径向基函数(RBF)的方法进行比较,对比结果如图3所示。比较中选取的性能指标是常用的归一化均方误差(NMSE),即
从图3中可以看出,相较于基于RBF的方法,本发明方法能显著提高频谱地图构建精度,在传感器个数为50时,构建精度提高约2dB,且随着传感器个数的增加,构建精度的提升更加明显。性能的提升主要得益于以下两个方面原因:首先,本发明方法中幂函数的线性组合更适合于对接收信号功率的无线电传播效应进行建模;其次,本发明方法中权重、中心点位置和衰减指数的自适应优化通常起着重要的作用。
采用本实施例的技术方案,效果是:(1)本发明方法无需利用关于发射机和传播环境的相关信息,从而应用场景方面更具备普适性;(2)本发明方法在建模过程中考虑到了电波传播的知识,在相同传感器数量和部署情况下能得到更高精度的频谱地图构建结果,或者在相同的频谱地图构建精度条件下,所需的观测点数量更少,从而对传感器数量和部署要求更低,更经济;(3)本发明方法无需进行参数调节,是一种非参数的频谱地图构建方法,更利于实现;(4)本发明使用幂函数的线性组合进行建模,结合使用迭代交替最小化方法对模型中权重、中心点位置和衰减指数的优化求解方法,得到更为精准的频谱地图。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种频谱地图的构建方法,其特征在于:包括以下步骤:
其中:观测值估计值权重α=[α0,α1,…,αI]T;Ξ=[1M×1|Ξ′],Ξ′∈RM ×I,且矩阵Ξ′的第(m,i)个元素Ξ′mi=(||xm-ci||κ+τ)-1;d取值为2;中心点数量I选取1-100个;M个频谱传感器的位置记为d为空间的维数;
其中:测量值y=[y1,y2,…,yM]T,||·||表示向量的2范数;
αj=(ΞTΞ)-1ΞTy 5);
其中:δ表示步长;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810245329.2A CN108566256B (zh) | 2018-03-23 | 2018-03-23 | 一种频谱地图的构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810245329.2A CN108566256B (zh) | 2018-03-23 | 2018-03-23 | 一种频谱地图的构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108566256A CN108566256A (zh) | 2018-09-21 |
CN108566256B true CN108566256B (zh) | 2021-02-02 |
Family
ID=63531937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810245329.2A Active CN108566256B (zh) | 2018-03-23 | 2018-03-23 | 一种频谱地图的构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108566256B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110346654B (zh) * | 2019-07-19 | 2021-03-16 | 中国人民解放军国防科技大学 | 基于普通克里金插值的电磁频谱地图构建方法 |
CN110445567B (zh) * | 2019-08-06 | 2021-06-25 | 中国人民解放军国防科技大学 | 一种电磁频谱地图的构建方法 |
CN112084274A (zh) * | 2020-08-24 | 2020-12-15 | 浙江云合数据科技有限责任公司 | 一种基于数据的时空频谱地图构建和预测方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101572897A (zh) * | 2009-03-24 | 2009-11-04 | 复旦大学 | 分布式认知无线电信道探测方法和系统 |
CN102025427A (zh) * | 2009-09-23 | 2011-04-20 | 华为技术有限公司 | 频谱检测的方法及用户设备与融合设备 |
CN102255675A (zh) * | 2010-05-19 | 2011-11-23 | 索尼公司 | 基于认知无线电的频谱感知装置、方法、程序 |
TWI489897B (zh) * | 2010-12-22 | 2015-06-21 | Htc Corp | 使用壓縮感測建構頻譜圖的方法 |
CN104251991B (zh) * | 2014-09-25 | 2016-11-23 | 中国科学院电子学研究所 | 一种基于稀疏度估计的分维度阈值迭代稀疏微波成像方法 |
CN105281779B (zh) * | 2015-11-04 | 2018-12-14 | 哈尔滨工业大学 | 多观测值向量稀疏度自适应压缩采样匹配追踪方法 |
CN107356819B (zh) * | 2016-05-10 | 2019-11-15 | 中国人民解放军理工大学 | 一种智能家用电磁辐射信息服务系统 |
-
2018
- 2018-03-23 CN CN201810245329.2A patent/CN108566256B/zh active Active
Non-Patent Citations (1)
Title |
---|
Reconstruction of a Gaussian random field with application to spectrum cartography;Geethu Joseph;《2016 IEEE International Conference on Communications (ICC)》;20161231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108566256A (zh) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Deep neural networks for channel estimation in underwater acoustic OFDM systems | |
CN108322409B (zh) | 基于广义正交匹配追踪算法的稀疏ofdm信道估计方法 | |
CN108566256B (zh) | 一种频谱地图的构建方法 | |
US10439594B2 (en) | Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation | |
CN110361778B (zh) | 一种基于生成对抗网络的地震数据重建方法 | |
US20110146371A1 (en) | Method for calibrating a transducer array | |
CN105790813A (zh) | 一种大规模mimo下基于深度学习的码本选择方法 | |
KR102067997B1 (ko) | 딥러닝을 이용한 무선 측위 장치 및 방법 | |
CN108880557B (zh) | 基于压缩感知的稀疏度自适应变步长匹配追踪方法 | |
Lytaev | Nonlocal boundary conditions for split-step padé approximations of the helmholtz equation with modified refractive index | |
Logutov et al. | Inverse barotropic tidal estimation for regional ocean applications | |
Katzfuss et al. | Tutorial on fixed rank kriging (FRK) of CO2 data | |
CN104333424B (zh) | 一种频谱检测与未知噪声方差跟踪估计方法与装置 | |
CN111597753B (zh) | 数据深度变化特征自适应的二维电阻率反演方法及系统 | |
CN109388778A (zh) | 一种迭代容积点无迹卡尔曼滤波方法 | |
Luo et al. | Efficient particle filtering through residual nudging | |
CN116155412A (zh) | 无线信道评估方法及系统 | |
CN106485074B (zh) | 一种基于自适应采样率的海洋温盐场采样方法 | |
CN110686610B (zh) | 基于自适应网格的光学变形测量方法及电子设备 | |
CN114139448A (zh) | 海基观测网站位布局优化的方法、系统、介质、终端及应用 | |
CN113740802B (zh) | 以自适应噪声估计进行矩阵补全的信号源定位方法和系统 | |
CN116559579A (zh) | 基于改进的VMD和Teager能量算子故障定位方法 | |
CN116017280A (zh) | 一种目标免携带设备的快速室内路径跟踪方法 | |
CN113723006B (zh) | 一种基于ls-svm的单台站地球变化磁场建模预测方法及系统 | |
CN112434439B (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 |