DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 2. 空间离散格式 2.1. 有限体积离散 一维守恒律方程为 其中 f ( u ) 为流通量。用有限体积法对(1)进行离散，首先将计算区间 剖分为N等分， a = x 1 2 ≤ x 3 2 ≤ ⋯ ≤ x N + 1 2 = b 且每个控制单元为 I j = [ x j − 1 2 ， x j + 1 2 ] ，其长度为 Δ x = b − a N ，对(1)在区间 I j 上进行积分，则有： ∫ x j − 1 2 x j + 1 2 u t d x + ∫ x j − 1 2 x j + 1 2 f ( u ) x d x = 0 令 ∫ u ( x j ， t ) d x Δ x = u ¯ ( x j ， t ) ，称其为积分平均值，于是可以得到 d u ¯ ( x j ， t ) Δ x = − f ( u j + 1 2 ， t ) − f ( u j − 1 2 ， t ) Δ x 利用数值通量逼近上式如下： d u ¯ ( x j ， t ) Δ x = − f ^ j + 1 2 − f ^ j − 1 2 Δ x 其中 u ¯ j ( t ) 是 的数值近似，数值流通量 f ^ j + 1 2 来近似逼近 f ( u j + 1 2 ， t ) ，其中 f ^ j + 1 2 = f j + 1 2 ( u j + 1 2 − ， u j + 1 2 + ) Lax-Friedrichs型的数值流通量为： f ^ ( u + ， u − ) = 1 2 [ f ( u + ) + f ( u − ) − α ( u + − u − ) ] (1) 其中 α = max u ∈ ( u + ， u − ) | f ′ ( u ) | 。 2.2. 对流项离散 本文构造了的非线性权的QUICK格式来逼近对流项 ，有如下格式 u j + 1 2 − = ( − 1 6 + 1 6 ω ) u ¯ j − 1 + ( 5 6 + 1 6 ω ) u ¯ j + ( 2 6 − 2 6 ω ) u ¯ j + 1 (2) (3) 其中 ω 为权，当解在光滑区域时，此时 ，数值格式变为三阶QUICK格式；当解在间断区域时，此时 ω = 1 ，数值格式变为一阶迎风格式，一阶格式可以抑制振荡，于是如何判断光滑与间断成为了关键问题。 为了有效的判断光滑与间断，使得格式在光滑处利用三阶QUICK格式，在间断处利用一阶迎风格式。根据对文献 [ 8 ] [ 9 ] [ 10 ] 里k个模板下的光滑因子 β r 进行分析， β r = ∑ l = 1 k − 1   Δ x 2 l − 1 ∫ x j − 1 2 x j + 1 2 ( d l d x l p r ( x ) ) 2 d x (4) 其中 是在k个模板 s r ( i ) = { x i − r ， ⋯ ， x i − r + k − 1 } ， r = 0 ， 1 ， 2 ⋯ ， k − 1 下构建的多项式。(4)右侧的积分是在 I j 上插值多项式 的所有导数的在 L 2 范数下的平方和。可知当 k = 3 时，根据文献 [ 11 ] [ 12 ] [ 13 ] 中对光滑度的测量 β r ，我们引入如下形式的光滑因子， r ( j ) = ( u ¯ j + 1 − 2 u ¯ j + u ¯ j − 1 ) 2 + ( u ¯ j + 1 − u ¯ j − 1 ) 2 通过分析计算可知，在连续时 r ( j ) 之间几乎相同，而在间断处 r ( j ) 之间存在较大的变化， r ( j ) 对数值解的间断部分比较敏感，通过分析文献 [ 6 ] ，我们给出了的如下光滑参数： 可以很好的衡量数值解的陡度和光滑程度，在光滑处通过泰勒展开可知 lim Δ x → 0 ∅ = 1 ，而在间断处 ∅ 的值变化比较大，由于 ∅ 值在间断处的变化不能够明确的界定间断与光滑，为了更加明确的判断，分析了以上的光滑因子与光滑参数的特性，我们给出 ∅ 的界定值，构造出了如下间断阈值： r a t i o = ( max { u ¯ j } − min { u ¯ j } ) 2 max { r ( j − 1 ) ， r ( j ) ， r ( j + 1 ) } + d d r max ， j = 1 ， 2 ， ⋯ ， N 其中 d d r max = max { ( u ¯ j + 1 − 2 u ¯ j + u ¯ j − 1 ) 2 } ，通过大量实验发现在间断处有 ∅ > r a t i o ，而在光滑处有 ，根据算例发现通过 与 的比较确实能够有效的判断间断与光滑区域。
system: [SYSTEM]

user: 3. 时间离散 对空间进行离散后得到关于时间的常微分方程组，为了保持整体格式的三阶精度，对时间利用三阶Runge-Kutta方法来进行离散，并且可以有效避免振荡，三阶Runge-Kutta格式如下： u ( 1 ) = u n + Δ t L ( u n ) u ( 2 ) = u n + 1 4 [ u ( 1 ) + Δ t L ( u ( 1 ) ) ] u n + 1 = 1 3 u n + 2 3 [ u ( 2 ) + Δ t L ( u ( 2 ) ) ]
system: [SYSTEM]

