運(yùn)用本底趨勢(shì)線模型評(píng)估2010年上海世博會(huì)對(duì)上海旅游

1、運(yùn)用本底趨勢(shì)線模型評(píng)估2010年上海世博會(huì)對(duì)上海旅游經(jīng)濟(jì)的影響力作者：信陽(yáng)師院數(shù)學(xué)07級(jí)4班三位同學(xué)摘要 本文主要采用本底趨勢(shì)線模型討論了世博會(huì)對(duì)上海旅游業(yè)的影響，本文建立了兩個(gè)模型。模型一首先對(duì)上海市2000 - 2009年旅游業(yè)的客流量、旅游業(yè)的收入和上海市的GDP等統(tǒng)計(jì)數(shù)據(jù)進(jìn)行內(nèi)插處理，采用本地趨勢(shì)線方法分別建立了上海市旅游業(yè)客流總量、旅游收入和GDP的預(yù)測(cè)模型，誤差分析表明該本底趨勢(shì)線模型得到的本底值能夠以非常高的精度和實(shí)際數(shù)據(jù)相吻合，這表明模型很好的反映了上海市旅游業(yè)的長(zhǎng)期發(fā)展趨勢(shì)和周期性，然后用該模型預(yù)測(cè)的2010年上海市旅游業(yè)客流量、收入和GDP的本底值，把2010年的預(yù)測(cè)值和實(shí)

2、際值做差，把差值看作世博會(huì)對(duì)上海市2010年旅游業(yè)的影響力，結(jié)果顯示世博會(huì)對(duì)上海市2010年的旅游業(yè)的各項(xiàng)指標(biāo)都有很大影響。模型二把上海市2000-2010年旅游業(yè)的的數(shù)據(jù)進(jìn)行內(nèi)插處理，建立了本地趨勢(shì)線預(yù)測(cè)模型，用該模型了預(yù)測(cè)了2011、2012和2013年上海市旅游業(yè)各項(xiàng)指標(biāo)的值，再把這些預(yù)測(cè)值與用模型一預(yù)測(cè)得到的2011、2012、2013的數(shù)值做差，來(lái)評(píng)價(jià)世博會(huì)對(duì)上海市2010年后旅游業(yè)的影響，結(jié)果顯示世博會(huì)對(duì)上海市2010后的旅游業(yè)仍有較大影響，但是比對(duì)2010年上海市旅游業(yè)的影響力要小。關(guān)鍵詞 本底趨勢(shì)線；旅游；世博會(huì)；預(yù)測(cè)1 問(wèn)題重述2010年上海世博會(huì)是首次在中國(guó)舉辦的世界博覽會(huì)

3、從1851年倫敦的“萬(wàn)國(guó)工業(yè)博覽會(huì)”開(kāi)始，世博會(huì)正日益成為各國(guó)人民交流歷史文化、展示科技成果、體現(xiàn)合作精神、展望未來(lái)發(fā)展等的重要舞臺(tái)請(qǐng)你們選擇感興趣的某個(gè)側(cè)面，建立數(shù)學(xué)模型，利用互聯(lián)網(wǎng)數(shù)據(jù)，定量評(píng)估2010年上海世博會(huì)的影響力。2 問(wèn)題提出與分析重大活動(dòng)對(duì)經(jīng)濟(jì)的特定區(qū)域，如旅游、運(yùn)輸?shù)扔兄卮蟮挠绊懀?dāng)這些重大活動(dòng)產(chǎn)生時(shí)，就能影響到特定時(shí)間內(nèi)（通常統(tǒng)計(jì)量是一年）經(jīng)濟(jì)活動(dòng)的總量，如國(guó)內(nèi)生產(chǎn)總值（GDP）。上海市是我國(guó)入境旅游業(yè)最發(fā)達(dá)的熱點(diǎn)城市之一,其接待的客流量和旅游外匯收入一直處于全國(guó)熱點(diǎn)城市排行榜的前幾位,旅游業(yè)在本市國(guó)民經(jīng)濟(jì)中占有極其重要的地位.旅游經(jīng)濟(jì)的發(fā)展主要有客流量決定，如何定量分析

