数学建模,教师的薪金

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

2、会评阅前进行编号）：评阅记录（供评阅时使用）：评阅人评分备注总分大学教师薪金模型一、摘要：模型一：利用Matlab建立X1-X7与Z（薪金）的线性关系，得到散点图与一般的线性回归模型（模型一），模型如下：模型二：由模型一整体效果与Stepwise Table图推断部分变量对Z（薪金）的影响并不显著，由残差分析法筛选出影响明显的变量X1、X4，将他们的平方项与交互项加入建立新的回归模型（模型二），模型如下：得到R2，F,P，与模型一相比较。模型一调用原始参数较全，但回归性差，模型二只启用影响明显的变量，回归模型更为显著，可靠度更高。于是得出：教师的薪金与工作时间、学历关系明显，与性别、是否受雇于

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

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

5、点图, 由此知Z与各变化量呈线性关系，于是可以建立线性回归模型： Z（薪金）为因变量,X1-X7分别表示对Z的值产生影响的各个变量,表示回归系数,表示随机变量. 用Matlab求解模型（见附录），得到的值与置信区间如下：参数参数估计值置信区间1.13111.0268 1.23530.00270.0023 0.0031-0.0229-0.1432 0.09740.0094-0.1005 0.11930.10890.0296 0.18820.0385-0.0670 0.14400.1817-0.0507 0.414200 0=0.7889F= 51.6934 P=0表一 由上表中=0.7889可以

6、知道薪金（）的78.89%可由模型确定,由P=0远小于0.05,模型从整体上看是成立的,但是还可以看到一个问题，即些部分的置信区间包含0点,因此我们可以知道这些变量对因变量的影响是不显著的.在Matlab中运行stepwise命令得到下图： 由图可知，在模型中对因变量的影响是不显著的.于是只保留，并将它们的交互项和平方项加入，建立逐步线性回归方程如下： 然后使用Matlab求解模型（程序见附录），得到值与置信区间如下：参数参数估计值置信区间6.90266.8557 6.94960.00430.0037 0.00490.17460.1073 0.2419-0.0001-0.0003 0.0001

7、-0.0000-0.0000 -0.0000-0.0228-0.0371 -0.0085= 0.9008F=152.6081P=0通过新建模型中得到的数据,可以看到明显提高,薪金的90.08%可由模型确定.远小于0.05,F远超过临界值，回归模型更为显著，可靠度增高.然后进行残差分析，在Matlab中运行命令rcoplot得到残差图如下：由图可知,除个别数据外,其他数据的残差离零点均较近,且残差的置信区间都包含零点.这说明回归模型能较好地符合原始数据，而个别异常点可以忽略.六、模型的评价优点:1. 该方案实用简单，可行性强，模型简单，易于理解。2. 模型一首先用简单的线性规划进行分析.结构简单

8、，计算方便，有利于对相似问题进行求解和对模型进行扩充。3. 模型二的建立是从一般问题到特殊问题的发展过程. 根据已知的数据，从常识和经验进行初步分析，并运用了逐步线性回归方法以及辅作散点图，决定取那几个回归变量及它们的函数形式.把对影响不显著的变量（）予以排除，又运用残值分析法建立新的回归模型.使得精确值增高，模型更合理.缺点： 1 该模型在处理此问题时有假设与理想化的思想，与实际问题的求解还有所差距.比如所求模型结果只达到了模型设想的80%左右.七、参考文献【01】赵静，数学建模与数学实验，北京，高等教育出版社，2003【02】苏彦华，MATLAB 7.0从入门到精通，北京，人民邮电出版社，

9、2010八、附录：1、 薪金模型数据表：编号ZX1X2X3X4X5X6X719987000000210151411000031028181101004125019110000510281901010061028190000007101827000001810723000000091290301100001012043001000011135231012010121204310001001311043800000014111841110000151127420000001612594211010017112742110000181127420001001910954700000120111352

10、000001211462520120102211825411000023140454000100241182540000002515945511211026145966000100271237671101002812376701010029149675010000301424781101003114247901000032134791110100331342920000013413109400010035181410300211036153410300000037143010311000038143911111010039194611411311040221611411411041183411

11、411411142141611700000143205213911010044208714000211145226415400211146220115811401147299215911511148169516201000049179216711010050169017300000151182717400000152260417511211053172019901000054172020900000055215920901410056185221001000057210421311010058185222000000159185222200000060221022211000061226622

12、301000062202722311000063185222700010064185223200000165199523500000166261624511311067232425311010068185225701000169205426000000070261728411311071194828711000072172029001000173260430811211074185230911010175194231900010076202732511000077194232611010078172032911010079204833700000080233434611211181172035

