差分方程建?;騙第1頁
差分方程建?;騙第2頁
差分方程建?;騙第3頁
差分方程建模或_第4頁
差分方程建?;騙第5頁
已閱讀5頁,還剩20頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)

文檔簡介

差分方程建?;虻?頁,課件共25頁,創(chuàng)作于2023年2月滿足一差分方程的序列yt稱為此差分方程的解。類似于微分方程情況,若解中含有的獨(dú)立常數(shù)的個數(shù)等于差分方程的階數(shù)時,稱此解為該差分方程的通解。若解中不含任意常數(shù),則稱此解為滿足某些初值條件的特解,例如,考察兩階差分方程易見與均是它的特解,而則為它的通解,其中c1,c2為兩個任意常數(shù)。類似于微分方程,稱差分方程

為n階線性差分方程,當(dāng)≠0時稱其為n階非齊次線性差分方程,而第2頁,課件共25頁,創(chuàng)作于2023年2月則被稱為方程對應(yīng)的齊次線性差分方程。若所有的ai(t)均為與t無關(guān)的常數(shù),則稱其為常系數(shù)差分方程,即n階常系數(shù)線性差分方程可分成(4.15)

的形式,其對應(yīng)的齊次方程為(4.16)

容易證明,若序列與均為方程(4.16)的解,則也是方程(4.16)的解,其中c1、c2為任意常數(shù),這說明,齊次方程的解構(gòu)成一個線性空間(解空間)。

此規(guī)律對于(4.15)也成立。第3頁,課件共25頁,創(chuàng)作于2023年2月方程(4.15)可用如下的代數(shù)方法求其通解:(步一)先求解對應(yīng)的特征方程

(4.17)

(步二)根據(jù)特征根的不同情況,求齊次方程(4.16)的通解

情況1

若特征方程(4.17)有n個互不相同的實(shí)根,…,,則齊次方程(4.16)的通解為

(C1,…,Cn為任意常數(shù)),情況2

若λ

是特征方程(4.17)的k重根,通解中對應(yīng)于λ的項(xiàng)為為任意常數(shù),i=1,…,k。情況3

若特征方程(4.17)有單重復(fù)根通解中對應(yīng)它們的項(xiàng)為為λ的模,為λ的幅角。

第4頁,課件共25頁,創(chuàng)作于2023年2月情況4

若為特征方程(4.17)的k重復(fù)根,則通解對應(yīng)于它們的項(xiàng)為為任意常數(shù),i=1,…,2k。

.若yt為方程(4.16)的通解,則非齊次方程(4.15)的通解為(步三)求非齊次方程(4.15)的一個特解

求非齊次方程(4.15)的特解一般要用到常數(shù)變易法,計(jì)算較繁。對特殊形式的b(t)也可使用待定系數(shù)法。第5頁,課件共25頁,創(chuàng)作于2023年2月對特殊形式的b(t)也可使用待定系數(shù)法

如果為t的多項(xiàng)式,形式的特解,其中為m次多項(xiàng)式;,將其代入(4.15)中(1)當(dāng)t不是特征根時,可設(shè)成形如(2)如果b是r重根時,可設(shè)特解:確定出系數(shù)即可。第6頁,課件共25頁,創(chuàng)作于2023年2月例4.13

求解兩階差分方程解對應(yīng)齊次方程的特征方程為,其特征根為,對應(yīng)齊次方程的通解為

