数学建模-大学生就业问题

1、'.2010-2011 第二学期数学建模课程设计2011年6月27日7月1日题目大学生就业问题第 11组组员 1组员 2组员 3组员 4姓名学号 08080602170808060218 0808060219 0808060220专业信计 0802信计 0802信计 0802信计 0802成绩;.'.论文摘要本文讨论了在新的形势下大学生的就业问题。20 世纪 90 年代以来，我国出现了一种前所未有的现象， 有着“天之骄子” 美誉的大学生也开始面临失业问题。大学生就业难问题已受到普遍关注。 大学生毕业失业群体正在不断扩大， 已成为我国扩大社会就业，构建和谐稳定社会的急需解决的社会

2、问题。本文针对我国现有的国情， 综合考虑了高校毕业生的就业率和高校招生规模的扩大之间的关系， 建立了定量分析的微分方程模型，随后又建立了了离散正交曲线拟合模型对得出的结果进行了检验， 并分析模型得出的结果得合理性。 最终得到生源数量与失业率之间的拟合多项式和拟合曲线， 并预测出了未来高校招生规模的变化趋势。在找到大学生失业规律以后， 本文还具体的对毕业生的性别、 出生地对失业的影响做出了定量分析。关键词： 大学生就业微分方程模型多项式曲线拟合MATLAB 软件1、问题重述大学生就业问题 ：如果我们将每年毕业的大学生中既没有找到工作又没有继续深造的情况视为失业， 就可以用失业率来反映大学生就业的

3、状况。 下面的表中给出了某城市的大学生失业数占城市总失业人数的比率， 比率的计算是按照国际劳工组织的定义，对 16 岁以上失业人员进行统计的结果。表 1年份失业率（ %）男性女性出生于城市出生于农村出生于城市出生于农村19894.46.05.710.819905.113.24.71119919.717.28.715.819921020817199312221019199411.32211.2211995112210.419199611.12110.220199710.7189.91619988.6148.213;.'.19996.9136.71120006.7116.31120018.

4、7166.913200211207.014200311.7247.218200410.6256.91620059.4206.31720069.1195.91520077.6165.41420086.314.14.912.620096.413.34.712.8请建立相应的模型对大学生就业状况进行分析找出其中的规律并讨论下面两个问题：(1) 、就业中是否存在性别歧视；(2) 、学生的出生对就业是否有影响。2、模型假设2.1 在本次研究中做出以下假设：(1)、假设毕业生求职时竞争是公平的；(2)、假设考研等继续深造的毕业生属于已就业人群；(3)、假设每个毕业生都有就业或者继续深造的意图(4)、假设就

5、业率和失业率之和为1；(5)、假设本文搜集的数据全部真实可靠；2.2 在定量分析性别、出生地对失业的影响时还要做以下假设：(1) 、假设毕业生就业情况只受性别、出生地等因素的影响；(2) 、假设具有上述同等条件的毕业生间就业机会相同(3) 、假设附件中的数据信息均合理；3、问题分析3.1 对问题的分析若要分析新失业群体产生的主要原因，并就其重要性给出各种因素的排序，就需要对搜集的数据进行整理，并进行系统的分析，划分为不同的体系和矛盾，然后我们考虑用Logistic模型分析。为了得到新失业群体对高校招生生源的影响和预测未来高校招生规模的变;.'.化趋势，我们考虑建立基于微分方程模型和离散

6、正交曲线拟合模型来进行求解，并将结果进行比较。3.2对大学生失业群体产生的客观原因分析及其重要性排序影响高校毕业生就业的主要因素的选取的基础是因素指标体系。 由于影响高校毕业生就业的因素有很多方面， 而且有些因素具有多方面表现特征， 因此对其进行描述，必须借助因素和因素体系。 根据科学性原则和定量与定性相结合的原则，我们将影响高校毕业生就业的因素分为高校扩招造成大学生就业市场供需失衡因素、教育结构与产业结构不协调因素体系以及高校毕业生的自身因素体系。教育结构与产业结构不协调因素体系是由于我国高等教育生产规模不断扩大，但不符合教育规律的压缩式追求低成本造成教育资源结构失调，表现在两方面：(1)、

