CN101631317B - 一种移动通信话务量的多步预测方法 - Google Patents
一种移动通信话务量的多步预测方法 Download PDFInfo
- Publication number
- CN101631317B CN101631317B CN2009100727034A CN200910072703A CN101631317B CN 101631317 B CN101631317 B CN 101631317B CN 2009100727034 A CN2009100727034 A CN 2009100727034A CN 200910072703 A CN200910072703 A CN 200910072703A CN 101631317 B CN101631317 B CN 101631317B
- Authority
- CN
- China
- Prior art keywords
- formula
- wavelet
- telephone traffic
- mobile communication
- communication telephone
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Telephonic Communication Services (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
一种移动通信话务量的多步预测方法,涉及移动通信话务量预测领域,它解决了现有移动通信话务量使用自回归滑动平均系列模型预测精度低以及使用小波分解过程中的盲目性。它的具体实现过程为:采用傅里叶谱分析的结果作为先验知识,对基于极大重叠离散小波变换的分解过程加以引导,提取出与傅里叶谱成分对应的各细节项和趋势项。针对小波分解带来的算法复杂度的增加,利用傅里叶谱先验知识对小波子层进行精简整合,对整合后的趋势项和周期项分别利用乘积季节求和自回归滑动平均模型预测,将预测结果相加,获得预测值。本发明所述方法还适用于居民自来水流量、城市公交车流量、电梯人流量和网络流量等实际多周期时间序列的预测。
Description
技术领域
本发明涉及移动通信话务量预测领域,具体涉及ARMA系列模型预测算法和基于先验知识的极大重叠离散小波变换分解整合算法。
背景技术:
现在移动通信的用户数和话务量保持着高速增长的势头,移动网络长久、稳定的运行,有赖于及时、有效的网络规划与优化。移动通信话务量超过一定容量时,极易造成交换系统过载,出现网络拥塞,给移动通信运营商和用户造成不可弥补的损失。因此,根据话务统计资料和其他业务信息对移动通信话务量的变化趋势进行预测,可以为通讯过程中峰值预警、基站配置、信道优化利用等问题提供决策支持。
话务量作为时间序列的一种,其预测方法可以追溯到时间序列预测的一般方法。时间序列预测就是根据已有的观察记录数据,在一定精度范围内预报未来行为,主要分为线性和非线性方法两类。事实上,基于线性模型的预测方法一般算法实现比较容易,速度快,且可实现多步预测,但缺点是预测精度较低;非线性预测方法预测精度较高,但缺点是算法复杂,且一般只能实现单步或少量多步预测。针对话务量的多步预测需求,需采用线性模型方法——ARMA系列模型方法进行预测。然而,话务量序列作为一种与人们生活作息联系紧密的多周期时间序列,同时又具有多尺度和非平稳特性,线性模型对其适应性较差,预测精度较低。
小波分析在计算机应用、信号处理、图像分析、非线性科学等领域中有着广泛的应用。考虑到小波分析良好的时频局部性和多分辨分析能力,采用小波方法对话务量进行分解,以期降低子层复杂度。
发明内容
本发明为解决现有的移动通信话务量使用ARMA系列模型预测精度低、以及小波分解过程中的盲目性,提出一种移动通信话务量的多步预测方法,该方法由以下步骤完成:
步骤一:对移动通信话务量数据进行傅里叶变换,将所述的傅里叶谱分析结果作为小波分解整合的先验知识;
步骤二:根据步骤一获得的傅里叶谱先验知识,采用极大重叠离散小波变换算法,选取Haar小波基和分解层数,对话务量序列进行小波分解,获得与傅里叶谱成分对应的趋势项与各细节项;
步骤三:根据步骤一获得的先验知识,将步骤二中的各细节项整合为周期项;
步骤四:采用乘积季节ARIMA模型,对步骤二获得的趋势项和步骤三获得的周期项分别进行预测,将两项预测结果相加,获得总的预测值。
本发明的优点:一、本发明采用极大重叠离散小波变换分解算法和乘积季节ARIMA模型相结合,实现了对移动通信话务量序列的多步预测,提高了预测精度;二、以傅里叶谱分析的结果作为话务量序列极大重叠离散小波变换分解的先验知识,解决了小波分解过程的盲目性;三、采用傅里叶谱先验知识作为小波子层精简整合的依据,降低了小波分解的算法复杂度。
附图说明
图1是哈尔滨市某小区话务量序列示意图,图2是哈尔滨市某小区话务量频谱周期成分示意图,图3是哈尔滨市某小区话务量序列小波分解示意图,图4是哈尔滨市某小区话务量序列小波分解精简整合的示意图,图5是采用本发明所述方法对精简整合后的话务量序列多步预测的效果图,图6是直接采用乘积季节ARIMA模型对原始话务量序列多步预测的效果图。
具体实施方式
具体实施方式一:一种移动通信话务量的多步预测方法,该方法由以下步骤完成:
步骤一:对移动通信话务量数据进行傅里叶变换,将所述的傅里叶谱分析结果作为小波分解整合的先验知识;
步骤二:根据步骤一获得的傅里叶谱先验知识,采用极大重叠离散小波变换算法,选取Haar小波基和分解层数,对话务量序列进行小波分解,获得与傅里叶谱成分对应的趋势项与各细节项;
步骤三:根据步骤一获得的先验知识,将步骤二中的各细节项整合为周期项;
步骤四:采用乘积季节ARIMA模型,对步骤二获得的趋势项和步骤三获得的周期项分别进行预测,将两项预测结果相加,获得总的预测值。
本实施方式中步骤一所述的对话务量数据进行频谱分析的过程中,将所述话务量频谱的频率幅度特性转化为周期幅度特性,作为小波分解的先验知识和小波分解后的模型整合的依据。
本实施方式中步骤二所述的采用极大重叠离散小波变换算法实现小波分解的过程为:
设话务量X为N点实序列,X={Xt,t=0,1,...N-1},其中N为话务量建模数据点数。话务量数据X以1小时为采样间隔,采用和对X进行循环滤波,通过公式一和公式二获得J0个基于极大重叠离散小波变换的小波系数和一个尺度系数;
公式一:
公式二:
式中J0为小波分解的层数,为第j层小波滤波器,为第j层尺度滤波器,j=1,2,...,J0,t=0,1,2,...,N-1,Lj=(2j+1)(L-1)+1, Lj表示第j层小波滤波器宽度,t表示时间,表示第J0层尺度滤波器宽度,L为首层尺度滤波器宽度;
公式三:
公式四:
公式五:
公式六:
通过上述公式实现了对话务量数据的极大重叠离散的小波分解,在分解过程中,选取Haar小波基,增大分解层数J0,直到恰好分解出傅里叶谱显示的周期最长的成分,而把其中周期更长、更平稳的成分作为趋势项。
本实施方式所述的步骤四中采用的乘积季节ARIMA模型是ARMA系列模型的一种。对具有趋势特性的时间序列,差分处理后,符合ARMA模型;对具有周期特性的时间序列,季节差分运算后,符合ARMA模型。针对既具有趋势特性,又具有周期特性的话务量序列,选用乘积季节ARIMA模型。
具体实施方式二:本实施方式为具体实施方式一的实施例:
现以黑龙江省哈尔滨市某小区话务量序列为例,并结合图1、图2、图3、图4、图5和图6说明本具体实施方式:
一、对黑龙江省哈尔滨市某小区话务量数据进行傅里叶谱分析,参见图1,所述分析过程为:首先对所述的话务量序列进行等间隔采样,采样周期为1小时,话务量的单位为爱尔兰(Erl);傅里叶谱分析的结果显示话务量时间序列数据在6、8、12、24、84、168小时等频谱成分上具有明显的周期特性,将其周期成分按振幅由大到小的顺序排列,顺次为:T=24小时、T=12小时、T=6小时、T=8小时、T=168小时,T=84小时等频谱成分,参见图2。将所述的频谱成分作为小波分解的先验知识;
二、根据上述傅里叶谱先验知识,采用具体实施方式一中步骤二所述的极大重叠离散小波变换算法,进行小波分解:设X为所述黑龙江省哈尔滨市某小区的话务量序列,将所述黑龙江省哈尔滨市某小区的话务量序列分解为X=S+D1+D2+D3+D4+D5+D6+D7;其中D1、D2、D3、D4、D5、D6、D7为频率较高的细节项,变化剧烈,且频率顺次由高到低;S为频率最低的趋势项,变化较为缓和;参见图3。对各子层分别进行傅里叶谱分析并加以验证,结果显示D7为T=168小时成分;D6为T=84小时成分;D4为T=24小时成分;D3为T=12小时成分;D2为T=6小时成分;
三、对小波分解得到的细节项进行整合,将所有细节项整合为一个周期项D,即:D=D1+D2+D3+D4+D5+D6+D7;此时确定T=168小时为公共周期,总的分解效果为X=D+S,参见图4;
四、对周期项D采用公共周期T=168小时作为季节差分的周期,对趋势项S采取T=1作为季节差分的周期,分别采用乘积季节ARIMA模型进行预测,将两项的预测结果等权相加,获得总的预测值,总的预测效果参见图5。
图6是直接采用乘积季节ARIMA模型的预测效果图;采用三种评价标准对图5和图6进行预测性能的评价:所述三种评价标准分别为:平均绝对误差MAE(Mean Absolute Error)、规范化均方根误差NRMSE(Normalized Root MeanSquare Error)和标准化误差NE(Normalized Error),所述三种评价标准的预测结果参见表1。
表1
结合图5和图6,本发明所述的方法可实现对话务量趋势更为准确的多步预测;由表1可知,采用本发明所述方法比直接采用乘积季节ARIMA模型在上述三种标准下预测误差均降低25%。
本发明针对小波分解的盲目性,提出了以傅里叶谱为先验知识;针对小波分解带来的算法复杂度的增加,提出了基于先验知识的模型精简整合方法。对整合后的各子层分别采用乘积季节ARIMA模型进行预测,从而实现了较高精度范围内的多步预测需求。本发明不仅适用于移动通信话务量预测问题,也适用于和人们生活联系紧密的实际多周期时间序列预测问题。
Claims (3)
1.一种移动通信话务量的多步预测方法,其特征是一种移动通信话务量的多步预测方法由以下步骤完成:
步骤一:对移动通信话务量数据进行傅里叶变换,将傅里叶谱分析结果作为小波分解整合的先验知识;
步骤二:根据步骤一获得的傅里叶谱先验知识,采用极大重叠离散小波变换算法,选取Haar小波基和分解层数,对话务量序列进行小波分解,获得与傅里叶谱成分对应的趋势项与各细节项;采用极大重叠离散小波变换算法实现小波分解的过程为:
话务量X为N点实序列,X={Xt,t=0,1,...N-1},其中N为话务量建模数据点数;话务量数据X以1小时为采样间隔,采用和对X进行循环滤波,通过公式一和公式二获得J0个基于极大重叠离散小波变换的小波系数和一个尺度系数;
公式一:
公式二:
式中J0为小波分解的层数,为第j层小波滤波器,为第j层尺度滤波器,j=1,2,...,J0,t=0,1,2,...,N-1,Lj=(2j+1)(L-1)+1,Lj表示第j层小波滤波器宽度,表示第J0层尺度滤波器宽度,L为首层尺度滤波器宽度;
公式三:
公式四:
公式五:
公式六:
步骤三:根据步骤一获得的先验知识,将步骤二中的各细节项整合为周期项;
步骤四:采用乘积季节ARIMA模型,对步骤二获得的趋势项和步骤三获得的周期项分别进行预测,将两项预测结果相加,获得总的预测值。
2.根据权利要求1所述的一种移动通信话务量的多步预测方法,其特征在于所述的步骤一中对话务量数据进行傅里叶变换的过程中,将话务量频谱的频率幅度特性转化为周期幅度特性。
3.根据权利要求1所述的一种移动通信话务量的多步预测方法,其特征在于步骤四中采用的乘积季节ARIMA模型是ARMA系列模型的一种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100727034A CN101631317B (zh) | 2009-08-19 | 2009-08-19 | 一种移动通信话务量的多步预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100727034A CN101631317B (zh) | 2009-08-19 | 2009-08-19 | 一种移动通信话务量的多步预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101631317A CN101631317A (zh) | 2010-01-20 |
CN101631317B true CN101631317B (zh) | 2011-12-28 |
Family
ID=41576210
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100727034A Expired - Fee Related CN101631317B (zh) | 2009-08-19 | 2009-08-19 | 一种移动通信话务量的多步预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101631317B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102088709A (zh) * | 2010-11-30 | 2011-06-08 | 哈尔滨工业大学 | 基于聚类和arima模型的话务量预测方法 |
CN103678514B (zh) * | 2013-11-26 | 2017-01-18 | 科大讯飞股份有限公司 | 一种业务趋势预测方法及系统 |
CN105577402B (zh) * | 2014-10-11 | 2019-05-24 | 北京通达无限科技有限公司 | 基于历史数据的业务异常监控方法及设备 |
DE102017221576A1 (de) * | 2017-11-30 | 2019-06-06 | Robert Bosch Gmbh | Verfahren zur Mittelung von pulsierenden Messgrößen |
CN109951243A (zh) * | 2017-12-20 | 2019-06-28 | 中国科学院深圳先进技术研究院 | 一种频谱预测方法、系统及电子设备 |
CN108173686A (zh) * | 2017-12-26 | 2018-06-15 | 北京工业大学 | 一种基于arima和小波变换的云数据中心请求流预测方法 |
CN108111353B (zh) * | 2017-12-26 | 2021-10-15 | 深圳广联赛讯股份有限公司 | 预付卡剩余流量预测方法、网络终端和存储介质 |
CN108389631A (zh) * | 2018-02-07 | 2018-08-10 | 平安科技(深圳)有限公司 | 水痘发病预警方法、服务器及计算机可读存储介质 |
CN115292304A (zh) * | 2022-10-09 | 2022-11-04 | 深圳市明源云科技有限公司 | 时序数据的异常检测方法、设备及计算机可读存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101247612A (zh) * | 2007-02-13 | 2008-08-20 | 中兴通讯股份有限公司 | 可视电话话务量预测方法及容量规模估计方法 |
CN101304590A (zh) * | 2008-04-14 | 2008-11-12 | 中国网络通信集团公司 | 移动通信网络的无线网络容量确定装置及方法 |
-
2009
- 2009-08-19 CN CN2009100727034A patent/CN101631317B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101247612A (zh) * | 2007-02-13 | 2008-08-20 | 中兴通讯股份有限公司 | 可视电话话务量预测方法及容量规模估计方法 |
CN101304590A (zh) * | 2008-04-14 | 2008-11-12 | 中国网络通信集团公司 | 移动通信网络的无线网络容量确定装置及方法 |
Non-Patent Citations (2)
Title |
---|
冯遂亮 等.《利用小波方差进行原子钟频率稳定度的估计》.《宇航计测技术》.2009,第29卷(第1期),全文. * |
刘大同 等.《基于残差预测修正的局部在线时间序列预测方法》.《电子学报》.2008,第36卷(第12A期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN101631317A (zh) | 2010-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101631317B (zh) | 一种移动通信话务量的多步预测方法 | |
Xu et al. | Understanding mobile traffic patterns of large scale cellular towers in urban environment | |
Diao et al. | A hybrid model for short-term traffic volume prediction in massive transportation systems | |
CN101604356B (zh) | 一种不确定性中长期水文预报模型的建立方法 | |
CN109587350A (zh) | 一种基于滑动时间窗口聚合的电信诈骗电话的序列异常检测方法 | |
CN107832709A (zh) | 基于经验小波变换进行特征提取的脑电情感识别方法 | |
Zhong et al. | A frequency domain approach to characterize and analyze load profiles | |
CN103646354A (zh) | 基于有效指标fcm和rbf神经网络的变电站负荷特性分类方法 | |
Mosavat-Jahromi et al. | Prediction and modeling of spectrum occupancy for dynamic spectrum access systems | |
CN102869033A (zh) | 基于arma-arch模型族的gsm通信系统gprs数据业务的预测方法 | |
CN111062539B (zh) | 基于次级用电量特性聚类分析的总电量预测方法 | |
CN109004997A (zh) | 一种频谱感知方法、装置和计算机可读存储介质 | |
Shi et al. | Discovering periodic patterns for large scale mobile traffic data: Method and applications | |
CN104268651A (zh) | 基于小波多尺度三次指数平滑模型的季节性能耗数据预测方法 | |
Scellato et al. | On nonstationarity of human contact networks | |
CN117676385A (zh) | 横向聚类lstm联邦学习下充电站负荷预测系统及方法 | |
CN109325635B (zh) | 一种基于自动补全的位置预测方法 | |
CN102065449A (zh) | 一种基于聚类ls-svm的移动通信话务量预测方法 | |
CN104881704A (zh) | 一种话务指标的预测方法、装置及电子设备 | |
CN103428320B (zh) | 基于模型拟合的智能手机功耗评估方法 | |
CN115130788B (zh) | 基于大数据分析的用电预测系统 | |
CN111194045A (zh) | 一种基于用户群体聚集行为模型进行节能的方法 | |
Wilson et al. | Periodicity in the autocorrelation function as a mechanism for regularly occurring zero crossings or extreme values of a Gaussian process | |
CN103441806B (zh) | 认知无线电的纯不连续马尔可夫过程频谱感知方法 | |
CN101630962B (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111228 Termination date: 20120819 |