數(shù)學(xué)建模,教師的薪金

1、昆明理工大學(xué)第六屆大學(xué)生數(shù)學(xué)建模競(jìng)賽承諾書(shū)我們仔細(xì)閱讀了昆明理工大學(xué)大學(xué)生數(shù)學(xué)建模競(jìng)賽的競(jìng)賽規(guī)則。我們完全明白，在競(jìng)賽開(kāi)始后參賽隊(duì)員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢(xún)等）與隊(duì)外的任何人研究、討論與賽題有關(guān)的問(wèn)題。我們知道，抄襲別人的成果是違反競(jìng)賽規(guī)則的。如果引用別人的成果或其他公開(kāi)的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻(xiàn)的表述方式在正文引用處和參考文獻(xiàn)中明確列出。我們鄭重承諾，嚴(yán)格遵守競(jìng)賽規(guī)則，以保證競(jìng)賽的公正、公平性。如有違反競(jìng)賽規(guī)則的行為，我們將受到嚴(yán)肅處理。評(píng)閱編號(hào)（由組委會(huì)評(píng)閱前進(jìn)行編號(hào)）：昆明理工大學(xué)第六屆大學(xué)生數(shù)學(xué)建模競(jìng)賽評(píng)閱專(zhuān)用頁(yè)評(píng)閱編號(hào)（由組委會(huì)評(píng)閱前進(jìn)行

2、編號(hào)）：評(píng)閱記錄（供評(píng)閱時(shí)使用）：評(píng)閱人評(píng)分備注總分大學(xué)教師薪金模型一、摘要：模型一：利用Matlab建立X1-X7與Z（新金）的線性關(guān)系，得到散點(diǎn)圖與一般的線性回歸模型（模型一），模型如下：模型二：由模型一整體效果與StepwiseTable圖推斷部分變量對(duì)Z（薪金）的影響并不顯著，由殘差分析法篩選出影響明顯的變量X1、X4，將他們的平方項(xiàng)與交互項(xiàng)加入建立新的回歸模型（模型二），模型如下：z6.90260.0043x10.1746x40.0001x1*x40*x1A20.0228x4A2得到R2,F,P,與模型一相比較。模型一調(diào)用原始參數(shù)較全，但回歸性差，模型二只啟用影響明顯的變量，回歸模型

3、更為顯著，可靠度更高。于是得出：教師的薪金與工作時(shí)間、學(xué)歷關(guān)系明顯，與性別、是否受雇于重點(diǎn)大學(xué)、是否接受過(guò)培訓(xùn)的關(guān)系較小，即女教師沒(méi)有受到不公正的待遇，婚姻狀況也不會(huì)影響收入。關(guān)鍵詞：回歸分析，互交作用，圖形結(jié)合，殘值分析法。二、問(wèn)題重述：某地人事部門(mén)為研究中學(xué)教師的薪金與他們的資歷、性別、教育程度及培訓(xùn)情況等因素之間的關(guān)系，要建立一個(gè)數(shù)學(xué)模型，分析人事策略的合理性，特別是考察女教師是否受到不公平的待遇，以及她們的婚姻狀況是否會(huì)影響收入。要求從當(dāng)?shù)亟處熤须S機(jī)選了3414位進(jìn)行觀察后，從中所保留的90個(gè)觀察對(duì)象的數(shù)據(jù)進(jìn)行分析。（1）進(jìn)行變量，建立變量與的回歸模型的關(guān)系，說(shuō)明教師薪金與哪些變量的關(guān)

4、系密切，是否存在性別和婚姻狀況上的差異。（2）除了變量本身之外，嘗試將它們的平方項(xiàng)或交叉項(xiàng)加到模型中，建立更好的模型三、模型假設(shè)：1、該地區(qū)的人事部門(mén)對(duì)中學(xué)教師的薪金調(diào)查是可信的；2、各參數(shù)對(duì)薪金的影響呈線性關(guān)系；3、工作時(shí)間、性別、教育程度及培訓(xùn)情況之間相互獨(dú)立，沒(méi)有交互作用；四、符號(hào)說(shuō)明：Z:月薪（元）；X1:工作時(shí)間（月）；X2:1男性，女性；X3:1男性或單身女性，0已婚女性；X4:學(xué)歷（取值0-6,值越大表示學(xué)歷越高）；X5:1受雇于重點(diǎn)大學(xué)，0其它；X6:1受過(guò)培訓(xùn)的畢業(yè)生，0未受過(guò)培訓(xùn)的畢業(yè)生或受過(guò)培訓(xùn)的肄業(yè)生；X7:1已兩年以上未從事教學(xué)工作，0其它。五、分析與建立模型：首先，

