數(shù)學(xué)建模題目

1、1、山區(qū)地貌：在某山區(qū)測得一些地點的高程如下表：（平面區(qū)域1200=x=4000,1200=y|d.-1-fl5005IntelLtwl26758RM5I：-3239340.981322業(yè)由飛印:3.570440T=455.337ModelHisto-y323 uoi/i-c:17311圖1逐步回歸分析結(jié)果由上圖可以看出：紅色說明從模型中移去的變量為X2、X4、X6,移除這三個變量后，再利用最小二乘法擬合一個多元回歸模型，有y276.76840.6077x10.7501x30.0041x52這個回歸模型的復(fù)判定系數(shù)R0.9813,調(diào)整復(fù)判定系數(shù)R0.9792。模型的剩余標(biāo)準(zhǔn)差為32.3934。

2、對模型進行F檢驗：F455.3368兩個回歸模型相比較，得到：后者的復(fù)判定系數(shù)與調(diào)整復(fù)判定系數(shù)的差距更小，與實際更加符合，因此所做的調(diào)整是有意義的，對于預(yù)測更加有利。數(shù)據(jù)沒有我全，沒法預(yù)測檢驗?zāi)攴輫袷杖雰|元工業(yè)總產(chǎn)值（億元）農(nóng)業(yè)總產(chǎn)值億元總?cè)丝谌f人就業(yè)人口萬人固定資產(chǎn)投資億元政入X)財明”198253312383177710165412001212：19835986P2646197811030081369P1367198472443106231610435718331643國859041P38671256411058512543P2005198610274449327891075073121

3、21221987120515252323310930037922199198815037658738651110264754235719891700172784266112704441026651990187187717506211433345172937表219821990年財政收入數(shù)據(jù)年份198219831984198519861987198819891990實際財政收入121213671643200521222199235726652937預(yù)測財政收入預(yù)測與實際之差平均值力差1多元回歸模型建立的程序：clc,clearloaddata.txt%表中的數(shù)據(jù)按照原來的排列存放在純文本文件da