13、50000018219423571100008321173801100018427423871121118527404031121118619424061101008722664370100008824364530100008920674580100009020004641121102、 与Z的关系及散点图：>> x1=7 14 18 19 19 19 27 30 30 30 31 31 38 41 42 42 42 42 47 52 52 54 54 54 55 66 67 67 75 78 79 91 92 94 103 103 103 111 114 114 114 117

14、139 140 154 158 159 162 167 173 174 175 199 209 209 210 213 220 222 222 223 223 227 232 235 245 253 257 260 284 287 290 308 309 319 325 326 329 337 346 355 357 380 387 403 406 437 453 458 464' |>> X1=ones(90,1) x1;>> Z=998 1015 1028 1250 1028 1028 1018 1072 1290 1204 1352 1204 1104 1

15、118 1127 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 1347 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 1720 1720 2159 1852 2104 1852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 1948 1720 2604 1

16、852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'>> x2=0 1 1 1 0 0 0 0 1 0 0 0 0 1 0 1 1 0 0 0 0 1 0 0 1 0 1 0 0 1 0 1 0 0 0 0 1 1 1 1 1 0 1 0 0 1 1 0 1 0 0 1 0 0 0 0 1 0 0 1 0 1 0 0 0 1 1 0 0 1 1 0 1 1 0 1 1 1 0 1 0 1 1 1 1 1 0 0 0 1'>> X2=ones

17、(90,1) x2;>> Z=998 1015 1028 1250 1028 1028 1018 1072 1290 1204 1352 1204 1104 1118 1127 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 1347 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 1720 1720 2159 1852

18、2104 1852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 1948 1720 2604 1852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'x3=0 1 1 1 1 0 0 0 1 1 1 0 0 1 0 1 1 0 0 0 1 1 0 0 1 0 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 0 1 0 0 1 1 1 1 0 0 1 1 0 1 1 1 0 0 1 1 1

19、 0 0 0 1 1 1 0 1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 1 1 1 1 1 1x3=0 1 1 1 1 0 0 0 1 1 1 0 0 1 0 1 1 0 0 0 1 1 0 0 1 0 1 1 1 1 1 1 0 0 0 0 1 1 1 1 1 0 1 0 0 1 1 1 1 0 0 1 1 0 1 1 1 0 0 1 1 1 0 0 0 1 1 1 0 1 1 1 1 1 0 1 1 1 0 1 0 1 1 1 1 1 1 1 1 1'>> X3=ones(90,1) x3;>> Z=998 1015 1028 1250 1

20、028 1028 1018 1072 1290 1204 1352 1204 1104 1118 1127 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 1347 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 1720 1720 2159 1852 2104 1852 1852 2210 2266 2027 1852 1852 1

21、995 2616 2324 1852 2054 2617 1948 1720 2604 1852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'>> x4=0 0 0 0 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 0 2 0 0 0 2 0 0 0 0 0 0 0 0 0 2 0 0 0 3 4 4 0 0 2 2 4 5 0 0 0 0 2 0 0 4 0 0 0 0 0 0 0 0 0 0 3 0 0 0 3 0 0 2 0 0 0 0 0

22、 0 2 0 0 0 2 2 0 0 0 0 2'>> X4=ones(90,1) x4;>> Z=998 1015 1028 1250 1028 1028 1018 1072 1290 1204 1352 1204 1104 1118 1127 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 1347 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992

23、1695 1792 1690 1827 2604 1720 1720 2159 1852 2104 1852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 1948 1720 2604 1852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'x5=0 0 1 0 1 0 0 0 0 0 0 1 0 0 0 1 0 1 0 0 0 0 1 0 1 1 1 1 0 1 0 1 0 1 1 0 0 1 1

24、1 1 0 1 1 1 0 1 0 1 0 0 1 0 0 1 0 1 0 0 0 0 0 1 0 0 1 1 0 0 1 0 0 1 1 1 0 1 1 0 1 0 0 0 1 1 1 0 0 0 1'>> X5=ones(90,1) x5;>> Z=998 1015 1028 1250 1028 1028 1018 1072 1290 1204 1352 1204 1104 1118 1127 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 134

25、7 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 1720 1720 2159 1852 2104 1852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 1948 1720 2604 1852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'>

26、> x6=0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 1 0 0 1 1 1 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 1 0 0 0 0 0 0 1 0 0 0 1 1 0 0 0 0 1'>> X6=ones(90,1) x6;>> Z=998 1015 1028 1250 1028 1028 1018 1072 1290 1204 1352 1204 1104 1118 11

27、27 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 1347 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 1720 1720 2159 1852 2104 1852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 1948 1720 2604 1852 19

28、42 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'>> x7=0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 1 0 1 1 1 1 0 0 1 1 0 0 0 0 0 0 1 0 0 0 0 0 1 1 0 0 1 0 0 0 1 0 1 0 0 0 0 0 1 1 0 1 1 1 0 0 0 0 0'>> X7=ones(90,1)

