数学建模,教师的薪金

1、昆明理工大学第六届大学生数学建模竞赛承诺书我们仔细阅读了昆明理工大学大学生数学建模竞赛的竞赛规则。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的。如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。评阅编号（由组委会评阅前进行编号）：昆明理工大学第六届大学生数学建模竞赛评阅专用页评阅编号（由组委会评阅前进行

2、编号）：评阅记录（供评阅时使用）：评阅人评分备注总分大学教师薪金模型一、摘要：模型一：利用Matlab建立X1-X7与Z（新金）的线性关系，得到散点图与一般的线性回归模型（模型一），模型如下：模型二：由模型一整体效果与StepwiseTable图推断部分变量对Z（薪金）的影响并不显著，由残差分析法筛选出影响明显的变量X1、X4，将他们的平方项与交互项加入建立新的回归模型（模型二），模型如下：z6.90260.0043x10.1746x40.0001x1*x40*x1A20.0228x4A2得到R2,F,P,与模型一相比较。模型一调用原始参数较全，但回归性差，模型二只启用影响明显的变量，回归模型

3、更为显著，可靠度更高。于是得出：教师的薪金与工作时间、学历关系明显，与性别、是否受雇于重点大学、是否接受过培训的关系较小，即女教师没有受到不公正的待遇，婚姻状况也不会影响收入。关键词：回归分析，互交作用，图形结合，残值分析法。二、问题重述：某地人事部门为研究中学教师的薪金与他们的资历、性别、教育程度及培训情况等因素之间的关系，要建立一个数学模型，分析人事策略的合理性，特别是考察女教师是否受到不公平的待遇，以及她们的婚姻状况是否会影响收入。要求从当地教师中随机选了3414位进行观察后，从中所保留的90个观察对象的数据进行分析。（1）进行变量，建立变量与的回归模型的关系，说明教师薪金与哪些变量的关

4、系密切，是否存在性别和婚姻状况上的差异。（2）除了变量本身之外，尝试将它们的平方项或交叉项加到模型中，建立更好的模型三、模型假设：1、该地区的人事部门对中学教师的薪金调查是可信的；2、各参数对薪金的影响呈线性关系；3、工作时间、性别、教育程度及培训情况之间相互独立，没有交互作用；四、符号说明：Z:月薪（元）；X1:工作时间（月）；X2:1男性，女性；X3:1男性或单身女性，0已婚女性；X4:学历（取值0-6,值越大表示学历越高）；X5:1受雇于重点大学，0其它；X6:1受过培训的毕业生，0未受过培训的毕业生或受过培训的肄业生；X7:1已两年以上未从事教学工作，0其它。五、分析与建立模型：首先，

5、调用所有相关变量，运用Matlab分别得到，Z与X1-X7之间的关系及散点图，由此知Z与各变化量呈线性关系，于是可以建立线性回归模型：Z（薪金）为因变量,X1-X7分别表示对Z的值产生影响的各个变量，表示回归系数，表示随机变量.用Matlab求解模型（见附录），得到的值与置信区间如下:参数参数估计值置信区怛01234567000r2=F=P=0表2由上表中R二可以知道薪金（z）的%河由模型确定,由P=0远小于，模型从整体上看是成立的，但是还可以看到一个问题，即些部分的置信区间包含0点，因此我们可以知道这些变量对因变量的影响是不显著的.在Matlab中运行stepwise命令得到下图：由图可知，

6、在模型中对因变量的影响是不显著的.于是只保留x1和x4,并将它们的交互项和平方项加入，建立逐步线性回归方程如下：ln（z）01*x12*x43*x1*x44*x1A25*x4A2然后使用Matlab求解模型（程序见附录），得到值与置信区间如下:参数参数估计值置信区1可012345r2=F=P=0通过新建模型中得到的数据,可以看到R2明显提高,薪金z的可由模型确定.P远小于,F远超过临界值，回归模型更为显著，可靠度增高.然后进行残差分析，在Matlab中运行命令rcoplot得到残差图如下：由图可知,除个别数据外,其他数据的残差离零点均较近,且残差的置信区间都包含零点这说明回归模型能较好地符合原

7、始数据，而个别异常点可以忽略.六、模型的评价优点:1. 该方案实用简单，可行性强，模型简单，易于理解。2. 模型一首先用简单的线性规划进行分析.结构简单，计算方便，有利于对相似问题进行求解和对模型进行扩充。3. 模型二的建立是从一般问题到特殊问题的发展过程.根据已知的数据，从常识和经验进行初步分析，并运用了逐步线性回归方法以及辅作散点图，决定取那几个回归变量及它们的函数形式.把对Z影响不显著的变量（Xi）予以排除，又运用残值分析法建立新的回归模型.使得精确值增高，模型更合理.缺点：1该模型在处理此问题时有假设与理想化的思想，与实际问题的求解还有所差距.比如所求模型结果只达到了模型设想的80%左

8、右.七、参考文献01】赵静，数学建模与数学实验，北京，高等教育出版社，2003【02苏彦华，MATLABA入门到精通，北京，人民邮电出版社，2010八、附录:1、薪金模型数据表:编号ZX1X2X3X4X5X6X71：998：7000001012国151411000031028181101004：1250：1911000015102819010100610281900000071018270000018国72：3000000091290301100001011204(30010000111352、1012010121204310001001311104；380000001411184111000

9、0151127420000001612594211010017M127421100001811274200010019109547000001201113国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的关系及散点图：> &

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