- 1、本文档共10页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
.
.
雨量预报方法的评价模型
李泳,易勋,贺望香
1. 问题分析
雨量预报对农业生产和城市工作和生活有重要作用,但准确、及时地对雨量作出预报是一个十分困难的问题。
针对问题一,气象部门提供了41天实测数据并且希望建立数学模型来评价两种雨量预报方法的准确性。通过分析数据可知,网格点数据代表的是53*47的网格点的预测雨量,而实测数据是以经纬度定位的观测点的实测雨量。很显然无法对两个定位方式不同的数据进行插值。所以我们选取与对应观测站点最近的5个网格点,通过计算网格点对对应观测站点的权重和网格点的预测雨量的乘积的和求出观测站点的预测雨量。然后就可以通过对两种方法的预测雨量与实测雨量进行比较得出哪种方法更准确。
由于两种方法与实测雨量很相近,很难对比出两种方法的准确性,所以我们采用预报偏差率来比较两种6小时雨量预报方法的准确性。
2. 模型假设
(1)观测站点之间距离的设置是不同的。
(2)雨量用毫米做单位,小于0.1毫米视为无雨
3. 符号说明
:表示第n个网格点到第i个观测站点的距离;(其中n=1,2,3,4,5; i=1-91)
:表示与第i个观测站点的距离最小的前5个网格点的权重;
:表示第i个网格点分别在某月某日某个时段的雨量值;(其中n=1,2,3,4,5)
:表示第i个观测点第j个时段的预测雨量值;
:表示第i个观测点第j个时段的实测雨量值。
4. 模型的建立与求解
4.1.1 首先,找出网格点与观测站点的散点图.从图中可以看出,观测站(红点)与网格点(蓝点)分布不均匀。且观测站点分布在网格点中间部分。
4.1.2 筛选出每一个站点周围的5个距离观测站点最近的网格点。5个网格点的选取通过先给定最大与最小值的范围,然后利用matlab的find()函数找出符合筛选条件的全部点。(见附录2)
4.1.3 运用欧拉公式求出符合条件的5个网格点到观测站点的距离的集合,并对其进行从小到大排序。
欧拉距离公式:
4.1.4 运用欧拉的倒数加权法计算出5个网格点到观测站的权重[1]
计算公式:
4.1.5 计算观测点的雨量预测值,取得每一个5个网格某月某日某时段对应的预测雨量值乘以其对应的权重即等于各个观测点各个时段的预测雨量值。
计算公式:
4.1.6 已经得到了两种方法的观测点的预测雨量,接下来则与实测数据比较,分析哪种方法的准确性高。
? 通过比较方法1与方法2的4个时段的预测雨量的平均值(见附录四)
预测雨量的平均值
第一时段
第二时段
第三时段
第四时段
方法一
1.2112
1.3588
3.9047
2.5623
方法二
2.4127
3.1605
2.7383
2.6722
实测雨量的平均值
第一时段
第二时段
第三时段
第四时段
实测
1.4827
2.1179
1.8007
1.6362
? 通过对上面两个表格的分析可得平均值无法得到准确评估,故采用求预报偏差率的方法来判断两种方法的准确性。
计算公式:
? 由于最终结果依然是一个164*91矩阵,所以对其求和之后再对他进行求列平均得到一个行向量,然后再对行向量平均得到一个均值
通过matlab程序,最终得到:
第一种方法预报偏离的平均值为 124.3207
第二种方法预报偏离的平均值为 124.4289
5. 参考文献
[1]《雨量预报数学模型》论文
/link?url=5VUoL8uZ-59EGbjdykKMlqXszLvVTv2ykGqwiEYdidwnmLFY-xZDOFIHGWr05jkT1YD4BGb5ENqkB6tTeKStB-NmelJmwvsqiLb-oqC0pKaqq-pf-to=pcqq.discussion
[2] 欧拉, 《基于欧式距离的最近邻改进算法》,知网,/Article/CJFDTOTAL-GXKX201004006.html
[3] H. P. Williams, 《数学规划模型建立与计算机应用》,北京,国防工业出版社,1991年;
六 附录
附录1
x=lat;
y=lon;
xi=A020618(:,2); %%为了不让编译器报错,故将文件名第一个改为A开头
yi=A020618(:,3);
plot(x,y,.b,xi,yi,.r)
axis([27.3 35.3 116.8 125]);
xlabel(经度),ylabel(纬度),title(观测站点与网格点分布图)
附录2
x=lat;
y=lon;
xi=A020618(:,2); %为了不让编译报错,故将文件改名为A字母开头
yi=A020618(:,3);
max_x1=max(xi);
min_x2=min(xi);
max_y1=max(yi);
min_y2=m
文档评论(0)