4、ta.txt中n,m=size(data);m=m-1;x=ones(30,1),data(:,1:6);y=data(:,7);b,bint,r,rint,stats=regress(y,x)%stats(4贖回的是殘差的樣本方差r2=stats(1)%提出復(fù)判定系數(shù)ad_r2=1-(1-r2)*(n-1)/(n-m-1)%計算調(diào)整復(fù)判斷系數(shù)f=stats(2)%提出F統(tǒng)計量tm=inv(x*x);%計算X*X的逆矩陣tm=diag(tm);%提出逆矩陣的對角線元素rmse=sqrt(stats(4)%計舁剩余標(biāo)準(zhǔn)差殘差的樣本標(biāo)準(zhǔn)差t=b./sqrt(tm)/rmse%求t統(tǒng)計量的值2逐步回

5、歸分析程序:clc,clearx0=195259834946157482207294418419535864554755879621364892161954707520491602662183297248195573755852961465223289825419568257155566282823018150268195783779857564653237111392861958102812355986599426600256357195911141681509672072617333844419601079187044466207258803805061961757115643465859

6、255901382711962677964461672952511066230196377910465146917226640852661964943125058470499277361293231965115215816327253828670175393196613221911687745422980521246619671249164769776368308141563521968118715656807853431915127303196913722101688806713322520744719701638274776782992344323125641971178031567908

7、522935620355638197218333365789871773585435465819731978368485589211366523746911974199336968919085937369393655197521214254932924213816846269219762052430995593717388344436571977218949259719497439377454723197824755590105896259398565509221979270260651150975424058156489019802791659211949870541896568826198

8、129276862127310007273280496810;x=x0(:,2:7);y=x0(:,8);stepwise(x,y,1:6)4、某?；饡幸还P數(shù)額為M元的基金，打算將其存入銀行或購買國庫卷。當(dāng)前銀行存款及各期國庫卷的利率見下表。假設(shè)國庫卷每年至少發(fā)行一次，發(fā)行時間不定。銀行存款稅后年利率(%)國修年利率(%)活期0.792半年期1.664一年期1.800二年期1.9442.55三年期2.1602.89五年期2.3043.14?；饡媱澰趎年內(nèi)每年用部分本息獎勵優(yōu)秀師生，要求每年的獎金額大致相同，且在n年末仍保留原基金數(shù)額。?；饡Ｍ@得最正確的基金使用計劃，以提高每年的獎

9、金額。請你幫助?；饡谌缦虑闆r下設(shè)計基金使用方案，并對M=5000萬元，n=10年及n=12年給出具體結(jié)果：1只存款不購國庫卷；2可存款也可購國庫卷；3學(xué)校在基金到位后的第3年要舉行百年校慶，基金會希望這一年的獎金比其他年度多20%以Xj(i0n1,j19)代表第i年第j種投資方式投資的金額數(shù)目以萬元計，其中i16代表存款數(shù)目，i7,8,9代表國庫券數(shù)目，yi(i1n)代表第i年末得到的本息和，z(i1n)代表第i年發(fā)的獎金額，由于每一年的獎金額大致相同，在此取zz,z為常數(shù)。假設(shè)該校每年的1月1日發(fā)獎金，且存款等都是1月1日進行的。1只存款不購國庫卷以獎金額最大為目標(biāo)函數(shù)，其約束條件有：6

10、?起初各種投資方式之和為M,即XojMj16?每一年得到的本息和，全部用于再次投資和獎金，即yzXj0j1?第n年末仍保留原基金數(shù)額M,即ynznM。所以可以得到目標(biāo)規(guī)劃為:目標(biāo)函數(shù)max6xojMji6s.tyizXij0jiynZnM當(dāng)M=5000萬元，n10年時，用lingo數(shù)學(xué)軟件編程計算得到如下結(jié)果:z108.787萬元，具體每一年各種投資方式的投資金額如下表所示：1活期2半年期3一年期4二年期5三年期6五年期第0年00396.7621200.4946195.61404207.129第1年0000195.614098.47287第2年0000098.47287第3年0000098.4

11、7287第4年0000098.47287第5年000004581.974第6年000000第7年000000第8年000000第9年000000表1n=10時，存款方式及金額當(dāng)n12時，z178.9138萬元，具體的投資方式及金額為:1活期2半年期3一年期4二年期5三年期6五年期第0年00474.56794150.110191.9997183.3225第1年0000191.9997183.3225第2年000004203.678第3年0000096.65341第4年0000096.65341第5年000096.65341第6年0000096.65341第7年000004580.154第8年00

12、0000第9年000000第10年000000第11年000000表2n=12時，存款方式及金額2可存款也可購國庫券根據(jù)基本原則，由于國庫券的利率更大，我們應(yīng)優(yōu)先考慮買國庫券。與1類似，進行目標(biāo)規(guī)劃，建立數(shù)學(xué)模型，lingo求解。n=10時，z148.857萬元，具體的投資金額為：01234567891活期00000000002半年期00000000003一年期277.01260000000004兩年期00000000005三年期00000000006五年期00000000007兩年國庫券260.50200000000008三年國庫券490.7406135.14110259.500200135

13、.14110009五年國庫券3971.7450126.9298126.929804448.4510000表3n=10時，存款及國庫券各年投資金額n12時，z280.5156萬元，具體的投資金額為:P012345678910111活期0000000000002半年期0000000000003一年期1095.00600000000005481.8424兩年000000000000期5三年期0000000000006五年期0000000000007兩年國庫券1029.741000000000008三年國庫券1939.853135.141101025.78100534.2004000009五年國庫券9

14、35.39990501.7421501.74210501.74215239.72200000表4n=12時，存款及國庫券各年投資金額3第三年獎金高于其他年20%類似于2，只是在第三年的獎金額為1.2Z,其他條件不變，得到每年獎金額Z為：135.7854萬元。投資方案如下表：年份01234567891活期00000000002半年期00000000003一年期271.2172000000271.2172004兩年期00000000005三年期00000000006五年期00000000007兩年國庫券255.0520000000255.0520008三年國庫券506.9365132.313802

15、54.071200132.3138506.9365132.313809五年國庫券3966.7940124.2743124.274304445.79503966.7940124.2743表5第三年獎金高于其他年份，n=10時的投資方案當(dāng)n=12時，得到z268.3706萬元，具體的投資方案如下表所示:012345678910111活期0000000000002半年期0000000000003一年1072.0000000000546期0979.9124兩年期0000000000005三年期0000000000006五年期0000000000007兩年國庫券1008.198000000000008三

16、5國庫學(xué)2003.875523.024401004.32000523.0244000009五年國庫券915.83040491.2452491.24520491.24525218.91300000表6第三年獎金高于其他年份，n=12時的投資方案lingo編寫的程序為：.model:sets:Var/z,x01,x02,x03,x05,x06,x11,x12,x13,x15,x16,x21,x22,x23,x25,x26,x31,x32,x33,x35,x36,x41,x42,x43,x45,x46,x51,x52,x53,x55,x56,x61,x62,x63,x64,x65,x71,x72,x

17、73,x74,x75,x81,x82,x83,x84,x91,x92,x93/;endsetsmax=z;x01+x02+x03+x05+x06=5000;x11+x12+x13+x15+x16-1.00792*x01-1.01664*x02-1.018*x03+z=0;x21+x22+x23+x25+x26-1.01944*x04-1.018*x13-1.00792*x11-1.01664*x12+z=0;x31+x32+x33+x35+x36-1.02160*x05-1.01994*x14-1.00792*x21-1.01664*x22-1.018*x23+z=0;x41+x42+x43+

18、x45+x46-1.00792*x31-1.01664*x32-1.018*x33-1.01994*x24-1.02160*x15+z=0;x51+x52+x53+x55+x56-1.00792*x41-1.01664*x42-1.018*x43-1.02304*x06-1.01944*x34-1.02160*x25+z=0;x61+x62+x63+x64+x65-1.00792*x51-1.01664*x52-1.018*x53-1.02304*x16-1.01944*x44-1.02160*x35+z=0;x71+x72+x73+x74+x75-1.00792*x61-1.01664*x6

19、2-1.018*x63-1.02304*x26-1.01944*x54-1.02160*x45+zx81+x82+x83+x84-1.00792*x71-1.01664*x72-1.018*x73-1.02304*x36-1.01944*x64-1.02160*x55+z=0;x91+x92+x93-1.00792*x81-1.01664*x82-1.018*x83-1.02304*x46-1.01944*x74-1.02160*x65+z=5000;z5000;end(2)model:max=z;X11+Y12+Y13+Y15=5000;X21+Y22+Y23+Y25=1.018*X11-z

20、;X31+Y32+Y33+Y35=1.051*Y12+1.018*X21-z;X41+Y42+Y43+Y45=1.0867*Y13+1.051*Y22+1.018*X31-z;X51+Y52+Y53+Y55=1.0867*Y23+1.051*Y32+1.018*X41-z;X61+Y62+Y63+Y65=1.157*Y15+1.0867*Y33+1.051*Y42+1.018*X51-z;X71+Y72+Y73=1.157*Y25+1.0867*Y43+1.051*Y52+1.018*X61-z;X81+Y82+Y83=1.157*Y35+1.0867*Y53+1.051*Y62+1.018*

21、X71-z;X91+Y92=1.157*Y45+1.0867*Y63+1.051*Y72+1.018*X81-z;X101=1.157*Y55+1.0867*Y73+1.051*Y82+1.018*X91-z;5000=1.157*Y65+1.0867*Y83+1.051*Y92+1.018*X101-z;endmodel:max=z;X11+X12+X13+X15+Y12+Y13+Y15=5000;X21+X22+X23+X25+Y22+Y23+Y25=1.018*X11-z;X31+X32+X33+X35+Y32+Y33+Y35=1.03888*X12+1.018*X21-z;X41+X4

22、2+X43+X45+Y42+Y43+Y45=1.0648*X13+1.03888*X22+1.018*X31+1.06394*Y12-z;X51+X52+X53+X55+Y52+Y53+Y55=1.0648*X23+1.03888*X32+1.018*X41+1.06394*Y22+1.10008*Y13-z;X61+X62+X63+X65+Y62+Y63=1.1152*X15+1.0648*X33+1.03888*X42+1.018*X51+1.06394*Y32+1.10008*Y23-z;X71+X72+X73+Y72+Y73=1.1152*X25+1.0648*X43+1.03888*

23、X52+1.018*X61+1.06394*Y42+1.10008*Y33+1.17125*Y15-z;X81+X82+X83+Y82=1.1152*X35+1.0648*X53+1.03888*X62+1.018*X71+1.06394*Y52+1.10008*Y43+1.17125*Y25-z;X91+X92=1.1152*X45+1.0648*X63+1.03888*X72+1.018*X81+1.06394*Y62+1.10008*Y53+1.17125*Y35-z;X101=1.1152*X55+1.0648*X73+1.03888*X82+1.018*X91+1.06394*Y72

24、+1.10008*Y63+1.17125*Y45-z;5000=1.1152*X65+1.0648*X83+1.03888*X92+1.018*X101+1.06394*Y82+1.10008*Y73+1.17125*Y55-z;endmodel:max=z;X11+X12+X13+X15+Y12+Y13+Y15=5000;X21+X22+X23+X25+Y22+Y23+Y25=1.018*X11-z;X31+X32+X33+X35+Y32+Y33+Y35=1.03888*X12+1.018*X21-z;X41+X42+X43+X45+Y42+Y43+Y45=1.0648*X13+1.0388

25、8*X22+1.018*X31+1.06856*Y12-z;X51+X52+X53+X55+Y52+Y53+Y55=1.0648*X23+1.03888*X32+1.018*X41+1.06856*Y22+1.10486*Y13-z;X61+X62+X63+X65+Y62+Y63=1.1152*X15+1.0648*X33+1.03888*X42+1.018*X51+1.06856*Y32+1.10486*Y23-z;X71+X72+X73+Y72+Y73=1.1152*X25+1.0648*X43+1.03888*X52+1.018*X61+1.06856*Y42+1.10486*Y33+1

26、.17633*Y15-z;X81+X82+X83+Y82=1.1152*X35+1.0648*X53+1.03888*X62+1.018*X71+1.06856*Y52+1.10486*Y43+1.17633*Y25-z;X91+X92=1.1152*X45+1.0648*X63+1.03888*X72+1.018*X81+1.06856*Y62+1.10486*Y53+1.17633*Y35-z;X101=1.1152*X55+1.0648*X73+1.03888*X82+1.018*X91+1.06856*Y72+1.10486*Y63+1.17633*Y45-z;5000=1.1152*

27、X65+1.0648*X83+1.03888*X92+1.018*X101+1.06856*Y82+1.10486*Y73+1.17633*Y55-z;End(3)model:max=z;X11+X12+X13+X15+Y12+Y13+Y15=5000;X21+X22+X23+X25+Y22+Y23+Y25=1.018*X11-z;X31+X32+X33+X35+Y32+Y33+Y35=1.03888*X12+1.018*X21-z;X41+X42+X43+X45+Y42+Y43+Y45=1.0648*X13+1.03888*X22+1.018*X31+1.06856*Y12-1.2*z;X5

28、1+X52+X53+X55+Y52+Y53+Y55=1.0648*X23+1.03888*X32+1.018*X41+1.06856*Y22+1.10486*Y13-z;X61+X62+X63+X65+Y62+Y63=1.1152*X15+1.0648*X33+1.03888*X42+1.018*X51+1.06856*Y32+1.10486*Y23-z;X71+X72+X73+Y72+Y73=1.1152*X25+1.0648*X43+1.03888*X52+1.018*X61+1.06856*Y42+1.10486*Y33+1.17633*Y15-z;X81+X82+X83+Y82=1.1

29、152*X35+1.0648*X53+1.03888*X62+1.018*X71+1.06856*Y52+1.10486*Y43+1.17633*Y25-z;X91+X92=1.1152*X45+1.0648*X63+1.03888*X72+1.018*X81+1.06856*Y62+1.10486*Y53+1.17633*Y35-z;X101=1.1152*X55+1.0648*X73+1.03888*X82+1.018*X91+1.06856*Y72+1.10486*Y63+1.17633*Y45-z;5000=1.1152*X65+1.0648*X83+1.03888*X92+1.018

30、*X101+1.06856*Y82+1.10486*Y73+1.17633*Y55-z;endmodel:max=z;X11+X12+X13+X15+Y12+Y13+Y15=5000;X21+X22+X23+X25+Y22+Y23+Y25=1.018*X11-z;X31+X32+X33+X35+Y32+Y33+Y35=1.03888*X12+1.018*X21-z;X41+X42+X43+X45+Y42+Y43+Y45=1.0648*X13+1.03888*X22+1.018*X31+1.06394*Y12-1.2*z;X51+X52+X53+X55+Y52+Y53+Y55=1.0648*X2

31、3+1.03888*X32+1.018*X41+1.06394*Y22+1.10008*Y13-z;X61+X62+X63+X65+Y62+Y63=1.1152*X15+1.0648*X33+1.03888*X42+1.018*X51+1.06394*Y32+1.10008*Y23-z;X71+X72+X73+Y72+Y73=1.1152*X25+1.0648*X43+1.03888*X52+1.018*X61+1.06394*Y42+1.10008*Y33+1.17125*Y15-z;X81+X82+X83+Y82=1.1152*X35+1.0648*X53+1.03888*X62+1.01

32、8*X71+1.06394*Y52+1.10008*Y43+1.17125*Y25-z;X91+X92=1.1152*X45+1.0648*X63+1.03888*X72+1.018*X81+1.06394*Y62+1.10008*Y53+1.17125*Y35-z;X101=1.1152*X55+1.0648*X73+1.03888*X82+1.018*X91+1.06394*Y72+1.10008*Y63+1.17125*Y45-z;5000=1.1152*X65+1.0648*X83+1.03888*X92+1.018*X101+1.06394*Y82+1.10008*Y73+1.17125*Y55-z;endmodel:max=z;X11+Y12+Y13+Y15=5000;X21+Y22+Y23+Y25=1.018*X11-z;X31+Y32+Y33+Y35=1.051*Y12+1.018*X21-z;X41+Y42+Y43+Y45=1.0867*Y13+1.051*Y22+1.018*X31-1.2*z;X51+Y52+Y53+Y55=1.0867*Y23+1.051*Y32+1.018*X41-z;X61+Y62+Y63+Y65=1.157*Y15+1.0867*Y33+1.051*Y42+1.018*X51-z;X