user: 4. 数值算例 4.1. 线性对流方程 { ∂ u ∂ t + a ∂ u ∂ x = 0 u ( x ， 0 ) = u 0 ( x ) ， x ∈ [ a ， b ] ， t > 0 (4) 4.1.1. 情形1 在给定初始条件 u ( x ， 0 ) = sin ( 2 π x ) ， x ∈ [ 0 ， 1 ] ， ，利用新格式计算当 Δ x = 20 ， 40 ， 80 ， 160 时，时间t = 0.1时的数值解，并求出格式的 L 1 ， L ∞ 误差及其精度阶(如表1)，计算公式如下： E L 1 = ‖ u n u m − u e x a c t ‖ L i = Δ x ∑ j = 1 n | u j n u m − u j e x a c t | E L ∞ = ‖ u n u m − u e x a c t ‖ L ∞ = max 1 ≤ j ≤ N | u j n u m − u j e x a c t | 精度阶的计算公式： o r d e r = ln ( E N / E 2 N ) ln 2 Table 1 Errors and orders for scheme Scheme mesh E L 1 o r d e r L 1 E L ∞ o r d e r L ∞ New scheme 20 3.876e−002 —— 9.69e−002 —— 40 1.28e−004 8.24228 2.012e−004 8.91172 80 1.4367e−005 3.12845 2.2976e−005 3.13043 160 1.006e−006 3.86292 1.598e−006 3.84579 表1. 格式误差与数值精度阶对比 由表可知此格式对于光滑初值条件下，新格式的数值解能够达到三阶精度，具有较高的分辨率。 4.1.2. 情形2 在给定间断复合初始条件 u ( x ， 0 ) = { 1 ，                                                                         1 ≤ x ≤ 3 x − 4 ，                                                               4 ≤ x ≤ 5 − x + 6 ，                                                           5 ≤ x ≤ 6 cos ( 0.5 π ( x − 8 ) ) ，                               7 ≤ x ≤ 9 exp ( − 4.5 ( x − 11 ) 2 ) ，                     10 ≤ x ≤ 12 0 ，                                                                       elsewhere 利用新格式计算当网格数 N = 200 时， CFL = 0. 5 ， x ∈ [ 0 ， 20 ] ，时间 t = 0.5 时的数值解并与精确解在图1中进行比较，精确解与QUICK格式下的数值解在图2中进行比较。 图1. 间断复合初始条件下的非线性QUICK格式和QUICK格式的数值解与精确解 由图1可知，QUICK格式与非线性QUICK格式相比，QUICK格式在间断处要产生的振荡，而新格式在间断处容易避免振荡，在光滑处的逼近效果也比较好。 4.1.3. 情形3 在给定间断初始条件 u ( x ， 0 ) = { 1 ，                                       0 ≤ x ≤ 0.2 4 x − 3 5 ，                     0.2 ≤ x ≤ 0.4 − 4 x + 13 5 ，               0.4 < x ≤ 0.6 1 ，                                     0.6 < x ≤ 0.8 0 ，                                     elsewhere ， x ∈ [ − 1 ， 2 ] ， CFL = 0.1 ， T = 0.2 当网格数 ，新格式的数值解及QUICK格式与精确解图像如图2。 图2. 间断初始条件下的非线性QUICK格式和QUICK格式的数值解与精确解 由图2可知QUICK的数值解在间断处会产生振荡，而新格式的数值解与精确解相比较，在间断处没有振荡，相比较新格式达到了比较好的逼近效果。 4.1.4. 情形4 对于有间断及相邻近处有极值的初始条件而言，对于格式的要求比较高，对于前面提到的格式，会在间断处产生振荡，不能够很好的区分光滑极值和间断，容易混淆这两种情况，从而在极值处的逼近效果不是很理想，下面我们给出在新格式下的数值解。 在给定有间断和极值的初始条件如下 u ( x ， 0 ) = { 1 − ( 10 x 3 ) 2 ，         x ∈ [ − 0.3 ， 0.3 ] 0 ，                                               其 他 ， x ∈ [ − 1 ， 1 ] ， CFL = 0. 5 利用新格式计算当 N = 600 时，时间 时的数值解和QUICK格式的数值解与精确解的比较如图3进行比较。 由图3可知，对于有光滑极值和极值邻近处有间断的初值问题，由数值解图像可知新格式无论是在间断处还是在极值点附近其逼近效果都非常好，而QUICK格式在间断处却会产生振荡，可知非线性QUICK格式比QUICK格式的逼近效果要好。 4.2. 无粘Burgers方程 无粘burgers方程解会随着时间的推移产生间断，给数值求解带来了很大困难，本文利用新格式求解此方程如下， 图3. 复杂结构初始条件下的非线性QUICK格式和QUICK格式的数值解与精确解 ∂ u ∂ t + ∂ ( u 2 2 ) ∂ x = 0 ， x ∈ [ 0 ， 2 ] ， t = 1.5 π 其初始条件为 u ( x ， 0 ) = 0.5 + sin ( 2 π x ) ，新格式的数值解图像及QUICK格式的数值解的图像如图4： 图4. Burgers方程的非线性QUICK格式和QUICK格式的数值解与精确解 由新格式与QUICK格式的数值解与精确解的图像比较可发现，在随着时间的推移，在逼近解的间断处时，QUICK格式下的数值解会有明显的振荡，而在非线性QUICK格式下的数值解会避免振荡，有很好的逼近效果。
system: [SYSTEM]

user: 5. 结论 本文基于有限体积下对双曲守恒方程进行数值离散，利用三阶QUICK格式与一阶迎风格式(FOU)加权后得到新的非线性QUICK格式来逼近对流项，本文通过引入间断阈值使得对流项的离散格式在间断处利用一阶迎风格式并且在光滑处使用三阶QUICK格式，通过算例分析，对于光滑初值能够达到三阶以上精度，对于有间断和极值的初值问题也可以避免振荡，使得数值解达到了比较好的逼近效果。
system: [SYSTEM]

user: 摘 要 高校科研团队是我国科研事业的重要组成部分，其研究成果对于国家的进步与发展具有重要意义。本文归纳总结影响高校科研团队的知识共享水平的影响因素：1) 知识特质性；2) 知识提供方与接受方个体差异；3) 知识共享渠道，具体表现为缺乏技术支持与交流手段不足；4) 知识共享环境。具体表现为缺乏良好的文化氛围与激励机制不完善。并对此提出组织要改善激励机制、缩小团队成员间的知识距离、营造知识共享的团队氛围、建立合理的团队结构、举办知识交流活动，从而提高高校知识创新能力，为我国高校科研发展做出更大的贡献。 关键词 :高校科研团队，知识共享，影响因素，对策 Copyright © 2018 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 问题的提出 在知识创新和知识传递的过程中，高校在这其中无疑扮演着必不可少的重要角色，而高校的主要职责在于教育和科研两大方面，我国的高校越来越把提高竞争力和科研创新能力作为自己的发展目标，甚至很多高校承担了大部分的国家级和省级重大项目，但是现在的知识经济的时代，单靠个人的能力实现某个科研成果已经变得不切实际，所以建立科研团队是现在高校或者企业进行科学研究的普遍形式 [ 1 ] 。 知识共享在知识管理过程中的地位愈来愈重要，知识共享直接影响知识的交流以及知识的创新，也影响着知识管理的效果，知识共享衡量知识管理的效果，知识管理的目的就在于知识的利用与知识的创新，知识必须在进行交流的过程中才能实现充分利用和创新，所以知识共享是知识管理的关键部分。 在我国高校科研团队中，他们的科研水平良莠不齐，具有很大的差距，许多高校具有高素质的人才，但科研水平较低，所以针对这样的现状，笔者针对为什么高校科研团队中配备了很多人才，但是科研水平反而很低做出了思考。所以本文归纳总结影响高校科研团队的知识共享水平的影响因素，最后根据影响因素提出高校科研团队知识共享水平的对策和建议，从而提高高校知识创新能力，为我国高校科研发展做出更大的贡献。
system: [SYSTEM]

