數(shù)學(xué)模型報告新作_第1頁
數(shù)學(xué)模型報告新作_第2頁
數(shù)學(xué)模型報告新作_第3頁
數(shù)學(xué)模型報告新作_第4頁
數(shù)學(xué)模型報告新作_第5頁
已閱讀5頁,還剩4頁未讀 繼續(xù)免費閱讀

下載本文檔

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

文檔簡介

9 實 驗 報 告課 程 數(shù)學(xué)模型實驗 系 別 數(shù)統(tǒng)學(xué)院 班 級 2012級 學(xué) 號 2012325433 姓 名 王靜怡 實驗報告一:Matlab解決預(yù)測問題1實驗課程名稱 數(shù)學(xué)實驗2實驗項目名稱 Matlab解決預(yù)測問題3 實驗?zāi)康暮鸵?了解用Matlab軟件求解圖論模型及層次分析模型的方法。4 實驗內(nèi)容和原理內(nèi)容:某公司想用全行業(yè)的銷售額作為自變量來預(yù)測公司的銷售額,下表給出了19771981年的公司銷售額和行業(yè)銷售額的分季度數(shù)據(jù)(單位:百萬元).(1) 畫出數(shù)據(jù)的散點圖,觀察用線性回歸模型擬合是否合適.(2) 建立公司銷售額對全行業(yè)銷售額的回歸模型,并用DW檢驗診斷隨機誤差項的自相關(guān)性.(3) 建立消除了隨機誤差項自相關(guān)性后的回歸模型.年季t公司銷售額y行業(yè)銷售額x年季t公司銷售額y行業(yè)銷售額x19771120.96127.3197931124.54148.32221.40130.041224.30146.43321.96132.7198011325.00150.24421.52129.421425.64153.119781522.39135.031526.36157.32622.76137.141626.98160.73723.48141.2198111727.52164.24823.66142.821827.78165.619791924.10145.531928.24168.721024.01145.342028.78171.75主要儀器設(shè)備 計算機與Windows 2000/XP系統(tǒng);Matlab等軟件。6操作方法與實驗步驟步驟:1.(1)打開Matlab,新建file-M文件。 (2)在M文件編輯窗口輸入以下程序。命令如下:y=20.96;21.40;21.96;21.52;22.39;22.76;23.48;23.66;24.10;24.01;24.54;24.30;25.00;25.64;26.36;26.98;27.52;27.78;28.24;28.78;x=127.3;130.0;132.7;129.4;135.0;137.1;141.2;142.8;145.5;145.3;148.3;146.4;150.2;153.1;157.3;160.7;164.2;165.6;168.7;171.7;n=length(y);plot(x,y,o);b,bint,r,rint,stats=regress(y,ones(n,1) x);b bintstatspauserstool(x,y);pausee=y-b(1)-b(2)*x;et=zeros(n-1,1);et_1=zeros(n-1,1);for i=1:n-1 et(i)=e(i+1); et_1(i)=e(i);endplot(et_1,et,o)pauseDW=0;sum2=0;for i=2:n DW=DW+(e(i)-e(i-1)2; sum2=sum2+e(i)*e(i);endDW=DW/sum2pauserho=1-DW/2;yh=zeros(n-1,1);xh=zeros(n-1,1);for i=2:n yh(i-1)=y(i)-rho*y(i-1); xh(i-1)=x(i)-rho*x(i-1);endbh,bhint,rh,rhint,statsh=regress(yh,ones(n-1,1) xh);bh,bhintstatshpauserstool(xh,yh);pausee=yh-bh(1)-bh(2)*xh;et=zeros(n-1,1);et_1=zeros(n-1,1);for i=1:n-2 et(i)=e(i+1); et_1(i)=e(i);endpauseDW=0;sum2=0;for i=2:n-1 DW=DW+(e(i)-e(i-1)2; sum2=sum2+e(i)*e(i);endDW=DW/sum2圖像及結(jié)果,執(zhí)行代碼可以得出如下結(jié)果:ans = -1.4548 -1.9047 -1.0048 0.1763 0.1732 0.1793stats = 1.0e+004 * 0.0001 1.4888 0 0.0000(3)點擊執(zhí)行按鈕,運行其代碼。y對x的散點圖如上擬合圖如下所示:9實驗二 Lingo的簡單用法1. 實驗課程名稱 數(shù)學(xué)實驗2. 實驗項目名稱 利用Lingo解決簡單的線性問題3. 實驗?zāi)康暮鸵?了解Lingo的基本知識,熟悉上機環(huán)境,掌握Lingo的基本用法解決簡單的線性問題。4. 實驗內(nèi)容和原理內(nèi)容:某儲蓄所每天的營業(yè)時間是上午9:00到下午5:00。根據(jù)經(jīng)驗,每天不同時間段所需要的服務(wù)員數(shù)量如下:時間段9101011111212112233445服務(wù)員數(shù)量43465688儲蓄所可以雇傭全時和半時兩類服務(wù)員,全時服務(wù)員每天報酬100元,從上午9:00到下午5:00工作,但中午12:00到下午2:00之間必須安排1h的午餐時間。儲蓄所每天可以雇傭不超過三名的半時服務(wù)員,每個半時服務(wù)員必須連續(xù)工作4h,報酬40元。問該儲蓄所應(yīng)如何雇傭全時和半時兩類服務(wù)員?如果不能雇傭半時服務(wù)員,每天至少增加多少費用?如果雇傭半時服務(wù)員的數(shù)量沒有限制,每天可以減少多少費用?原理:利用Lingo編寫的程序得到最優(yōu)解,得出問題的答案。5. 主要儀器設(shè)備計算機與Windows 2000/XP系統(tǒng),Lingo軟件。6.操作方法與實驗步驟步驟:1:設(shè)儲蓄所每天雇傭的全時服務(wù)員中以12:001:00為午餐時間的有x1名,以1:002:00為午餐時間的有x2名;半時服務(wù)員中從9:00,10:00,11:00,12:00,1:00開始工作的分別為y1,y2,y3,y4,y5名。 2:打開Lingo,輸入以下程序Min=100*x1+100*x2+40*y1+40*y2+40*y3+40*y4+40*y5 ; ; ; ; ; ; ; ;且為整數(shù)。7. 實驗結(jié)果和分析實驗結(jié)果:1:經(jīng)過Lingo最

溫馨提示

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

評論

0/150

提交評論