4、客流量與旅游經(jīng)濟(jì)的關(guān)系？在舉行世博的情況下，將如何影響旅游業(yè)的發(fā)展？因此,本底趨勢(shì)線模型可以反映在沒(méi)有特殊事件和因素的影響下,旅游業(yè)發(fā)展的自然趨勢(shì)。就重大事件活動(dòng)而言,其性質(zhì)是一種超常規(guī)的特殊因素,對(duì)于上海世博會(huì)這類一次性的重大活動(dòng)來(lái)說(shuō),其特殊性就更加突出,因而,應(yīng)用孫根年12的本底趨勢(shì)線模型方法分析重大事件活動(dòng)對(duì)旅游業(yè)的影響是合適的.本文將利用本底趨勢(shì)線模型預(yù)測(cè)2010年旅游業(yè)的客流量、旅游收入和GDP，把這個(gè)預(yù)測(cè)值看作是在沒(méi)有世博會(huì)的情況下旅游業(yè)的自然增長(zhǎng)情況，然后把這個(gè)預(yù)測(cè)值和世博會(huì)帶來(lái)的旅游收入的值進(jìn)行做差，把這個(gè)差值和自然增長(zhǎng)值的比值看作世博會(huì)的影響力。3 符號(hào)說(shuō)明代表所要訂正的年份

5、；為第年客流量統(tǒng)計(jì)值 ；表示第年客流量?jī)?nèi)插值；表示第年旅游收入統(tǒng)計(jì)值；表示第年旅游收入內(nèi)插值；為進(jìn)行直線內(nèi)插的公差值；t 為時(shí)間變量，從2000年開(kāi)始依次取t=1，2，3；、表示直線內(nèi)插的起始點(diǎn)、終止點(diǎn)；表示最終本底趨勢(shì)線模型方程計(jì)算的本底值與內(nèi)插值的相關(guān)系數(shù)。4 模型建立與求解對(duì)本底趨勢(shì)線的自然延伸還可預(yù)測(cè)旅游業(yè)發(fā)展的基本趨勢(shì)，而建立旅游本底趨勢(shì)線需要以長(zhǎng)期的旅游統(tǒng)計(jì)數(shù)據(jù)為依據(jù). 本文所采集的統(tǒng)計(jì)數(shù)據(jù)主要取自上海年鑒,統(tǒng)計(jì)項(xiàng)目包括上海總客流量、上海旅游總收入,GDP。 下面建立的模型一是利用上海市2000-2009年旅游業(yè)的統(tǒng)計(jì)數(shù)據(jù)來(lái)預(yù)測(cè)2010年的旅游業(yè)的自然增長(zhǎng)值。4.1模型一 建立旅游

6、本底趨勢(shì)線,需要對(duì)因突發(fā)事件所造成統(tǒng)計(jì)數(shù)據(jù)的異常變化運(yùn)用內(nèi)插法進(jìn)行訂正。第一步用原始統(tǒng)計(jì)數(shù)據(jù)作出長(zhǎng)期變化的實(shí)際統(tǒng)計(jì)表；第二部根據(jù)目視觀察從統(tǒng)計(jì)表上找出需進(jìn)行訂正的時(shí)間區(qū)段，并確定適合直線內(nèi)插的起始點(diǎn) 和終止點(diǎn)；第三步利用內(nèi)插公式。訂正異常統(tǒng)計(jì)數(shù)據(jù). 在該公式中, d 為內(nèi)插的公差值,其計(jì)算公式為。2000-2009年上海市旅游業(yè)在客流總量、旅游收入方面部分年份出現(xiàn)下滑,在統(tǒng)計(jì)表上形成較明顯的凹形谷,所以需對(duì)其進(jìn)行內(nèi)插訂正后才適合建本地趨勢(shì)線模行進(jìn)行預(yù)測(cè)；上海生產(chǎn)總值無(wú)此明顯波動(dòng),不作內(nèi)插訂正. 表1 上海旅游經(jīng)濟(jì)實(shí)際值年份總客流量（萬(wàn)人次）旅游總收入（億元）上海市生產(chǎn)總值（億元）2000802