7、学科与专业的结构失衡，即学科结构和专业结构与社会经济结构不相适应。表现为：有些学科和专业人才的过度教育；有些学科和专业人才的教育资源配置不足目前我国高等教育的学科与专业设置基本上由各个学校自主决定， 学校出于经济效益和学校发展需要， 主要考虑的因素是市场行情， 而市场调节的盲目性和滞后性使得不同专业稀缺程度的变化与经济结构转换的需求变化不同步，致使许多热门应用学科类专业低水平重复现象严重， 并冲击专业教学的质量， 造成这些专业的毕业生结构性过剩。 从学科专业设置的历史沿革看， 许多高校的学科专业建设更多是以学科自身内在逻辑的发展为依据和基础， 而较少参照现实社会经济建设和发展领域的需要。 当社

8、会经济结构特别是产业结构的变化迅速， 高等教育结构未能迅速地做出相应的调整， 其培养的人才则不能满足社会和劳动力市场的需要，因此产生了知识性失业与职位空缺的矛盾。(2)、层次结构失衡，即高等学校的学历学位教育层次比例及构成与经济社会发展 的需求结构不匹配，出现了职业刚性失业和职位空缺并存的现象。在我国目前的高等人才培养中，社会需求存在对学历教育的“符号效应” ，即对能力的需求不如对学历的要求。 而从需求看， 社会对各级人才的需求结构呈 “金字塔型”。社会不仅需要从事高深学问研究和创造发明的学术性人才，更需要把现有科技转化为现实生产力和产品的大批熟练工人和技术人员。 因此，教育供给应该在专科、本

9、科和研究生教育的数量上形成合理的比例。 然而，整个教育层次的扩招和缺乏鲜明特色的学科内容， 不仅使得这种比例开始失调， 而且造成各层次的毕业生不能达到应有的质量要求。同时，替代性增加，引发了就业市场中的“挤占效应”。 即博士研究生挤占硕士研究生位置， 硕士研究生挤占本科生位置， 本科生挤占专科生位置，专科生又因为其质量不高、特色不鲜明而找不到位置。4、模型建立4.1 微分方程模型;.'.因为高校毕业生的就业率和招生人数都可视为随时间动态变化， 所以我们考虑通过建立微分方程模型去认识和解决有关毕业生就业和计划招生规模的实际问题。为此，我们做出如下的基本假设：(1)高校毕业生就业人数的变化

10、率与毕业生的综合素质 (如品学表现 )和社会的需求呈正的线性相关。(2)部分毕业生的主客观原因(比如，没有顺利完成学业， 或者想继续报考研究生，或者就业意识淡薄，就业观念差，对自己估计不足等 ) 影响了自身的就业，因而对毕业生的就业产生了阻滞作用。(3)高校当年的计划招生人数与毕业生总人数成正比，比例系数为c 。其中假设 (2)借鉴了人口增长阻滞模型 中的“阻滞”的思想。我们引入如下符号：N (t) ：时刻 t 高校毕业生的总人数，N (t )0 ；M (t ) ：时刻 t 高校计划招生的总人数，M (t)0 ；r (t ) ：时刻 t 毕业生的就业率 ( 即：就业人数 / 毕业总人数 ) ，

11、 0r (t)1 ；R(t ) ：时刻 t 社会对于毕业生的需求率( 即：需求人数 / 毕业生总人数 ) 。记 r (t 0)r 0 与 M (t0 )M 0 分别为时刻 t 0 的高校毕业生就业率与高校招生人数。很明显，对于需求率R(t ) 而言，我们有：当R(t )1 时， 毕业生供大于求；当 R(t )1 时，毕业生供求平衡；当R(t)1 时，毕业生供不应求。由于社会的劳动力需求是与国家的经济运行情况正相关的，故我们这里的需求率R(t ) 还反映了社会经济发展的GDP 速度。4.2基本微分方程模型的建立首先，根据模型的假设 (1)和(2)，我们有：drNR(t )N(1 r ) NrNd

