資源描述：
《20115394王福臨數(shù)學(xué)實(shí)驗作業(yè)(函數(shù)迭代)》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、重慶大學(xué)學(xué)生實(shí)驗報告實(shí)驗課程名稱數(shù)學(xué)實(shí)驗開課實(shí)驗室DS1401學(xué)院計算機(jī)學(xué)院年級大二專業(yè)計科六班學(xué)生姓名王福臨學(xué)號20115394開課時間至學(xué)年第學(xué)期總成績教師簽名數(shù)理學(xué)院制開課學(xué)院、實(shí)驗室：DS1401實(shí)驗時間：2013年4月13日課程名稱數(shù)學(xué)實(shí)驗實(shí)驗項目名稱函數(shù)迭代實(shí)驗項目類型驗證演示綜合設(shè)計其他指導(dǎo)教師何光輝成績實(shí)驗?zāi)康摹ぃ?）理解迭代的基本概念；·（2）掌握迭代數(shù)列的系列圖形表示方法；·（3）以一類特殊二次函數(shù)（Logistic函數(shù)）為例，掌握二次函數(shù)迭代序列的收斂性分析方法；·（4）熟悉編寫函數(shù)迭代的MATLAB
2、程序；·（5）了解二元函數(shù)迭代方法及其圖形特征基礎(chǔ)實(shí)驗一、實(shí)驗內(nèi)容1、迭代以下函數(shù)，分析其收斂性。任選一個完成。使用線性聯(lián)結(jié)圖、蛛網(wǎng)圖或費(fèi)根鮑姆圖對參數(shù)a進(jìn)行討論與觀察，會得到什么結(jié)論？2.生物種群的數(shù)量問題種群的數(shù)量（為方便起見以下指雌性）因繁殖而增加，因自然死亡和人工捕獲而減少。記xk(t)為第t年初k歲（指滿k-1歲，未滿k歲，下同）的種群數(shù)量，bk為k歲種群的繁殖率（1年內(nèi)每個個體繁殖的數(shù)量），dk為k歲種群的死亡率（1年內(nèi)死亡數(shù)量占總量的比例），hk為k歲種群的捕獲量（1年內(nèi)的捕獲量）。今設(shè)某種群最高年齡為5歲（
3、不妨認(rèn)為在年初將5歲個體全部捕獲），b1=b2=b5=0，b3=2，b4=4，d1=d2=0.3，d3=d4=0.2，h1=400，h2=200，h3=150，h4=100。A.??建立xk(t+1)與xk(t)的關(guān)系（k=1,2,...5,t=0,1,...），.如為簡單起見，繁殖量都按年初的種群數(shù)量xk(t)計算，不考慮死亡率。B.??用向量表示t年初的種群數(shù)量，用bk和dk定義適當(dāng)?shù)木仃嘗，用hk定義適當(dāng)?shù)南蛄縣，將上述關(guān)系表成的形式。C.?設(shè)t=0種群各年齡的數(shù)量均為1000，求t=1種群各年齡的數(shù)量。又問設(shè)定的捕
4、獲量能持續(xù)幾年。D.?種群各年齡的數(shù)量等于多少，種群數(shù)量x(t)才能不隨時間t改變。記D的結(jié)果為向量x*,給x*以小的擾動作為x(0)，觀察隨著t的增加x(t)是否趨于x*,分析這個現(xiàn)象的原因。二、實(shí)驗過程（一般應(yīng)包括實(shí)驗原理或問題分析，算法設(shè)計、程序、計算、圖表等，實(shí)驗結(jié)果及分析）1．開啟軟件平臺——MATLAB，開啟MATLAB編輯窗口；2．根據(jù)各種數(shù)值解法步驟編寫M文件3．保存文件并運(yùn)行；4．觀察運(yùn)行結(jié)果(數(shù)值或圖形)；5．根據(jù)觀察到的結(jié)果寫出實(shí)驗報告，并淺談學(xué)習(xí)心得體會。應(yīng)用實(shí)驗（或綜合實(shí)驗）1、迭代以下函數(shù)，分析
5、其收斂性。任選一個完成。使用線性聯(lián)結(jié)圖、蛛網(wǎng)圖或費(fèi)根鮑姆圖對參數(shù)a進(jìn)行討論與觀察，會得到什么結(jié)論？解：選擇第3)題完成（1）線性聯(lián)結(jié)圖解法:>>a=[0.1,0.3,0.5,0.8];>>x1=[];x2=[];x3=[];x4=[];x1(1)=0.5;x2(1)=0.5;x3(1)=0.5;x4(1)=0.5;>>fori=2:20x1(i)=x1(i-1)^4-a(1);x2(i)=x2(i-1)^4-a(2);x3(i)=x3(i-1)^4-a(3);x4(i)=x4(i-1)^4-a(4);end>>n=1:20
6、;subplot(2,2,1),plot(n,x1),title('a=0.1,x0=0.5')subplot(2,2,2),plot(n,x2),title('a=0.3,x0=0.5')subplot(2,2,3),plot(n,x3),title('a=0.5,x0=0.5')subplot(2,2,4),plot(n,x4),title('a=0.8,x0=0.5')由上圖可知：a=0.1,x0=0.5，a=0.3,x0=0.5和a=0.5，x0=0.5時，函數(shù)迭代收斂，a=0.8,x0=0.5時函數(shù)迭代不收斂。(
7、2)蛛網(wǎng)模型解法：a=0.5;ezplot('x^4-0.5',[-0.6,-0.2]),holdonezplot('x',[-0.6,-0.2]);x1=[];x1(1)=0.2;fori=2:50x1(i)=x1(i-1)^4-a;plot([x1(i-1),x1(i-1)],[x1(i-1),x1(i)]);plot([x1(i-1),x1(i)],[x1(i),x1(i)]);end由上圖可知，函數(shù)a=0.5,x0=0.2時，函數(shù)迭代收斂(3)費(fèi)根鮑姆圖解法：iter.m：functionroot=iter(x0,
8、a)x(1)=x0;fori=2:100x(i)=a*sin(pi*x(i-1));endroot=x;died.mclf;x=[];x0=0.2;holdon;fora=0:0.01:1root=iter(x0,a);plot(a.*ones(size(root(51:100))),root(51:1