user: 2. 高校科研团队知识共享的理论基础 近年来，知识至今没有一个明确而且统一的定义，学者们也在不同的维度对知识做出了自己理解的解释。越来越多人倾向于认为知识是一种组织资源。我们认为知识管理最重要的是知识的资源利用和价值转化，所以本文运用了达文波特和普鲁萨克对知识的解释：“知识是一种包含了结构化的经验、价值观、语境信息、专家见解和直觉等要素的动态的混合体，它为评估和利用新经验与信息提供了环境和框架”。 O’Dell和Grayson认为组织文化是推动知识管理成功实施的不可或缺的要素之一 [ 2 ] 。汤建影等对合作研发高校科研团队间的知识共享进行研究，认为重视人性关怀与情感交流、鼓励创新的组织文化能激发各成员共享知识，加速组织内部知识创造过程，提高合作研发高校科研团队间知识共享的质与量，相反，不合适的企业文化会阻碍知识在企业内部的流通与转换，隐性知识则不能转换为显性知识，扼杀团队成员的创新能力 [ 3 ] 。所以，建立包容、信任开放和鼓励创新的组织文化会对高校科研团队的知识管理起到巨大的推动作用。Teece认为企业通过利用信息技术将零散的信息和知识流结合在一起，可以实现对内外部知识的实时的监测 [ 4 ] 。Brown和Duguid认为信息技术可以支持不同团体间不同程度的正式和非正式的交流，实现知识在不同结构边界的跨越 [ 5 ] 。 尽管研究高校科研团队知识共享已经取得一定研究成果，但鲜有文献从影响知识共享意愿因素视角，探寻其对科研团队创新绩效的影响。李霞通过界定高校科研团队概念，分析知识共享行为与学习行为对团队绩效作用机制 [ 6 ] 。晋琳琳利用统计学的相关知识，探讨拥有不同学科背景下的科研团队成员之间知识共享和知识整合是否对团队创新绩效产生影响 [ 7 ] 。韩成标挖掘现有的高校科研绩效评价体系中存在的问题，在此基础上提出改善的意见 [ 8 ] 。
system: [SYSTEM]

user: 3. 高校科研团队知识共享水平的影响因素 (一) 知识的特性 OECD在《以知识为基础的经济》一书中，将知识按表达方法分为两类，一类是显性知识，一类是隐性知识。所谓显性知识是指可以通过正常的语言方式传播的知识，典型显性知识主要是以专利、科学发明和特殊技术等存在的知识，所以显性知识是可以表达的，也是有载体的 [ 9 ] 。然而隐性知识往往是个人或组织经过长期积累而拥有的知识，通常不易用言语表达，也不可能传播给别人或传播起来非常困难，所以隐性知识不易被认识到、不易衡量其价值也不易被其他人所理解和掌握。由于隐性知识的存在，就很造成在高校科研团队成员之间知识的不对称性，同时信息掌握者也会衡量隐性知识带来的价值，从而导致其不愿意进行知识共享 [ 10 ] 。 (二) 知识提供方与接受方个体差异 在高校科研团队中，一个完整的知识共享过程必须是知识提供方有意愿将知识分享给别人，在这个前提下也需要知识提供方能够准确的将知识向接收方表达出来，当知识提供方分享出自己的知识，也要考虑接受方的意愿以及其接受程度，如果太难或者太复杂，接受方对于知识的理解和消化就非常困难 [ 11 ] 。即使接受方理解了知识的内涵，其后期利用知识的能力也有待考量，所以这是高校科研团队知识共享水平的关键因素。 (三) 知识共享渠道 1) 缺乏技术支持 相关的技术是知识共享活动得以进行的重要手段，如数据采集手段、群件技术、知识地图、新型检索技术、计算机网络技术。如果缺乏相应的技术支持，会使得知识交流难以进行，不利于知识需求者的使用，另外，对于已有知识、信息，缺乏应有的信息技术对它们进行整理和综合，会使得大量的知识资源得不到有效的利用 [ 12 ] 。所以，没有技术的支持，对知识共享的水平会有很大的影响。 2) 交流手段不足 高校科研团队之间如果不进行知识的交流，对于科研成果会有很大的影响，所以，成员之间要进行紧密的交流，但是要进行知识他们要进行知识共享，就需要有必要的交流手段，或者交流方式或交流平台，比如经验交流会、讨论会、团队学习等等，离开了这些交流手段，自然就会影响知识共享活动的进行。但是从每个成员的角度出发，他们被组织定格在一个狭小的范围里面，按照管理者的安排去进行工作，无法突破岗位对于个人的束缚，所以彼此之间的交流机会无疑就会很少。所以，类似于经验交流会、讨论会之类的交流手段和交流平台就很难付诸现实。但是此类交流手段的不足，自然会妨碍知识共享活动的正常进行。 (四) 知识共享环境 要进行知识共享活动，首先就需要有一定的客观环境，如果环境合适，知识共享活动就会顺利进行，如果环境不合适，知识共享活动就很难在各个成员之间进行。 1) 缺乏良好的文化氛围 众所周知，良好的文化氛围可以对人们的行为起规范、凝聚、激励和强化作用，从而促进知识共享的顺利进行。但是，如果没有良好的文化环境，在科研团队中，就会被一片拒绝或消极对待知识共享的文化所围绕，这就很明显会对知识共享产生很强大的阻碍作用。 2) 激励机制不完善 良好的激励机制对于科研团队的创造力具有很好的促进作用。由于只是共享只是人们应该积极参与的行为，而并非“理所当然”的行为，同时这其中也包含着一定的利益风险，所以，只要有了一定的激励机制，人们才乐于参与知识共享的活动，而极力机制不完善，知识拥有者一般是不会主动去进行知识分享的。所以，没有完善的激励机制，就难以调动知识拥有者参与知识共享活动的积极性，从而影响科研团队的科研水平。
system: [SYSTEM]