7、.94908.634771.1720018458.26957.255210.1220029033.531182.65741.0320037922.481250.26694.2320049007.051472.118072.8320059573.291604.26915461731.110296.97200710865.62074.812001.16200811640.372060.3013698292407.9214900.93表2 數(shù)據(jù)序列直線內(nèi)插參數(shù)及內(nèi)插數(shù)值數(shù)據(jù)序列名稱內(nèi)插曲短段起始點(diǎn)內(nèi)插區(qū)段終至點(diǎn)公差值修正值公式2003年修正后的值2

8、004年修正后的值2008年修正后的值總客流量20012006368.279194.799562.56旅游總收入20002007166.62241.39表3內(nèi)插值年份總客流量（萬(wàn)人次）旅游總收入（億元）上海市生產(chǎn)總值（億元）2000802.94908.634771.1720018458.26957.255210.1220029033.531182.65741.0320039194.79）1250.26694.2320049562.561472.118072.8320059573.291604.26915461731.110296.97200710865.62074.8

9、12001.16200811640.372241.3913698200912218.292407.9214900本地趨勢(shì)線=長(zhǎng)期趨勢(shì)項(xiàng)+長(zhǎng)期波動(dòng)項(xiàng)+隨機(jī)波動(dòng)項(xiàng)。中國(guó)旅游業(yè)自然增長(zhǎng)是線性增長(zhǎng)或指數(shù)增長(zhǎng)，長(zhǎng)期波動(dòng)項(xiàng)選取為邏輯線，周期波動(dòng)選取正弦函數(shù)。直線模型： ；指數(shù)模型： ；邏輯線模型：；正弦線模型：；在此基礎(chǔ)上，將直線模型、指數(shù)模型與邏輯線模型、正弦線模型進(jìn)行組合，提出了六種新的數(shù)學(xué)模型以補(bǔ)充上述模型。這六種模型是：直線指數(shù)復(fù)合模型：直線邏輯線復(fù)合模型：直線正弦線復(fù)合模型：指數(shù)邏輯線復(fù)合模型： 指數(shù)正弦線復(fù)合模型：邏輯線正弦線復(fù)合模型： （1）分別用直線模型、指數(shù)模型、邏輯線模型、正弦線模型、

10、直線指數(shù)復(fù)合模型、直線邏輯線復(fù)合模型、直線正弦線復(fù)合模型、指數(shù)邏輯線復(fù)合模型、指數(shù)正弦線復(fù)合模型、邏輯線正弦線復(fù)合模型對(duì)表（3）各個(gè)指標(biāo)進(jìn)行擬和，確定有關(guān)參數(shù)，獲得各個(gè)指標(biāo)的擬和趨勢(shì)線模型方程。上海旅游總客流量的各個(gè)擬和趨勢(shì)線模型方程如下分別為：直線模型： 指數(shù)模型： 圖1邏輯線模型：正弦線模型： 圖2直線指數(shù)復(fù)合模型：直線邏輯線復(fù)合模型： 圖3直線正弦線復(fù)合模型：指數(shù)邏輯線復(fù)合模型： 圖4指數(shù)正弦線復(fù)合模型：邏輯線正弦線復(fù)合模型： 圖5上海旅游總收入的各個(gè)擬和趨勢(shì)線模型方程和擬合圖像（橫坐標(biāo)t 為時(shí)間變量，從2000 年開(kāi)始依次取t=1，2，3，縱坐標(biāo)表示對(duì)應(yīng)t的旅游總收入（億元）分別為：直

11、線模型： 指數(shù)模型： 圖6邏輯線模型：正弦線模型： 圖7直線指數(shù)復(fù)合模型：直線邏輯線復(fù)合模型：圖8直線正弦線復(fù)合模型：指數(shù)邏輯線復(fù)合模型： 圖9指數(shù)正弦線復(fù)合模型：邏輯線正弦線復(fù)合模型： 圖10.上海市生產(chǎn)總值的各個(gè)擬和趨勢(shì)線模型方程分別為直線模型： 指數(shù)模型： 圖11邏輯線模型：正弦線模型：、直線指數(shù)復(fù)合模型：直線邏輯線復(fù)合模型：直線正弦線復(fù)合模型：指數(shù)邏輯線復(fù)合模型：圖12指數(shù)正弦線復(fù)合模型：邏輯線正弦線復(fù)合模型： 圖13（2）利用得到的以上各個(gè)擬和趨勢(shì)線模型方程，計(jì)算得出2010年各個(gè)指標(biāo)本底值與各個(gè)指標(biāo)內(nèi)插值的相關(guān)系數(shù)，根據(jù)該相關(guān)系數(shù)確定最終的趨勢(shì)線模型和方程。確定步驟：選出最大的相關(guān)

