【正文】 r2(r2= amp。 r1(r1=12 amp。 r2=rqm_f2。 rqf2=dlmread(ff2)。 ff2=[39。)。F:\5\data\FORECAST\39。 r(r=6 amp。 lonm=tmp(:,3)。F:\5\data\MEASURING\39。 b=abs(rqm(j,k)rqm_f2(j))。V439。,39。data\FORECAST\39。%預(yù)測網(wǎng)格點(diǎn)經(jīng)度 latf=dlmread(39。%實(shí)測數(shù)據(jù) 41 latm=tmp(:,2)。 for i=[618:628,701:730] fm=[39。 把經(jīng)緯度分別設(shè) 為 x， y 軸，雨量設(shè)為 z 軸， 對方法一預(yù)測的數(shù)據(jù)進(jìn)行插值， 求出 由 方法一得到 的 91 個(gè) 站點(diǎn) 的 預(yù)測 雨量 _1rqm f ，同理求出 91 個(gè) 站點(diǎn) 由 方法二得到的 預(yù)測 雨量 _2rqm f ，再與 91 個(gè)站點(diǎn) 的 實(shí)測 雨量 rqm 進(jìn)行 作差得到某站點(diǎn)的誤差 ， 誤差較小者則相應(yīng)地讓 12nn或 自增 1,91 個(gè)站點(diǎn)比完之后再對 4個(gè)時(shí)段進(jìn)行比較， 12nn或 中較大者相應(yīng)地讓 12N或 N 自增 1，然后再對 41 天進(jìn)行比較，12N或 N 中較大者相應(yīng)地讓 12MM或 自增 1，最后可得到兩種方法在 41 天內(nèi)更準(zhǔn)確的天數(shù)； 針對問題二，循環(huán)分別讀取出某天某時(shí)段的 91 個(gè)站點(diǎn)的實(shí)測雨量，和網(wǎng)格點(diǎn)預(yù)測雨量，用散亂節(jié)點(diǎn)插值函數(shù) griddata 插值求得某天某時(shí)段 91 個(gè)站點(diǎn)的預(yù)測雨量，由此便得到某天某時(shí)段 91 個(gè)站點(diǎn)的實(shí)測雨量 rqm 、方法一預(yù)測雨量 _1rqm f 、方法二預(yù)測雨量 _2rqm f ，接著讓 _1rqm f 、 _2rqm f 和 rqm 按照題目給出的范圍進(jìn)行雨量的等級劃分，然后對方法一、二分別求出等級誤差 1a 、 2a 和雨量誤差 1c 、 2c ； 若等級誤差為 0，則滿意度為 1，否則，按照分類討論中的表格進(jìn)行等級誤差的確定，并代入函數(shù) x? ?求出某站點(diǎn)的滿意度 ；最后把 41 天的 4 個(gè)時(shí)段的 91 個(gè)站點(diǎn)的滿意度分別加起來得到兩種方法的總滿意度 1S 和 2S . 模型假設(shè) （ 1） 實(shí)際測量的雨量沒有誤差，不受風(fēng)速等影響； （ 2） 方 法一、二預(yù)測的雨量不會(huì)因?yàn)槿藶橐蛩囟e(cuò)誤。同時(shí)設(shè)立 91 個(gè)觀測站點(diǎn)實(shí)測這些時(shí)段的實(shí)際雨量，由于各種條件的限制，站點(diǎn)的設(shè)置是不均勻的。 問題一 分別對 每一天每一時(shí)段兩種方法預(yù)測的雨量 進(jìn)行 插值 ，得到兩種方法在 91 個(gè)站點(diǎn)上的預(yù)測雨量 ，并分別與實(shí)際雨量 對應(yīng) 進(jìn)行比較，哪種方法預(yù)測的誤差 較小的總的次數(shù)最多 ，即為預(yù)測的比較準(zhǔn)確 的 方法，由 Matlab 編程求解得知方法一預(yù)報(bào)的比方法二更準(zhǔn)確。預(yù)報(bào)數(shù)據(jù)在文件夾 FORECAST 中， 實(shí) 測數(shù)據(jù) 在文件夾 MEASURING 中，其中的文件都可以用Windows 系統(tǒng)的“寫字板”程序打開閱讀。 模型求解 （一） Matlab 求解 1Step :利用 dlmread 函數(shù) 從 6 月 18 日到 7 月 30 日循環(huán)讀取實(shí)測站點(diǎn)數(shù)據(jù)，并賦給矩陣 tmp ， 從矩陣 tmp 提取出第 4 到 7 列分別賦給矩陣 latm 、 lonm 、rqm ； 2Step ： 讀取網(wǎng)格點(diǎn)的經(jīng)、緯度賦給矩陣 lonf 、 latf ； 3Step ： 讓時(shí)段從 1 到 4 循環(huán)，分別讀取方法一、二預(yù)測的網(wǎng)格點(diǎn)的雨量賦給矩 陣 1rqf 、 2rqf ，并通過 griddata 函數(shù)插值得到方法一、二在 91 個(gè)站點(diǎn)上預(yù)測雨 量 _1rqm f 、 _2rqm f ； 4Step ： 站點(diǎn)從 1 到 91 循環(huán)，對應(yīng)站點(diǎn)的預(yù)測雨量和實(shí)測雨量相減取絕對值，分別得到方法一、二的 誤差 a 、 b ； 5Step ： a 與 b 中小者則對應(yīng)讓 1n 或 2n 自增 1，知直到 91 個(gè)站點(diǎn)循環(huán)完；然后再讓 1n 與 2n比較，較大者則對應(yīng)地讓1N或2N自增 1，直到 4 個(gè)時(shí)段循環(huán)完；接著讓1N與2N比較，兩者中較大者對應(yīng)地讓1M或2M自增 1，直到 41 天循環(huán)完 （二）求解結(jié)果與分析 最后求得 1M =17， 2M =10,；從結(jié)果可知，在 41 天中對兩種方法進(jìn)行 4 個(gè)時(shí)段的比較 ，其中有 17 天方法一更準(zhǔn)確，有 10 天方法二更準(zhǔn)確，還有 14 天，方法一、二在 4 個(gè)時(shí)段中分別都有兩個(gè)時(shí)段更準(zhǔn)確，即兩種方法在這 14 天的 4 個(gè)時(shí)段均一樣準(zhǔn)確。02039。%4 個(gè)時(shí)段 lonf=dlmread(39。39。,num2str(i),num2str(k),39。_dis239。V439。 end end end if n1n2 N1=N1+1。,num2str(i),39。 r(r)=0。r25)=4。)。,39。f39。)。 r1)=1。 r160)=5。 r212)=3。 c1=abs(rqm_f1(j)rqm(j,k))。 else if c1= amp。 c112 x1(j)=2。 else if c1=10 amp。 s1(j)=1./(1+*x1(j).^)。 else if rqm_f2(j)=0 amp。 rqm_f2(j)25 if c2= amp。 s2(j)=1./(1+*x2(j).^)。c230 x2(j)=3。 else if c2=60 x2(j)=6。 end S1 S2