。
[0017] (4)、返回步驟(1)，并使n=n+l，直至得出n+m年的糧食消費(fèi)量。
[0018] 實(shí)施例1，以糧食消費(fèi)量為口糧消費(fèi)量，n=2007，m=2為例進(jìn)行詳細(xì)介紹，若干年選 擇1978年~2007年這30年，這30年間的口糧消費(fèi)量和與口糧消費(fèi)量有關(guān)的變量的數(shù)值為已 知數(shù)據(jù)，記1978年為第1年，2007年為第30年，其中口糧消費(fèi)量和與口糧消費(fèi)量有關(guān)的變量 的數(shù)值可以從《中國農(nóng)村統(tǒng)計(jì)年鑒》查詢到。
[0019] 如圖1和圖2所示的糧食消費(fèi)量動(dòng)態(tài)預(yù)測方法，其中依次包括以下步驟： (1)確定口糧消費(fèi)量的影響因子:根據(jù)第1~30年口糧消費(fèi)量數(shù)值和第1~30年與口糧消 費(fèi)量變化有關(guān)的國內(nèi)人口數(shù)量、城鎮(zhèn)化水平、恩格爾系數(shù)以及農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)等變量 的數(shù)值，計(jì)算這些變量和口糧消費(fèi)量之間的關(guān)聯(lián)度。在判斷變量是否為影響因子的時(shí)候，需 要看該變量和口糧消費(fèi)量之間的關(guān)聯(lián)度是否大于〇. 7，大于0.7則該變量對口糧消費(fèi)量的變 化影響較大其為影響因子。在對該預(yù)測方法進(jìn)行驗(yàn)證的時(shí)候，發(fā)現(xiàn)大于0.7的關(guān)聯(lián)度的取值 對于提高預(yù)測的準(zhǔn)確度具有很大的影響作用。通過以上方法，最終確定國內(nèi)人口數(shù)量、城鎮(zhèn) 化水平、恩格爾系數(shù)以及農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)這4個(gè)變量為口糧消費(fèi)量的影響因子。其中變 量和口糧消費(fèi)量之間的關(guān)聯(lián)度的計(jì)算方法會(huì)在步驟(2)中給出。在此，選取變量和口糧消費(fèi) 量的關(guān)聯(lián)度作為判斷該變量是否為影響因子，可以提高判斷地準(zhǔn)確性，避免隨機(jī)年份的誤 差影響預(yù)測結(jié)果。
[0020] 人口數(shù)量、城鎮(zhèn)化水平、恩格爾系數(shù)以及農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)和口糧消費(fèi)量的關(guān) 聯(lián)度如表1所示： 表1影響因子和口糧消費(fèi)量之間的關(guān)聯(lián)度
同時(shí)從表1可以看出，農(nóng)村和城鎮(zhèn)有些影響因子和口糧消費(fèi)量之間的關(guān)聯(lián)度是不同的。 為了提高預(yù)測結(jié)果的準(zhǔn)確率，口糧消費(fèi)量為農(nóng)村口糧消費(fèi)量或城鎮(zhèn)口糧消費(fèi)量。即，在預(yù)測 的時(shí)候?qū)⑥r(nóng)村口糧和城鎮(zhèn)口糧區(qū)分開來預(yù)測，通過這種方法可以避免因兩者的影響因子和 口糧消費(fèi)量之間的關(guān)聯(lián)度的不同而造成的預(yù)測結(jié)果誤差較大的問題。
[0021] (2)計(jì)算第卜30年口糧消費(fèi)量和各個(gè)影響因子之間的關(guān)聯(lián)度，其中各個(gè)影響因子 僅以農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)這個(gè)影響因子為例進(jìn)行詳細(xì)描述，其它的影響因子和口糧消費(fèi)量 之間的關(guān)聯(lián)度的計(jì)算方法依次類推即可。第1~30年口糧消費(fèi)量和農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)之間 的關(guān)聯(lián)度的計(jì)算方法如下： 第一步:建立口糧消費(fèi)量時(shí)間數(shù)列和農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)時(shí)間數(shù)列： 對口糧消費(fèi)量建立時(shí)間數(shù)列，2,3, .-,30);其中kt為第t年的口 糧消費(fèi)量;ki為第一年的口糧消費(fèi)量;k2為第二年的口糧消費(fèi)量;k3Q為第30年的口糧消費(fèi) 量； 對影響因子的值建立時(shí)間數(shù)列6=(心汀2、'"43〇)&=1，2,3，'"，30);其中匕為第七年的 農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù);f\為第一年的農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù);f2為第二年的農(nóng)產(chǎn)品生產(chǎn)價(jià)格指 數(shù);f 3〇為第30年的農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)； 第二步:無量綱化:對口糧消費(fèi)量建立時(shí)間數(shù)列kt和影響因子建立時(shí)間數(shù)列ft均進(jìn)行 無量綱化處理后得出無量綱口糧消費(fèi)量時(shí)間數(shù)列PjP無量綱農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)時(shí)間數(shù) 列f\。其中，對口糧消費(fèi)量的時(shí)間數(shù)列k t和農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)的時(shí)間數(shù)列ft進(jìn)行無量綱 化處理的方法為:將口糧消費(fèi)量時(shí)間數(shù)列k t中的值分別與kt中的初始值（SPh)相除而得到 無量綱口糧消費(fèi)量時(shí)間數(shù)列k\;農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)時(shí)間數(shù)列f t中的值分別與ft中的初始 值(即f\)相除而得到無量綱農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)時(shí)間數(shù)列ft 1。
[0022] 第三步:求絕對差值:求無量綱口糧消費(fèi)量時(shí)間數(shù)列k\和無量綱農(nóng)產(chǎn)品生產(chǎn)價(jià)格 指數(shù)時(shí)間數(shù)列f\在第t年的絕對差值序列fk t，并從fkt中找出絕對差值的最大值fkmax和絕 對差值的最小值fk min; 第四步:求出關(guān)聯(lián)系數(shù):根據(jù)公式Lt=(fkmax+fkmin)/(fkt+fk max)得出農(nóng)產(chǎn)品生產(chǎn)價(jià)格指 數(shù)在對應(yīng)年份的關(guān)聯(lián)系數(shù)，其中Lt為第t年農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)和口糧消費(fèi)量的關(guān)聯(lián)系數(shù)； 第五步:求出關(guān)聯(lián)度:根據(jù)公式Κ=αΛ 2+···+Ι^)/30得出農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)與口糧 消費(fèi)量的關(guān)聯(lián)度，其中R為關(guān)聯(lián)度為第1年農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)和口糧消費(fèi)量的關(guān)聯(lián)系 數(shù)，L 2為第2年農(nóng)產(chǎn)品生產(chǎn)價(jià)格指數(shù)和口糧消費(fèi)量的關(guān)聯(lián)系數(shù)，L3Q為第30年農(nóng)產(chǎn)品生產(chǎn)價(jià)格 指數(shù)和口糧消費(fèi)量的關(guān)聯(lián)系數(shù)。
[0023] (3)預(yù)測第31年的各個(gè)影響因子的值:根據(jù)第1~30年影響因子的數(shù)值，預(yù)測第31年 各個(gè)影響因子的數(shù)值。其中預(yù)測第31年的各個(gè)影響因子的值的方法為：利用ARMA(p，d，q) 模型預(yù)測出第31年影響因子的值，方法為： 第一步，輸入第卜30年影響因子的數(shù)值，對其進(jìn)行平穩(wěn)性檢驗(yàn)，若平穩(wěn)，則可直接進(jìn)行 第三步;若非平穩(wěn)，則進(jìn)行第二步； 第二步，原始數(shù)據(jù)進(jìn)行差分處理，直至序列達(dá)到平穩(wěn)； 第三步，建立預(yù)測模型； 第四步，對預(yù)測模型進(jìn)行殘差檢驗(yàn)，確定最終模型； 第五步，根據(jù)第四步確定的最終模型預(yù)測第31年影響因子的值。
[0024] 其中，平穩(wěn)性檢驗(yàn)、對原始數(shù)據(jù)進(jìn)行差分處理、建立預(yù)測模型、對預(yù)測模型進(jìn)行殘 差檢驗(yàn)確定最終模型，根據(jù)最終模型預(yù)測第31年各個(gè)影響因子的值的方法為現(xiàn)有技術(shù)。
[0025] 在本方法中，通過加入殘差檢驗(yàn)可以提高預(yù)測出的影響因子的數(shù)值的準(zhǔn)確性，進(jìn) 而提高口糧消費(fèi)量的預(yù)測準(zhǔn)確率。
[0026]根據(jù)以上步驟，預(yù)測出的第31年即2008年城鎮(zhèn)各個(gè)影響因子的數(shù)值如表2所示： 圭9滿、、腦丨山的畝義木宮5
古 (4) 建立多元回歸模型：根據(jù)口糧消費(fèi)量和影響因子之間的關(guān)系，建立多元回歸模型 731=1^+1^(+11:^+11^2+1^13+ε，其中y3i為第31年口糧消費(fèi)量;X、Xi、X 2、χ3為預(yù)測出的第31年 各個(gè)影響因子的值;υ〇、ΙΛ、U2、U 3、U4和ε為根據(jù)第1~30年影響因子和口糧消費(fèi)量之間的關(guān)聯(lián) 度計(jì)算出的回歸參數(shù)和隨機(jī)誤差； (5) 計(jì)算回歸參數(shù)和隨機(jī)誤差：由第1~30年口糧消費(fèi)量和影響因子之間的關(guān)聯(lián)度得出 回歸參數(shù)11〇、111、1]2、1]3、1]4和隨機(jī)誤差￡，方法為:米用最小二乘法確定，通過最小二乘法確定 多元回歸方程的回歸參數(shù)和隨機(jī)誤差為現(xiàn)有技術(shù)。
[0027] 在預(yù)測的時(shí)候，回歸參數(shù)和隨機(jī)誤差根據(jù)影響因子和口糧消費(fèi)量之間的關(guān)聯(lián)度的 變化而變化，從而有別于現(xiàn)有的預(yù)測方法中始終利用一個(gè)回歸參數(shù)和一個(gè)隨機(jī)誤差進(jìn)行預(yù) 測，進(jìn)而提高了預(yù)測的準(zhǔn)確率。
[0028] (6)預(yù)測第31年的口糧消費(fèi)量：由公式ε，得出第31年 的口糧消費(fèi)量731441、辦、乂3為預(yù)測出的第31年各個(gè)影響因子的值； (7)預(yù)測第32年口糧消費(fèi)量:根據(jù)步驟(3)中預(yù)測的第31年的各個(gè)影響因子的值和步驟 (6 )中預(yù)測出的第31年即2008年度的口糧，返回步驟(2 )得出第1~31年影響因子和口糧消費(fèi) 量之間的關(guān)聯(lián)度，并根據(jù)第1~31年影響因子按步驟(3)預(yù)測出第32年即2009年度的各個(gè)影 響因子的值，根據(jù)第1~31年影響因子和口糧消費(fèi)量之間的關(guān)聯(lián)度按照步驟(4)、（5)計(jì)算出 新的回歸參數(shù)和隨機(jī)誤差，將預(yù)測出的第32年即2009年度的各個(gè)影響因子的值和新的回歸 參數(shù)和隨機(jī)誤差代入第(4)步建立的模型中得出第32年即2009年度的口糧消費(fèi)量。
[0029] 在進(jìn)行后面年份預(yù)測的時(shí)候，對步驟(2)、（3)、（4)、（5)、（6)依次重復(fù)即可，從而實(shí) 現(xiàn)口糧消費(fèi)量的長期預(yù)測。
[0030] 根據(jù)