雨量預(yù)報方法的評價論-wenkub
2023-01-23 11:38:54 本頁面
【正文】 格點降雨量 2rqf 某天某個時段方法二預(yù)測的網(wǎng)格點降雨量 r1 rqm_f 方法一插值得到的 91個站點預(yù)測雨量 r2 rqm_f 方法二插值得到的 91個站點預(yù)測雨量 1a 方法一的等級差 2a 方法二的等級差 1c 方法一的雨量差 2c 方法二的雨量差 12xx、 方法一、二對應(yīng)的等級差 1s 方法一某天某時段某站點的滿意度 2s 方法二某天某時段某站點的滿意度 12SS、 方法一、二的總滿意度 1ss 方法一某天某時段的滿意度 2ss 方法二某天 某時段的滿意度 1sss 方法一某天的滿意度 2sss 方法二某天的滿意度 分析與建立模型 問題一 模型建立 （一）目標(biāo)確定 比較方法一、方法二預(yù)測的準(zhǔn)確性，總體來看哪種方法預(yù)測的誤差較小。 MEASURING 中包含了 41 個名為 日期 .SIX 的文件，如 表示2022 年 6 月 18 日 晚上 21 點開始的連續(xù) 4 個時段 各站點的實測數(shù)據(jù) （雨量） ，這些文件的數(shù)據(jù)格式是： 站號 緯度 經(jīng)度 第 1 段 第 2 段 第 3 段 第 4 段 58138 58139 58141 58143 58146 …… 雨量用毫米做單位，小于 毫米視為無雨。 氣象部門希望建立一種 科學(xué)評價預(yù)報方法好壞的數(shù)學(xué)模型與方法 。 由 Matlab 編程解得 1S =14394 ， 2S =14923 ， 即方法二的滿意度高于方法一 。雨量預(yù)報方法的評價 摘要 本文研究了 雨量預(yù)報的準(zhǔn)確性問題 ,運用 Matlab 軟件編程， 通過插值和對比計數(shù) 的方法 ，比較方法一和方法二 預(yù)測 的準(zhǔn)確性；問題二研究的是公眾對雨量預(yù)報的滿意程度，運用模糊數(shù)學(xué)的知識，將模糊的評價等級數(shù)字化 ，用偏小型柯西隸屬函數(shù)來刻畫公眾對兩種方法的滿意度。 關(guān)鍵詞 : 插值 偏小型柯西隸屬函數(shù) 對比計數(shù) 問題重述 雨量預(yù)報對農(nóng)業(yè)生產(chǎn)和城市 工作和 生活有重要作用，但準(zhǔn)確、及時地對雨量作出預(yù)報是一個十分困難的問題，廣受世界各國 關(guān)注。氣象部門提供了 41 天的用 兩種不同方法 的 預(yù)報數(shù)據(jù) 和相應(yīng)的 實測數(shù)據(jù) 。 (1) 請建立數(shù)學(xué)模型來 評價 兩種 6 小時雨量預(yù)報方法的 準(zhǔn)確性 ； (2) 氣象部門將 6 小時降雨量分為 6 等： — 毫米為小雨， — 6 毫米為中雨， — 12 毫米為大雨， — 25 毫米為暴雨， — 60 毫米為大暴雨，大于 毫米為特大暴雨。 （二 ）模型的建立 分別以經(jīng)、緯度為 x、 y 坐標(biāo)建立平面直角坐標(biāo)系，得到網(wǎng)格點與實測點的相對位置，如下圖： 通過循環(huán)分別讀取出某天某時段的 91 個站點的實測雨量，和網(wǎng)格點預(yù)測雨量，用散亂節(jié)點插值函數(shù) griddata 插值求得某天某時段 91 個站點的預(yù)測雨量，由此便得到某天某時段 91 個站點的實測雨量 rqm 、方法一預(yù)測雨量 _1rqm f 、方法二預(yù)測雨量_2rqm f ，接著讓 _1rqm f 、 _2rqm f 分別和 rqm 進(jìn)行比較， 91 個站點循環(huán)完后，比較1n 與 2n ，即可得到某天某時段中更準(zhǔn)確的方法，再進(jìn)行 4 個時段的循環(huán)，比較 1N 和 2N 即可得到某天內(nèi)更準(zhǔn)確的方法，最后再對 41 天循環(huán)，比較 1M 和 2M ， 即可得出 41 天中更為準(zhǔn)確的方法；具體實施 步驟參考模型求解。 附錄 問題一程序 ： clear all clc M1=0。,39。]。%實測經(jīng)度 rqm=tmp(:,4:7)。39。,39。N2=0。f39。%第 1 方法預(yù)測雨量 41x4 ff2=[39。,num2str(i),num2str(k),39。%第 1 方法預(yù)測雨量 53x47 rqf2=dlmread(ff2)。%插值得 91 個站點的預(yù)測值 rqm_f2=griddata(lonf,latf,rqf2,lonm,latm,39。n2=0。 else if ab n2=n2+1。 else if N2N1 M2=M2+1。02039。 tmp=dlmread(fm)。 r=rqm。 r(r= amp。 r(r=12 amp。 r(r=60)=6。39。,39。F:\5\data\FORECAST\39。_dis139。,39。]。V439。)。 r1(r1= amp。 r1(r1=6 amp。 r1(r1=25 amp。 r2(r2= amp。 r2(r2=6 amp。 r2(r2=25 amp。 a2=abs( r2(j)r(j,k) )。 s1(j)=1./(1+*x1(j).^)。 s1(j)=1./(1+*x1(j).^)。 end end else if rqm_f1(j)=6 amp。 else if c1=6 amp。 c1= x1(j)=3。 s1(j)=1./(1+*x1(j).^)。 else if c1=20 amp。c140
公司管理相關(guān)推薦
鄂ICP備17016276號-1