【正文】 subplot(2,2,1),plot(x,y,39。y=p39。) y= *6^3+ *6^2+ *6+y = clear a=[ ]。) y= *6^3+ *6^*6+y = clear b=[2550 3065 3662 5024 7749100 100 100 132 4351788 1843 2725 2552 33782858 4010 5325 5112 71011498 1580 1772 2680 30711920 1950 2170 3520 4243]39。) y= *6^3+ *6^*6+y = clear a=[ ]。subplot(2,2,1),plot(x,y,39。p=b*a/10^4。 result(:,5)=temp*u39。 end temp(:,1)=(L./sum)39。 ]。M(5)1/2*h(5)。M(7)w(7)。附錄[1]生產(chǎn)計(jì)劃準(zhǔn)確度表學(xué)科類(lèi)別20012002200320042005數(shù)據(jù)偏差的平方和計(jì)劃可信度均值計(jì)劃可信度均值　計(jì)算機(jī)類(lèi)經(jīng)管類(lèi)數(shù)學(xué)類(lèi)英語(yǔ)類(lèi)兩課類(lèi)機(jī)械、能源類(lèi)化學(xué)、化工類(lèi)地理、地質(zhì)類(lèi)環(huán)境類(lèi)[2]讀者滿(mǎn)意度表　2001200220032004200520012005期望平均計(jì)算機(jī)類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差經(jīng)管類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差數(shù)學(xué)類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差外語(yǔ)類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差兩課類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差機(jī)械、能源類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差化學(xué)、化工類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差地理、地質(zhì)類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差環(huán)境類(lèi)滿(mǎn)意度期望滿(mǎn)意度方差[3]市場(chǎng)占有率分析 p3 = y=*6^*6^2+*6+y = p3 = y= *6^*6^2+ *6+y = p3 = y=*6^*6^2+ *6+y = p3 = y= *6^3+*6^2+*6+ y = y= *6^3+*6^2+*6+y = p3 = y=*6^*6^2+ *6+ y = p3 = y= *6^3+ *6^2 *6+y = p3 = y= *6^3+ *6^2 + *6+ y = p3 = y=*6^*6^2+ *6+y = [4]擬合函數(shù)的選取分析在對(duì)2006年各個(gè)出版分社的總銷(xiāo)售額進(jìn)行預(yù)測(cè)的過(guò)程中，為了提高計(jì)算的精確度，我們對(duì)自變量（即年份）進(jìn)行一定的縮放處理，即使用均值與標(biāo)準(zhǔn)差對(duì)其的數(shù)據(jù)進(jìn)行規(guī)范化：X1=(xmean(x))/std(x)，同時(shí)對(duì)縱坐標(biāo)（金額）做放縮處理，均縮小10^4,以萬(wàn)元標(biāo)記，然后用規(guī)范化后的數(shù)據(jù)進(jìn)行一次，二次與三次多項(xiàng)式函數(shù)擬和，其三次擬合圖象與殘差圖象如下：,擬和程度最好。缺點(diǎn)：模型的缺點(diǎn)主要有：層次分析法受主觀判斷的影響，可能會(huì)有偏差?？傻孟铝行畔ⅲ? ① 將第（2）步所得的合成權(quán)重排序結(jié)果，利用選優(yōu)排序法，可得 ② 另外各分社2006年申請(qǐng)書(shū)號(hào)數(shù)為 表2 2006年各分社申請(qǐng)書(shū)號(hào)數(shù)分社（）申請(qǐng)數(shù)量分社（）申請(qǐng)數(shù)量1計(jì)算機(jī)類(lèi)1106機(jī)械能源類(lèi)762經(jīng)管類(lèi)667化工化學(xué)類(lèi)403數(shù)學(xué)類(lèi)2228地理地質(zhì)類(lèi)404外語(yǔ)類(lèi)1189環(huán)境類(lèi)405兩課類(lèi)72—————————————③ 人力資源的限制表3 各分社人力資源限制表分社（）年最大工作量（書(shū)號(hào)數(shù)）分社（）年最大工作量（書(shū)號(hào)數(shù)）1計(jì)算機(jī)類(lèi)1146機(jī)械能源類(lèi)722經(jīng)管類(lèi)1147化工化學(xué)類(lèi)443數(shù)學(xué)類(lèi)1208地理地質(zhì)類(lèi)634外語(yǔ)類(lèi)1029環(huán)境類(lèi)725兩課類(lèi)111———————————————— 根據(jù)上面的約束條件以及組合權(quán)系數(shù)，利用lingo軟件求解整數(shù)規(guī)劃可以求得每個(gè)分社的書(shū)號(hào)分配量，繪成表格如表4，軟件求解過(guò)程詳見(jiàn)附錄6。 求解完第三層對(duì)于第二層所有準(zhǔn)則的最大特征值的歸一化特征向量后，就可得到方案層各元素關(guān)于目標(biāo)的合成權(quán)重向量為 這樣就可以對(duì)各分社的綜合實(shí)力進(jìn)行排序了，在整個(gè)模型中，由于第三層對(duì)于第二層的一致性指標(biāo)都是0，所以模型具有整體的一致性。 具體方法如下： 首先根據(jù)讀者對(duì)不同教材（）的打分，其中為分社所擁有的課程數(shù)，將四個(gè)評(píng)價(jià)方面的分值求平均值，不妨設(shè)為，然后按對(duì)該分社進(jìn)行調(diào)查的所有問(wèn)卷的項(xiàng)求均值，設(shè)為，并按照歷年的統(tǒng)計(jì)數(shù)據(jù)計(jì)算出該分社的年平均的滿(mǎn)意度（求均值），用它來(lái)概括表示讀者對(duì)該分社的滿(mǎn)意程度。這樣，如何合理量化以上這五個(gè)評(píng)價(jià)指標(biāo)就成了我們需要解決的問(wèn)題。這樣既可以滿(mǎn)足從整體上對(duì)各種因素的綜合考慮，又能夠體現(xiàn)增強(qiáng)強(qiáng)勢(shì)產(chǎn)品的支持力度原則。 本文還借助了統(tǒng)計(jì)分析法詳細(xì)分析了幾種擬合歷年相關(guān)統(tǒng)計(jì)數(shù)據(jù)的函數(shù)的優(yōu)劣性，最后選取了較好的擬合函數(shù)，從而使得2006年的相關(guān)指標(biāo)預(yù)測(cè)更加可信。二、模型假設(shè)（1） 每年的問(wèn)卷調(diào)查結(jié)果都是有效的，即調(diào)查對(duì)象都是如實(shí)作答，這就保證了所搜集信息的可靠性，事實(shí)上，這一點(diǎn)使很容易滿(mǎn)足的，對(duì)于這類(lèi)關(guān)系切身利益的調(diào)查（因?yàn)樽x者都希望獲得更好的教材），同學(xué)們一般都會(huì)積極相應(yīng)的；（2） 當(dāng)各分社需求書(shū)號(hào)總量大于總社書(shū)號(hào)總量時(shí)，總社會(huì)按照增加強(qiáng)勢(shì)產(chǎn)品支持力度的原則來(lái)進(jìn)行書(shū)號(hào)資源的配置，這是很符合常理的，因?yàn)閺?qiáng)勢(shì)產(chǎn)品擁有更大的市場(chǎng)份額、更高的讀者滿(mǎn)意度以及更多的收益等，所以分配書(shū)號(hào)應(yīng)優(yōu)先考慮強(qiáng)勢(shì)產(chǎn)品；（3） 不考慮各分社的人力資源的調(diào)整，而僅以平均人力資源數(shù)作為評(píng)價(jià)分社生產(chǎn)能力的尺度，這樣可以適當(dāng)?shù)睾?jiǎn)化模型的復(fù)雜度；（4） 與題目相一致，我們假定同一種課程的不同書(shū)目?jī)r(jià)格差別不大，同時(shí)銷(xiāo)售量也相近，這樣就可以獲得每一中課程的均價(jià)。滿(mǎn)意度直接影響著出版社的長(zhǎng)遠(yuǎn)發(fā)展，所以也要作為一個(gè)衡量各分社綜合實(shí)力的標(biāo)準(zhǔn)。這樣通過(guò)兩種數(shù)學(xué)模型的結(jié)合，就可以確定出合理的書(shū)號(hào)優(yōu)化配置方案，較好地滿(mǎn)足題目地要求。（2）利用層次分析法（AHP）確定各分社實(shí)力排序分析題目，我們將優(yōu)化書(shū)號(hào)資源配置作為目標(biāo)層、生產(chǎn)計(jì)劃可信度、讀者滿(mǎn)意度、市場(chǎng)占有率以及利潤(rùn)作為準(zhǔn)則層，各分社的綜合實(shí)力作為方案層，由此建立層次分析模型。細(xì)化滿(mǎn)意度表（2001年）分社（）細(xì)化滿(mǎn)意度指數(shù)計(jì)算機(jī)類(lèi)經(jīng)管類(lèi)數(shù)學(xué)類(lèi)外語(yǔ)類(lèi)兩課類(lèi)機(jī)械、能源類(lèi)化學(xué)、化工類(lèi)地理、地質(zhì)類(lèi)環(huán)境類(lèi) ③ 市場(chǎng)占有率在研究市場(chǎng)占有率的時(shí)候，我們可以通過(guò)各分社歷年市場(chǎng)占有率的走勢(shì)來(lái)確定擬合函數(shù)，如圖2，以下各項(xiàng)均以計(jì)算機(jī)類(lèi)分社為例說(shuō)明，其他各類(lèi)分社請(qǐng)見(jiàn)附錄3。 表6優(yōu)化分配方案及其與歷年收益比較20012002200320042005計(jì)算機(jī)4347574846經(jīng)管2526332730數(shù)學(xué)120120120120120外語(yǔ)8587756878兩課6170637268機(jī)械能源5757525557化工化學(xué)3635292626地質(zhì)地理3328334036環(huán)境4030374438實(shí)際總收益860．6122922．91081113．9061435．8381788．8優(yōu)化方案收益900．3011980．4321138．7401450．7541799．7 從表中，我們可以看出模型的合理性，它較好地提高了出版社地年收益。十、參考文獻(xiàn)[1] 阮曉青 周義倉(cāng)編，數(shù)學(xué)建模引論，北京：高等教育出版社，[2] 魏宗舒等編，概率論與數(shù)理統(tǒng)計(jì)教程，北京：高等教育出版社，(2004 重印)[3] 吳祈宗編，運(yùn)籌學(xué)與最優(yōu)化方法，北京：機(jī)械工業(yè)出版社，[4] 姜啟源 謝金星 葉俊編，數(shù)學(xué)模型，北京：高等教育出版社，[5] 清源計(jì)算機(jī)工作室編，MATLAB基礎(chǔ)及其應(yīng)用，北京：機(jī)械工業(yè)出版社，[6] 唐五湘，程桂枝編，Excel在預(yù)測(cè)中的應(yīng)用，北京：電子工業(yè)出版社，致出版社的一封信出版社的領(lǐng)導(dǎo):你們好!根據(jù)提供的數(shù)據(jù)資料，而且查閱了一些相關(guān)的信息，得知資源配置在出版社當(dāng)中占有重要地位，它直接影響到出版社當(dāng)年的經(jīng)濟(jì)效益和長(zhǎng)遠(yuǎn)發(fā)展戰(zhàn)略。SUM(mat(i):M(i))=500。M(7)h(7)。 。 sm=bijiao(satisfy(5,:))。 for j=1:9 sum=sum+L(j)。y1i=polyval(p1,xi)。,xi,y1i,39。b=[8767 9351 11393 13461 16694441 513 699 894 1314100 183 223 269 7331581 1965 2089 4248 4247672 626 989 1465 23111348 1617 1518 2240 2955100 177 467 480 273621 735 759 1088 2123687 657 903 1316 2106173 258 448 500 759]39。*39。*39。subplot(2,2,2),plot(x,y,39。,xi,y1i,39。39。subplot(2,2,2),plot(x,y,39。)p1=polyfit(x,y,3)p1 = xi=1::5。)p1=polyfit(x,y,3)p1 = x=1:1:5。)p1=polyfit(x,y,3)p1 = x=1:1:5。39。subplot(2,2,2),plot(x,y,39。*39。 L=M(:,2)39。result=zero