球面最短距离的求法公式及其推导.DOCVIP

2005高教社杯全国大学生数学建模竞赛 (C题) 雨量预报方法的评价模型 学校名称：江西财经大学 参赛队员：肖潇 杨靖 雅孺 指导老师： 数学建模小组 2005年9月19日  雨量预报方法的评价模型 肖潇 杨靖 雅孺 摘 要： 本文使用方差分析模型、相关性检验模型和创建效用模型分析了气象部门两种雨量预测方法的准确性，利用MATLAB软件强大的计算功能，对问题给出的大量数据进行分析，并得到了一致的结果。 首先，建立了坐标转换公式，将问题中给出的网格点和观测站点的经纬度转化为平面直角坐标，解决了网格点不等距问题。并利用已知的网格点上雨量的预测数据通过二维Lagrange插值多项式，求出了91个观测站点处雨量的预测值； 其次，使用方差分析的方法，，通过比较两种预测方法下所有预测数据与实际数据的总方差、不同观测站点处预测数据与实际数据的方差以及不同时段下预测数据与实际数据的方差，均得到预测方法一较预测方法二的准确性更高的一致结论，并对观测站点的位置进行分析，发现长江三角洲地区和沿海地区的观测站点更适用第二种预测方法；同时发现了五个数据异常的观测站点，并对其成因进行了分析； 第三，使用二分类方法，对两种预测方法的雨量预测数据和实际数据之间分别作了相关系数检验(r—检验)、独立性检验(—检验)和回归系数检验(Q—检验)，均同样得到了与方差分析方法一致的结果：预测方法一较预测方法二的准确性更高； 第四，对于问题2，通过将预测数据转化为分类等级，引入指数效用模型： ，， 分别针对两种预测方法，求出了所有预测等级的总效用、各观测站点上的预测等级的效用以及不同的四个时段上预测等级的效用，通过对两种预测方法效用的比较，依然得出了预测方法一比预测方法二效用更高、准确性更好的结论； 　　最后，作为文章的结束，阐述了本文模型的优点及缺陷，同时也指出了需要改进的地方。 关键字：雨量预报、预测方法、Lagrange插值、误差分析、相关性分析、效用分析 一、问题的重述 雨量预报对农业生产和城市工作和生活有重要作用，但准确、及时地对雨量作出预报是一个十分困难的问题，广受世界各国关注。我国某地气象台和气象研究所正在研究6小时雨量预报方法，即每天晚上20点预报从21点开始的4个时段（21点至次日3点，次日3点至9点，9点至15点，15点至21点）在某些位置的雨量，这些位置位于东经120度、北纬32度附近的53×47的等距网格点上。同时设立91个观测站点实测这些时段的实际雨量，由于各种条件的限制，站点的设置是不均匀的。 气象部门希望建立一种科学评价预报方法好坏的数学模型与方法。气象部门提供了41天的用两种不同方法的预报数据和相应的实测数据。小于0.1毫米视为无雨。 请建立数学模型来评价两种6小时雨量预报方法的准确性； 气象部门将6小时降雨量分为6等：0.1—2.5毫米为小雨，2.6—6毫米为中雨，6.1—12毫米为大雨，12.1—25毫米为暴雨，25.1—60毫米为大暴雨，大于60.1毫米为特大暴雨。若按此分级向公众预报，如何在评价方法中考虑公众的感受？ 二、模型的分析与假设 根据预测点（网格点）和实测点的经纬度，画出预测点和实测点的位置图（图1），但由于经纬度坐标并不对应于直角坐标系，得到的并不是问题中给出的规则网格图。同时发 图（一） 现91个观测点不但分布不均匀,而且基本上没有落在预测点上，因此，我们考虑将预测点和实测点的经纬度坐标转化为平面直角坐标。 坐标转换后91个实测点并没有在网格点上，因此我们将使用二维插值多项式计算出实测点处两种预测方法预报的降雨量的近似值，然后将其与实测点处降雨量的实际值进行比较，给出评价模型和问题的解答。 鉴于本问题中给出的数据十分庞大，（不完全统计，数据量达80万个），所以在分析过程中，本文大量使用了MATLAB软件中提取文件中数据的dlmread函数。 根据本文的需要，我们提出以下合理的假设： 1．雨量用毫米做单位，小于0.1毫米视为无雨； 2．地球为规则球体，地球半径取WGS84标准中的长半轴6378.137公里，短半轴6356.752公里的平均值6371公里； 3．预测点位置位于东经120度、北纬32度附近的53×47的等距网格点上； 4．不考虑人为因素所造成的测量值的误差； 5．题目中所给经纬度坐标值和雨量数据为真实值，或者与真实情况相符合。 三、符号说明 如无特别说明，本文中的符号特指以下意义： ：经纬度坐标； ：修正后欧氏坐标； ：球心角； ：经度差； ：两点的球面距离； ：降雨量； ：表示第i个观测点第j天第k时段雨量的实际值； ：表示第i个观测点第j天第k时段雨量的预测值； ：表示用雨量的预测值与实际值的方差； ：