【正文】 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。%預測網(wǎng)格點經(jīng)度 latf=dlmread(39。%實測數(shù)據(jù) 41 latm=tmp(:,2)。 for i=[618:628,701:730] fm=[39。 把經(jīng)緯度分別設(shè) 為 x， y 軸，雨量設(shè)為 z 軸， 對方法一預測的數(shù)據(jù)進行插值， 求出 由 方法一得到 的 91 個 站點 的 預測 雨量 _1rqm f ，同理求出 91 個 站點 由 方法二得到的 預測 雨量 _2rqm f ，再與 91 個站點 的 實測 雨量 rqm 進行 作差得到某站點的誤差 ， 誤差較小者則相應地讓 12nn或 自增 1,91 個站點比完之后再對 4個時段進行比較， 12nn或 中較大者相應地讓 12N或 N 自增 1，然后再對 41 天進行比較，12N或 N 中較大者相應地讓 12MM或 自增 1，最后可得到兩種方法在 41 天內(nèi)更準確的天數(shù)； 針對問題二，循環(huán)分別讀取出某天某時段的 91 個站點的實測雨量，和網(wǎng)格點預測雨量，用散亂節(jié)點插值函數(shù) griddata 插值求得某天某時段 91 個站點的預測雨量，由此便得到某天某時段 91 個站點的實測雨量 rqm 、方法一預測雨量 _1rqm f 、方法二預測雨量 _2rqm f ，接著讓 _1rqm f 、 _2rqm f 和 rqm 按照題目給出的范圍進行雨量的等級劃分，然后對方法一、二分別求出等級誤差 1a 、 2a 和雨量誤差 1c 、 2c ； 若等級誤差為 0，則滿意度為 1，否則，按照分類討論中的表格進行等級誤差的確定，并代入函數(shù) x? ?求出某站點的滿意度 ；最后把 41 天的 4 個時段的 91 個站點的滿意度分別加起來得到兩種方法的總滿意度 1S 和 2S . 模型假設(shè) （ 1） 實際測量的雨量沒有誤差，不受風速等影響； （ 2） 方 法一、二預測的雨量不會因為人為因素而錯誤。同時設(shè)立 91 個觀測站點實測這些時段的實際雨量，由于各種條件的限制，站點的設(shè)置是不均勻的。 問題一 分別對 每一天每一時段兩種方法預測的雨量 進行 插值 ，得到兩種方法在 91 個站點上的預測雨量 ，并分別與實際雨量 對應 進行比較，哪種方法預測的誤差 較小的總的次數(shù)最多 ，即為預測的比較準確 的 方法，由 Matlab 編程求解得知方法一預報的比方法二更準確。預報數(shù)據(jù)在文件夾 FORECAST 中， 實 測數(shù)據(jù) 在文件夾 MEASURING 中，其中的文件都可以用Windows 系統(tǒng)的“寫字板”程序打開閱讀。 模型求解 （一） Matlab 求解 1Step :利用 dlmread 函數(shù) 從 6 月 18 日到 7 月 30 日循環(huán)讀取實測站點數(shù)據(jù)，并賦給矩陣 tmp ， 從矩陣 tmp 提取出第 4 到 7 列分別賦給矩陣 latm 、 lonm 、rqm ； 2Step ： 讀取網(wǎng)格點的經(jīng)、緯度賦給矩陣 lonf 、 latf ； 3Step ： 讓時段從 1 到 4 循環(huán)，分別讀取方法一、二預測的網(wǎng)格點的雨量賦給矩 陣 1rqf 、 2rqf ，并通過 griddata 函數(shù)插值得到方法一、二在 91 個站點上預測雨 量 _1rqm f 、 _2rqm f ； 4Step ： 站點從 1 到 91 循環(huán)，對應站點的預測雨量和實測雨量相減取絕對值，分別得到方法一、二的 誤差 a 、 b ； 5Step ： a 與 b 中小者則對應讓 1n 或 2n 自增 1，知直到 91 個站點循環(huán)完；然后再讓 1n 與 2n比較，較大者則對應地讓1N或2N自增 1，直到 4 個時段循環(huán)完；接著讓1N與2N比較，兩者中較大者對應地讓1M或2M自增 1，直到 41 天循環(huán)完 （二）求解結(jié)果與分析 最后求得 1M =17， 2M =10,；從結(jié)果可知，在 41 天中對兩種方法進行 4 個時段的比較 ，其中有 17 天方法一更準確，有 10 天方法二更準確，還有 14 天，方法一、二在 4 個時段中分別都有兩個時段更準確，即兩種方法在這 14 天的 4 個時段均一樣準確。02039。%4 個時段 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