原方程有形如的特解。代入原方程求得,,故原方程的通解為在應(yīng)用差分方程研究問題時,一般不需要求出方程的通解,在給定初值后,通??捎糜?jì)算機(jī)迭代求解,但我們常常需要討論解的穩(wěn)定性。對差分方程(4.15),若不論其對應(yīng)齊次方程的通解中任意常數(shù)C1,…,Cn如何取值,在時總有,則稱方程(7.14)的解是穩(wěn)定的,否則稱其解為不穩(wěn)定的.根據(jù)通解的結(jié)構(gòu)不難看出,非齊次方程(4.15)穩(wěn)定的充要條件為其所有特征根的模均小于1。第7頁,課件共25頁,創(chuàng)作于2023年2月例4.14(市場經(jīng)濟(jì)的蛛網(wǎng)模型)在自由競爭的市場經(jīng)濟(jì)中,商品的價格是由市場上該商品的供應(yīng)量決定的,供應(yīng)量越大,價格就越低。另一方面,生產(chǎn)者提供的商品數(shù)量又是由該商品的價格決定的,價格上升將刺激生產(chǎn)者的生產(chǎn)積極性,導(dǎo)致商品生產(chǎn)量的增加。反之,價格降低會影響生產(chǎn)者的積極性,導(dǎo)致商品生產(chǎn)量的下降。在市場經(jīng)濟(jì)中,對每一商品事實(shí)上存在著兩個不同的函數(shù):(1)供應(yīng)函數(shù)x=f(P),它是價格P的單增函數(shù),其曲線稱為供應(yīng)曲線。(2)需求函數(shù)x=g(P),它是價格P的單降函數(shù),其曲線稱為需求曲線,供應(yīng)曲線與需求曲線的形狀如圖所示。第8頁,課件共25頁,創(chuàng)作于2023年2月記t時段初市場上的供應(yīng)量(即上一時段的生產(chǎn)量)為xt,市場上該商品的價格為Pt。商品成交的價格是由需求曲線決定的,即隨著

,Mt將趨于平衡點(diǎn)M*,即商品量將趨于平衡量x*,價格將趨于平衡價格P*。圖中的箭線反映了在市場經(jīng)濟(jì)下該商品的供應(yīng)量與價格的發(fā)展趨勢。xoPP0P2P*P1xx1x2x0x*需求曲線供應(yīng)曲線M0M2M1M*①PoM3M2M1②但是,如果供應(yīng)曲線和需求曲線呈圖①中的形狀,則平衡點(diǎn)M*是不穩(wěn)定的,Mt將越來越遠(yuǎn)離平衡點(diǎn)。第9頁,課件共25頁,創(chuàng)作于2023年2月圖①和圖②的區(qū)別在哪里,如何判定平衡點(diǎn)的穩(wěn)定性呢?但是,如果供應(yīng)曲線和需求曲線呈圖②中的形狀,則平衡點(diǎn)M*是不穩(wěn)定的,Mt將越來越遠(yuǎn)離平衡點(diǎn)。即使初始時刻的供應(yīng)量和價格對應(yīng)于平衡點(diǎn),一點(diǎn)微小的波動也會導(dǎo)致市場供求出現(xiàn)越來越大的混亂。上述用圖示法分析市場經(jīng)濟(jì)穩(wěn)定性的討論在經(jīng)濟(jì)學(xué)中被稱為市場經(jīng)濟(jì)的蛛網(wǎng)模型。不難看出,在圖①中平衡點(diǎn)M*處供應(yīng)曲線的切線斜率大于需求曲線切線斜率的絕對值,而在圖②中情況恰好相反。第10頁,課件共25頁,創(chuàng)作于2023年2月現(xiàn)在利用差分方程方法來研究蛛網(wǎng)模型,以驗(yàn)證上述猜測是否正確。我們知道,平衡點(diǎn)M*是否穩(wěn)定取決于在M*附近供、需曲線的局部性態(tài)。為此,用M*處供、需曲線的線性近似來代替它們,并討論此線性近似模型中M*的穩(wěn)定性。設(shè)供應(yīng)曲線與需求曲線的線性近似分別為

和式中,a、b分別為供應(yīng)曲線在M*處的切線斜率與需求曲線在M*處切線斜率的絕對值。

根據(jù)市場經(jīng)濟(jì)的規(guī)律,當(dāng)供應(yīng)量為xt時,現(xiàn)時段的價格,又對價格,由供應(yīng)曲線解得下一時段的商品量

