数学建模竞赛2005年C题论文1[1].doc

雨量预报方法的评价 摘要 雨量预报对农业生产和城市工作和生活有重要作用，但准确、及时地对雨量做出预报是一个很困难的问题，广受世界各国关注。我国某地气象台和气象研究所正在研究6小时雨量预报方法，由于受到学科发展水平的限制，目前国内外降雨数值预报水平还不高.为了使预报方法更为准确，使天气预报更好的服务于公众生活，我们用数学模型来分析研究这一问题. 文中我们建立了比较两种预测降雨量方法优劣的数学模型.即根据2491个网格点的纬度、经度和降水量的预测值，采用二维插值的方式，分别对91个观测站点的降雨量进行预测，利用Matlab软件中的griddata函数： =griddata =griddata 然后将其与实测值对比，求出预测值与实际值之间的误差，利用Matlab软件中的矩阵范数函数norm N1==norm()=405.3782，N2==norm()=416.1976 根据范数的含义,所得范数越小，即误差越小.因为有N1N2，故可得出结论： 第一种方法比第二种方法预测雨量的准确性更高. 为了解决如何在评价方法中考虑公众的感受的问题，我们将第一题中通过二维插值得到的91个气象站41天的预测值用分级形式输出，即无雨、小雨、中雨、大雨、暴雨、大暴雨、特大暴雨.将两种方法输出的雨量预报情况与实际降雨量情况进行比较， 统计出每种方法准确预报、空报、漏报的次数，误报次数越少的，对应的方法准确性应越高，公众对其可信度越高.程序运算结果得到： 降雨预测结果 两种预测方法的预测结果 第一种预测方法 第二种预测方法 预测值等于实测值的元素个数 12053 11921 预测值大于实测值的元素个数 2476 2552 预测值小于实测值的元素个数 395 451 预测值等于实测值代表观测站点预报准确，预测值大于实测值代表观测站点空报的次数或对天气状况预测过于恶劣，预测值小于实测值代表观测站点漏报的次数或对天气恶劣状况估计不足. 得到两种预测方法的准确率分别为80.7625%，79.8780%. 可见运用第一种方法时，误报的次数较少，准确率较高，故第一种方法较好. 为了使雨量预报方法准确性更高，适用范围更广，我们给出了改进建议.一、问题重述 雨量预报对农业生产和城市工作和生活有重要作用，但准确、及时地对雨量做出预报是一个很困难的问题，广受世界各国关注.我国某地气象台和气象研究所正在研究6小时雨量预报方法，即每天晚上20点预报从21点开始的4个时段（21点至次日3点，次日3点至9点，9点至15点，15点至21点）在某些位置的雨量，这些位置位于东经120度、北纬32度附近的53×47的等距网格点上.同时设立91个观测站点实测这些时段的实际雨量，由于各种条件的限制，站点的设置是不均匀的. 气象部门希望建立一种科学评价预报方法好坏的数学模型与方法.气象部门提供了41天的用两种不同方法的预报数据和相应的实测数据. 雨量用毫米做单位，小于0.1毫米视为无雨. (1)请建立数学模型来评价两种6小时雨量预报方法的准确性； (2)气象部门将6小时降雨量分为6等：0.1-2.5毫米为小雨，2.6-6毫米为中雨，6.1-12毫米为大雨，12.1-25毫米为暴雨，25.1-60毫米为大暴雨，大于60.1毫米为特大暴雨.若按此分级向公众预报，如何在评价方法中考虑公众的感受？ 二、模型假设 1.天气变化状况是局部连续的. 2.各个观测站点设备及测量水平相同，不存在技术上的误差. 三、符号约定 ——网格点的纬度构成的矩阵 ——网格点的经度构成的矩阵 ——采用第一种预测方法时，网格点处的降雨量预测值 ——采用第二种预测方法时，网格点处的降雨量预测值 ——按照第一种预测方法，第天，第个时段的预测结果,是一个91维的列向量. (=1,2……,41;=1,2,3,4) ——按照第二种预测方法，第天，第个时段的预测结果，是一个91维的列向量. (=1,2……,41;=1,2,3,4) ——第个气象站点在第个时段降雨量的实测值 (=1,2……,91;=1,2,3,4) 四、模型的建立与求解 两种预测方法的优劣比较 衡量一种降水量预测方法的优劣，依据就是由这种方法预报的天气状况能够准确的反映实际的天气变化.因此我们可以这样建立模型：将题目中给出的预测和实测两种数据导入Matlab软件.lat，lon数据导入后作为两个矩阵的形式，代表网格点的相应位置；其余数据为相应网格点处降雨量的预测值. 根据上述对应关系，我们可以对已经给出的预测值采用二维插值的方式，找出它们之间的关系:,分别对91个观测站点的降雨量进行预测，然后将预测值与实测值对比；利用矩阵范数，得到预测值与实际值之间的误差，将这两个误差相比，误差小的，相应的预测方法就比较准确. 算法步骤： 以20