user: 4. 提高高校科研水平的对策及建议 (一) 改善激励机制 因为高校科研团队是一个文化水平相对较高的团队，而成员自身获取知识的过程非常困难，要付出时间成本和金钱成才能形成一个自己独特的知识体系所以有利己主义的成员自身会在做知识共享前衡量知识共享会不会给自己带来相应的报酬，这种报酬是不是大于自身失去竞争力的风险，只有收益大于风险，才会有知识共享的意愿，因此，建立一个良好的激励机制是提高知识共享水平的必然要求，在团队成员进行知识共享以后，希望得到相应的报酬(表扬、奖金、津贴、福利等)，这种激励属于保健因素，虽然不能提高知识共享的质量，但是可以相应提高知识共享的意愿和知识共享的数量 [ 13 ] 。同时可以给进行知识共享的成员一定的晋升机会，让有知识共享成员得到认同感，这样既满足外在保健因素，也能满足知识共享成员的精神需求，让知识共享成员得到成就感，从而达到另外一种层次，即激励因素，这种因素虽然不能提高知识共享水平的数量，但是可以提高知识共享的质量，使得知识共享成员之间有一种和平且互惠互利的感觉。 (二) 缩小团队成员间的知识距离 知识共享的目的是实现知识利用甚至实现知识创新，知识的不对称性可以说是影响知识共享水平的基本根源所在，如果只是共享主体之间的知识距离太大，那么其他多么完善也是徒劳，因为在前期的知识交流过程就无法实现知识共享的实现，那么讨论只是共享水平也就没有意义。 缩短知识共享主体之间的知识距离可以从两个方面考虑，一是建立一个知识交换平台，个体成员对自身知识体系的积累是存在很大的差异的，知识共享主体之间相互融合彼此的知识资源，才能有效缩短知识距离，二是探索知识共享主体之间知识共享模式，在平时进行知识交流的时候多组织相关会议和讲座，缩短知识共享主体之间文化差异和习惯差异，进而实现知识共享。 (三) 营造知识共享的团队氛围 科研团队建设初期，团队成员之间来自于不同的学术背景，有着不同的研究习惯，有各自的研究领域，为了一个共同的科研目的组成一个团队，需要一定的时间进行磨合才能够习惯团队里的每一个人，如果团队内部气氛不够活跃，团队成员不能够畅所欲言，高校科研团队是一个学术性团队，在一个没有互相信任的不和谐的环境下，知识共享是无法很成功地实现的，因此在高校科研团队中，团队带头人的影响力是非常重要的，团队带头人应该带领团队成员为团队创造一个畅所欲言、言论自由、观点自由，可以进行头脑风暴进行知识交流，并且互相信任、易于沟通的氛围，在这种知识共享环境下，团队成员可以专心努力做研究，在一个和谐的学习型团队下，和团队共同发展，共同为团队实现知识创新做出贡献。 (四) 建立合理的团队结构 团队结构有垂直结构和扁平结构，垂直结构的团队类型在知识交流时受到繁琐的上下级限制，机构森严，层级与层级之间关系疏远，很大程度上影响了知识共享的效率和效果，成员之间没有很好的沟通机制，而对于扁平化的团队结构，团队成员之间层级设计较少，有利于扁平化的团队结构，团队结构中团队带头人管理权力最大，层级与层级之间没有复杂的组织权利机制，成员和成员之间能够很容易地直接进行联系和沟通，不会有层级压力导致知识交流壁垒的情况。垂直化的团队结构主要遵循上级安排下级执行的模式进行知识共享，这样往往迫于权力的压力，成员很难言论和观点自由，而扁平化的团队结构使成员与成员之间关系相对较为平等，能够有更多的机会参与到科研项目的决策或意见层次，在科研项目中有问题也可以观点自由，在扁平化结构的科研团队中可以最大程度地实现知识的流动，最大程度地发挥团队的优点，从而使得团队知识共享效率和效果显著，为团队实现进一步的知识创新。 (五) 举办知识交流活动 在知识经济时代，知识的更新速度越来越快，阶段性的知识交流显得更加重要，前面四点建议中提到改善激励机制，缩小团队成员之间的知识距离，营造良好的团队氛围，建立合理的团队结构，但是在这些基础上，进行物理性的定期知识交流还是有必要的，虽然网络平台的构建能够很好地实现知识共享，但是网络平台是对于已有知识的存储和再次利用，而且只是类型大多是面对显性知识，对于更新的知识如果上传到网络平台，让团队成员各自进行账号登录学习，成员自身可能难免漏掉一些重要的知识没有学习到，而进行物理的知识交流(会议、讲座、头脑风暴等)，这种形式的知识交流可以消除团队成员之间知识的不对称性，提高后期进行科研项目的效率。团队成员在良好和谐、互相信任的团队氛围下，在完善的激励制度下，做到以传统的物理知识交流和网络平台交流相结合的模式进行知识共享。高校科研团队不仅要发扬和创造传统物理知识交流的机会，也不能忽略基于互联网技术知识的网络平台交流相结合，及时更新知识，为科研团队积累竞争力实现更多的知识创新。
system: [SYSTEM]

user: 摘 要 随着工程项目的大型化和复杂化，消防工程项目的进度管理越来越复杂。如何提高施工进度生成效率是理论研究者和工程管理人员面临的难题。本文基于BIM技术，采用改进遗传算法提高消防工程施工进度生成效率，使施工能力最大化。考虑到构件之间的约束依赖关系，将这些约束规则与施工进程进行关联，从而为施工进度的自动生成提供可能。此外，此方法可从项目的BIM模型中获取足够的信息，也可以在3D CAD系统中得到有效的可视化信息。最后，通过案例验证了消防工程施工进度生成效率的可行性和有效性。 关键词 :消防工程，项目进度管理，进度生成效率，BIM，GA算法 Copyright © 2019 by author and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着工程项目的大型化和复杂化，作为一项系统工程，消防工程项目的进度管理越来越复杂 [ 1 ]。根据计划规定的期限合理地规划成本或根据最低成本的要求寻求最佳生产周期是项目管理者追求的重要目标 [ 2 ] [ 3 ]。如何寻求一个任务的调度序列与资源分配方案，使得消防工程施工进度计划总体用时最短是管理者面临的重要难题。现代项目所具有的技术难度大、结构形式复杂、信息量多、项目管理难度大等特点，使得传统项目管理方法越来越不适用建筑业的发展形势。 BIM技术顺应建筑业发展对信息化技术的需求而产生并受到广泛关注。BIM技术作为一共享的知识资源，包含了建设项目全生命周期所有物理和功能特性的数字化表达。该技术的出现为项目管理理论和技术的发展提供了新的思路 [ 4 ]。当前基于BIM的各项研究和应用不断发展，各项应用研究之间也出现了较大的交叉性和互用性，比如基于BIM的5D模拟技术、基于BIM的进度风险管理等，因此基于BIM的施工进度计划自动生成模型的研究不仅能够解决施工进度计划的编制问题，也能够促进其他技术的深入研究和应用的推广。 近年来大量研究多将智能算法应用于资源配置及资源受限的进度编制问题。褚晓剑等对ZGC项目的进度管理的改进进行了研究 [ 5 ]。董兆佳等对赢得值法在海洋工程项目进度控制中的应用进行了分析 [ 6 ]。许靖对钢结构进度管理进行了分析 [ 7 ]。许海燕应用BIM对项目进度管理进行了研究 [ 8 ]。Hegazy等则将遗传算法与过程模拟相结合进行资源优化等工作 [ 9 ]。Symeon等提出了一个方法来安排资源受限工程利用基于蚁群优化算法，该方法应用于资源受限的网络，用于检查资源可用性约束影响关键路径计算和项目完成时间 [ 10 ]。然而，对如何提高消防工程施工进度生成效率等问题的研究还较少。并且，输入的3D模型越复杂，处理整个3D模型元素进度的算法越困难。因此，研究生成施工进度的效率问题是非常有必要的。 针对此问题，本文提出了改进遗传算法并采用BIM技术解决现代消防工程量庞大、结构与技术复杂化等问题，通过对三维设计的工业基础类文件进行解析，得到实体构件的位置和几何信息，结合使用本体描述施工过程中的各种约束规则，从而得到构件间的约束依赖关系；然后根据构件间的约束依赖关系生成约束矩阵；最后利用算法生成消防工程项目进度计划。该成果为进度管理提供了一种新的思路和方法。
system: [SYSTEM]