12、系數(shù)模型；若按模型計(jì)算的2010年本底值，若出現(xiàn)負(fù)值，則去掉該模型，接著再選模型計(jì)算本底值與內(nèi)插值相關(guān)系數(shù)平方和大的模型，以此類推，直到得出2010年本地值為正值的相應(yīng)相關(guān)系數(shù)表4三種趨勢(shì)線模型的本底值與內(nèi)插植的相關(guān)系數(shù)的總客流量旅游總收入上海市生產(chǎn)總值直線模型0.7961 0.9913 0.9879 指數(shù)模型0.7597 0.9951 0.9977 邏輯線模型0.9456 0.9962 0.9984 正弦線模型0.2572 0.1674 0.9877 直線指數(shù)復(fù)合模型0.8563 0.9961 0.9983 直線邏輯線復(fù)合模型0.9966 0.9913 0.9878 直線正弦線復(fù)合模型0.9

13、435 0.9965 0.9993 指數(shù)邏輯線復(fù)合模型0.7957 0.996 0.998 指數(shù)正弦線復(fù)合模型0.9069 0.9969 0.9980 邏輯線正弦線復(fù)合模型0.5470 0.2756 0.2396 表5選定的本底趨勢(shì)線最終模型表總客流量旅游總收入上海生產(chǎn)總值選定的本底趨勢(shì)線最終模型直線-邏輯線復(fù)合模型指數(shù)-正弦線復(fù)合模型直線-正弦線復(fù)合模型相關(guān)系數(shù)的平方0.99660.99690.9993所以得出上海旅游總客流量的最終的本地趨勢(shì)線模型方程： .(1) 相關(guān)系數(shù)=0.9966 ，t 為時(shí)間變量，從2000 年開(kāi)始依次取t=1，2，3上海旅游總收入的最終本底趨勢(shì)線模型方程為：(2)

14、相關(guān)系數(shù)=0.9969，t 為時(shí)間變量，從2000 年開(kāi)始依次取t=1，2，3上海市國(guó)內(nèi)生產(chǎn)總值的最終本底趨勢(shì)線模型方程為： (3)相關(guān)系數(shù)=0.9993，t 為時(shí)間變量，從2000 年開(kāi)始依次取t=1，2，3將各個(gè)指標(biāo)本底趨勢(shì)線與統(tǒng)計(jì)量相比較如下： 圖14用2000-2009年的統(tǒng)計(jì)旅游總?cè)舜螖?shù)據(jù)得到的2010年的本底線、統(tǒng)計(jì)量圖15 用2000-2009年的上海統(tǒng)計(jì)的旅游總收入、年生產(chǎn)總值得到的2010年的本底線、統(tǒng)計(jì)量模型解釋 由圖一可看出2010年總客流量的影響值明顯高于往年客流量的影響值（各個(gè)指標(biāo)的實(shí)際值與本底值之差定義為上海世博對(duì)各個(gè)指標(biāo)的影響值）, 由最終擬合函數(shù)(1)式可得到2

15、010年各個(gè)指標(biāo)的本底值為12339萬(wàn)人次，總客流量的影響值為7031.61萬(wàn).在理論上,這些數(shù)值是在不舉辦世博會(huì)的情況下,上海在旅游客流量上損失的數(shù)額，由客流量的影響值得到世博會(huì)對(duì)2010年客流量的影響率（定義影響值與實(shí)際值的比率為世博會(huì)對(duì)各個(gè)指標(biāo)的影響率）為，類似分析，如圖二分別將旅游總收入和生產(chǎn)總值的本底線趨勢(shì)線與統(tǒng)計(jì)量相比較，得到上海2010年旅游總收入影響值為780.41億元，上海生產(chǎn)總值的影響值2000.89億元，同樣可以得到上海世博會(huì)對(duì)上海當(dāng)年旅游總收入和生產(chǎn)總值GDP的貢獻(xiàn)率分別為和，而以旅游業(yè)占GDP的百分比也就是說(shuō)2010年上海世博會(huì)對(duì)上海當(dāng)年旅游經(jīng)濟(jì)影響還是非常大的。下面