5、調(diào)用所有相關(guān)變量，運(yùn)用Matlab分別得到，Z與X1-X7之間的關(guān)系及散點(diǎn)圖，由此知Z與各變化量呈線性關(guān)系，于是可以建立線性回歸模型：Z（薪金）為因變量,X1-X7分別表示對(duì)Z的值產(chǎn)生影響的各個(gè)變量，表示回歸系數(shù)，表示隨機(jī)變量.用Matlab求解模型（見(jiàn)附錄），得到的值與置信區(qū)間如下:參數(shù)參數(shù)估計(jì)值置信區(qū)怛01234567000r2=F=P=0表2由上表中R二可以知道薪金（z）的%河由模型確定,由P=0遠(yuǎn)小于，模型從整體上看是成立的，但是還可以看到一個(gè)問(wèn)題，即些部分的置信區(qū)間包含0點(diǎn)，因此我們可以知道這些變量對(duì)因變量的影響是不顯著的.在Matlab中運(yùn)行stepwise命令得到下圖：由圖可知，

6、在模型中對(duì)因變量的影響是不顯著的.于是只保留x1和x4,并將它們的交互項(xiàng)和平方項(xiàng)加入，建立逐步線性回歸方程如下：ln（z）01*x12*x43*x1*x44*x1A25*x4A2然后使用Matlab求解模型（程序見(jiàn)附錄），得到值與置信區(qū)間如下:參數(shù)參數(shù)估計(jì)值置信區(qū)1可012345r2=F=P=0通過(guò)新建模型中得到的數(shù)據(jù),可以看到R2明顯提高,薪金z的可由模型確定.P遠(yuǎn)小于,F遠(yuǎn)超過(guò)臨界值，回歸模型更為顯著，可靠度增高.然后進(jìn)行殘差分析，在Matlab中運(yùn)行命令rcoplot得到殘差圖如下：由圖可知,除個(gè)別數(shù)據(jù)外,其他數(shù)據(jù)的殘差離零點(diǎn)均較近,且殘差的置信區(qū)間都包含零點(diǎn)這說(shuō)明回歸模型能較好地符合原

7、始數(shù)據(jù)，而個(gè)別異常點(diǎn)可以忽略.六、模型的評(píng)價(jià)優(yōu)點(diǎn):1. 該方案實(shí)用簡(jiǎn)單，可行性強(qiáng)，模型簡(jiǎn)單，易于理解。2. 模型一首先用簡(jiǎn)單的線性規(guī)劃進(jìn)行分析.結(jié)構(gòu)簡(jiǎn)單，計(jì)算方便，有利于對(duì)相似問(wèn)題進(jìn)行求解和對(duì)模型進(jìn)行擴(kuò)充。3. 模型二的建立是從一般問(wèn)題到特殊問(wèn)題的發(fā)展過(guò)程.根據(jù)已知的數(shù)據(jù)，從常識(shí)和經(jīng)驗(yàn)進(jìn)行初步分析，并運(yùn)用了逐步線性回歸方法以及輔作散點(diǎn)圖，決定取那幾個(gè)回歸變量及它們的函數(shù)形式.把對(duì)Z影響不顯著的變量（Xi）予以排除，又運(yùn)用殘值分析法建立新的回歸模型.使得精確值增高，模型更合理.缺點(diǎn)：1該模型在處理此問(wèn)題時(shí)有假設(shè)與理想化的思想，與實(shí)際問(wèn)題的求解還有所差距.比如所求模型結(jié)果只達(dá)到了模型設(shè)想的80%左