29、 x7;>> Z=998 1015 1028 1250 1028 1028 1018 1072 1290 1204 1352 1204 1104 1118 1127 1259 1127 1127 1095 1113 1462 1182 1404 1182 1594 1459 1237 1237 1496 1424 1424 1347 1342 1310 1814 1534 1430 1439 1946 2216 1834 1416 2052 2087 2264 2201 2992 1695 1792 1690 1827 2604 1720 1720 2159 1852 2104 1

30、852 1852 2210 2266 2027 1852 1852 1995 2616 2324 1852 2054 2617 1948 1720 2604 1852 1942 2027 1942 1720 2048 2334 1720 1942 2117 2742 2740 1942 2266 2436 2067 2000'3、 模型一程序：>> x1;>> x2;>> x3;>> x4;>> x5;>> x6;>> x7;>> x=x1,x2,x3,x4,x5,x6,x7;>>

31、; stepwise(x,z)>> b,bint,r,rint,stats=regress(z,x)Warning: R-square and the F statistic are not well-defined unless X has acolumn of ones.Type "help regress" for more information.> In regress at 162b = 1.0e+003 * 1.1311 0.0027 -0.0229 -0.0229 0.1089 0.0385 0.1817bint = 1.0e+003 *

32、1.0268 1.2353 0.0023 0.0031 -0.1432 0.0974-0.1005 0.1193 0.0296 0.1882-0.0670 0.1440-0.0507 0.4142 0.0000 0.0000r = 1.0e+003 * 1.6545 0.8444 0.1855 0.8335 -0.1303 1.6284 0.8731 1.6045 0.8096 0.4959 0.7509 0.9522 1.5871 0.7857 1.5784 0.1334 0.7835 0.9283 0.8296 0.8187 0.7053 0.7575 0.9022 1.5524 0.36

33、23 0.8762 0.0791 -0.2346 0.3981 0.0552 0.3894 0.0269 0.7318 0.8153 1.0530 1.4459 0.6510 -0.0165 0.3433 0.4525 -0.2855 0.6775 -0.0773 0.2346 0.2042 0.2691 -0.2741 0.2091 -0.1382 0.5559 0.5537 0.1016 0.1287 1.2156 -0.1065 0.1048 -0.2381 0.4538 1.1874 0.3925 0.0766 0.3903 0.5264 0.4277 0.4212 0.0587 -0

34、.3250 -0.7352 1.1048 -0.0260 0.2513 -0.8069 -0.1873 -1.1846 0.3265 0.1687 -0.4836 -0.4901 0.9376 -1.0078 0.1605 0.0992 -0.6887 -1.0969 -1.1316 -0.6574 -0.3883 -0.4231 -0.4339 -0.5262rint = 1.0e+003 * 0.1748 3.1341 -0.6232 2.3119 -1.2830 1.6541 -0.6350 2.3020 -1.5468 1.2861 0.1476 3.1092 -0.5915 2.33

35、76 0.1227 3.0863 -0.6609 2.2801 -0.9648 1.9566 -0.6306 2.1323 -0.5146 2.4190 0.1047 3.0695 -0.6867 2.2581 0.0957 3.0611 -1.3399 1.6067 -0.6890 2.2561 -0.5395 2.3961 -0.6383 2.2975 -0.6500 2.2874 -0.6803 2.0908 -0.7170 2.2319 -0.5665 2.3710 0.0690 3.0358 -1.0889 1.8135 -0.5934 2.3457 -1.3984 1.5565 -

36、1.6607 1.1915 -1.0699 1.8662 -1.4238 1.5342 -1.0792 1.8580 -1.4538 1.5077 -0.7419 2.2056 -0.6554 2.2860 -0.3692 2.4751 -0.0385 2.9304 -0.8291 2.1311 -1.4994 1.4664 -1.1231 1.8097 -0.9679 1.8728 -1.6954 1.1244 -0.7984 2.1535 -1.5624 1.4078 -1.1914 1.6607 -1.2227 1.6311 -1.1037 1.6418 -1.5958 1.0476 -

37、1.2660 1.6842 -1.6245 1.3482 -0.9221 2.0338 -0.9243 2.0317 -1.3617 1.5650 -1.3465 1.6040 -0.2610 2.6923 -1.1202 0.9073 -1.3701 1.5798 -1.7244 1.2482 -1.0229 1.9305 -0.2873 2.6621 -1.0892 1.8742 -1.3979 1.5511 -1.0913 1.8719 -0.9364 1.9892 -1.0482 1.9036 -1.0545 1.8969 -1.4065 1.5240 -1.8092 1.1592 -2.1647 0.6943 -0.3630 2.5727 -1.4867 1.4347 -1.2239 1.7264 -2.2359 0.622