16、進(jìn)一步評(píng)估上海世博會(huì)對(duì)上海旅游經(jīng)濟(jì)未來(lái)3年會(huì)有怎樣的影響4.2 模型二在2000-2009年的統(tǒng)計(jì)數(shù)據(jù)基礎(chǔ)上，用(1)(2)(3)三個(gè)擬合函數(shù)分別得出在沒(méi)有舉辦世博會(huì)的情況下，上海市2010-2013年旅游業(yè)自然增長(zhǎng)的的本底值如表6。年份總客流量（萬(wàn)人次）旅游總收入（億元）上海市生產(chǎn)總值（億元）2010123392731.6169282011127883109.1176602012132373471.8187742013136863813.719683 表6 用2000-2009年的統(tǒng)計(jì)數(shù)據(jù)預(yù)測(cè)得到的2010-2013年的本底值考慮上海世博會(huì)，即在表1的基礎(chǔ)上加上2010年上?？偪土髁康慕y(tǒng)計(jì)預(yù)

17、測(cè)值（用據(jù)上海統(tǒng)計(jì)資料信息得知的上海2010年一月到四月的累計(jì)客流量的1.5倍作為2010年上海1-4月和11-12月的累計(jì)客流量，同樣將已統(tǒng)計(jì)到的2010年5-8月上海的累計(jì)客流量乘以1.5所得的結(jié)果作為2010年上海5-10月的累計(jì)客流量.最后得出就可以作為2010年上海全年總客流量）為19370.61萬(wàn)人次。用同樣的方法可以估算出2010年上海的旅游總收入為3512.01億元和年生產(chǎn)總值18928.89億元。下面對(duì)2000年和2010年的值用內(nèi)插法分別得出2001-2009年的內(nèi)插值如表7。年份客流量（萬(wàn)人次）旅游總收入（億元）上海市生產(chǎn)總值（億元）2000820.940908.6347

18、71.1720012659.7071168.9686186.94220024516.4741429.3067602.71420036373.2411689.6449018.48620048230.0081949.98210434.258200510086.7752210.3211850.03200611943.5422470.65813265.802200713800.3092730.99614681.574200815657.0762991.33416097.346200917513.8433251.67217513.118201019370.6103512.01018928.890 表7

19、由2000年和2010年的統(tǒng)計(jì)值用內(nèi)插法分別得出2001-2009年的內(nèi)插值用模型一中的擬合方法(選出最大的相關(guān)系數(shù)模型,按模型計(jì)算2011-2013年本底值，若出現(xiàn)負(fù)值，則去掉該模型，接著再選模型計(jì)算本底值與內(nèi)插值相關(guān)系數(shù)平方次大的模型，以此類推，直到得出2011-2013年本地值為正值的相應(yīng)相關(guān)系數(shù))得到相關(guān)系數(shù)為的直線模型方程：由直線模型預(yù)測(cè)得出在舉行世博的情況下2011-2013年的各項(xiàng)本底值表8。年份客流量（萬(wàn)人次）旅游總收入（億元）上海市生產(chǎn)總值（億元）201121225.5773772.34820344.662201223080.5444032.68621760.43420132

20、4935.5114293.02423176.206表8 用2000年和2010年的統(tǒng)計(jì)數(shù)據(jù)預(yù)測(cè)得到的2010-2013年的本底值對(duì)在不舉行世博的情況下得到的2011-2013年的各項(xiàng)指標(biāo)的本底值表(6)，和對(duì)舉行世博的情況下得出的2011-2013年各項(xiàng)指標(biāo)的本底值表(8)，進(jìn)行做差得出影響值表9。年份客流量（萬(wàn)人次）旅游總收入（億元）上海市生產(chǎn)總值（億元）20118438663.28162684.66220129844560.92502986.434201311250479.33413493.206合計(jì)295321703.54079164.302表9世博會(huì)各年各個(gè)指標(biāo)的影響值（表8減去表6的

