第三次交卷我們仔細(xì)閱讀了生數(shù)學(xué)建模競賽的規(guī)則

大學(xué)生數(shù)學(xué)建模賽區(qū)評(píng)閱編號(hào)（由賽區(qū)評(píng)閱前進(jìn)行編號(hào)賽區(qū)評(píng)閱記錄（可供賽區(qū)評(píng)閱時(shí)使用統(tǒng)一編號(hào)（由賽區(qū)送交前編號(hào)評(píng)閱編號(hào)（由評(píng)閱前進(jìn)行編號(hào)本文根據(jù)近十年（1999.1-2008.6）震級(jí)在4.0級(jí)以上的統(tǒng)計(jì)數(shù)據(jù)（見附錄中表1）和全球各帶分布情況（見圖1，利用均生函數(shù)和方差分析的數(shù)學(xué)建模方法“分級(jí)擬合方法”提高預(yù)測準(zhǔn)確率，對(duì)各帶中短期趨勢進(jìn)行了預(yù)測。1帶分布情況進(jìn)行數(shù)據(jù)整理，得到三其次，考慮到預(yù)測需要對(duì)趨勢預(yù)測基本準(zhǔn)確，我們引入“雙評(píng)分準(zhǔn)則，,即在擬合時(shí)考慮了數(shù)量預(yù)報(bào)和趨勢預(yù)報(bào)的平衡問題,選取合適的調(diào)整系數(shù),從而提高了地震中期趨勢預(yù)測的準(zhǔn)確率.另外采用了“分級(jí)擬合”的思路,利用殘差序列逐級(jí)提高擬CSC[1

N(k)

]*Q(kCSCm短 預(yù)測趨勢，其中短 趨勢預(yù)測模型為：X(Na)=X(i1)(Nm最后我們對(duì)2008年下半年和2009年的趨勢進(jìn)行了預(yù)測并對(duì)預(yù)測結(jié)果進(jìn)行了評(píng)價(jià)，結(jié)果顯示預(yù)測結(jié)果基本符合實(shí)際情況，準(zhǔn)確率在72%以上。三個(gè)帶的預(yù)測準(zhǔn) 預(yù) 趨 均生函數(shù)雙評(píng)分準(zhǔn) 分級(jí)擬

一、問題的重，對(duì)人類的生命財(cái)產(chǎn)具有相當(dāng)大的危害性，5.12汶川造成了巨大的損失截至6月23日12時(shí)已造成69181人遇難，374171人受傷，18498人對(duì)以往活動(dòng)進(jìn)行分析可以幫助人類更清楚的了解的規(guī)律，見附錄表一：國家臺(tái)網(wǎng)中心發(fā)布的震級(jí)4.0以上的記錄數(shù)據(jù)二、問題的分周期分析是預(yù)測的法。從序列本身出發(fā),我們可以利用各種中我們可以借鑒氣象學(xué)中用于中期預(yù)測的數(shù)學(xué)方法,建立中短期趨勢預(yù)測模結(jié)合預(yù)測在某些情況下只要求趨勢準(zhǔn)確而非數(shù)值精確的特點(diǎn)提出一種FAICBIC分準(zhǔn)則CSC(CoupleScoreCriterion)是從數(shù)量預(yù)報(bào)評(píng)分和趨勢預(yù)報(bào)評(píng)分兩個(gè)示數(shù)量評(píng)分,稱之為精評(píng)分;S2三、模型準(zhǔn)通過對(duì)附錄表一的數(shù)據(jù)處理，我們得 帶主要分布情況然后帶劃分三個(gè)區(qū)，分別為帶一-----東亞-南亞帶二-----西太平洋帶三----- 帶最大震級(jí)-時(shí)間統(tǒng)計(jì)結(jié)

四、模型假

洋上下上7下上6下6上下上下8上下上7下上8下8上下上8下1、0.50.52、假設(shè)不同趨勢的判定標(biāo)準(zhǔn)為0.5，即當(dāng)后一次的震級(jí)預(yù)測值、實(shí)際值分別五、符號(hào)說L：：調(diào)整系數(shù)，表示數(shù)值預(yù)報(bào)與趨勢預(yù)報(bào)的權(quán)重模型分

六、時(shí)間與震級(jí)的中短期預(yù)測模第k級(jí)擬合(k>=2)是對(duì)第k-1級(jí)的殘差序列進(jìn)行擬合。在對(duì)序列進(jìn)行分級(jí)擬合 模型的建N設(shè)原始震級(jí)時(shí)間序列為XXtNXX1,X

Xj

XN

jjwi+1(i=0,1...m)L(i+1),i+1XXN殘差數(shù)據(jù)序列 N

XX

t,XX

XX

(

XX

(N)j1,

X值 X

(NaX

(N

X

(b

NXXN

XX1,X

Xj

XN

1)1)X1,X1X1(X(NN++X2+X2+X2(+X(NN+++++++i1 Xi1(1),Xi1 X+( Xi1++++++++++

(N

XX

X

N,N

X(XX(X(NN

X(j)X

mXXtm

XX

NNNXN+a(a>=1)個(gè)時(shí)間段（a）的預(yù)X(NaN模型——時(shí)間與震級(jí)的中短期預(yù)測模型m未來第a個(gè)時(shí)間段 最大震級(jí)預(yù)測值為：X(Na)=Xm

(NXX1,X

Xj

XN

X

X(i1)(1),X(i1)

X

(

X

(N

X(i1)(Na)X(i1)(bjL(i

mXXt{X(i1)

模型的求NXXtNXX1,X

Xj

XN

Xjj，N

1

1LX =

(jvL) MnL

為滿足nNMN L1MH

XX

X對(duì)于的中長期預(yù)測要求是趨勢基本正確，因此我們利用了雙評(píng)分準(zhǔn)則2用隱式表達(dá)：CSC=S1+S2,S1:數(shù)值預(yù)報(bào)評(píng)分（精評(píng)分)，S2:CSC[1NN(k

N(k)

]*Q(k)Q(k2 2NQ(k) [(XNt

(t)X3N(k){11C[X(t)]C[

t

(t)X(t1)

(t)，XL(t)XL(t1)XL(t)X(t)X(t)

X(t)U,UCU

X(t)（其中：UU=0.5。為調(diào)整系數(shù)，表示數(shù)值預(yù)報(bào)與趨勢預(yù)報(bào)的權(quán)重，我們假設(shè)兩項(xiàng)權(quán)重相等，即=1.0）計(jì)算CSC值，取CSC值最小的周期為一級(jí)擬合的最優(yōu)周期。采用分級(jí)擬合提高擬合的精度，即在第一級(jí)擬合中選取CSC值最小的LDD

=N

， NN

表示在第一級(jí)擬合中擬合誤差<=d的元素?cái)?shù)量X1X

(t)=

(t)

(t)；（其中：L(1)為第一級(jí)擬合最優(yōu)周期D對(duì)殘差序列重復(fù)步驟a選取第二級(jí)擬合最優(yōu)周期求 D

DD

，則繼續(xù)進(jìn)行下一級(jí)的殘差擬合； DD

d.d.

(t

X

Xm

(t 中的元素值相加，作為原始序 的擬合序

...

存在長L期，那么由此周期構(gòu)成的新序列是X(1X(2X(L)

(L1)

(L2)...

(2LX(1)X(L1)和

(2L

、

(L2)X(2L2FLF(F

S(L)/(L S/(NS(

L L2nL[XL(j)X]2jLnLSj1

[jvL]

為滿足

NFL)

LLL 序列而言可能存在多個(gè)周期計(jì)算不同周期下的F(L)值將通過檢驗(yàn)的F(L)值從大到小排序，根據(jù)實(shí)際需要選取所需的主要優(yōu)勢周期。5.3.5以1999- -2008的數(shù)據(jù)資料作為原始數(shù)據(jù)對(duì)未來三年進(jìn)行預(yù)測并檢，通過與已有數(shù)據(jù)的比較‘東亞—南亞’大陸斷裂帶和西太平洋帶半年段預(yù)測震72%83%80%。根據(jù)上表預(yù)測的未來一年半東亞-南亞、西太平洋最大震級(jí)強(qiáng)度及趨勢，未來三年東上下上下東亞—降平平西太平洋平平平上下上下上下東洋平降平

七、模型的優(yōu)缺點(diǎn)分1、我們在的預(yù)測模型中，用到了均值生成函數(shù)法（等效于移動(dòng)平均數(shù)法）去構(gòu)造新序列，為趨勢的分析和預(yù)測提供了一種很好的方法23預(yù)報(bào)和趨勢預(yù)報(bào)的平衡問題，從而提高了中期趨勢預(yù)測的準(zhǔn)確率4、用F八、參考文董肇君，系統(tǒng)工程與運(yùn)籌學(xué)（第2版 ：國防工 朝倫巴根，賈德彬，數(shù)值計(jì)算方法 ：中國水利水 姜啟源，謝金星，葉俊，數(shù)學(xué)模型（第三版）[M].：高等教育。[4]馬文靜，張曉清，張，均生函數(shù)在預(yù)測模型中的應(yīng)用，西北學(xué)報(bào)，221；2000。

九、附緯 經(jīng)深 震內(nèi)鄂倫春自治旗與阿榮旗--云南盈江縣與交界地洛浦與日土交界地中海----蘇門答臘西--福建省市永春---青海海西、交內(nèi)巴林右-8蘇門答臘西--河北河北內(nèi)牙克石俄羅斯堪察加地堪察加?xùn)|北地青海治多與班戈青海治多與交界---內(nèi)陳巴爾虎旗牙克石市阿瓦提與云南會(huì)澤與會(huì)東交-內(nèi)烏拉特中--重慶市石柱土自治---內(nèi)東烏珠穆沁島花蓮東俄、蒙、俄、蒙、俄、蒙、-Ms5.9內(nèi)巴林、阿魯旗0Ms5.0花-0中國-0內(nèi)托克00000伽師0青海交0000甘省000內(nèi)西烏珠穆沁0000000--0000000000000-0琉球群島西-澳大利亞以00青海交界（境內(nèi)若羌00-000000---000內(nèi)五-000000000州00000000000000000000000000本州以0000000000-0洋0000000-000000000000-00-000000-0000000洋0000000000000000000000000000000000000000000000000中國與交000000花蓮東-000內(nèi)錫林浩特 0Ms5.0 105.1105759249 335538 1041055425 964105464-21331054445 33 3310510565445104121105105105104105531 531 5 31 15 15631531 531 63124 -7- - 38102333 555535 - 7- - 100- 1002.89633 - 0.9 3333-6-7--7--415- 150- -- -- 72.6 - 7- - - 101- - - - - 3- 7 - - - - - 6 7 4 23 10133336--22 38 1.2 2458120-151712028 -61-120105 28 33333315- 1055 0.1 - - 5-20-174.2 49 24 61 8 23 - 120- -20-178.2540---7--2440 33 15 1258 2.1 -- - 1222.197330.6 7 793.1 72.29733 5.2 7 6.8 3.9 - - - - 4.6- - 242424 151397995123--75--71374---755-122571555343333-12223121107237526 29128105014486-7505--26 -3333337--72702.6 7 7 5 05 05 0 0 05 5 0 -- 0 05 77 5 0 7 -03.1 0 117- 0 0 0 0525 16025 160 54024-5024 050--7--051000424 94-0000000076.8 16010 -16810-1607 123-12536009410941010-076-4 705--0700750980-----16- 01210--744 7008700005241210 -71.9 6-070-0-82457-07.0070050-0-0-15400-70500500524 005000-07701407005000600-07-0-70-1027277104-12124 -0000-19 700508905070005000--070-07-050001230-0505000121121-0-024 0-01215001215050512150012150024 50012124005000507000000509.1- 00000-0042 542 05 029050-7-707-05050006-70000-16700 1010fori=1:702 forholdplot(dian1(:,3)*750/360+375,195-fori=1:ls 1999)+floor((mod(floor(dian(i,1)/100),100)- forasd=1:4forl=2:floor(n/2)fori=1:lforj=0:floor(n/l)-1fori=2:nifyue(i)-yue(i-1)>0.5ifyue(i)-yue(i-1)<-0.5ifxl(mod(i,l))-