【正文】 4 到 7 列分別賦給矩陣 latm 、 lonm 、rqm ,讓矩陣 r =rqm ，并讓 r 按照題目給出的雨量范圍進行降雨量等級劃分 ； 2Step ：讀取網(wǎng)格點的經(jīng)、緯度賦給矩陣 lonf 、 latf ； 3Step ：讓時段從 1 到 4 循環(huán)，分別讀取方法一、二預測的網(wǎng)格點的雨量賦給矩 陣 1rqf 、 2rqf ，并通過 griddata 函數(shù)插值得到方法一、二在 91 個站點上預測雨 量 _1rqm f 、 _2rqm f ，分別讓 1r = _1rqm f 、 2r = _2rqm f ，并讓 1r 、 2r 按照 題目給出的雨量范圍進行 降雨量等級劃分 ； 4Step ：站點從 1 到 91 循環(huán)， 對應站點的余量等級對應相減取絕對值 ，分別得到方法一、二的 等級 誤差1a、2a， 對應站點的預測雨量和實測雨量相減取絕對值，分別得到方法一、二的雨量誤差 1c 、 2c ； 5Step ： 對兩種方法，若等級差為 0，則滿 意度 s=0； 否則 ， 則按分類討論中的情況編程即可 ， 對應的等級差代入隸屬函數(shù) x? ?；最終把 41 天的 4 個時段的 91個站點的滿意度加起來得到方法一、二的總滿意度 1S 和 2S （二）求解結果與分析 最終 求得 ， 方法一的總滿意度 1S =14394， 方法二的總滿意度 2S =14923。,num2str(i),39。data\FORECAST\39。)。_dis139。]。)。 else if n1n2 N2=N2+1。.SIX39。 r(r= amp。 r(r=25 amp。 latf=dlmread(39。f39。,num2str(i),num2str(k),39。 rqm_f2=griddata(lonf,latf,rqf2,lonm,latm,39。 r1(r1= amp。 r1(r1=60)=6。 r2(r2=12 amp。 c2=abs( rqm_f2rqm(j,k) )。 c1= x1(j)=2。 s1(j)=1./(1+*x1(j).^)。 c120 x1(j)=2。 else if c1=40 amp。 rqm_f2(j)6 if c2= amp。c26 x2(j)=1。 end end end else if rqm_f2(j)=25 if c2= amp。 s2(j)=1./(1+*x2(j).^)。 s2(j)=1./(1+*x2(j).^)。 S2=sum(sss2)。 s2(j)=1./(1+*x2(j).^)。 else if c2=20 amp。 c2= x2(j)=3。 end end else if rqm_f2(j)=6 amp。 s2(j)=1./(1+*x2(j).^)。c140 x1(j)=4。 s1(j)=1./(1+*x1(j).^)。 else if c1=6 amp。 s1(j)=1./(1+*x1(j).^)。 a2=abs( r2(j)r(j,k) )。 r2(r2=6 amp。 r1(r1=25 amp。 r1(r1= amp。V439。,39。F:\5\data\FORECAST\39。39。 r(r=12 amp。 r=rqm。02039。 else if ab n2=n2+1。%插值得 91 個站點的預測值 rqm_f2=griddata(lonf,latf,rqf2,lonm,latm,39。,num2str(i),num2str(k),39。f39。,39。%實測經(jīng)度 rqm=tmp(:,4:7)。,39。 （二 ）模型的建立 分別以經(jīng)、緯度為 x、 y 坐標建立平面直角坐標系，得到網(wǎng)格點與實測點的相對位置，如下圖： 通過循環(huán)分別讀取出某天某時段的 91 個站點的實測雨量，和網(wǎng)格點預測雨量，用散亂節(jié)點插值函數(shù) griddata 插值求得某天某時段 91 個站點的預測雨量，由此便得到某天某時段 91 個站點的實測雨量 rqm 、方法一預測雨量 _1rqm f 、方法二預測雨量_2rqm f ，接著讓 _1rqm f 、 _2rqm f 分別和 rqm 進行比較， 91 個站點循環(huán)完后，比較1n 與 2n ，即可得到某天某時段中更準確的方法，再進行 4 個時段的循環(huán)，比較 1N 和 2N 即可得到某天內(nèi)更準確的方法，最后再對 41 天循環(huán)，比較 1M 和 2M ， 即可得出 41 天中更為準確的方法；具體實施 步驟參考模型求解。氣象部門提供了 41 天的用