21、對(duì)應(yīng)值）定義影響值與實(shí)際值的比率為世博會(huì)對(duì)各個(gè)指標(biāo)的影響率，由表9各項(xiàng)指標(biāo)的影響值與本底值表8的比率為世博會(huì)對(duì)各個(gè)指標(biāo)的影響率，得到世博會(huì)各年各項(xiàng)指標(biāo)的影響率表10年份客流量旅游總收入上海市生產(chǎn)總值201165.9821.3313.20201274.3616.1613.72201382.2012.5715.07 表10世博會(huì)對(duì)2011-2013年各個(gè)指標(biāo)的影響率 單位：%模型解釋 根據(jù)表9的計(jì)算結(jié)果：2011-2013 年，上海世博會(huì)對(duì)上海旅游總?cè)舜蔚挠绊懼禐?9532萬(wàn)人次，旅游總收入的影響值為1703.5407億元；生產(chǎn)總值的影響值為1954.302億元。在理論上，表10的數(shù)值可以表示在不

22、舉辦世博會(huì)的情況下,上海旅游總客流量、旅游總收入和上海生產(chǎn)總值上將損失的效應(yīng)，可以得出世博對(duì)2011-2013年上?？土髁康挠绊懧史謩e為、。對(duì)三年的旅游總收入影響率為、。有此數(shù)據(jù)可以看出世博會(huì)對(duì)上海這三年旅游經(jīng)濟(jì)的影響是很大。5 模型評(píng)價(jià)從上面預(yù)測(cè)值與實(shí)際值的比較可以發(fā)現(xiàn)所以，用本底趨勢(shì)線是能很好地反映世博會(huì)對(duì)旅游業(yè)的影響。因此模型得到的數(shù)據(jù)能夠很好的來(lái)預(yù)測(cè)世博會(huì)對(duì)上海市旅游經(jīng)濟(jì)的影響。但因?yàn)閺幕ヂ?lián)網(wǎng)上不容易得到相關(guān)的統(tǒng)計(jì)數(shù)據(jù)，因此模型只取了上海市過(guò)去十年的統(tǒng)計(jì)數(shù)據(jù)進(jìn)行預(yù)測(cè)，這里統(tǒng)計(jì)數(shù)據(jù)較少，可能影響模型反映世博會(huì)對(duì)上海市經(jīng)濟(jì)的影響力評(píng)估的精度。參考文獻(xiàn)：1 孫根年.我國(guó)境外旅游本底趨勢(shì)線的建

23、立及科學(xué)意義J.地理科學(xué),1998. 8 (5):442448.2 孫根年.我國(guó)6大境外客源市場(chǎng)旅游本底趨勢(shì)線的建立J.系統(tǒng)工程論與實(shí)踐,2000,20(1):140143.3孫根年,茍小東.上海市入境旅游本底趨勢(shì)線的建立及其應(yīng)用.陜西師范大學(xué)學(xué)報(bào),2006,6:版次:104諸大建,姜富明.世博會(huì)對(duì)上海的影響和對(duì)策M(jìn).上海:同濟(jì)大學(xué)出版社,2004.5孫根年.西安市境外旅游業(yè)發(fā)展的動(dòng)態(tài)趨勢(shì)及波動(dòng)周期J.西北大學(xué)學(xué)報(bào)(自然科學(xué)版),2001,31: 514 517.附錄：客流量的模型擬合：subplot(1,2,1);x=1:10;y=802.9 8458.26 9033.53 9194.79

24、9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;p=polyfit(x,y,1);x=linspace(1,10,10);z=polyval(p,x);plot(x,z,'k',x,y,'O');plegend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('直線模型')subplot(1,2,2);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56

25、9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('beta(1)*exp(beta(2)*x)','beta','x');beta=nlinfit(x,y,fun,50 1)g=beta(1)*exp(beta(2)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('指數(shù)模型')f

26、igure;subplot(1,2,1);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x)','beta','x');beta=nlinfit(x,y,fun,50 10 10)g=beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x);plot(x,y,'O&#

27、39;,x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('邏輯線模型')subplot(1,2,2);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('beta(1)*sin(beta(2)*x+beta(3)','beta','x')

