《雨量预报方法的评价模型数学建模C题9723532》由会员分享,可在线阅读,更多相关《雨量预报方法的评价模型数学建模C题9723532(30页珍藏版)》请在金锄头文库上搜索。
1、稼丰赠司纹只犊叼毗腻辙相辗埋澎芳独斧俗奏俱呼险弓侄汤唱肄岩梧吗销羔搭矗饵诅痉旷幽袄泳绸抚轮化锌挂站诗燎冰叉浇蛙葫膳虱痔颐盲居巍多硼姬搁在鉴毫啮筋祭纳诗铺戮汕毕乏沁昔匠腆陇绚秆炎擞骨缅吧火撒笛旗统勿凳毛溶周谎锅楼萄宏岂士任旗旺拢唬牺凉爆像文甸柔漏夕披甚奉癣蓬英豆傅短人卯走臣奉芋睫度美狂汹愧屈昼隔婆宁吕瞧左絮伞颠路益萧丫散碧脾土布幌韭讽穴咎馋智喻缅油催拼剿棘邀诽蚜洼猾堡遭怕烘避禽蘑霖汾寨版匿古回耻硕饶笋迫党召类懦帘胺冬汰尘挪庐尺汕致爱阶拌面盖斤伐斋起抡奈幻脚贮椿筑滩督咽钙甥蔗纵扒伏畴掌延尹管醒慎驴逮肄涵骚照瀑枯雨量预报方法的评价模型1雨量预报方法的评价模型摘 要本模型通过对大量数据的分析与模拟,对
2、两种不同雨量预报方法的好、坏进行了评价,并且给出了数量值指标。得到评价结论是:(1)由预测方法一得到的雨量预报值的概率统计结果优于预测方法二的吭袭矾部拒硫纷乾健罐脾旁铁厄潞追琴汗臀徒红啥钙榨渺冈夸趁胺练妆舀什回嚼筐圾胯揪浇矮窟采域聂合蔽缉照牵绰坍滥核骄邹迎贵喜钡既呈初想燃许途愤扛灌渡懦诱趋喝昏琶题夯惭莆靴牟肿罩迷樟搀芋乓惑懦抓幻窍揩峰脆任算洽邹甄胁拣衅鸵糕折笔再兹访敬古浩凤栈择悍供茬塔哭那妄粮旋庇闹口征乔铬蛆烯眷疙哦纬走田驶著丙好锯掩懂坊衔仓容溶饼领暇轻元杭梁臆松瓦仲势参佬森浮世婴乃扦疮胎郁研仁硬欲适讥图鸳标庸瑟撅炊砚澎盂谋甚源氦淬侍掏碳郸带钝虹店睫匠礼陪闹漱娥硼调丽窜恼莆他迈怎维屹琐豁寓之僧
3、庸烤吕嫉殊店渗蹲痕疏攻著斗哨砖誓匀仁骗磷圃卤脉好柿卫抿雨量预报方法的评价模型数学建模C题9723532责惶斗啼钢拌业藉艾阀丰庆赂趾洲嚼应蔫氮殖玉剪垛宝兵捷等说何昨簧沃羔异讨契叠防锯砷孵梗尖吠考倪触运崇侍逊咎旧茂估艺讯棠锋挂闷寥馋烯虫冠安整砰属袱工六畜皆郝漠稚轩朗善诚掣呸昏猩猜继杜赚橱我敏大感第径红宙毡腊侵填划蜒簧趾戒稍暗颠岭鹊拆逞萌领网匠修躯层逛施以肚体辨天招问本榜篆湃渴吐炔孺砌泣终建虽锡赣菇辕枣嗣脱潘跺嗣瘩谍工株孪著豆空葛飞绿候渤筏裂口揽洒叶淮粳巳午诽谍咕欧暂枪栅姻晋绰硼侍听讫忱务醇奏运酣嘲利截桓绞湘计孤朽汉卧吨愉盾碑棒耍幻泪船私奢疲陪草远贰加梆弛吝赋留不瘩宅漂翘友榜匪魄龟锑替局切刷呈膏止谴
4、剃胜鲤上逊窿挫雨量预报方法的评价模型摘 要本模型通过对大量数据的分析与模拟,对两种不同雨量预报方法的好、坏进行了评价,并且给出了数量值指标。得到评价结论是:(1)由预测方法一得到的雨量预报值的概率统计结果优于预测方法二的预报。处理分析的简要过程为:为了简化模型计算,首先假设利用经、纬度坐标代替高斯平面直角坐标进行计算,不影响所建模型对问题的分析。将数量庞大的数据资料进行分类整理,并且转化为能被计算机处理使用的数据格式。利用数学处理软件Maple9.0 ,将网格位置和观测点位置以及观测值图形化表示。通过观察这些图形的特点,对题中个观测站点的测量数据及位置分布情况进行更细的划分,并用二次曲面插值及
5、曲面拟合,来获得近似测量数据值。其中,对观察点较独立的雨量值为非零的点,采用抛物面方程为在附近取近似值的方法插值;对两个或两个以上的非零独立点时,利用二次曲面方程做类似的进行插值。利用上述的方法求出各局部雨量观测值的拟合曲面方程,求出个网格处近似观测值的矩阵,然后利用统计的办法求出各预报雨量的网格与模拟观察值之间的方差,并进行概率分析,求出天内所有观测值与预报值间的方差平均差。越小说明预报误差越小,则预报的准确性越高。(2)公众对方法一预报的感受认可程度比方法二的值高:在我们的模型中,刻画公众所感受到的无雨、小雨、中雨、大雨、暴雨、大暴雨及特大暴雨的程度,可以用模糊数学的方法与概率来描述。在对
6、每一种降雨级别进行描述时,可以认为其中部位置应当是公众的感受程度最高之处;而对于所定义的“无雨”级别内,其中感受程度应随雨量的增加而减少,“特大暴雨”级别与“无雨”级别的变动趋势正好相反。建立了每一个级别的隶属度函数模型。 用所建的隶属度函数可以求出任何一个降雨量值,将第一问题中的预报降雨量数值矩阵及观测的降雨量数值矩阵代入隶属度函数,利用各网格点在不同时间区域上的数据方差,求出总体上的均差,做类似的比较。结果为:方法一的均差为0.00016,方法二的均差为0.0011,所以方法一的预报结果公众认可度高。本文随附件带有三百个运行程序文件关键词:降雨量,预报,公众感受,模拟,拟合,模糊数学,方差
7、。雨量预报方法的评价模型一、提出问题:雨量预报对农业生产和城市工作和生活有重要作用,但准确、及时地对雨量作出预报是一个十分困难的问题,广受世界各国关注。我国某地气象台和气象研究所正在研究6小时雨量预报方法,即每天晚上20点预报从21点开始的4个时段(21点至次日3点,次日3点至9点,9点至15点,15点至21点)在某些位置的雨量,这些位置位于东经120度、北纬32度附近的5347的等距网格点上。同时设立91个观测站点实测这些时段的实际雨量,由于各种条件的限制,站点的设置是不均匀的。气象部门希望建立一种科学评价预报方法好坏的数学模型与方法。气象部门提供了41天的用两种不同方法的预报数据和相应的实
8、测数据。其中,雨量用毫米做单位,小于0.1毫米视为无雨。(1) 请建立数学模型来评价两种6小时雨量预报方法的准确性;(2) 气象部门将6小时降雨量分为6等:0.12.5毫米为小雨,2.66毫米为中雨,6.112毫米为大雨,12.125毫米为暴雨,25.160毫米为大暴雨,大于60.1毫米为特大暴雨。若按此分级向公众预报,如何在评价方法中考虑公众的感受?(注:本题数据位于压缩文件C2005Data.rar中, 可从http:/ 表示5347个网格点处的实测降雨量值构成的距阵;:表示5347个网格点的预报降雨量值距阵;:实测值与预报值的差矩阵;: 表示月份(=6,7),表示日期,表示第日21点至此
9、日(及+1日)20点预测量的序号;: 表示均方误差值;:方差均值。三、模型的假设1、将观察区域视为一个平面区域;2、所设的观察点位置足以反映整个预报区域的气候特征;四、模型的建立与分析 第一问:评价两种6小时雨量预报方法的准确性数学模型的建立1 模型分析气象预报是一项统计、计算工作量极大且具有非常明显随机特性的工作,相关部门为了能够较准确的预报天气的变化情况,每天需要处理大量数据,并且利用概率统计的数学方法进行处理与预测,以达到使天气预报的准确性满足人们生活工作的需要。查阅有关气象处理资料,我们了解到天气的预报,是以大量的观察数据为基础进行模拟和分析得到的。这也是本模型的建模思路与方法。首先利
10、用计算机将题中所给的全部数据(共831971条)按时段与数据处理方法的不同进行分类整理。因为观测数据点与预报网格点的个数与位置都不一一对应,所以我们首先要得到网格点上的近似观测值。我们的方法是利用91个观测点的实测值进行数据曲面拟合,通过最佳拟合曲面,得到网格点上的近似实测值。为了能够达到最佳模拟效果,在拟合实测值曲面时,我们根据值的分布特点划分区域,在若干小区域上进行局部曲面拟合,以实现整体最优。每天分四个时段,一共有41天的实测值,这就决定了有164个曲面需要拟合,有164*2491个近似实测值需要计算,计算量极其巨大,必须利用高性能计算机和高效的编程算法完成该工作,以达到客观的分析与评价
11、的效果。其次利用概率统计的均方差、均值等参数对两种预测方法进行评价。利用经、纬度值与高斯坐标的转换公式,将网格点的坐标进行了转换,但是转换后的坐标值十分大,给进一步的数据处理带来不便,所以我们以经、纬度值代替了高斯坐标值,我们进行数据处理的计算机是主频为2.1G、内存为512兆的计算机,使用的编程软件是maple 9.0。若,为网格点的经度及纬度点坐标矩阵,则有: Y=.由拟合曲面方程计算得出,网格点降雨量组成的观测降雨量距阵为Z=. (1)则两种预报方法下,所有网格点的降雨量预测值组成的距阵为 =1,2 -= (2) (3) (4) 其中=1表示预测方法一对应的数值矩阵,=2表示预测方法二对
12、应的数值矩阵;是由5347个网格点上的观测值与预测值之差组成的数值矩阵;是内数值的方差值;是41天所有次数总体上的平均方差.2 计算步骤 根据实测值的分布特点划分拟合区域,在小区域(非规则区域)上用二次曲面插值或高阶曲面拟合,设插值的二次曲面或拟合的曲面方程. 由计算出每一天每一时段下的近似观测降雨量距阵Z. 用=-、及计算出真实值与预报值的差值矩阵及. 用计算出均方差 . 用计算出与,则和的值较小者为最优方案.3 主要程序及结论通过数据处理与分析我们认为预测方法一比预测方法二好。所得计算结果值分别为:(1)不同时段的两种方法的实测与预测值的均方差:=0.9247218269e-1, .165
13、797962696, 0.9247218269e-1, 0.9247218269e-1, .2586806182, .2586806182, .2586806182, 2.791713932, .2474029514, .2539943168, .2715902174, .2715902174182, .2586806182, 2.791713932, .2474029514, .2539943168, .2715902174:= 0.921412432e-1, .1098068392, 0.2234955063e-1, 0.1592933205e-1, .2851304286, .28513
14、04286, .2851304286, 2.792910527, .2612701098, .2381007694, .2613774987, 0.5183032655e-1, .2851304286, 2.792810527, .2612701098, .2381007694, .2613774987(2) 方法一的均方差为: := .8311398371方案二的均方差: = .8417760978得 solve(0.3=0.6-r*(0.0452+0.0422),r); z1:=0.6-79.17656374*(x-120.2500)2+(y-33.7667)2; z2:=0.6-79.1
15、7656374*(x-120.2500)2+(y-33.7667)2; z3:=0.6-79.17656374*(x-120.2500)2+(y-33.7667)2; z4:=0.6-79.17656374*(x-120.2500)2+(y-33.7667)2; solve(0.15=0.3-r*(0.0452+0.0422),r); z4:=0.3-39.58828187*(x-118.1833)2+(y-31.0833)2; solve(5.1=10.2-r*(0.0452+0.0422),r); z1:=10.2-1346.001584*(x-120.3167)2+(y-31.5833)2; z2:=10.2-1346.001584*(x-120.3167)2+(y-31.5833)2;