【文章內(nèi)容簡(jiǎn)介】 第一年的投資資金(2)到第二年的時(shí)候，該年的投資資金=第一年投資剩余的資金+第一年投資項(xiàng)目2收回的本金和捐贈(zèng)資金，即：(3)到第三年的時(shí)候，該年的投資金額=第一年投資剩余的資金+項(xiàng)目4收回的本金，即：(4)到第四年的時(shí)候，投資資金需滿足的條件為：(5)到第五年的時(shí)候，投資資金需滿足的條件為：確定模型表達(dá)式：綜上可以得到，公司收益最大的單目標(biāo)優(yōu)化模型為：根據(jù)上述建立的單目標(biāo)最優(yōu)化模型，運(yùn)用MATLAB求解(求解代碼見附錄七)，得到該公司的決策如下表(14)：表(14)：?jiǎn)栴}五的投資方案(億元)第一年第二年第三年第四年第五年總利潤(rùn)項(xiàng)目1066項(xiàng)目260066項(xiàng)目30項(xiàng)目43330項(xiàng)目500項(xiàng)目644400項(xiàng)目703000項(xiàng)目800000貸款0200存款0000、改進(jìn)與推廣優(yōu)點(diǎn)：(1) 模型一中運(yùn)用單目標(biāo)優(yōu)化模型進(jìn)行求解，對(duì)列出的表達(dá)式通過Lingo軟件編程進(jìn)行計(jì)算，計(jì)算過程嚴(yán)密，模型易懂，算法不太難，且思路清晰，便于理解。(2) 對(duì)于問題二中預(yù)測(cè)今后五年各項(xiàng)目獨(dú)立投資及項(xiàng)目之間相互影響下的投資的到期利潤(rùn)率、風(fēng)險(xiǎn)損失率的問題運(yùn)用了灰色預(yù)測(cè)模型。不需要大量的樣本，也不需要考慮序列的分布規(guī)律，能很好的對(duì)其進(jìn)行預(yù)測(cè)，并引入了關(guān)聯(lián)度對(duì)其進(jìn)行驗(yàn)證得到比較滿意的結(jié)果。(3) 模型建立過程中，數(shù)據(jù)量較大，利用計(jì)算機(jī)軟件對(duì)數(shù)據(jù)進(jìn)行處理，使模型建立時(shí)的計(jì)算量大大降低。(4) 模型的適用性強(qiáng)，可以解決類似的投資求最大收益的問題。缺點(diǎn)：(1)模型只考慮了項(xiàng)目4，6，8同時(shí)投資時(shí)的相互影響，沒有考慮其他項(xiàng)目同時(shí)投資時(shí)的相互影響，也沒有考慮項(xiàng)目不同時(shí)投資但相互影響的情況。(2) GM(1，1)模型不是通用的，即不是對(duì)所有的預(yù)測(cè)都有效，只有對(duì)臨近時(shí)間段的預(yù)測(cè)的準(zhǔn)確率才較高，隨著時(shí)間的發(fā)展，其誤差率越來越大。模型中只考慮了項(xiàng)目4，項(xiàng)目6以及項(xiàng)目8同時(shí)投資時(shí)的影響，可以考慮其他項(xiàng)目同時(shí)投資時(shí)的影響，使模型更準(zhǔn)確。所給出的數(shù)據(jù)有限，可以調(diào)查更多的數(shù)據(jù)，提高預(yù)測(cè)模型的準(zhǔn)確程度。單目標(biāo)規(guī)劃模型不僅可以解決公司投資決策問題，還可以解決其他求最大，最小值問題，如：油田開發(fā)的產(chǎn)量?jī)?yōu)化，城市道路網(wǎng)雙層規(guī)劃問題。[1] 薛秀謙 朱開永，.[2] 肖偉 劉忠，.[3] 盛驟 謝式千,.[4] .[5] .[6] (中國(guó)人民銀行，貨幣政策司)附錄一：?jiǎn)栴}一中運(yùn)用MATLAB的求解代碼sets:a/1..8/:costrate。y/1..5/:amoney。d/1..8/:mmoney。yeartime/1..5/ 。task/1..8/ 。x(yeartime,task):money。endsetsdata:costrate=,,。mmoney=6,3,4,3,3,2,4,3。enddatamax=costrate(1)*money(5,1)+costrate(2)*money(5,2)+costrate(3)*money(4,3)+costrate(4)*money(4,4)+costrate(5)*money(3,5)+costrate(6)*money(3,6)+costrate(7)*money(2,7)+costrate(8)*money(3,8)amoney(1)。@for(x(i,j):money(i,j)=mmoney(j))。amoney(1)=20。@for(y(i)|ine2:money(i,7)=0)。@for(y(i)|ine3:money(i,8)=0)。@for(x(i,j)|(ige4)and(jge5):money(i,j)=0)。money(5,3)=0。money(5,4)=0。@sum(task(j)|jle6:money(1,j))=amoney(1)。amoney(2)=costrate(1)*money(1,1)+costrate(2)*money(1,2)。@sum(task(j)|jle7:money(2,j))=amoney(2)。amoney(3)=costrate(1)*money(2,1)+costrate(2)*money(2,2)+costrate(3)*money(1,3)+costrate(4)*money(1,4)。money(3,8)+@sum(task(j)|jle6:money(3,j))=amoney(3)。amoney(4)=costrate(1)*money(3,1)+costrate(2)*money(3,2)+costrate(3)*money(2,3)+costrate(4)*money(2,4)+costrate(5)*money(1,5)+costrate(6)*money(1,6)。@sum(task(j)|jle4:money(4,j))=amoney(4)。amoney(5)=costrate(1)*money(4,1)+costrate(2)*money(4,2)+costrate(3)*money(3,3)+costrate(4)*money(3,4)+costrate(5)*money(2,5)+costrate(6)*money(2,6)。money(5,1)+money(5,2)=amoney(5)。附錄二%獨(dú)立投資加權(quán)平均移動(dòng)預(yù)測(cè)第七個(gè)項(xiàng)目預(yù)測(cè)clear allclcdata=importdata(39。H:\數(shù)學(xué)建模作業(yè)\第四輪作業(yè)\39。)。for i=1:20 for j=1:8 A(i,j)=data(2*i,j)/data(2*i1,j)。 endendw=[,]。for i=1:3 A(i+17,7)=A(i+173:i+171,7)39。*w39。endA(1:20,1)%%%混合投資第七個(gè)項(xiàng)目加權(quán)平均移動(dòng)預(yù)測(cè)代碼clear allclcdata=importdata(39。H:\數(shù)學(xué)建模作業(yè)\第四輪作業(yè)\39。)。for i=1:19 for j=1:7 A(i,j)=data(2*i,j)/data(2*i1,j)。 endendw=[,]。for i=1:2 A(i+18,7)=A(i+194:i+192,7)39。*w39。endA(1:20,7) 附錄三%獨(dú)立投資第一個(gè)項(xiàng)目灰色預(yù)測(cè)代碼clear allclcdata=importdata(39。H:\數(shù)學(xué)建模作業(yè)\第四輪作業(yè)\39。)。for i=1:20 for j=1:8 A(i,j)=data(2*i,j)/data(2*i1,j)。 endendx0=A(1:20,1)。n=length(x0)。x1(1)=x0(1)。for i=1:n1 x1(i+1)=x1(i)+x0(i+1)。 endfor k=2:n z(k)=*x1(k)+*x1(k1)。endz(1)=x0(1)。B=[z(2:n)39。,ones(size(z(2:n)))39。]。Y=x0(2:n)。E=inv(B39。*B)*B39。*Y。for i=1:n+4 X(i+1)=(x0(1)E(2)/E(1))*exp(E(1)*i)+E(2)/E(1)。endX(1)=x0(1)。for i=2:n+5 X0(i)=X(i)X(i1)。endX0(1)=X(1)。X0%%%混合投資第一個(gè)項(xiàng)目灰色預(yù)測(cè)代碼clear allclcdata=importdata(39。C:\Documents and Settings\Administrator\桌面\39。)。for i=1:19 for j=1:7 A(i,j)=data(2*i,j)/data(2*i1,j)。 endendx0=A(1:19,1)。n=length(x0)。x1(1)=x0(1)。for i=1:n1 x1(i+1)=x1(i)+x0(i+1)。 endx1for k=2:n z(k)=*x1(k)+*x1(k1)。endz(1)=x0(1)。B=[z(2:n)39。,ones(size(z(2:n)))39。]。Y=x0(2:n)。E=inv(B39。*B)*B39。*Y。for i=1:n+5 X(i+1)=(x0(1)E(2)/E(1))*exp(E(1)*i)+E(2)/E(1)。endX(1)=x0(1)。for i=2:n+6 X0(i)=X(i)X(i1)。endX0(1)=X(1)。X0X0附錄四%該程序主要采用逐步外推預(yù)測(cè)的方法，能夠較為精準(zhǔn)的預(yù)測(cè)將來的發(fā)展趨勢(shì)。%clcclearclose all%% 原始數(shù)據(jù)輸入%%pp為項(xiàng)目一1986至2005年利潤(rùn)率 pp=[。。。]39。Sr_1=zeros(3,17)。Sc_1=zeros(1,17)。 for k=1:4 if k==4 Sc_1=pp(4:20)。 else Sr_1(k,:)=pp(k:k+16)。 endend%處理原始數(shù)據(jù)將值轉(zhuǎn)到[0 1]區(qū)間，pn未處理后數(shù)據(jù) Sr_2=(Sr_1)/。 Sc_2=(Sc_1)/。 NodeNum1 =5。 % 隱層第一層節(jié)點(diǎn)數(shù)%NodeNum2=40。 % 隱層第二層節(jié)點(diǎn)數(shù)TypeNum = 1。 % 輸出維數(shù)TF1 = 39。tansig39。TF2 = 39。logsig39。 TF3 = 39。tansig39。net=newff(minmax(Sr_2),[NodeNum1,TypeNum],...{TF1 TF3},39。traingdx39。)。%網(wǎng)絡(luò)創(chuàng)建%traingdm=50。=50000。%訓(xùn)練次數(shù)設(shè)置=1e50。%訓(xùn)練目標(biāo)設(shè)置=。 net=train(net,Sr_2,Sc_2)。an=sim(net,Sr_2)。a=an*。[m,b,r]=postreg(a,Sc_1)。 %輸出值與原值的線性擬合save mynetdate nety=an39。 plot(1:length(Sc_1),Sc_1,39。*39。,1:length(Sc_1),a,39。o39。)。title(39。o表示預(yù)測(cè)值 *表示實(shí)際值39。)grid onerror=(Sc_1a)./Sc_1。figure plot(1:length(error),error,39。.39。)title(39。相對(duì)誤差變化圖39。)grid on%================================================================ %以下為未來利潤(rùn)率預(yù)測(cè)程序 yc(1)=pp(20)。 test=pp(17:19)。for m=1:6test=[test(2:3) yc(m)]。sdtest=test/。testres=sim(net,sdtest39。)。res=testres*。yc(m+1)=res。end yc %yc(2:6)是今后五年的預(yù)測(cè)銷售量附錄五問題三運(yùn)用Lingo的求解代碼sets:y/1..5/:amoney。d/1..8/:mmoney。yeartime/1..5/ 。yeartime2/1..4/。task/1..8/ 。x(yeartime,task):money。b(yeartime,task):costrate。c(yeartime2,task):costratexh。endsetsdata:mmoney=6,6,3,3,4,3,3。costrate=