差分方程matlab解法ppt課件_第1頁
差分方程matlab解法ppt課件_第2頁
差分方程matlab解法ppt課件_第3頁
差分方程matlab解法ppt課件_第4頁
差分方程matlab解法ppt課件_第5頁
已閱讀5頁,還剩5頁未讀 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

1、用用MatlabMatlab求解差分方程問題求解差分方程問題一階線性常系數(shù)差分方程高階線性常系數(shù)差分方程線性常系數(shù)差分方程組一階線性常系數(shù)差分方程 瀕危物種的自然演化和人工孵化 問題 Florida沙丘鶴屬于瀕危物種,它在較好自然環(huán)境下,年均增長率僅為1.94%,而在中等和較差環(huán)境下年均增長率分別為 -3.24% 和 -3.82%,假設(shè)在某自然維護區(qū)內(nèi)開場有100只鶴,建立描畫其數(shù)量變化規(guī)律的模型,并作 數(shù)值計算。模型建立 記第k年沙丘鶴的數(shù)量為xk,年均增長率為r,那么第k+1年鶴的數(shù)量為 xk+1=(1+r)xk k=0,1,2 知x0=100, 在較好,中等和較差的自然環(huán)境下 r=0.0

2、194, -0.0324,和-0.0382 我們利用Matlab編程,遞推20年后察看沙丘鶴的數(shù)量變化情況Matlab實現(xiàn) 首先建立一個關(guān)于變量n ,r的函數(shù) function x=sqh(n,r) a=1+r; x=100; for k=1:n x(k+1)=a*x(k); end 在command窗口里調(diào)用sqh函數(shù) k=(0:20); y1=sqh(20,0.0194); y2=sqh(20,-0.0324); y3=sqh(20,-0.0382); round(k,y1,y2,y3)利用plot 繪圖察看數(shù)量變化趨勢 可以用不同線型和顏色繪圖 r g b c m y k w 分別表示

3、紅綠蘭蘭綠洋紅黃黑白色 : + o * . X s d 表示不同的線型 plot(k,y1,k,y2,k,y3) 在同一坐標(biāo)系下畫圖 plot(k,y2,:) plot(k,y2,-) plot(k,y2,r) plot(k,y2,y) plot(k,y2,y,k,y1,:) plot(k,y2,k,y1,:) plot(k,y2,oy,k,y1,:)用gtext(r=0.0194),gtext(r=-0.0324),gtext(r=-0.0382)在圖上做標(biāo)志。 function x=fhsqh(n,r,b) a=1+r; X=100; For k=1:n X(k+1)=a*x(k)+b; end k=(0:20) ; %一個行向量 y1=(20,-0.0324,5); 也是一個行向量 round( k , y 1 ) 對k,y1四舍五入,但 是 不改動變量的值 plot( k , y1) k y1 是行向量列向量都可以 也可以察看200年的開展趨勢,以及在較差條件下的開展趨勢,也可以調(diào)查每年孵化數(shù)量變化的影響。一階線性常系數(shù)差分方程的解、平衡點及其穩(wěn)定性一階線性常系數(shù)差分方程的解、平衡點及其穩(wěn)定性 自然環(huán)境下,b=0 人工孵化條件下 令xk=xk+1=x得 差分方程的平衡

溫馨提示

  • 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)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論