8、右.七、參考文獻(xiàn)01】趙靜，數(shù)學(xué)建模與數(shù)學(xué)實(shí)驗(yàn)，北京，高等教育出版社，2003【02蘇彥華，MATLABA入門(mén)到精通，北京，人民郵電出版社，2010八、附錄:1、薪金模型數(shù)據(jù)表:編號(hào)ZX1X2X3X4X5X6X71：998：7000001012國(guó)151411000031028181101004：1250：1911000015102819010100610281900000071018270000018國(guó)72：3000000091290301100001011204(30010000111352、1012010121204310001001311104；380000001411184111000

9、0151127420000001612594211010017M127421100001811274200010019109547000001201113國(guó)000001211462520120102211825411000023140415400010101241182540000002515945511211026【1459166000101012712376711010028123767010100291496；750100010；30：14247811010031142479010000321347I9111010:0:3313421920000013413109400010035181

10、410300211036M534M0300000037143010311000038143911111010039194611411311040221611411411041183411411411142p416：117000001143205213911010044208714000211145226415400211146：2201；15811401147299215911511148：1695：1620100010149M792M6711010050169017300000151：1827：1740000011152260417511211053172019901000054M72020

11、9000001015521591209014100561852210010000571210412131101010158185222000000159185222200000060221022211000016122662230100006220272231100006318522270001010164185223200000165199523500000166,6161245113110672324125311010068185225701000169；2054I260000000170【26171284113110711948287110000721720I29001000173260

12、4130811211074185230911010175194231900010076【20271325110000771942326110100781720329110100792048133700000080233434611211181172035500000182：1942：35711000101831211713801100018427423871121118512740：4031121111861942：40611010087226643701000088：243645301000:01892067p580100009020004641121102、與Z的關(guān)系及散點(diǎn)圖：> &

13、gt;x1=7141819191927303030313138414242424247525254545455666767757879919294103103103111114114114117139140154158159162167173174175199209209210213220222222223223227232235245253257260284287290308309319325326329337346355357380387403406437453458464'|> >X1=ones(90,1)x1;> >Z=99810151028125010

14、281028101810721290120413521204110411181127125911271127109511131462118214041182159414591237123714961424142413471342131018141534143014391946221618341416205220872264220129921695179216901827260417201720215918522104185218522210226620271852185219952616232418522054261719481720260418521942202719421720204823

15、341720194221172742274019422266243620672000'> >x2=011100001000010110000100101001010000111110100110100100001001010001100110110111010111110001'>>X2=ones(90,1)x2;>>Z=9981015102812501028102810181072129012041352120411041118112712591127112710951113146211821404118215941459123712371

16、4961424142413471342131018141534143014391946221618341416205220872264220129921695179216901827260417201720215918522104185218522210226620271852185219952616232418522054261719481720260418521942202719421720204823341720194221172742274019422266243620672000'x3=011110001110010110001100101111110000111110100

17、111100110111001110001110111110111010111111111x3=011110001110010110001100101111110000111110100111100110111001110001110111110111010111111111'>>X3=ones(90,1)x3;>>Z=998101510281250102810281018107212901204135212041104111811271259112711271095111314621182140411821594145912371237149614241424

18、13471342131018141534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 17201720 2159 1852 2104 1852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 19481720 2604 1852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 20672000'> >

19、;x4=000000000020000000002000200000000020003440022450000200400000000003000300200000020002200002'> >X4=ones(90,1)x4;> >Z=998101510281250102810281018107212901204135212041104111811271259112711271095111314621182140411821594145912371237149614241424134713421310181415341430143919462216183414

20、16205220872264220129921695179216901827260417201720215918522104185218522210226620271852185219952616232418522054261719481720260418521942202719421720204823341720194221172742274019422266243620672000'x5=001010000001000101000010111101010110011110111010100100101000001001100100111011010001110001'>

21、; >X5=ones(90,1)x5;> >Z=998101510281250102810281018107212901204135212041104111811271259112711271095111314621182140411821594145912371237149614241424134713421310181415341430143919462216183414162052208722642201299216951792169018272604172017202159185221041852185222102266202718521852199526162324

22、18522054261719481720260418521942202719421720204823341720194221172742274019422266243620672000'> >x6=000000000010000000001000100000000010001110011110000100000000000001000100100000010001100001'> >X6=ones(90,1)x6;> >Z=998101510281250102810281018107212901204135212041104111811271

23、25911271127109511131462118214041182159414591237123714961424142413471342131018141534143014391946221618341416205220872264220129921695179216901827260417201720215918522104185218522210226620271852185219952616232418522054261719481720260418521942202719421720204823341720194221172742274019422266243620672000'> >x7=000000100000000000110000000000001000000011011110011000000100000110010001010000011011100000'> >X7=ones(90,1)x7;> >Z=99810151028125010281028101810721290120413521204110