user: 2. 改进GA算法 遗传算法(genetic algorithm， GA)具有快速的全局搜索能力，在遗传算法中，因偏重的差异，交叉的概率和变异的概率各不相同，通过交叉和变异这两个操作使其具备平衡整体和局部的搜索能力 [ 11 ]。这两个操作的高效性使遗传算法具有两大优点：一是使遗传算法具有一定程度的局部随机搜索能力，当迭代到一定次数已接近较优解时，通过变异操作可以加快向较优解的收敛速度；二是能维持可行解的多样性，防止呈现早熟的现象。在不同行业的调度和旅行商问题(TSP)的解决方案中，遗传算法作为一种全局优化方法得到了广泛的应用。许多领域的研究人员长时间的研究验证了GA方法稳定可靠的性能和调度优化方面的可靠性。此外，该方法有不同的反应机制，还能避免局部最优。但遗传算法并没有利用迭代过程中的反馈信息，导致大量重复迭代，求准确解的效率不高。为了克服基本遗传算法的缺陷，本文借鉴蚁群算法中信息素反馈思想，提出基于信息素的改进遗传算法。 蚁群算法(ant colony algorithm， ACO)是一种求解组合优化问题的群体寻优算法 [ 12 ] [ 13 ]。其工作原理是通过信息素的正反馈，实现算法的寻优。因此，本文采用蚁群算法中信息素正反馈的思想，对遗传算法进行改进。首先在遗传算法初始化过程中，利用蚁群算法的全局搜索能力，在进行遍历的过程中引入变异和交叉操作，提高可行解的多样化程度，避免算法的早熟，提高算法的鲁棒性和获得较优解的效率。具体改进策略为： 1) 信息素的优化。在遗传算法初始化时，将蚁群算法中各路径的信息素初值设为 τ max ，初始时通过获得不同路径中的信息素。信息素的更新机制为所有蚂蚁都遍历完所有的任务后，只有相对路径短的蚁群才进行信息素更新，即只取最短路径长度在前90%的蚁群进行信息素更新，而所有路径的更新方程均采用： τ i j ( t + 1 ) = ρ ⋅ τ i j ( t ) + ∑ Δ τ i j k ( t ) (1) 其中， τ i j ( t ) 为路径 ( i ， j ) 在 t 时刻的信息素强度； Δ τ i j k ( t ) 为蚁群 k 在路径 ( i ， j ) 上单位路程的信息素； ρ 为信息素的保留率； Q u 为常数，是单位路程的信息素总量； L k 表示第 k 次蚁群在本次循环中所走过的路径的长度。 2) 变异操作的优化。为优化算法在变异操作中的随机性，提高算法的效率。根据变异概率 P m 随机选择变异点，并且依据信息素浓度 f i j ( t ) 决定是否互换变异位置，确保变异位置中交换后信息素浓度比交换前的信息素浓度高。 3) 改进遗传算法步骤 步骤1：通过蚁群算法得到初始值，并根据初始值完成遗传算法的初始化； 步骤2： i = i + 1 ； 步骤3： L k 出发点随机定位于当前解的可行域中，根据交叉概率对目标进行交叉操作，对比后若信息素浓度提高，则接受新值，否则拒绝； 步骤3：根据变异概率对目标进行变异操作，若信息素浓度提高，则接受新值，否则拒绝； 步骤4：对适应度函数 L k 小于给定值的路径，按式(1)进行信息素的更新，并进行轨迹强度的修改； 步骤5：若 i < 给定的迭代次数，则转步骤2，否则，则输出目前较优解。
system: [SYSTEM]

user: 3. 案例分析 选择某工程钢结构施工计划作为对象，研究建立的仿真分析模块依托BIM技术，基于该平台的数据接口导入并集成过程施工计划与各项任务的资源数据。 采用同一个工程，分别用传统的蚁群算法和遗传算法进行计算对比。在本研究中，只关注结构模型和内部元素，而管道、设备、暖通空调等并不包括在这些模型中。BIM模型如图1所示。图1中的模型是一个包括36根柱子和72根梁，一共108个元素的结构模型。从3D模型的详细输入情况也表明了3D钢结构模型输入的复杂性。 图1. BIM 3D模型 1) IFC处理 在采用BIM技术过程中，IFC主要描述的是把实体建筑中的信息描述成模型形式，这也表明IFC用来描述实体建筑中的信息是非常适用的，具体表现为： a) IFC定义中的实体类层次结构与实体建筑中概念之间的层次结构非常相似； b) IFC中每个实体构件都具备相关的属性，如几何位置信息、形状信息、材料信息等，这与实体建筑中概念和实体的属性机制相似； c) IFC中的项目的组成关系，例如一个项目中包含多个场地、一个场地则包含多个建筑、一个建筑又包括多个楼层，这与实体建筑中的组合关系非常类似。 2) 约束处理 在消防工程中，存在着各种各样的约束。然而，在标准中，并没有针对这些约束规则提供有效的表达方式。因此需要使用一种方式来描述这些约束规则并与施工进程进行关联，从而为施工进度的自动生成提供可能。 参照历史相似项目的进度计划以及施工过程中的约束规则，为本项目的建设提供参考，通过记录同一类型的约束规则，可辅助进度计划的制定。使用规范的语言对不同类型的建筑历史信息进行记录，可以有效发挥历史数据的作用，为未来的项目提供参考，这正体现了大数据的用途，可以很大程度地降低进度计划制定的难度。 3) 改进GA的应用 本文采用改进GA算法导出消防工程施工进度计划。算法操作及其分析过程中的适用性也是非常匹配。算法整体流程为： 第一步：解析IFC文件。 第二步：构建描述实体构件信息和施工约束的信息库，实体构件的信息可以通过IFC文件解析结果自动转换，施工约束规则主要通过人工方式进行输入。 第三步：依据实体的位置、几何形状信息，结合本体中的规则，构建MoCC约束矩阵。 第四步：使用GA-CRO混合算法进行分析，得到实体构件之间的依赖顺序关系以及进度计划表。 算法运行的预期结果是实现满足施工能力和稳定性模型约束下的多个完整的施工顺序。如果它可以表明本次设计研究满足了既定的目标，那么可以证明该算法适用于自动稳定的发展建设项目调度并且具有一定程度的效率保证。 4) 结果分析 优化对比结果见表1。 Table 1 Comparison of the result 解的优化程度 生成时间减少率 Improved GA algorithm 15% 10% Ant colony algorithm 13% 12% 表1. 优化结果对比 表1表明采用遗传算法，对于给定的3D模型，能生成稳定、可构成的施工项目进度，较好的解决了进度生产问题。试验结果验证了所采用的约束条件策略应用有效性，可保证消防工程施工计划的固定时性，与兼顾全局的调度视角，使GA有效地利用资源实现对原计划的优化，并具有良好的鲁棒性。GA优势明显，这是因为其具有更强的解空间搜索能力。在充分迭代执行的情况下，GA有更好的整体优化能力。因此，从实验结果分析来看，能证明生成施工进度的效率提高。 这种方法不仅为项目能生成稳定的施工进度，这也显示了消防工程序列的三维模型，说明完成三维模型的时间跨度的可行性。
system: [SYSTEM]