12、t（ 1） R(t)() rN其中的比例系数，与分别与需求人数，就业人数和未就业人数有关，故分别称为需求因子，就业因子与阻滞因子。在本文中，我们均假设0 。其次，根据模型的假设 (3)，将方程 (1)的两边同乘系数 c ，我们得到;.'.drNR(t )M(1 r )MrMdt（ 2） R(t) ()r M于是，根据方程 (1) ，当毕业生的总人数为常数时，我们得到高校毕业生的就业率满足一阶线性微分方程模型：dr)rR(t)(dt（ 3）r (t 0)r 0最后，又根据方程 (2)，当毕业生的就业率为常数时，我们得到高校计划招生的总人数满足一阶线性微分方程模型：dM1 R(t )dtr

13、M (t0 )M 05、模型求解为方便模型 (3) 和 (4) 的求解，我们假设在模型 (3) 和(4) 中社会对于毕业生的需求率是常数 ( 此时记为 R)。5.1研究毕业生的失业率模型方程 (3) 有解 ( 即失业率 ) 为：atr t 1 er 0 r c r c其中 r caRa是模型 (3) 的不稳定平衡点 ( 失业率 ) 。我们有以下结论：当 rc0 时，这表明：只要影响毕业生失业的因素较大( 或者社会对毕业生的需求量较小 ) ，就存在着不稳定的毕业生失业率；当 r c 1 时，这表明：只要影响毕业生失 2业的因素非常大，就会出现不稳定的低失业率。当 r c r 0 时，这表明：只要

14、不稳定的失业率低于初始的失业率， 就有毕业生的失业率超过不稳定的失业率r c 。进一步可知，毕业生的失业率达到0所需要的时间为：t1ln r car c r 0(5)(6);.'.故当 R( 供不应求 ) 时，有 t0，这表明：需求率越大，达到低失业率的时间越短； 当0时，有 t1lnRaR r 0a ，由此可见，即使阻滞因子很小，达到低失业率也需要一定的时间。5.2研究高校的招生规模我们有方程 (4) 的解 ( 即高校计划招生的总人数) 为：aM t e r r r c t M 0因此，我们有以下结论：当 rrc 时，有 M (t )M (t1)，表明：高校招生总人数规模宜降低；当

15、rrc 时，有 M (t )M (t1)，表明：高校招生总人数规模宜保持不变；当 rrc 时，有 M (t)M (t1) ，表明：高校招生总人数规模宜扩大。 当 R时，有 r c0 和 M (t ) M (t 1)，此时高校招生总人数规模宜扩大。总之，高校应当按照毕业生的失业率或者社会对于毕业生的失业率去确定其计划招生的规模。6、模型分析及改进6.1模型分析如果在模型的假设 (3)的基础上，将高校招生人数的相对变化率按照毕业生的失业率去进行调整，即dMr M ，其中的比例系数r 为 r 的函数，那么dt就有dNr N ，再由 (1)知高校毕业生的失业率满足一阶非线性微分方程模dt型：dr ar

16、 rRt dtrt 0r 0特别地，当取 R tR0 和r0时，只要0a， R0，方程 (8)就有惟一稳定的平衡点 ( 失业率 )0aR01，且 rs0 的充要条件是：rs0aR00a 。这表明：在招生规模扩大和需求率较大的条件下，将会得到稳定的失业率。(7)(8);.'.如果进一步将模型的假设（3）改进为高校当年计划招生的人数相对于毕业生人数的变化率，按照其毕业生的就业率去进行调整，则有dM1 r N ，其中dt的比例系数1r为 r的函数，并且还注意到毕业生的人数为年前 ( 通常取4 的招生人数，即 NtMt，从而就有高校计划招生的总人数满足一阶线性时滞微分方程模型：dM1 rM t

