CN108282805A - 一种cdma2000基站电磁辐射预测方法 - Google Patents
一种cdma2000基站电磁辐射预测方法 Download PDFInfo
- Publication number
- CN108282805A CN108282805A CN201810071024.4A CN201810071024A CN108282805A CN 108282805 A CN108282805 A CN 108282805A CN 201810071024 A CN201810071024 A CN 201810071024A CN 108282805 A CN108282805 A CN 108282805A
- Authority
- CN
- China
- Prior art keywords
- models
- prediction
- model
- cdma2000
- predicted 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W24/00—Supervisory, monitoring or testing arrangements
- H04W24/06—Testing, supervising or monitoring using simulated traffic
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/147—Network analysis or design for predicting network behaviour
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种CDMA2000基站电磁辐射预测方法,其步骤如下:取CDMA2000基站电磁辐射历史数据作为模型训练数据,将训练数据通过预测模型进行训练来确定三个预测模型的参数,再将训练好的三个预测模型的预测值的变化趋势与对应时段话务量变化趋势进行相关性分析,得出相关系数,再根据相关系数确定三个预测模型的组合权重,得到权重系数,再用权重系数与三个预测模型相结合对CDMA2000基站电磁辐射做组合预测。本发明分析了CDMA2000基站电磁辐射预测的方式,该方法能够准确有效的进行长时间的对CDMA2000基站电磁辐射进行预测,具有一定的社会效益。
Description
技术领域
本发明涉及一种CDMA2000基站电磁辐射预测方法。
背景技术
目前针对通信基站周围的电磁辐射预测时,常运用训练好的预测模型进行预测。例如,申请号2018100095052公布了一种GSM基站电磁辐射预测方法,在基站电磁辐射的预测中,先对预测模型进行训练,然后使用训练好的模型进行预测,但在这种预测中,若进行长时间预测时,其预测精度会不断下降,基站电磁辐射变化与话务量的变化息息相关,目前在已公开的文献中,没有考虑通过话务量变化趋势与预测变化趋势的关系来提升预测效果。
针对现有技术中存在的不足,本专利提出将训练数据通过预测模型进行训练来确定三个预测模型的参数,再将训练好的三个预测模型的预测值的变化趋势与对应时段话务量变化趋势进行相关性分析,得出相关系数,再根据相关系数确定三个预测模型的组合权重,得到权重系数,再用权重系数与三个预测模型相结合对CDMA2000基站电磁辐射做组合预测。通过实验表明,本专利提出的预测方法能准确有效的进行长时间的对CDMA2000基站电磁辐射进行预测。
发明内容
为了解决上述技术问题,本发明提供一种CDMA2000基站电磁辐射预测方法,能准确有效的进行长时间的对CDMA2000基站电磁辐射进行预测。
本发明解决上述技术问题的技术方案包括以下步骤:
1)、取CDMA2000基站电磁辐射历史数据作为训练数据,将训练数据通过三个预测模型进行训练,三个预测模型分别为:A模型、P模型、M模型,并确定三个预测模型的参数,其中A模型表达式为:
yA(t)=d1y(t-1)+d2y(t-2)+…+dny(t-n) (1)
在上式(1)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数d1,d2,...,dn;
P模型表达式为:
在上式(2)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数ε,θ;
M模型表达式为:
在上式(3)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数c;
2)、将预测时的输入数据,分别输入步骤1训练确定的A模型、P模型、M模型,获得预测值yA(t),yA(t+1),...,yA(t+k)、yP(t),yP(t+1),...,yP(t+k)、yM(t),yM(t+1),...,yM(t+k);
3)、根据步骤2得到的单支预测值yA(t),yA(t+1),...,yA(t+k)、yP(t),yP(t+1),...,yP(t+k)、yM(t),yM(t+1),...,yM(t+k)的变化趋势分别与对应时段话务量变化趋势做相关性分析,求得三组预测模型的预测值变化趋势与对应时段话务量变化趋势的相关系数lA,lp,lM相关系数表达式为:
其中,lA、lp、lM为A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数,h(t+k)为预测点对应的话务量,k=0,1,…,n,单位为Erl,为n+1个预测点对应的话务量平均值,单位为Erl,yA(t+k)为通过A模型做的单支预测值,k=0,1,…,n,单位为V/m,为A模型所做n+1个预测值的平均值,单位为V/m,yp(t+k)为通过P模型做的单支预测值,k=0,1,…,n,单位为V/m,为P模型所做n+1个预测值的平均值,单位为V/m,yM(t+k)为通过M模型做的单支预测值,k=0,1,…,n,单位为V/m,为M模型所做n+1个预测值的平均值,单位为V/m;
4)、根据步骤3得到相关系数求得A模型、P模型、M模型的权重系数QA,QP,QM,根据权重系数对A模型、P模型、M模型的预测值做组合预测,表达式为:
y(ttk)=QAyA(t+k)+QPyP(t+k)+QMyM(t+k) (7)
其中,y(t+k)为模型的组合预测值,k=0,1,…,n,单位为V/m。
上述的一种CDMA2000基站电磁辐射预测方法,所述步骤4)中,其特征在于,A模型、P模型、M模型的权重系数QA,QP,QM的计算表达式为:
其中,QA,QP,QM分别为A模型、P模型、M模型做组合预测的权重系数,lA、lp、lM为A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数。
本发明的有益效果在于:通过CDMA2000基站电磁辐射训练数据,对预测模型进行训练来确定三个预测模型的参数,再将训练好的三个预测模型的预测值的变化趋势与对应时段话务量变化趋势进行相关性分析,得出相关系数,再根据相关系数确定三个预测模型的组合权重,得到权重系数,再用权重系数与三个预测模型相结合对CDMA2000基站电磁辐射做组合预测。所建立的模型能准确有效的进行长时间的对CDMA2000基站电磁辐射进行预测。该方法对基站建设和环境保护有极较大的参考价值,具有一定的社会效益。
附图说明
图1为本发明的流程图。
具体实施方式
本实施例在以本发明内容为前提下进行,给出了详细的实施步骤,但本发明的保护范围不限于下述的实施例。
本次实施的实验地在空旷平坦的区域,实施对象为校教学楼屋顶基站,实验仪器为频谱分析仪AT6030D,测量对象为CDMA2000基站。
下面结合附图和实施例对本发明做进一步的说明。
图1为本发明技术方案的框图,具体步骤为:
步骤一:取CDMA2000基站电磁辐射历史数据作为训练数据,将训练数据通过三个预测模型进行训练,三个预测模型分别为:A模型、P模型、M模型,并确定三个预测模型的参数,其中A模型表达式为:
yA(t)=d1y(t-1)+d2y(t-2)+…+dny(t-n) (11)
在上式(11)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数d1,d2,...,dn;
P模型表达式为:
在上式(12)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数ε,θ;
M模型表达式为:
在上式(13)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数c。
在本次实施例中,用CDMA2000的历史数据对A模型进行训练,在模型中n取值为9,用历史数据A1,f=[0.1292,0.1325,0.1362,0.1388,0.1365,0.1314,0.1267,0.1292,0.1295],A2,f=[0.1325,0.1362,0.1388,0.1365,0.1314,0.1267,0.1292,0.1295,0.1304],...,A9,f=[0.1295,0.1304,0.1358,0.1279,0.1311,0.1264,0.1278,0.1235,0.1268],f=1,2,...,9,作为输入,a1=[0.1304],a2=[0.1358],...,a9=[0.1302],作为对应的输出,对模型A进行训练,获得d1,d2,...,d9的值,其中d1=0.1357,d2=0.0985,d3=0.1217,d4=0.1112,d5=0.1001,d6=0.0963,d7=0.1537,d8=0.1221,d9=0.08751;用CDMA2000的历史数据对P模型进行训练,在模型中n取值为8,用历史数据P1,9=[0.1292,0.1304,0.1358,0.1279,0.1311,0.1264,0.1278,0.1235,0.1268],P2,9=[0.1295,0.1304,0.1358,0.1279,0.1311,0.1264,0.1278,0.1235,0.1268],作为输入,p1=[0.1235],p2=[0.1302],作为对应的输出,对模型P进行训练,获得ε,θ的值,其中ε=6.091,θ=0.10889;用CDMA2000的历史数据对M模型进行训练,在模型中n取值为9,用历史数据M1,9=[0.1292,0.1304,0.1358,0.1279,0.1311,0.1264,0.1278,0.1235,0.1268],M2,9=[0.1295,0.1304,0.1358,0.1279,0.1311,0.1264,0.1278,0.1235,0.1268],作为输入,m1=[0.1235],m2=[0.1302],作为对应的输出,对模型M进行训练,获得c的值,其中c=970.8180414。
步骤二:将预测时的输入数据,分别输入步骤1训练确定的A模型、P模型、M模型,获得预测值yA(t),yA(t+1),...,yA(t+k)、yP(t),yP(t+1),...,yP(t+k)、yM(t),yM(t+1),...,yM(t+k)。
在本次实施例中,将预测点前的历史数据作为模型预测时的输入数据,分别输入步骤1)中训练好的A模型,P模型,M模型进行预测获得预测值,然后再以上一次的预测值与之前的历史数据作为输入,获得下一次的预测值,以此循环获得后续14个时间点的预测值,yA(t),yA(t+1),...,yA(t+13)、yP(t),yP(t+1),...,yP(t+13)、yM(t),yM(t+1),...,yM(t+13),其中预测值如下表所示:
表1单支预测值
步骤三:根据步骤2得到的单支预测值yA(t),yA(t+1),...,yA(t+k)、yP(t),yP(t+1),...,yP(t+k)、yM(t),yM(t+1),...,yM(t+k)的变化趋势分别与对应时段话务量变化趋势做相关性分析,求得三组预测模型的预测值变化趋势与对应时段话务量变化趋势的相关系数lA,lp,lM相关系数表达式为:
其中,lA、lp、lM为A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数,h(t+k)为预测点对应的话务量,k=0,1,…,n,单位为Erl,为n+1个预测点对应的话务量平均值,单位为Erl,yA(t+k)为通过A模型做的单支预测值,k=0,1,…,n,单位为V/m,为A模型所做n+1个预测值的平均值,单位为V/m,yP(t+k)为通过P模型做的单支预测值,k=0,1,…,n,单位为V/m,为P模型所做n+1个预测值的平均值,单位为V/m,yM(t+k)为通过M模型做的单支预测值,k=0,1,…,n,单位为V/m,为M模型所做n+1个预测值的平均值,单位为V/m。
在本次实施例中,将单支预测值yA(t+k),yP(t+k),yM(t+k)与对应时段话务量的值h(t+k),输入上式获得相关系数lA,lp,lM其中对应时段的话务量的值如下表所示;
表二对应时段的话务量
k | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 |
h(t+k) | 6.01236 | 6.110205 | 6.17664 | 6.31557 | 6.16482 | 6.00678 | 6.049575 | 6.186435 | 6.239175 | 6.15342 | 5.96187 | 5.90589 | 6.01125 | 6.0548 |
最终求得A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数lA=0.8325,lp=0.8266,lM=0.7811。
步骤四:根据步骤3得到相关系数求得A模型、P模型、M模型的权重系数QA,QP,QM,根据权重系数对A模型、P模型、M模型的预测值做组合预测,表达式为:
y(t+k)=QAyA(t+k)+QPyP(t+k)+QMyM(t+k) (17)
其中,y(t+k)为模型的组合预测值,k=0,1,…,n,单位为V/m。
在本次实施例中,先将步骤三求得的相关系数lA,lp,lM代入QA,QP,QM的计算表达式求得三个模型的组合权重:
其中,QA,QP,QM分别为A模型、P模型、M模型做组合预测的权重系数,lA、lp、lM为A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数,最终求得组合权重QA=0.3412,QP=0.3387,QM=0.3201,再将求得的组合权重对三个预测模型做组合预测,其中预测值如下表所示:
表3预测值y(t+k)与测量值的展示
K | 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 |
y(t+k) | 0.133557 | 0.135735 | 0.137233 | 0.140389 | 0.136984 | 0.133448 | 0.134379 | 0.137328 | 0.138558 | 0.136751 | 0.132547 | 0.1313 | 0.133553 | 0.13452 |
测量值 | 0.13386 | 0.1351 | 0.137275 | 0.139935 | 0.13718 | 0.133995 | 0.134425 | 0.13709 | 0.138175 | 0.136315 | 0.13285 | 0.131245 | 0.13396 | 0.13513 |
从实验结果可以看出,预测值与测量值在预测时间k不断增加时,预测值与测量值的都比较相近,说明利用此方法能准确有效的进行长时间的对CDMA2000基站电磁辐射进行预测,同时实验结果验证了本发明所使用方法的有效性。
Claims (2)
1.一种CDMA2000基站电磁辐射预测方法,其特征在于,包括以下步骤:
1)、取CDMA2000基站电磁辐射历史数据作为训练数据,将训练数据通过三个预测模型进行训练,三个预测模型分别为:A模型、P模型、M模型,并确定三个预测模型的参数,其中A模型表达式为:
yA(t)=d1y(t-1)+d2y(t-2)+…+dny(t-n) (1)
在上式(1)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数d1,d2,...,dn;
P模型表达式为:
在上式(2)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数ε,θ;
M模型表达式为:
在上式(3)中,将历史数据输入预测模型进行训练,获得最佳预测模型参数c;
2)、将预测时的输入数据,分别输入步骤1训练确定的A模型、P模型、M模型,获得预测值yA(t),yA(t+1),...,yA(t+k)、yP(t),yP(t+1),...,yP(t+k)、yM(t),yM(t+1),...,yM(t+k);
3)、根据步骤2得到的单支预测值yA(t),yA(t+1),...,yA(t+k)、yP(t),yP(t+1),...,yP(t+k)、yM(t),yM(t+1),...,yM(t+k)的变化趋势分别与对应时段话务量变化趋势做相关性分析,求得三组预测模型的预测值变化趋势与对应时段话务量变化趋势的相关系数lA,lp,lM相关系数表达式为:
其中,lA、lp、lM为A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数,h(t+k)为预测点对应的话务量,k=0,1,…,n,单位为Erl,为n+1个预测点对应的话务量平均值,单位为Erl,yA(t+k)为通过A模型做的单支预测值,k=0,1,…,n,单位为V/m,为A模型所做n+1个预测值的平均值,单位为V/m,yP(t+k)为通过P模型做的单支预测值,k=0,1,…,n,单位为V/m,为P模型所做n+1个预测值的平均值,单位为V/m,yM(t+k)为通过M模型做的单支预测值,k=0,1,…,n,单位为V/m,为M模型所做n+1个预测值的平均值,单位为V/m;
4)、根据步骤3得到相关系数求得A模型、P模型、M模型的权重系数QA,QP,QM,根据权重系数对A模型、P模型、M模型的预测值做组合预测,表达式为:
y(t+k)=QAyA(t+k)+QPyP(t+k)+QMyM(t+k) (7)
其中,y(t+k)为模型的组合预测值,k=0,1,…,n,单位为V/m。
2.如权利要求1所述的一种CDMA2000基站电磁辐射预测方法,所述步骤4)中,其特征在于,A模型、P模型、M模型的权重系数QA,QP,QM的计算表达式为:
其中,QA,QP,QM分别为A模型、P模型、M模型做组合预测的权重系数,lA、lp、lM为A模型、P模型、M模型预测值的变化趋势与对应时段话务量变化趋势的相关系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810071024.4A CN108282805B (zh) | 2018-01-25 | 2018-01-25 | 一种cdma2000基站电磁辐射预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810071024.4A CN108282805B (zh) | 2018-01-25 | 2018-01-25 | 一种cdma2000基站电磁辐射预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108282805A true CN108282805A (zh) | 2018-07-13 |
CN108282805B CN108282805B (zh) | 2021-04-13 |
Family
ID=62805020
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810071024.4A Active CN108282805B (zh) | 2018-01-25 | 2018-01-25 | 一种cdma2000基站电磁辐射预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108282805B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101588401A (zh) * | 2008-05-22 | 2009-11-25 | 中兴通讯股份有限公司 | 一种电磁辐射量的统计计算方法 |
CN103076505A (zh) * | 2012-12-27 | 2013-05-01 | 广东省辐射防护协会 | Td-scdma移动通信基站环境电磁辐射三维空间预测方法 |
CN103179574A (zh) * | 2011-12-20 | 2013-06-26 | 上海无线通信研究中心 | 基于地理位置与定位精度的信道分类感知方法 |
EP1946574B1 (en) * | 2005-09-27 | 2014-06-04 | Telecom Italia S.p.A. | Method and system for estimating traffic distribution in a cellular mobile radio communications network |
CN103874090A (zh) * | 2014-03-31 | 2014-06-18 | 湘潭大学 | 一种gsm通信基站电磁辐射预测方法 |
CN105592167A (zh) * | 2016-03-11 | 2016-05-18 | 北京波尔通信技术股份有限公司 | 一种基站电磁辐射监测*** |
-
2018
- 2018-01-25 CN CN201810071024.4A patent/CN108282805B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1946574B1 (en) * | 2005-09-27 | 2014-06-04 | Telecom Italia S.p.A. | Method and system for estimating traffic distribution in a cellular mobile radio communications network |
CN101588401A (zh) * | 2008-05-22 | 2009-11-25 | 中兴通讯股份有限公司 | 一种电磁辐射量的统计计算方法 |
CN103179574A (zh) * | 2011-12-20 | 2013-06-26 | 上海无线通信研究中心 | 基于地理位置与定位精度的信道分类感知方法 |
CN103076505A (zh) * | 2012-12-27 | 2013-05-01 | 广东省辐射防护协会 | Td-scdma移动通信基站环境电磁辐射三维空间预测方法 |
CN103874090A (zh) * | 2014-03-31 | 2014-06-18 | 湘潭大学 | 一种gsm通信基站电磁辐射预测方法 |
CN105592167A (zh) * | 2016-03-11 | 2016-05-18 | 北京波尔通信技术股份有限公司 | 一种基站电磁辐射监测*** |
Non-Patent Citations (1)
Title |
---|
何玮珊: "基于小波变换组合模型的忙时话务量预测", 《中国优秀硕士学位论文全文数据库信息科技辑(月刊 )》 * |
Also Published As
Publication number | Publication date |
---|---|
CN108282805B (zh) | 2021-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111651941B (zh) | 一种全球电离层电子总含量预测的算法 | |
CN109580003B (zh) | 一种静止气象卫星热红外数据估算近地面大气温度方法 | |
CN107590565A (zh) | 一种构建建筑能耗预测模型的方法及装置 | |
CN110555561A (zh) | 一种中长期径流集合预报方法 | |
CN102914623B (zh) | 一种湿度传感器温度补偿的融合方法 | |
CN105259318B (zh) | 一种基于气象参数的恶臭ou值预测方法和*** | |
CN111737912B (zh) | 一种基于深度神经网络的mwhts模拟亮温计算方法 | |
CN104778298A (zh) | 基于egmm的高斯过程回归软测量建模方法 | |
CN103679210B (zh) | 基于高光谱图像解混的地物识别方法 | |
CN102628938B (zh) | 基于噪声先验的联合高斯模型雷达目标稳健识别方法 | |
CN105653502A (zh) | 一种基于遗传算法的通信基站电磁辐射相关性分析方法 | |
CN104048676B (zh) | 基于改进粒子滤波的mems陀螺随机误差补偿方法 | |
CN104504475A (zh) | 基于ar*-svm混合建模的雾霾时间序列预测方法 | |
CN110595960B (zh) | 一种基于机器学习的pm2.5浓度遥感估算方法 | |
CN107688863A (zh) | 自适应迭代加强的短期风速高精度组合预测方法 | |
CN110595968A (zh) | 一种基于静止轨道卫星的pm2.5浓度估算方法 | |
CN111105068A (zh) | 基于序回归学习的数值模式订正的方法 | |
CN103235890A (zh) | 卫星短时临近降水预报***及降水预报方法 | |
CN117933095B (zh) | 一种基于机器学习的地表发射率实时反演及同化方法 | |
CN108183754B (zh) | 一种gsm基站电磁辐射预测方法 | |
CN109992748B (zh) | 一种河口潮位调和分析方法 | |
CN103093425A (zh) | 一种气动热辐射指纹库的建立方法及其应用 | |
Zheng et al. | Improved iterative prediction for multiple stop arrival time using a support vector machine | |
CN108282805A (zh) | 一种cdma2000基站电磁辐射预测方法 | |
CN117391139A (zh) | 基于改进UNet神经网络的天气现象预报订正方法 |
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 |