user: 4. 结论 本论文以消防工程施工进度生成效率优化为研究对象，基于GA算法和BIM技术，将三维设计作为输入信息，通过对三维设计的工业基础类文件进行解析，得到实体构件的位置和几何信息，结合施工过程中的各种约束规则，从而得到构件间的约束依赖关系。然后，根据构件间的约束依赖关系生成约束矩阵。最后，选择算法生成建筑项目进度计划表。针对BIM技术，本文提出一种提高消防工程施工进度生成效率的优化方法。此外，计算结果能在3D CAD系统中进行有效的可视化模拟。 但是本研究在某些方面也有一些缺陷和不足，可从以下几个方面进一步深入研究： 1) 在一些进度计划中，在资源的有效性方面可能会出现进度计划数据的重叠。因此，在后续的研究中可针对此难点进行深入分析，进一步提高效率。 2) 本文算法只考虑了构件和构件间的约束关系，生成了初始进度计划，也可以进一步研究将机械、人力等因素加入到模型中，对初始进度计划进行优化和完善，与已有的基于混合算法和BIM的施工进度优化的成果相结合，生成更加全面和准确的施工进度计划。这方面的研究也极具价值。
system: [SYSTEM]

user: 摘 要 针对王家岭煤矿软岩顶板沿空巷道围岩稳定性控制中存在的问题，本文通过试验研究、理论分析、数值计算等研究方法，确定了厚泥岩顶板物理力学性质，在此基础上，分析了现有支护中存在的问题，给出在支护优化设计中应考虑的因素，通过数值模拟分析得到了沿空巷道的合理支护方案及参数，并进行了工程应用。结果表明，采用新支护方案的300 m试验段支护效果明显改善，巷道开掘支护3个月内断面顶底板平均移近量为160 mm，顶板完整性保持较好。 关键词 :软岩顶板，沿空掘巷，数值模拟，支护设计 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 沿空掘巷作为工作面回采巷道的一种布置及掘进方式，因其具有减少煤炭资源损失、简化接续工艺等优点，在越来越多的矿井中得到应用，国内外许多学者对沿空巷道在应用中存在的问题进行了大量研究，特别是沿空巷道的围岩控制问题。柏建彪 [ 1 ]、朱学富 [ 2 ]、尤跟林 [ 3 ]、戴磊 [ 4 ]、郝金鹏 [ 5 ]、宗彩建 [ 6 ] 等针对不同情况的沿空巷道给出了相应的支护控制措施；李磊 [ 7 ] 探讨了复合顶板沿空巷道的支护控制问题；诸葛长华 [ 8 ] 提出了耦合让压的沿空巷道支护方法；张农 [ 9 ]、张进鹏 [ 10 ] 研究了相向掘进巷道围岩控制技术，这些成果推动了沿空掘巷技术的应用。 然而采矿地质条件复杂多样，特别是在大采高及采场动压条件下，覆岩切顶冒落影响范围大，沿空回采巷道多次出现大变形、顶板离层、碎涨冒落等难以控制的失稳破坏情况。为提高此种回采巷道围岩支护的稳定性，基于煤岩工程性质和锚(索)支护原理，对沿空巷道支护参数进行设计优化，为巷道围岩的稳定性控制提供依据。
system: [SYSTEM]

user: 2. 工作面沿空掘巷支护理论分析 王家岭矿采区接续工作面回风巷均采用沿空掘巷，巷道断面初步设计为5.2 × 3.6 m。通过对巷道顶板岩层的工程性质研究，针对巷道围岩特点，巷道围岩稳定性与控制技术需考虑三个因素：围岩的应力状态、围岩的赋存状态及力学性能、支护技术。采用以下技术思路： 2.1. 重视对巷道预留煤层的保护 2#煤预留煤层是直接顶的保护层，煤层有效锚固是顶板安全的基础，顶板支护应遵循“一次性支护和高预应力与预应力扩散”的支护原则。 2.2. 注意直接顶的防水 2#煤直接顶属亲水性岩石，且有膨胀性，瞬时的膨胀率和膨胀力均较大，这也是部分巷道顶板离层破碎的原因。作为2#煤巷道顶板支护的锚固层，其工程性质应特别重视。支护设计中兼顾煤层与直接顶，重视锚杆(索)支护与煤岩层的相互作用关系。 2.3. 避免直接顶裸露，尽量沿底掘进 沿2#煤掘进巷道，应尽量沿底板掘进，掘进支护注意预留煤层保护，避免直接顶裸露。
system: [SYSTEM]

user: 3. 锚杆支护参数数值模拟动态分析 采用FLAC3D软件进行数值计算，沿空掘巷煤柱宽度为5 m，建立数值计算模型；根据设计思路及已知断面尺寸，确定巷道顶板锚杆支护间距取900 mm，确定以下4组模拟方案(表1)。 Table 1 Numerical simulation calculation schem 模拟参数 方案编号 锚杆间、排距/mm 锚杆长度/mm 锚杆排距 帮锚杆间距 顶板 两帮 帮锚杆间距 1-1 1000 850 2400 1800 1-2 1000 900 2400 1800 1-3 1000 950 2400 1800 锚杆排距 2-1 900 900 2400 1800 2-2 1000 900 2400 1800 2-3 1100 900 2400 1800 顶锚杆长度 3-1 1000 900 2400 1800 3-2 1000 900 2500 1800 3-3 1000 900 2600 1800 帮锚杆长度 4-1 1000 900 2400 1800 4-2 1000 900 2400 2000 4-3 1000 900 2400 2200 表1. 数值模拟计算方案 3.1. 帮锚杆间距 巷道两帮围岩变形及锚杆受力受帮部锚杆间距影响较大，而巷道顶板锚杆受力及顶板围岩变形受其影响较小，如图1所示。两帮锚杆间距为900 mm和950 mm时巷道围岩变形及锚杆受力变化不大，结合巷道断面高度，两帮锚杆间距为900 mm布置4根锚杆间距较为适宜。 图1. 帮锚杆间距模拟 3.2. 锚杆排距 巷道锚杆受力及断面变形受锚杆排距影响较大，特别是顶板围岩变形和顶板锚杆受力，对比锚杆排距为1100 mm和900 mm两种模拟方案结果，前者顶板下沉量约为230 mm，相比后者高近一半；锚杆排距为1000 mm和900 mm两方案下巷道围岩变形及锚杆受力差距较小，如图2所示。综合工程效益考虑，锚杆排距为1000 mm较为适宜。 图2. 锚杆排距模拟 3.3. 顶锚杆长度 巷道锚固受力和断面变形受顶板锚杆间距影响较大，特别是巷道顶板下沉量与顶板锚杆长度的变化联系紧密，巷道顶板锚杆长度增加提高了其在巷道支护中承担的围岩压力，能够一定程度减小两帮锚杆支护的压力，如图3所示。采用长2600 mm的顶板锚杆能够提高巷道断面结构稳定性，但是考虑到顶板岩石的弱水特性，顶板锚杆长度选用2500 mm较符合岩石工程条件。 图3. 顶板锚杆长度模拟 3.4. 帮锚杆长度 巷道锚固受力和断面变形受帮锚杆长度影响较大，特别是巷道两帮移近量与帮锚杆长度的联系密切，巷道帮锚杆长度增加提高了其对巷道两帮尤其是上帮煤柱的变形控制，两帮承载力的提高减小了顶锚杆承受的巷道支护围岩压力，对比帮锚杆长为2000 mm和2200 mm两种方案模拟结果表明，两者巷道围岩变形和锚杆受力相差不大，帮锚杆长度可选用2000 mm，满足工程需求，如图4所示。 图4. 帮锚杆长度模拟
system: [SYSTEM]

