雨量預(yù)報(bào)方法的評(píng)價(jià)論-預(yù)覽頁(yè)
【正文】 報(bào)方法，即 每天晚上 20 點(diǎn)預(yù)報(bào) 從 21 點(diǎn)開始的 4 個(gè)時(shí)段（ 21 點(diǎn)至次日 3 點(diǎn)，次日 3 點(diǎn)至 9 點(diǎn)， 9 點(diǎn)至 15 點(diǎn)， 15 點(diǎn)至 21 點(diǎn)）在某些位置的雨量，這些位置位于東經(jīng) 120 度、北緯 32 度附近的 53 47（緯度） 的等距 網(wǎng)格點(diǎn)上。預(yù)報(bào)數(shù)據(jù)在文件夾 FORECAST 中， 實(shí) 測(cè)數(shù)據(jù) 在文件夾 MEASURING 中，其中的文件都可以用Windows 系統(tǒng)的“寫字板”程序打開閱讀。若按此分級(jí)向公眾預(yù)報(bào)，如何在評(píng)價(jià)方法中考慮公眾感受 . 問(wèn)題分析 針對(duì)問(wèn)題一，在 53 47 的 經(jīng)緯度 等距 網(wǎng)格點(diǎn)上 有方法一、二的預(yù)測(cè)雨量 ，在 91 個(gè)觀測(cè)站點(diǎn) 上有 41 天每天 4 個(gè)時(shí)段的 實(shí)測(cè) 雨量 ，而且站點(diǎn)的設(shè)置 不與網(wǎng)格格點(diǎn)重合。 模型求解 （一） Matlab 求解 1Step :利用 dlmread 函數(shù) 從 6 月 18 日到 7 月 30 日循環(huán)讀取實(shí)測(cè)站點(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ù)測(cè)的網(wǎng)格點(diǎn)的雨量賦給矩 陣 1rqf 、 2rqf ，并通過(guò) griddata 函數(shù)插值得到方法一、二在 91 個(gè)站點(diǎn)上預(yù)測(cè)雨 量 _1rqm f 、 _2rqm f ； 4Step ： 站點(diǎn)從 1 到 91 循環(huán)，對(duì)應(yīng)站點(diǎn)的預(yù)測(cè)雨量和實(shí)測(cè)雨量相減取絕對(duì)值，分別得到方法一、二的 誤差 a 、 b ； 5Step ： a 與 b 中小者則對(duì)應(yīng)讓 1n 或 2n 自增 1，知直到 91 個(gè)站點(diǎn)循環(huán)完；然后再讓 1n 與 2n比較，較大者則對(duì)應(yīng)地讓1N或2N自增 1，直到 4 個(gè)時(shí)段循環(huán)完；接著讓1N與2N比較，兩者中較大者對(duì)應(yīng)地讓1M或2M自增 1，直到 41 天循環(huán)完 （二）求解結(jié)果與分析 最后求得 1M =17， 2M =10,；從結(jié)果可知，在 41 天中對(duì)兩種方法進(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)確。M2=0。02039。%實(shí)測(cè)字符串 41 tmp=dlmread(fm)。%4 個(gè)時(shí)段 lonf=dlmread(39。)。39。 for k=1:4 ff1=[39。,num2str(i),num2str(k),39。data\FORECAST\39。_dis239。 rqm_f1=griddata(lonf,latf,rqf1,lonm,latm,39。V439。 for j=1:91 a=abs(rqm(j,k)rqm_f1(j))。 end end end if n1n2 N1=N1+1。 end end end M1 M2 問(wèn)題二程序 ： clear all clc for i=[618:628,701:730] fm=[39。,num2str(i),39。 latm=tmp(:,2)。 r(r)=0。r6)=2。r25)=4。 lonf=dlmread(39。)。39。,39。]。f39。 rqf1=dlmread(ff1)。)。 r1=rqm_f1。 r1)=1。 r112)=3。 r160)=5。 r2)=1。 r212)=3。 r260)=5。 c1=abs(rqm_f1(j)rqm(j,k))。 else if rqm_f1(j)=0 amp。 else if c1= amp。 rqm_f1(j)25 if c1= amp。 c112 x1(j)=2。 s1(j)=1./(1+*x1(j).^)。 else if c1=10 amp。c130 x1(j)=3。 s1(j)=1./(1+*x1(j).^)。 else if c1=60 x1(j)=6。 else if rqm_f2(j)=0 amp。 else if c2= amp。 rqm_f2(j)25 if c2= amp。 c212 x2(j)=2。 s2(j)=1./(1+*x2(j).^)。 else if c2=10 amp。c230 x2(j)=3。 s2(j)=1./(1+*x2(j).^)。 else if c2=60 x2(j)=6。 end sss1(i)=sum(ss1)。 end S1 S2
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1