雨量預(yù)報方法的評價論-預(yù)覽頁
【正文】 報方法，即 每天晚上 20 點預(yù)報 從 21 點開始的 4 個時段（ 21 點至次日 3 點，次日 3 點至 9 點， 9 點至 15 點， 15 點至 21 點）在某些位置的雨量，這些位置位于東經(jīng) 120 度、北緯 32 度附近的 53 47（緯度） 的等距 網(wǎng)格點上。預(yù)報數(shù)據(jù)在文件夾 FORECAST 中， 實 測數(shù)據(jù) 在文件夾 MEASURING 中，其中的文件都可以用Windows 系統(tǒng)的“寫字板”程序打開閱讀。若按此分級向公眾預(yù)報，如何在評價方法中考慮公眾感受 . 問題分析 針對問題一，在 53 47 的 經(jīng)緯度 等距 網(wǎng)格點上 有方法一、二的預(yù)測雨量 ，在 91 個觀測站點 上有 41 天每天 4 個時段的 實測 雨量 ，而且站點的設(shè)置 不與網(wǎ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)，分別讀取方法一、二預(yù)測的網(wǎng)格點的雨量賦給矩 陣 1rqf 、 2rqf ，并通過 griddata 函數(shù)插值得到方法一、二在 91 個站點上預(yù)測雨 量 _1rqm f 、 _2rqm f ； 4Step ： 站點從 1 到 91 循環(huán)，對應(yīng)站點的預(yù)測雨量和實測雨量相減取絕對值，分別得到方法一、二的 誤差 a 、 b ； 5Step ： a 與 b 中小者則對應(yīng)讓 1n 或 2n 自增 1，知直到 91 個站點循環(huán)完；然后再讓 1n 與 2n比較，較大者則對應(yīng)地讓1N或2N自增 1，直到 4 個時段循環(huán)完；接著讓1N與2N比較，兩者中較大者對應(yīng)地讓1M或2M自增 1，直到 41 天循環(huán)完 （二）求解結(jié)果與分析 最后求得 1M =17， 2M =10,；從結(jié)果可知，在 41 天中對兩種方法進行 4 個時段的比較 ，其中有 17 天方法一更準確，有 10 天方法二更準確，還有 14 天，方法一、二在 4 個時段中分別都有兩個時段更準確，即兩種方法在這 14 天的 4 個時段均一樣準確。M2=0。02039。%實測字符串 41 tmp=dlmread(fm)。%4 個時段 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 問題二程序 ： 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號-1