第11頁,課件共25頁,創(chuàng)作于2023年2月由此導(dǎo)出一階差分方程:(4.18)此差分方程的解在(b/a)<1時是穩(wěn)定的,從而證實(shí)了我們的猜測。注意到a和b的實(shí)際含義,上述結(jié)果在經(jīng)濟(jì)學(xué)上可作如下解釋:當(dāng)a>b時,顧客需求對價格的敏感度較?。ㄐ∮谏a(chǎn)者的敏感程度),商品供應(yīng)量和價格會自行調(diào)節(jié)而逐步趨于穩(wěn)定;反之,若a<b(商品緊缺易引起顧客搶購),該商品供售市場易造成混亂.如果生產(chǎn)者對市場經(jīng)濟(jì)的蛛網(wǎng)模型有所了解,為了減少因價格波動而造成的經(jīng)濟(jì)損失,他應(yīng)當(dāng)提高自己的經(jīng)營水平,不應(yīng)當(dāng)僅根據(jù)上一周期的價格來決定現(xiàn)階段的生產(chǎn)量。例如可以根據(jù)本時段與前一時段價格的平均值來確定生產(chǎn)量。此時,若t時段的商品量為xt

時,仍有第12頁,課件共25頁,創(chuàng)作于2023年2月(4.21)將(4.19)式、(4.21)式代入(4.20)式,整理得(4.19)但t+1時段的商品量則不再為而被修正為(4.20)由(4.19)式得(4.22)(4.22)式是一個常系數(shù)二階線性差分方程,特征方程為其特征根為第13頁,課件共25頁,創(chuàng)作于2023年2月記。若,則此時差分方程(4.22)是不穩(wěn)定的。

,若,此時特征根為一對共軛復(fù)數(shù),。

由線性差分方程穩(wěn)定的條件,當(dāng)r<2即b<2a時(4.22)式是穩(wěn)定的,從而M*是穩(wěn)定的平衡點(diǎn)。不難發(fā)現(xiàn),生產(chǎn)者管理方式的這一更動不僅使自己減少了因價格波動而帶來的損失,而且大大消除了市場的不穩(wěn)定性。生產(chǎn)者在采取上述方式來確定各時段的生產(chǎn)量后,如發(fā)現(xiàn)市場仍不穩(wěn)定(b≥2a),可按類似方法試圖再改變確定生產(chǎn)量的方式,此時可得到更高階的差分方程。對這些方程穩(wěn)定性條件的研究很可能會導(dǎo)出進(jìn)一步穩(wěn)定市場經(jīng)濟(jì)的新措施。第14頁,課件共25頁,創(chuàng)作于2023年2月例4.15

國民經(jīng)濟(jì)的穩(wěn)定性國民收入的主要來源是生產(chǎn),國民收入的開支主要用于消費(fèi)資金、投入再生產(chǎn)的積累資金及政府用于公共設(shè)施的開支。現(xiàn)在我們用差分方程方法建立一個簡略的模型,粗略地分析一下國民經(jīng)濟(jì)的穩(wěn)定性問題。再生產(chǎn)的投資水平It取決于消費(fèi)水平的變化量,設(shè)政府用于公共設(shè)施的開支在一個不太大的時期內(nèi)變動不大,設(shè)為常數(shù)G。故由可得出。將及代入。

記yt為第t周期的國民收入,Ct為第t周期的消費(fèi)資金。Ct的值決定于前一周期的國民收入,設(shè)第15頁,課件共25頁,創(chuàng)作于2023年2月(4.23)(4.23)式是一個二階常系數(shù)差分方程,其特征方程為,相應(yīng)特征根為(4.24)成立時才是穩(wěn)定的。(4.24)式可用于預(yù)報經(jīng)濟(jì)發(fā)展趨勢。現(xiàn)用待定系數(shù)法求方程(4.23)的一個特解令代入(4.23)式,得故當(dāng)(4.24)式成立時,差分方程(4.23)的通解為其中ρ為的模,ω為其幅角。第16頁,課件共25頁,創(chuàng)作于2023年2月例如,若取,易見,此時關(guān)系式(4.12)成立,又若取y0=1600,y1=1700,

G=550,則由迭代公式求得y2=1862.5,y3=2007.8,y4=2110.3,y5=2171.2,y6=2201.2,y7=2212.15,y8=2213.22,y9=2210.3,…。易見第17頁,課件共25頁,創(chuàng)作于2023年2月例4.16

商品銷售量預(yù)測

