數(shù)學(xué)模型報(bào)告新作

9 實(shí) 驗(yàn) 報(bào) 告課 程 數(shù)學(xué)模型實(shí)驗(yàn) 系 別 數(shù)統(tǒng)學(xué)院 班 級(jí) 2012級(jí) 學(xué) 號(hào) 2012325433 姓 名 王靜怡 實(shí)驗(yàn)報(bào)告一：Matlab解決預(yù)測(cè)問題1實(shí)驗(yàn)課程名稱 數(shù)學(xué)實(shí)驗(yàn)2實(shí)驗(yàn)項(xiàng)目名稱 Matlab解決預(yù)測(cè)問題3 實(shí)驗(yàn)?zāi)康暮鸵?了解用Matlab軟件求解圖論模型及層次分析模型的方法。4 實(shí)驗(yàn)內(nèi)容和原理內(nèi)容：某公司想用全行業(yè)的銷售額作為自變量來預(yù)測(cè)公司的銷售額，下表給出了19771981年的公司銷售額和行業(yè)銷售額的分季度數(shù)據(jù)（單位：百萬元）.(1) 畫出數(shù)據(jù)的散點(diǎn)圖，觀察用線性回歸模型擬合是否合適.(2) 建立公司銷售額對(duì)全行業(yè)銷售額的回歸模型，并用DW檢驗(yàn)診斷隨機(jī)誤差項(xiàng)的自相關(guān)性.(3) 建立消除了隨機(jī)誤差項(xiàng)自相關(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è)備 計(jì)算機(jī)與Windows 2000/XP系統(tǒng)；Matlab等軟件。6操作方法與實(shí)驗(yàn)步驟步驟：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）點(diǎn)擊執(zhí)行按鈕，運(yùn)行其代碼。y對(duì)x的散點(diǎn)圖如上擬合圖如下所示：9實(shí)驗(yàn)二 Lingo的簡(jiǎn)單用法1. 實(shí)驗(yàn)課程名稱 數(shù)學(xué)實(shí)驗(yàn)2. 實(shí)驗(yàn)項(xiàng)目名稱 利用Lingo解決簡(jiǎn)單的線性問題3. 實(shí)驗(yàn)?zāi)康暮鸵?了解Lingo的基本知識(shí)，熟悉上機(jī)環(huán)境，掌握Lingo的基本用法解決簡(jiǎn)單的線性問題。4. 實(shí)驗(yàn)內(nèi)容和原理內(nèi)容：某儲(chǔ)蓄所每天的營(yíng)業(yè)時(shí)間是上午9:00到下午5:00。根據(jù)經(jīng)驗(yàn)，每天不同時(shí)間段所需要的服務(wù)員數(shù)量如下：時(shí)間段9101011111212112233445服務(wù)員數(shù)量43465688儲(chǔ)蓄所可以雇傭全時(shí)和半時(shí)兩類服務(wù)員，全時(shí)服務(wù)員每天報(bào)酬100元，從上午9:00到下午5:00工作，但中午12:00到下午2:00之間必須安排1h的午餐時(shí)間。儲(chǔ)蓄所每天可以雇傭不超過三名的半時(shí)服務(wù)員，每個(gè)半時(shí)服務(wù)員必須連續(xù)工作4h，報(bào)酬40元。問該儲(chǔ)蓄所應(yīng)如何雇傭全時(shí)和半時(shí)兩類服務(wù)員？如果不能雇傭半時(shí)服務(wù)員，每天至少增加多少費(fèi)用？如果雇傭半時(shí)服務(wù)員的數(shù)量沒有限制，每天可以減少多少費(fèi)用？原理：利用Lingo編寫的程序得到最優(yōu)解，得出問題的答案。5. 主要儀器設(shè)備計(jì)算機(jī)與Windows 2000/XP系統(tǒng)，Lingo軟件。6.操作方法與實(shí)驗(yàn)步驟步驟：1：設(shè)儲(chǔ)蓄所每天雇傭的全時(shí)服務(wù)員中以12:001:00為午餐時(shí)間的有x1名，以1:002:00為午餐時(shí)間的有x2名；半時(shí)服務(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. 實(shí)驗(yàn)結(jié)果和分析實(shí)驗(yàn)結(jié)果：1:經(jīng)過Lingo最