28、;beta=nlinfit(x,y,fun,900 100 100)g=beta(1)*sin(beta(2)*x+beta(3);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('正弦線模型')figure;subplot(1,2,1);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.

29、37 12218.29;fun=inline('beta(1)+beta(2)*x+beta(3)*exp(beta(4)*x)','beta','x');beta=nlinfit(x,y,fun,1 10 10 1)g=beta(1)+beta(2)*x+beta(3)*exp(beta(4)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('

30、直線指數(shù)復(fù)合模型')subplot(1,2,2);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('beta(1)+beta(2)*x+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x)','beta','x');beta=nlinfit(x,y,fun,10 10 20 30 0)g=beta(1)+beta(2)*x+beta(3)*ones(1

31、,10)./(1+exp(beta(4)-beta(5)*x);plot(x,y,'O',x,g,'K')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('直線邏輯線復(fù)合模型')figure;subplot(1,2,1);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('be

32、ta(1)+beta(2)*x+beta(3)*sin(beta(4)*x+beta(5)','beta','x');beta=nlinfit(x,y,fun,500 100 900 100 100)g=beta(1)+beta(2)*x+beta(3)*sin(beta(4)*x+beta(5);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('直線正弦線復(fù)

33、合模型')subplot(1,2,2);y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('beta(1)*exp(beta(2)*x)+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x)','beta','x');beta=nlinfit(x,y,fun,500 1 10 100 100)g=beta(1)*exp(beta(2)*x)+beta(3)*ones(

34、1,10)./(1+exp(beta(4)-beta(5)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('指數(shù)邏輯線復(fù)合模型')figure;subplot(1,2,1);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('b

35、eta(1)*exp(beta(2)*x)+beta(3)*sin(beta(4)*x+beta(5)','beta','x');beta=nlinfit(x,y,fun,50 1 100 100 500)g=beta(1)*exp(beta(2)*x)+beta(3)*sin(beta(4)*x+beta(5);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title(

36、9;指數(shù)正弦線復(fù)合模型')subplot(1,2,2);x=1:10;y=802.9 8458.26 9033.53 9194.79 9562.56 9573.29 10299.6 10865.6 11640.37 12218.29;fun=inline('beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x)+beta(4)*sin(beta(5)*x+beta(6)','beta','x');beta=nlinfit(x,y,fun,10 10 20 30 10 20)g=beta(1)*ones(

37、1,10)./(1+exp(beta(2)-beta(3)*x)+beta(4)*sin(beta(5)*x+beta(6);plot(x,y,'O',x,g,'K')legend('統(tǒng)計(jì)量','本底線')ylabel('萬(wàn)人次')xlabel('年份')title('邏輯線正弦線復(fù)合模型')旅游總收入的模型擬合：subplot(1,2,1);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241

38、.39 2407.92;p=polyfit(x,y,1);x=linspace(1,10,10);z=polyval(p,x);plot(x,z,'k',x,y,'O');plegend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線模型')subplot(1,2,2);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun

39、=inline('beta(1)*exp(beta(2)*x)','beta','x');beta=nlinfit(x,y,fun,50 1)g=beta(1)*exp(beta(2)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('指數(shù)模型')figure;subplot(1,2,1);x=1:10;y=908.63 957.25 11

40、82.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun=inline('beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x)','beta','x');beta=nlinfit(x,y,fun,50 10 10)g=beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線

41、9;)ylabel('億元')xlabel('年份')title('邏輯線模型')subplot(1,2,2);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun=inline('beta(1)*sin(beta(2)*x+beta(3)','beta','x');beta=nlinfit(x,y,fun,900 100 100)g=beta(1)*sin(beta(2)*x+b

42、eta(3);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('正弦線模型')figure;subplot(1,2,1);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun=inline('beta(1)+beta(2)*x+beta(3)*exp(beta(4)*

43、x)','beta','x');beta=nlinfit(x,y,fun,1 10 10 1)g=beta(1)+beta(2)*x+beta(3)*exp(beta(4)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線指數(shù)復(fù)合模型')subplot(1,2,2);x=1:10;y=908.63 957.25 1182.6 1250.2 1

44、472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun=inline('beta(1)+beta(2)*x+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x)','beta','x');beta=nlinfit(x,y,fun,10 10 20 30 0)g=beta(1)+beta(2)*x+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x);plot(x,y,'O',x,g,'K')lege

45、nd('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線邏輯線復(fù)合模型')figure;subplot(1,2,1);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun=inline('beta(1)+beta(2)*x+beta(3)*sin(beta(4)*x+beta(5)','beta','x');b

46、eta=nlinfit(x,y,fun,500 100 900 100 100)g=beta(1)+beta(2)*x+beta(3)*sin(beta(4)*x+beta(5);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線正弦線復(fù)合模型')subplot(1,2,2);y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074

47、.8 2241.39 2407.92;fun=inline('beta(1)*exp(beta(2)*x)+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x)','beta','x');beta=nlinfit(x,y,fun,500 1 10 100 100)g=beta(1)*exp(beta(2)*x)+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量&#

48、39;,'本底線')ylabel('億元')xlabel('年份')title('指數(shù)邏輯線復(fù)合模型')figure;subplot(1,2,1);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 2074.8 2241.39 2407.92;fun=inline('beta(1)*exp(beta(2)*x)+beta(3)*sin(beta(4)*x+beta(5)','beta','x');beta=nlin

49、fit(x,y,fun,50 1 100 100 500)g=beta(1)*exp(beta(2)*x)+beta(3)*sin(beta(4)*x+beta(5);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('指數(shù)正弦線復(fù)合模型')subplot(1,2,2);x=1:10;y=908.63 957.25 1182.6 1250.2 1472.11 1604.26 1731.1 207

50、4.8 2241.39 2407.92;fun=inline('beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x)+beta(4)*sin(beta(5)*x+beta(6)','beta','x');beta=nlinfit(x,y,fun,10 10 20 30 10 20)g=beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x)+beta(4)*sin(beta(5)*x+beta(6);plot(x,y,'O',x,g,'K')l

51、egend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('邏輯線正弦線復(fù)合模型')生產(chǎn)總值的模型擬合：subplot(1,2,1);x=1:10;y=4771.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;p=polyfit(x,y,1);x=linspace(1,10,10);z=polyval(p,x);plot(x,z,'k',x,y,'O'

52、;);plegend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線模型')subplot(1,2,2);x=1:10;y=4771.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;fun=inline('beta(1)*exp(beta(2)*x)','beta','x');beta=nlinfit(x,y,fun,50 1)

53、g=beta(1)*exp(beta(2)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('指數(shù)模型')figure;subplot(1,2,1);x=1:10;y=4771.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;fun=inline('beta(1)*ones(1,10

54、)./(1+exp(beta(2)-beta(3)*x)','beta','x');beta=nlinfit(x,y,fun,50 10 10)g=beta(1)*ones(1,10)./(1+exp(beta(2)-beta(3)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('邏輯線模型')subplot(1,2,2);x=1:10;y=47

55、71.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;fun=inline('beta(1)*sin(beta(2)*x+beta(3)','beta','x');beta=nlinfit(x,y,fun,900 100 100)g=beta(1)*sin(beta(2)*x+beta(3);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel(&

56、#39;億元')xlabel('年份')title('正弦線模型')figure;subplot(1,2,1);x=1:10;y=4771.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;fun=inline('beta(1)+beta(2)*x+beta(3)*exp(beta(4)*x)','beta','x');beta=nlinfit(x,y,fun,1 10 10 1)g=beta(1)+beta(2

57、)*x+beta(3)*exp(beta(4)*x);plot(x,y,'O',x,g,'k')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線指數(shù)復(fù)合模型')subplot(1,2,2);x=1:10;y=4771.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;fun=inline('beta(1)+beta(2)*x+

58、beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x)','beta','x');beta=nlinfit(x,y,fun,10 10 20 30 0)g=beta(1)+beta(2)*x+beta(3)*ones(1,10)./(1+exp(beta(4)-beta(5)*x);plot(x,y,'O',x,g,'K')legend('統(tǒng)計(jì)量','本底線')ylabel('億元')xlabel('年份')title('直線邏輯線復(fù)合模型')figure;subplot(1,2,1);x=1:10;y=4771.17 5210.12 5741.03 6694.23 8072.83 9154.18 10296.97 12001.16 13698 14900;fun=inline('beta(1)+beta(2)*x+