user: 4. 锚杆支护参数确定及应用 综上所以模拟方案结果，20109工作面材料道锚杆支护参数：顶锚杆间排距为900 × 1000 mm，锚杆长2500 mm；帮锚杆间排距850 × 1000 mm，锚杆长2000 mm。图5为材料道两种支护方案效果的数值模拟结果对比，其中(a)为数值计算采用的初始方案，(b)为通过动态设计对比分析调整优化后确定的支护方案。 图5. 巷道围岩塑性区 对比两者围岩塑性区分布模拟结果可看出，模拟分析所得支护方案对塑性区应变控制特别是对于煤柱受剪应力作用塑性区发育的控制作用效果提升明显，提高断面整体支护强度增强了围岩自身结构的稳定性，同时保护了煤柱。 改进后的支护方案现场应用效果表明：300 m试验段支护效果明显改善，3个月内顶底板移进量为约160 mm，顶板完整性较好。
system: [SYSTEM]

user: 5. 结论 1) 通过煤岩物理力学试验，表明巷道顶板为膨胀型泥岩，易离层、潮解碎涨，为巷道支护优化设计提供计算依据； 2) 理论分析认为，此类巷道支护设计时需考虑围岩的应力状态、围岩的赋存状态及力学性能、支护技术等因素，并给出了设计优化思路； 3) 采用数值模拟方法对巷道锚杆支护参数进行动态计算分析，优选出合理的支护方案及参数，并进行了现场应用，结果表明，改进后的方案支护效果良好，可满足安全生产要求。
system: [SYSTEM]

user: 摘 要 大体积混凝土在施工过程中，由于混凝土的量大，水泥的水化热热量大，在混凝土内外散热不均匀以及受到内外约束的情况下，混凝土内部易产生较大的温度应力，导致混凝土发生开裂。以南盘江特大桥工程的劲性骨架钢筋混凝土拱桥为依托，为防止拱座基坑大体积混凝土温度裂缝的产生，通过计算进行模拟混凝土浇注和水化的放热过程并布置散热管控制混凝土的温度。 关键词 :大体积混凝土，防裂，温控，水化热 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 从以往大体积混凝土的施工经验中可知，大体积混凝土在施工过程中，由于混凝土的量大，水泥的水化热热量大，在混凝土内外散热不均匀以及受到内外约束的情况下，混凝土内部易产生较大的温度应力，导致混凝土发生开裂 [ 1 ] [ 2 ] [ 3 ] [ 4 ] 。因此，大体积混凝土施工中的温度控制是防止混凝土开裂的关键，如何进行混凝土的温度控制即采取那些技术措施防止混凝土开裂，是桥拱座大体积混凝土施工前必须研究解决的重要课题 [ 5 ] [ 6 ] [ 7 ] 。
system: [SYSTEM]

user: 2. 混凝土产生开裂的原因分析 2.1. 水泥水化热的影响 混凝土浇注后，水泥在水化过程中要释放大量的水化热，而大体积混凝土结构断面较厚，热量聚集在结构内部不易散发，使结构内部温度升高。由于混凝土的导热性能较差，浇注初期混凝土强度及弹性模量都很低，对水化热引起的急剧温升约束不大。随着混凝土龄期的增长，弹性模量的增高，对混凝土内部降温收缩的约束也就越来越大，以致产生很大的拉应力，当混凝土的抗拉强度不足以抵抗这种拉应力时，便开始出现温度裂缝。 2.2. 外界气候的影响 大体积混凝土浇注时间较长，外界气温昼夜变化大。气温越高，混凝土浇注温度越高，同时在高温条件下，混凝土不易散热。昼夜温差或混凝土浇注后气温突然下降，使混凝土内外产生较大温差(>25℃)，而温度应力则是由温差所引起的温度变形造成的，因此温差越大，温度应力也越大，就越容易出现裂缝。 2.3. 混凝土的收缩变形影响 混凝土中的水分要蒸发，会引起混凝土体积的干缩。在高温下，混凝土中水分蒸发快，很容易引起干缩变形。如拌合、浇注、振捣等方法不当，混凝土内粗细骨料分配不匀，密实度不等，收缩系数不一，在约束不同的情况下，即可在脆弱部位引起开裂。
system: [SYSTEM]

user: 4. 结论 通过以上研究，可以得出如下结论： 1) “内排外保”的温度控制方式，可以减少混凝土内外温差。 2) 通水后降温可以通过循环水管的进水流量控制，实际施工时可根据实测温度值来控制循环水冷却的周期；按计算结果分析要将各龄期内外温差、环境和表面温差控制在20℃以内完全可以做到。 3) 采取“内排外保”措施后，大体积混凝土温度应力不大于抗拉应力且抗裂安全系数均大于1.15，满足要求，混凝土不会产生裂缝。
system: [SYSTEM]

user: 摘 要 中国是一个农业大国，几千年来传统的原始落后的农耕社会使中国农业发展滞后于全社会经济的发展 [ 1 ] 。新世纪中国发展的关键在于解决农民的发展问题，其实质就在于提高农民的实际收入。本文为研究中国农村居民收入的影响因素，通过建立模型并加以检验，从而得出农作物总播种面积与人均国内总产值对农村居民收入的线性关系显著存在且在经济角度存在正相关。 关键词 :农民年人均纯收入，国内总产值，模型回归，影响因素 Copyright © 2018 by author and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 自改革开放以来，中国经济迅速发展，GDP逐年攀升，人民的生活水平总体上有很大的改善 [ 2 ] 。但农村经济发展和农民生活水平的提高却面临着严重的困难，成了制约中国经济发展的重大因素。伴随着人民总体生活水平的提高，中国的城乡基尼系数却在不断的扩大，这不仅影响了中国经济的健康正常发展，还有可能引发一系列的社会问题。因此，“三农”问题成为中国经济研究的热点问题，提高农村居民收入迫在眉睫。
system: [SYSTEM]

