雨量預(yù)報(bào)方法的評(píng)價(jià)論-免費(fèi)閱讀
【正文】 end sss1(i)=sum(ss1)。 s2(j)=1./(1+*x2(j).^)。 else if c2=10 amp。 c212 x2(j)=2。 else if c2= amp。 else if c1=60 x1(j)=6。c130 x1(j)=3。 s1(j)=1./(1+*x1(j).^)。 rqm_f1(j)25 if c1= amp。 else if rqm_f1(j)=0 amp。 r260)=5。 r2)=1。 r112)=3。 r1=rqm_f1。 rqf1=dlmread(ff1)。]。39。 lonf=dlmread(39。r6)=2。 latm=tmp(:,2)。 end end end M1 M2 問(wèn)題二程序 ： clear all clc for i=[618:628,701:730] fm=[39。 for j=1:91 a=abs(rqm(j,k)rqm_f1(j))。 rqm_f1=griddata(lonf,latf,rqf1,lonm,latm,39。data\FORECAST\39。 for k=1:4 ff1=[39。)。%實(shí)測(cè)字符串 41 tmp=dlmread(fm)。M2=0。若按此分級(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)重合。我國(guó)某地氣象臺(tái)和氣象研究所正在研究 6 小時(shí) 雨量預(yù)報(bào)方法，即 每天晚上 20 點(diǎn)預(yù)報(bào) 從 21 點(diǎn)開(kāi)始的 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)上。 問(wèn)題二 為了刻畫(huà)公眾對(duì)兩種預(yù)報(bào)方法的滿意度，我們引進(jìn)了偏小型柯西隸屬函數(shù) 1(1 ( ) )S xa??? ??， 同樣地對(duì) 每一天每一時(shí)段兩種方法預(yù)測(cè)的雨量進(jìn)行插值，得到兩種方法在 91 個(gè)站點(diǎn)上的預(yù)測(cè)雨量并按照給定的范圍等級(jí)化，求出兩種方法與實(shí)測(cè)數(shù)據(jù)的等級(jí)差 x ， 并代入確定的隸屬函數(shù) x? ?， 最終將各天 各 時(shí)段各站點(diǎn)的滿意度加起來(lái)，即求得方法一、方法二的總滿意度 1S 、 2S 。 FORECAST 中的文件 和 分別包含網(wǎng)格點(diǎn)的經(jīng)緯度，其余文件名為 f 日期 i_dis1 和 f 日期 i_dis2，例如 f6181_dis1 中包含 2022 年 6 月 18日 晚上 20 點(diǎn) 采用第一種方法預(yù)報(bào)的第一時(shí)段數(shù)據(jù)（其 2491 個(gè)數(shù)據(jù)為該時(shí)段各網(wǎng)格點(diǎn)的雨量），而 f6183_dis2 中包含 2022 年 6 月 18 日 晚上 20 點(diǎn) 采用第二種方法 預(yù)報(bào)的 第三時(shí)段數(shù)據(jù)。綜上所述，方法一比方法二更為準(zhǔn)確； 問(wèn)題二 模型建立 （一）目標(biāo)確定 通過(guò) 引入偏小型柯西隸屬分布函數(shù)， 求出某天某時(shí)段某個(gè)站點(diǎn)的滿意度，然后對(duì) 41 天的 4 個(gè)時(shí)段的 91 個(gè)站點(diǎn)的滿意度進(jìn)行相加， 求出兩種方法的總滿意度 1S 和 2S ； （二） 模型準(zhǔn)備 引入偏小型柯西隸屬分布函數(shù) 為 了 描 述 滿 意 度 ， 我 們 引 入 了 偏 小 型 柯 西 隸 屬 分 布 函 數(shù) 1(1 ( ) )S xa??? ??，其中橫坐標(biāo) x 表示降雨量等級(jí)差： 當(dāng) x =1 時(shí)， S =,； 當(dāng) x =6 時(shí)， S =； 令 a =0 把以上兩點(diǎn)代入函數(shù)中科解得 ? =， ? =，即柯西隸屬分布函數(shù)為 x? ? 函數(shù)圖像如下圖： 分類討論 為了防止誤差太大，我們細(xì)化了幾種情況： 方法一、二的站點(diǎn)預(yù)測(cè)降雨量 ① 雨 量 差 和 等 級(jí) 差 ② 雨 量 差 和 等 級(jí) 差 ③ 雨 量 差 和 等 級(jí) 差 ~6 ~ 1 ~6 2 ~25 ~6 1 6~12 2 12~25 3 25以 上 ~10 1 10~20 2 20~30 4 40~50 4 50~60 5 60~ 6 （三）模型的建立 通過(guò)循環(huán)分別讀取出某天某時(shí)段的 91 個(gè)站點(diǎn)的實(shí)測(cè)雨量，和網(wǎng)格點(diǎn)預(yù)測(cè)雨量，用散亂節(jié)點(diǎn)插值函數(shù) griddata 插值求得某天某時(shí)段 91 個(gè)站點(diǎn)的預(yù)測(cè)雨量，由此便得到某天某時(shí)段 91 個(gè)站點(diǎn)的實(shí)測(cè)雨量 rqm 、方法一預(yù)測(cè)雨量 _1rqm f 、方法二預(yù)測(cè)雨量_2rqm f ，接著讓 _1rqm f 、 _2rqm f 和 rqm 按照題目給出的范圍進(jìn)行雨量的等級(jí)劃分，然后對(duì)方法一、二分別求出等級(jí)誤差 1a 、 2a 和雨量誤差 1c 、 2c ； 若等級(jí)誤差為 0，則滿意度為 1，否則，按照分類討論中的表格進(jìn)行等級(jí)誤差的確定，并代入函數(shù) x? ?求出某站點(diǎn)的滿意度 ；最后把 41 天的 4 個(gè)時(shí)段的 91 個(gè)站點(diǎn)的滿意度分別加起來(lái)得到兩種方法的總滿意度 1S 和 2S 模型求解 （一） Matlab 求解 1Step :利用 dlmread 函數(shù) 從 6 月 18 日到 7 月 30 日循環(huán)讀取實(shí)測(cè)站點(diǎn)數(shù)據(jù)，并賦給矩陣 tmp ， 從矩陣 tmp 提取出第
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1