(實(shí)例)某商品前5年的銷售量見表。現(xiàn)希望根據(jù)前5年的統(tǒng)計(jì)數(shù)據(jù)預(yù)測第6年起該商品在各季度中的銷售量。從表中可以看出,該商品在前5年相同季節(jié)里的銷售量呈增長趨勢,而在同一年中銷售量先增后減,第一季度的銷售量最小而第三季度的銷售量最大。預(yù)測該商品以后的銷售情況,一種辦法是應(yīng)用最小二乘法建立經(jīng)驗(yàn)?zāi)P?。即根?jù)本例中數(shù)據(jù)的特征,可以按季度建立四個經(jīng)驗(yàn)公式,分別用來預(yù)測以后各年同一季度的銷售量。例如,如認(rèn)為第一季度的銷售量大體按線性增長,可設(shè)銷售量由15253217152430151320271512182614111625121234第五年第四年第三年第二年第一年銷售量季度

年份第18頁,課件共25頁,創(chuàng)作于2023年2月求得a=1.3,b=9.5。根據(jù)預(yù)測第六年起第一季度的銷售量為

=17.3,=18.6,…如認(rèn)為銷售量并非逐年等量增長而是按前一年或前幾年同期銷售量的一定比例增長的,則可建立相應(yīng)的差分方程模型。仍以第一季度為例,為簡便起見不再引入上標(biāo),以表示第t年第一節(jié)季度的銷售量,建立形式如下的差分方程:或等等。上述差分方程中的系數(shù)不一定能使所有統(tǒng)計(jì)數(shù)據(jù)吻合,較為合理的辦法是用最小二乘法求一組總體吻合較好的數(shù)據(jù)。以建立二階差分方程為例,為選取a0,a1,a2使最小,解線性方程組:第19頁,課件共25頁,創(chuàng)作于2023年2月即求解得a0=-8,a1=-1,a2=3。即所求二階差分方程為

第20頁,課件共25頁,創(chuàng)作于2023年2月雖然這一差分方程恰好使所有統(tǒng)計(jì)數(shù)據(jù)吻合,但這只是一個巧合。根據(jù)這一方程,可迭代求出以后各年第一季度銷售量的預(yù)測值y6=21,y7=19,…等。上述為預(yù)測各年第一季度銷售量而建立的二階差分方程,雖然其系數(shù)與前5年第一季度的統(tǒng)計(jì)數(shù)據(jù)完全吻合,但用于預(yù)測時預(yù)測值與事實(shí)不符。憑直覺,第六年估計(jì)值明顯偏高,第七年銷售量預(yù)測值甚至小于第六年。稍作分析,不難看出,如分別對每一季度建立一差分方程,則根據(jù)統(tǒng)計(jì)數(shù)據(jù)擬合出的系數(shù)可能會相差甚大,但對同一種商品,這種差異應(yīng)當(dāng)是微小的,故應(yīng)根據(jù)統(tǒng)計(jì)數(shù)據(jù)建立一個共用于各個季度的差分方程。為此,將季度編號為t=1,2,…,20,令或等,利用全體數(shù)據(jù)來擬合,求擬合得最好的系數(shù)。以二階差分方程為例,為求a0、a1、a2使得

最小第21頁,課件共25頁,創(chuàng)作于2023年2月求解線性方程組即求解三元一次方程組解得a0=0.6937,a1=0.8737,a2=0.1941,故求得二階差分方程(t≥21)

根據(jù)此式迭代,可求得第六年和第七年第一季度銷售量的預(yù)測值為y21=17.58,y25=19.16還是較為可信的。第22頁,課件共25頁,創(chuàng)作于2023年2月例4.16

人口問題的差分方程模型在§3.2中,我們已經(jīng)討論了人口問題的兩個常微分方程模型——Malthus模型和Verhulst模型(又稱Logistic模型)。前者可用于人口增長的短期預(yù)測,后者在作中、長期預(yù)測時效果較好。1、離散時間的Logistic模型在研究人口或種群數(shù)量的實(shí)際增長情況時,有時采用離散化的時間變量更為方便。例如,有些種群具有相對較為固定的繁殖期,按時段統(tǒng)計(jì)種群數(shù)量更接近種群的實(shí)際增長方式。人口增長雖無這種特征,但人口普查不可能連續(xù)統(tǒng)計(jì),任何方式的普查都只能得到一些離散時刻的人口總量(指較大范圍的普查)。這樣,如何

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論