17、dt(9)M t0M 0d NdM又由 dtdt 和 (1)知，高校毕业生的失业率仍满足模型(8) ，只不过r1 r 。对于1r10 和 R tR0 均为常数的情形。由于微分方程(6-2-9) 的 特征 方程e1 有无 穷 多 个根 ，故 我们 可 利用 Taylor 展开e 1将其改写成近似的二次方程11 ，在舍去负根后，得到正根121，从而就有：14111min1 ,11(10)以及方程 (9)的特解 M tM 0e tM 0 e 1t，这表明高校扩大招生人数规模的实际增长率 应低于相对变化率1。又由于有R0R01aa故如果按照实际的去扩大招生规模 ( 招生人数为： M 0 e t) 将会

18、有利于降低模型 (8)的稳定失业率。随着我国高等教育质量的不断提高， 经济发展新情况的出现和国家就业政策的调整，社会对毕业生的需求将发生较大的波动。如果在模型的假设（3）的基础上，将高校招生人数的相对变化率按照社会对毕业生的需求率去进行调整，那么就有：dM2 R M ，其中的比例系数2 R 为 R 的函数。再由 (1) 知高校毕dt业生的就业率满足一阶非线性微分方程模型：;.'.dra2R rRdtrt0r0特别地，当取 R tR0，2RR00时，只要 R,R0 a，方程(11) 就有稳定的失业率rsR00，且出现低失业率 rs0的充要条件是：R0a需求率 R0a ，。而当时，只要 R

19、0，就有 rs1 ，这表明当需求率充分大时，可得到稳定的低失业率。如果将模型假设 （3） 改进为高校当年计划招生人数相对于毕业生人数的变化率按照社会对毕业生的需求率去进行调整，即有dM3R N ，其中比例系数dt3 R 为 R 的函数，且注意到毕业生的人数为年前 (4) 的招生人数，即N t M t，那么可得到高校计划招生的总人数满足一阶线性时滞微分方程模型：dM3R MtdtM t 0M 0对方程 (12) 的讨论类似于 (9) ，结果是：对于 R tR0和 3R3 为常数的情形，模型方程 (12) 的特征根满足(11)(12)1 min ,31(13)这表明扩大招生人数规模的实际增长率应低

20、于相对变化率3 ，而且按照实际增长率去扩大招生规模 ( 招生人数为： M 0 e t ) 也必将会有利于降低 (11) 的稳定失业率。最后，我们研究高校毕业生就业率和招生人数规模数学模型的一般化。如果将模型的假设（ 3）改进为高校当年计划招生人数相对于毕业生人数M t 的变化率综合按照当年的就业率 r t 和未来年( 如取4 ) 社会对于毕业生的需求率 R t去进行调整， 即： dMrt , R tN，其中比例系数r , R 为 rdt;.'.和 R的函数，且注意到毕业生的人数为年前的招生人数，即N tMt，那么我们有高校计划招生的总人数满足一阶线性时滞微分方程模型：dMM tr t

21、, R tdt(14)M t0 M 0又由 dNdM 和 (1)知，高校毕业生的就业率满足一阶非线性时滞微分方程模dtdtdrr , R trR ta型：dt(15)r r0 r06.2 定量分析在上文中已经建立和定性分析了关于高校毕业生就业率的如下微分方程模型drR(t )()r（ 1）dtr (t0 )r 0其中 r (t ) 表示 t 时刻毕业生的就业率 (即:就业人数 /毕业生人数 , 0r (t) 1 )；R(t) 表示 t 时刻社会对于毕业生的需求率(即:需求人数 /毕业生人数 )。因为在实际调查中需求率是难以完整统计的 , 所以我们应根据这种实际情况将模型 (1) 进行相应修改