user: 2. 相关理论文献与数据来源 2.1. 文献综述 2.1.1. 人均国内总产值 人均国内生产总值，即“人均GDP”，常作为发展经济学中衡量经济发展状况的指标，是最重要的宏观经济指标之一，它是人们了解和把握一个国家或地区的宏观经济运行状况的有效工具 [ 3 ] 。 2.1.2. 农产品生产价格指数 农产品生产价格指数是反映一定时期内，农产品生产者出售农产品价格水平变动趋势及幅度的相对数。该指数可以客观反映全国农产品生产价格水平和结构变动情况，满足农业与国民经济核算需要。农业收入是农民收入的重要来源，农产品生产价格指数也直接关系到农民的切实利益。 2.1.3. 化肥施用量 化肥施用量对农产品的产量有较大的影响，化肥施用的合理水平对于提高农产品产量有着举足轻重的作用，化肥施用量与农村居民的收入成正相关关系。 2.1.4. 农作物总播种面积 农作物播种面积与我国农村居民收入之间存在正相关，农作物播种面积越大，总产值就会增加，从而使农村居民收入提高。 2.2. 数据来源 本文所采用的人均国内总产值、农产品生产价格指数、乡村从业人员、化肥施用量、国家财政支农支出、农作物总播种面积的数据均来自中国国家统计局官方网站，即《2017年国家统计年鉴》中2017年的数据见表1。 Table 1 Relevant years and dat 年份 X 1 人均国内总产值(元) X 2 农产品生产价格指数(上年 = 100) X 3 化肥施用量(万吨) X 4 农作物总播种面积 Y 我国农民年人均纯收入(元) 1996 5074 119.9 3593.7 149，879.3 1577.7 1997 5845.89 104.2 3827.9 152，380.6 1926.1 1998 6420.18 95.5 3980.7 153，969.2 2090.1 1999 6796.03 92 4083.7 155，705.7 2162 2000 7158.5 87.8 4124.32 156，372.8 2210.3 2001 7857.68 96.4 4146.41 156，299.9 2253.4 2002 8621.71 103.1 4253.76 155，707.9 2366.4 2003 9398.05 99.7 4339.39 154，635.5 2475.6 2004 10，541.97 104.4 4411.56 152，415 2622.2 2005 12，335.58 113.1 4636.58 153，552.6 2936.4 2006 14，185.36 101.4 4766.22 155，487.7 3254.9 2007 16，499.7 101.2 4927.69 152，149.5 3587 2008 20，169.46 118.49 5107.83 153，463.9 4140.4 2009 23，707.71 114.06 5239.02 156，265.7 4760.6 2010 25，607.53 97.6 5404.4 158，613.6 5153.2 2011 30，015.05 110.94 5561.68 160，674.8 5919 2012 35，181.24 116.45 5704.24 162，283.2 6977.3 2013 39，544 102.7 5838.85 163，415.7 7916.6 2014 43，269.64 103.2 5911.86 164，626.9 8926.7 2015 47，203 99.8 5995.94 165，446.2 10，489 2016 50，251 101.7 6022.60 166，373.81 11，422 表1. 相关年份及数据 被解释变量：Y——我国农民年人均纯收入(元)。解释变量：X 1 ——人均国内总产值(元)；X 2 ——农产品生产价格指数(上年 = 100)；X 3 ——化肥施用量(万吨)；X 4 ——农作物总播种面积。
system: [SYSTEM]

user: 3. 模型的建立 ADF检验及协协整检验 确定模型： Y = β 0 + β 1 X 1 + β 2 X 2 + β 3 X 3 + β 4 X 4 + μ 其中 β 0 、 β 1 、 β 2 、 β 3 、 β 4 为待估参数，而 μ 为随机误差项。 由于所用数据为时间序列数据，需要检验其平稳性。 在命令窗口genr输入YT = log(Y)，XiT = log(Xi)，从而得到新的形式，再进行图示检验检验。通过该时间序列的时间路径图来粗略地判断其是否存在平稳性。 对被解释变量和各解释变量进行图示检验见图1，图2，图3，图4，图5。可以看出，每个变量在图中均表现出在不同的时间段具有持续上升的趋势，因此，该时间序列是非平稳序列。 图1. 图示检验 图2. 图示检验 图3. 图示检验 图4. 图示检验 图5. 图示检验 并用EG两步法考察它们之间是否存在协整关系见图6。现在对YT，X1T，X2T，X3T，X4T分别考察其单阶整数。 YT： 图6. 协整关系 从检验结果看，在1%，5%，10%三个显著水平下，单位根检验的Mackinnon临界值分别为−3.831511，−3.02997，−2.655194，t检验统计值3.324114大于相应临界值，从而不能拒绝H 0 ，表明YT序列存在单位根，是非平稳序列。 为了得到YT的单整阶数，我们选取带截距项，滞后差分项选2阶见图7。 图7. 单位根检验 从检验结果看，在1%，5%，10%三个显著水平下，单位根检验的Mackinnon临界值分别为−3.886751，−3.052169，−2.666593，t检验统计值−4.580844小于相应临界值，从而拒绝H 0 ，表明y序列不存在单位根，是平稳序列，即y序列是二阶单整的。Y~II(2) 采用同样的方法，可检验到X1T，X2T，X3T，X4T序列也是二阶单整的。见图8~图12。 X1T： 图8. 单位根检验 X2T： 图9. 单位根检验 X3T 图10. 单位根检验 X4T： 图11. 单位根检验 为了检验回归残差的平稳性，令et = Resid，选择无截距项，无趋势项的ADF检验： 图12. 单位根检验 在5%的显著性水平下，t检验统计量值为−3.323011，小于相应临界值，从而拒绝H 0 ，表明残差序列不存在单位根，是平稳序列，说明YT与X1T，X2T，X3T，X4T之间存在协整关系(修正放在所有检验后面)。 在经典线性回归的五个基本假定(零均值，同方差，无自相关，解释变量与扰动项不相关，无多重共线性)得到满足的情况下，可以使用普通最小二乘法(以下简称OLS)求出 β 0 、 β 1 、 β 2 、 β 3 、 β 4 的估计量。利用Eviews软件中的最小二乘法计算。 估计模型的参数得下结果见图13： 图13. 模型回归结果参数 YT = − 8 . 799122 + 1.319021 * X1T − 0.44159 6 * X2T − 2.8640 19 * X3T + 2 . 56 0 192 * X4T SE = 13.59739             0.180664                       0.185549                         0.758267                         0.843425 T = ( − 0. 647118 )         ( 7 . 3 00 953 )                 ( − 2 . 379937 )               ( − 3 . 777 0 57 )                 ( 3 .0 3547 0 ) R 2 = 0.995998                   R 2 ¯ = 0.994997                 F = 995.3959                     D . W . = 1 . 4 0 618