江西师范大学:熊军军、许谞、许盛敏指导教师:温利民 1、附录一、等距网格点和观测站点的散点分布图程序 X=[纬度] x是由文件夹 FORECAST中lat文件的数据生成的一个纬度值矩阵,x 大小为53*4 y=经度]:%y是由文件夹 FORECAST中lon文件的数据生成的一个经度值矩阵,y 大小为53*47 lat=32983333300033666733.800033.483333.033333.2333 33.766733.383333.200032.100032.300032.000032.683332.8000 32933332416732.333332.200032.866732.183332533332.3833 32.333332.066731.800031950031.333331.566731.700031.0833 31983331.750031766731950031.433331.366731266731.8833 31650031.583331416731.066731.150031.900031.100031.4000 31.366731.616731.200031.050031.233331466730.300030.8500 30.683330.933330.300030.616730.083330.883331.1333309833 0.966730.233330.050030.850030.850030.783330.000030.5333 30.516730.633330.233330.200030.883331.133331.000030.9333 30.616730.266730.066730.733330.033330.250029.866729.7167 29.783329816729.7000299667]; lon=18.51671188500119266711980001198167119.033 19.3000 120.2500120.150012048331182667118.3000118.80001190167119.4500 119833311941671199333120.0000120.316711946671204500120.5667 121.1833121.6000121.6667118.8500118.3833118.5000118.5167 1195833119.55001199333119166711948331198167120.6333120.2667 120.7333120.3167120.95001204333120.63331212000121.36671214833 121.2500121.4500121.4333121.7833121.5333121.1000118.1333118.3167 1184000118.7500118.5333118.9833118.5833119416711918331198833 19.6833119700011995001200833120.9000120.7333120.6333120.0667 120.6833120.5333120.1667120.3167121.1667121.1167121.2500121.4833 121.0833121.2167121.15001224500122.1000122.18331184333118.2833 l18.18331196833120.2500121.7500]; y,g); xlabel(纬度", Font Weight,bold); ylabel(经度, Font Weight,bold); legend(实测站点,网格结点) xis([2736117125]) 3、附录二搜索观察站点周围最近的结点以及根据距离长短排序 a=[经度] %a是由文件夹 FORECAST中lon文件的数据生成的一个经度值矩 阵a大小为53*47 b=[纬度] %b是由文件夹 FORECAST中lat文件的数据生成的一个纬度值 矩阵,b大小为53*47 c=[经度,纬度] %c为9个实测站点对应的经度和纬度,c大小为91*2 =0.25; %e为可接受度数差 for k=1: 91 m=1 for i=1: 53 for j=1: 47 if abs((b(ij)c(k,1)<=e&abs(a(ij)-c(k,2)<=e%在可接受度数差范 围内,搜索几个网格点 2005年全国大学生数学建模竞赛全国一等奖江西师范大学：熊军军、许谞、许盛敏 指导教师：温利民 2005 年全国大学生数学建模竞赛全国一等奖 1、附录一、等距网格点和观测站点的散点分布图程序 x=[纬度]; % x 是由文件夹 FORECAST 中 lat 文件的数据生成的一个纬度值矩阵,x 大小为 53*47 y=[经度]； % y 是由文件夹 FORECAST 中 lon 文件的数据生成的一个经度值矩阵,y 大小为 53*47 lat=[32.9833 33.3000 33.6667 33.8000 33.4833 33.0333 33.2333 33.7667 33.3833 33.2000 32.1000 32.3000 32.0000 32.6833 32.8000 32.9333 32.4167 32.3333 32.2000 32.8667 32.1833 32.5333 32.3833 32.3333 32.0667 31.8000 31.9500 31.3333 31.5667 31.7000 31.0833 31.9833 31.7500 31.7667 31.9500 31.4333 31.3667 31.2667 31.8833 31.6500 31.5833 31.4167 31.0667 31.1500 31.9000 31.1000 31.4000 31.3667 31.6167 31.2000 31.0500 31.2333 31.4667 30.3000 30.8500 30.6833 30.9333 30.3000 30.6167 30.0833 30.8833 31.1333 30.9833 30.9667 30.2333 30.0500 30.8500 30.8500 30.7833 30.0000 30.5333 30.5167 30.6333 30.2333 30.2000 30.8833 31.1333 31.0000 30.9333 30.6167 30.2667 30.0667 30.7333 30.0333 30.2500 29.8667 29.7167 29.7833 29.8167 29.7000 29.9667]; lon=[118.5167 118.8500 119.2667 119.8000 119.8167 119.0333 119.3000 120.2500 120.1500 120.4833 118.2667 118.3000 118.8000 119.0167 119.4500 119.8333 119.4167 119.9333 120.0000 120.3167 119.4667 120.4500 120.5667 121.1833 121.6000 121.6667 118.8500 118.3833 118.5000 118.5167 118.1833 119.5833 119.5500 119.9333 119.1667 119.4833 119.8167 120.6333 120.2667 120.7333 120.3167 120.9500 120.4333 120.6333 121.2000 121.3667 121.4833 121.2500 121.4500 121.4333 121.7833 121.5333 121.1000 118.1333 118.3167 118.4000 118.7500 118.5333 118.9833 118.5833 119.4167 119.1833 119.8833 119.6833 119.7000 119.9500 120.0833 120.9000 120.7333 120.6333 120.0667 120.6833 120.5333 120.1667 120.3167 121.1667 121.1167 121.2500 121.4833 121.0833 121.2167 121.1500 122.4500 122.1000 122.1833 118.4333 118.2833 118.1833 119.6833 120.2500 121.7500]; plot(lat,lon,'dk',x,y,'.g'); xlabel('纬度','FontWeight','bold'); ylabel('经度','FontWeight','bold'); legend('实测站点','网格结点') axis([27 36 117 125 ]) 3、附录二 搜索观察站点周围最近的结点以及根据距离长短排序 a=[经度]； % a 是由文件夹 FORECAST 中 lon 文件的数据生成的一个经度值矩 阵,a 大小为 53*47 b=[纬度]； % b 是由文件夹 FORECAST 中 lat 文件的数据生成的一个纬度值 矩阵,b 大小为 53*47 c=[经度，纬度]; % c 为 91 个实测站点对应的经度和纬度,c 大小为 91*2 g=zeros(91,20); e=0.25; % e 为可接受度数差 for k=1:91 m=1; for i=1:53 for j=1:47 if abs( b(i,j)-c(k,1))<=e&abs(a(i,j)-c(k,2))<=e % 在可接受度数差范 围内，搜索几个网格点