22、, 不把需求率考虑进去 , 于是我们有下面的微分方程模型：dr)r u(dt（ 2）r (t0 )r 0其中比例系数与分别与就业人数和未就业人数有关。以下均假设0 。在本文中 ,我们的主要目的是通过调查实际数据计算出参数、，从而给出与（ 2）相应的微分方程定量模型，并做出定量的分析。参数、的计算方法如下:设高校招生人数，毕业生人数与就业率的数据如表所示。;.'.表 6.2.1 历年全国高校招生人数和失业率年份（ t ）招生人数 / 万人 (M）毕业人数 / 万人 (N）就业率 /%(r ）2001268.281140.92002320.511450.82003382.17212.20.

23、72004447.342540.732005504.463260.7262006548.584130.722007567.364790.7120085995590.720096296110.68（该数据来自中华人民共和国人力资源和社会保障部）将 dr()ru 进行差商近似处理，再根据最优平方逼近法我们得到dtr i1r i)r iti(1ti为确定参数、，不妨设x,y 则方程（ ）化为：3yr i1r i0rixt iti11n( y rixr i1r i) ，将数据代入s0,s记 s(x, y)titix0 可得方程ni 11y1n( y ri 1 xri1ri)0titin i 111nr

24、r i1r i)(r)0( yi1xintiti1i 11由方程（ 5）便可解出 x, y ,从而确定出参数、，得到相应的微分方程定量模型。我们根据参考文献 13得到参数为：0.25,0.1进而微分方程化为dr(0.250.35)r0.25dtr (2001)0.7（ 3）（ 4）（ 5）;.'.dM0.25(0.35) MdtrM (t0)480求的解析解为：r (t)0.012e0.35 t0.1740.350.714) t(M (t) eM 06.3 离散正交曲线拟合模型离散正交多项式定义 1 如果两个多项式 P( x) ， Q(x) 满足mP(xi)Q (xi )0i0则称 P

25、( x) 与 Q ( x) 在点集 x , x , x 上是离散正交的。设 P( x), P ( x), P ( x)0 1m01n为多项式， Pk 为 k 次多项式，如果满足m0( jk )Pj (xi ) Pk (xi)（5-15）0( jk )i 0Ak则称 P0 (x), P1( x), , Pn( x) 为点集 x0, x1, , xm 上的离散正交多项式系。对于给定的节点 x , x , x ，可以按下列公式构造离散正交多项式系01m P( x), P ( x), P ( x)(nm)01nP0( x)1P ( x) ( x）（5-16）) P ( x110Pk 1( x)( xk

26、 1)Pk( x)kPk1(x)(k1,2,3, n 1)其中mk1xiP2k (xi )i0kk(k0,1,2, n1)kk1mk P2 k( xi)i 0(k0,1,2, n1)（5-17）(k0,1,2, n);.'.这样的P(x) 的首项系数为 1 的 k 次多项式。定理 2由式（ 5-16），（5-17）构造的多项式系 P( x), P ( x), P ( x)为点集01n x0, x1, , xm 上的离散正交多项式系。用离散正交多项式作曲线拟合设 ( xi, yi ) (i0,1, m) 为给定数据 P0( x), P1( x), Pn( x) 为点集 xi 上的离散正交

27、系，为由其所有线性组合生成的多项式集合，Span P0( x), P1(x), Pn( x) 。用正交多项式进行最小二乘曲线拟合，亦即求nmmn2yi2P( x)akPk( x)使其满足 Q P(xi )ak k( xi)yi min 利k0i 0i0k 0用多项式系 P0 (x), P1( x), Pn( x) 的离散正交性易知，此时法方程组为：mmP2 0(xi )P2 0( xi) yii 0a0i 0mm22P1 ia1P1 ii( x )( x) yi0i 0manmP2n (xi )P2n (xi )yii 0i 0其解显然为mmPk (xi) yiPk (xi) yiaki0i

28、0(k0,1, n)（5-21）mkP2k( xi )i0所以，容易得到拟合多项式nP( x)akPk (x)k0且其偏差平方和mmnQ P( xi ) yi 2y2ika2 ki0i 0k 0最后可得所求的拟合多项式yP(x)a0P0( x)a1P1( x)anPn( x)同样我们利用表的数据，采用上述原理和算法进行拟合，得到拟合多;.'.项式的系数如下：表离散正交曲线拟合系数系数0123求得结果67109500-5.9E+0717186619-1667018并且得到离散正交拟合曲线图6.3.2 如下：离散正交拟合曲线150010005000-500源生-1000-1500-2000-2500-30000.260