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

user: 5. 经验和体会 通过本次为时连续9天不中断的细致调查，第一调查队取得了预期成效，累计采集各类样本163份。与此同时，考察队员们有不少的收获和体会，也积累了一些宝贵的经验。 一是项目开展的及时而有效。“第三次全国农作物种质资源普查与收集行动”项目的实施，湖北有幸作为全国第一批启动的省市，钟祥市的普查与收集行动取得较好成效，抢救性地收集到了一批有价值的种质资源。 二是强有力的领导为项目实施提供了强有力的组织保障。在湖北省种子管理局和湖北省农科院的统一领导和协调下，第一调查队组成(由资深的水稻资源专家张再君博士任组长，队员分别来自湖北省农科院和武汉市农业科学院且都是各领域的优秀人才)后人员没有变动，途中保持了队伍的稳定，确保了考察、收集活动时判断的一致性。 三是细致入微的普查工作提高了项目实施的效率。2015年，钟祥市种子管理局组织专班对全境内农作物种质资源进行了系统普查，为第一调查队提供了全面而详实的第一手资料，使得队员们能直接奔赴目的地开展针对性的工作，到任何一地都不会空手而归。 四是强烈的事业心和高度的使命感为项目的实施提供了质量保证。第一调查队的5位队员都是工作在一线的优秀的农业科研工作者，长期从事着水稻、玉米、果树和蔬菜等作物的种质资源收集、创新、育种和推广，都有强烈的事业心和高度的责任感。能参加“第三次全国农作物种质资源普查与收集行动”湖北调查队，大家都感到使命光荣、责任重大。队员们都有各自的专业工作、繁琐的日常事务要处理，都是有家有口的需要照顾，接到任务后都没有丝毫犹豫，义无反顾地奔赴钟祥。在钟祥的9天里，队员们工作认真、负责，每天回到住地已近很辛苦了，但坚持及时整理资料，确保了收集到的基本不出现重复。 五是“懂农业、爱农村、爱农民”的亲情思维保证了项目的顺利实施。第一调查队在钟祥期间，虽然每天都有当地专家、基层干部做向导，但进村入户后，队员们还要直接面对一个个农户。每到一处，面对淳朴的大爷大妈、大哥大嫂，队员们都能端正态度、放下身段，诚恳地向他们嘘寒问暖、认真地向他们讲解政策、虚心地向他们请教、仔细地倾听他们的介绍，如此下来，考察和收集行动就非常顺利。 六是良好的科学素养和扎实的专业知识为项目的实施提供了技术支撑。第一调查队的5位队员不是取得学士学位就是获得博士学位，都是各自学科的优秀人才，日常工种中都体现出了良好的科学素养、有丰富的专业知识储备。在每天的考察、收集行动中，队员们利用专业知识和专业技能进行科学判断，才能去伪存真，才能尽量避免重复，也才能尽量不出错。
system: [SYSTEM]

user: 6. 建议 农业主管部门、科研院所应该加大农作物种质资源的保护力度，尤其是果树资源和野菜资源还应加强原产地保护。要组织科研院所加快已收集种质资源主要生物学特性的鉴定。组织有关高校、科研机构对收集的种质资源进行提纯复壮，开展遗传特性研究和优异基因挖掘，进而为地方种质资源的开发利用打下基础。
system: [SYSTEM]

user: 摘 要 为探讨多效唑在水稻生产中进行喷施的最佳浓度，本研究以2个杂交水稻品种为材料，在田间自然条件下，研究了秧苗期和拔节期叶面喷施不同浓度多效唑对水稻生长及其产量的影响。结果表明，秧苗期叶面喷施450 mg/L多效唑溶液可将水稻移栽时苗高降低30%，单蔸分蘖数增加19.4%以上，且该浓度下单蔸有效穗数、穗总粒数和结实率较大，产量较高。拔节期叶面喷施360~450 mg/L多效唑溶液可使水稻成熟期株高降低4.2%~9.6%，且该浓度下单蔸有效穗数和穗总粒数较大，产量最高。水稻生产上，秧苗期和拔节期叶面喷施多效唑的最佳浓度分别为450 mg/L和360~450 mg/L。 关键词 :水稻，多效唑，产量，生长发育 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. 引言 水稻是世界主要粮食作物之一，是全球近50%人口的主要食物来源，其生产对保障全球粮食安全具有至关重要的意义 [ 1 ]。多效唑是一种高效低毒的植物生长延缓剂和杀菌剂，具有调节植物生长、抑制茎秆伸长、促进分蘖、提高产量和增强抗倒性和抗逆性等多种作用 [ 2 ]。 随着我国人口老龄化到来，农村劳动力人口减少的现象日趋严重，在水稻生产过程中，大力发展轻简化栽培技术是解决农村劳动力减少的重要措施 [ 3 ]。直播、机插等栽培技术是水稻轻简化栽培体系的重要内容，具有省时、省工、省本等诸多优点 [ 4 ]。我国早在20世纪50年代初就开始研究和应用直播稻，近年来，越来越受到稻农青睐，不推自广，面积逐年扩大 [ 5 ] ；虽然我国机插水稻的研究起步较晚，但近年发展较快，至2014年，全国水稻机插秧面积已发展到播种面积的38%，并有迅速扩大发展的趋势 [ 6 ]。倒伏是水稻，尤其是直播稻生产上面临的一大难题。倒伏后，水稻叶片的光合作用、茎秆的运输功能都会受阻，导致结实率下降，空瘪粒增多，品质降低，产量潜力难以发挥，并且增加收获的难度 [ 7 ]。研究表明，喷施多效唑可降低作物株高，提高抗倒伏程度，其作用机理是多效唑抑制赤霉素的合成，提高IAA氧化酶的活性，调节植株体内IAA和ABA的含量，从而达到矮化植株的效果 [ 8 ]。何冲霄等研究表明，在拔节期喷施多效唑可显著降低株高，增加水稻抗倒伏能力 [ 9 ] ；闫凯莉等研究也发现，多效唑对水稻秧苗生长抑制和产量的提高具有明显的作用效果，以400 mg/L的促分蘖和抑制茎秆伸长效果最好 [ 10 ]。机插秧对秧苗的形态及秧苗素质要求较高，而机插秧苗具有播种密度高、生长空间小等特点，这就使得秧苗适栽秧龄期短，秧龄弹性小 [ 11 ]。王长新认为，在机插育秧上运用多效唑可有效控制株高和叶龄，培育出适合机插的30 d秧龄的长龄秧 [ 12 ]。 目前，多效唑在水稻生产中的应用研究主要是关于授粉期的调控和试管苗的离体培养等方面 [ 13 ] [ 14 ] ，而在田间条件下，多效唑在水稻生产中进行喷施的最佳浓度尚无明确定论。本研究在大田自然条件下，以2个不同杂交水稻品种为材料，探讨秧苗期和拔节期叶面喷施不同浓度多效唑对水稻生长发育和产量的影响，以期为水稻轻简化栽培和抗倒栽培技术提供试验依据。
system: [SYSTEM]

user: 2. 材料与方法 2.1. 试点概况 试验于2016年在黄冈市农业科学院梅家墩农场(30˚27'N， 114˚5'E)进行，气候为典型的亚热带大陆性季风气候。试验地前茬空闲，土壤为砂壤土，地势平坦，排灌方便，肥力中等偏上。 2.2. 材料 选用杂交晚稻新品种28优158和巨风优441为材料，均由黄冈市农业科学院水稻所中稻室提供。 多效唑为江苏七洲绿色化工股份有限公司生产的15%可湿性多效唑粉剂。 2.3. 试验设计与方法 试验在大田自然条件下，采用微/小区试验法，在当地习惯种植模式下进行，试验1以28优158为试验材料，试验2以28优158和巨风优441为试验材料。 秧苗期试验：以秧苗期(一叶一心期)叶面喷施不同浓度多效唑(0 (CK)， 300， 450， 600， 700 mg/L)为试验处理，共5个处理，3次重复，完全随机排列。每个微/小区播种面积约0.6 m 2 (0.6 m × 1.0 m)，于播后第7天按设计叶面喷施不同浓度多效唑溶液，喷施量为75 mL/m 2 ，之后，待秧龄达到25天时移栽，每个微/小区移栽面积约6.67 m 2 (约0.01亩，2 m × 3.33 m)。 拔节期试验：以拔节期叶面喷施不同浓度多效唑(试验1：0 (CK)、180、360、540 mg/L；试验2：0 (CK)、450、540、630 mg/L)为试验处理，共4个处理，3次重复，完全随机排列。每个微/小区面积约6.67 m 2 (约0.01亩，2 m × 3.33 m)，于移栽后第22天按试验设计叶面喷施不同浓度多效唑溶液，喷施量为75 mL/m 2 。 秧田管理：试验1和试验2分别于2016年5月10日和2016年6月23日播种育秧。播种前1天，施洋丰牌水稻专用复合肥(含N量15%，含P 2 O 5 量7%，含K 2 O量8%) 300 kg/hm 2 作底肥，播种后第10天和第20天分别追施尿素(含N量46%) 112.5 kg/hm 2 ，其余管理采用当地最优方式进行。 移栽后管理：试验1和试验2分别于2016年6月4日和2016年7月17日移栽，株行距16.7 cm × 20 cm (5 × 6寸)，密度约300，000穴/hm 2 ，每穴1粒谷苗。移栽前1天，施洋丰牌水稻专用复合肥(含N量15%，含P 2 O 5 量7%，含K 2 O量8%) 900 kg/hm 2 作底肥，移栽后第7天追施尿素(含N量46%) 112.5 kg/hm 2 ，移栽后第15天追施尿素(含N量46%) 75 kg/hm 2 和氯化钾(含60% K 2 O) 112.5 kg/hm 2 。水分管理按照深水返青，浅水分蘖，够苗晒田，足水保胎，薄水壮籽，干干湿湿灌浆的管理原则进行。试验期间适时进行人工除草和防治病虫害。 2.4. 数据收集 移栽前1天，每微/小区选取代表性秧苗15株，对苗高(秧苗基部到最长叶叶尖的长度)进行测量；移栽后分别于第15天(试验1)和第20天(试验2)，对每个微/小区内10蔸的分蘖数进行计数；收获前1天，每微/小区选取代表性植株15蔸，对株高进行测量，之后取5蔸分穗进行考种，包括穗总粒数、实粒数和千粒重等，并计算结实率。收获时，各微/小区去边行单收单打，晒干后称取稻谷产量，折算得到单位面积产量。 2.5. 数据分析 应用SPSS 20.0统计软件进行方差分析，采用Duncan’s多重比较检验差异显著性，数据以平均值 ±标准误表示。
system: [SYSTEM]

user: 4. 讨论与结论 秧龄弹性过小是水稻轻简化栽培–机插育秧的一个关键问题，特别是在一年多熟的地区，茬口比较紧张，育秧期间温度高，秧苗生长快，秧苗适宜机插的天数短，林昌明等研究表明，采取化控与控水相结合的半旱育秧，能有效地增加秧龄弹性 [ 15 ] ，王长新研究也发现，在机插育秧上运用多效唑可有效控制株高和叶龄 [ 12 ]。本研究表明，秧苗期叶面喷施450 mg/L多效唑溶液可将水稻移栽时苗高降低30%，这和闫凯莉等的研究结论 [ 10 ] 相近。据试验和大田观察，喷施多效唑可使水稻分蘖数增加，本研究发现，秧苗期叶面喷施450 mg/L多效唑溶液可将水稻单蔸分蘖数增加19.4%以上。秧苗期喷施多效唑对后期水稻产量形成及最终产量影响不显著 [ 11 ] ，本研究发现，秧苗期叶面喷施多效唑对水稻成熟期株高、穗总粒数、结实率、千粒重和产量影响均不显著，喷施450 mg/L多效唑溶液时，单蔸有效穗数、穗总粒数和结实率较大，产量较高。以上说明，秧苗期叶面喷施多效唑对水稻秧苗的生长发育有明显的调节作用，能够显著地抑制秧苗的生长和促进分蘖，而对成熟期株高、穗总粒数、结实率、千粒重和产量影响不显著，且以叶面喷施450 mg/L多效唑的效果最好。 最近数年来，新品种的不断出现，特别是穗大、粒多、粒重的两系杂交品种推广以来，每年水稻抽穗到灌浆及蜡熟期，倒伏率都在8%~20%，严重影响农机收割，水稻产量大幅减少 [ 16 ] [ 17 ]。汤日圣等 [ 18 ] 指出，水稻拔节前后施用多效唑能使节间缩短，茎杆高度降低，从而明显提高防倒效果。本研究表明，拔节期叶面喷施360~450 mg/L多效唑可使水稻成熟期株高降低4.2%~9.6%，且据本试验观察使用过多效唑的田地没有出现水稻倒伏的现象。本研究还发现，拔节期叶面喷施多效唑对水稻成熟期产量构成和产量影响均不显著，当叶面喷施360~450 mg/L多效唑时，单蔸有效穗数和穗总粒数较大，产量最高。以上说明，该试验条件下，水稻抗倒生产上拔节期叶面喷施多效唑的最佳浓度为360~450 mg/L。
system: [SYSTEM]

user: 摘 要 本文从理论上分析了利用分层MC抽样技术可以提高模拟的计算精度，以及推导了分层抽样方差缩减的效果(通过VRER S 值或粗糙估计量方差的比较)，给出了分层MC模拟的具体算法；且根据推导发现分层概率加权方法、分层匹配样本方法与最优分层方法这三种方法中，最优分层方法方差缩减效果比分层匹配样本方法方差缩减效果好。并把所研究的成果应用于欧式看涨期权的定价之中，对于期权理论的研究和完善具有重要的基础作用。 关键词 :Monte-Carlo模拟，分层MC，方差缩减 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. 引言 期权既是高效的风险对冲工具又是高效的投资工具，一直被投资者青睐。1973年著名学者Black与Scholes [ 1 ] 基于一些假定的条件下开创性的推导出标准欧式期权定价的模型。国内外众多学者以及金融业界人士就期权定价问题开展了广泛的研究，但随着研究的深入人们发现只有在特定的条件下一些期权才具有解析解，期权都没有解析解或是解析解求解非常困难，如：美式期权、亚式期权等这类特殊的期权是没有解析解的。因此。在大多数的情形下，期权定价一般需要寻求数值解法；目前，期权定价的数值求解有：二叉树 [ 2 ] 、二项式模型 [ 3 ] 、有限差分法(隐式或显式)以及随机模拟 [ 4 ] 等方法。尤其是其中的蒙特卡洛模拟法由于其具有良好的适用性引起人们特别的重视。尤其是在实际应用中，如多标的资产、复合基础资产、标的资产既带跳跃又有随机波动率的条件且期权模型又要考虑与标的资产路径相关 [ 5 ] [ 6 ] ，蒙特卡罗方法则可以大显身手。 但是蒙特卡罗模拟方法自身也有不足之处，就是模拟结果的精度、波动性等受到模拟的路径数量的重大影响 [ 7 ] [ 8 ] ，所以纵观蒙特卡罗模拟方法在金融衍生物定价的完善与发展主要体现在模拟估计误差减少技术与扩大其应用范围。模拟估计误差减少技术主要是针对模拟过程所产生的离散化误差与统计误差，离散化误差即由样本路径模拟过程所导致，而统计误差即由于有限样本及不同的抽样方法所引起的估计误差。其中，统计误差可以通过方差缩减技术来解决，离散化误差则通过采取有效的高阶收敛性的样本路径模拟的方式来逼近 [ 9 ] [ 10 ] 。下文将重点介绍分层抽样法结合Monte-Carlo模拟，在提高模拟精度技术的数学原理基础上重点说明其在期权模拟定价上的使用，然后给出精确的模拟误差与精度值。
system: [SYSTEM]

user: 2. 分层Monte-Carlo方法 2.1. 分层Monte-Carlo方法的产生 在我们以下的讨论模型中，均以欧式看涨期权 S ( T ) = S ( t ) exp { ( r − 1 2 σ 2 ) ( T − t ) + σ ( B ( T ) − B ( t ) ) } 为例，其中T风险价格到期日， S ( t ) 表示当前时刻风险资产价格， B ( T ) − B ( t ) ∼ N ( 0 ， T − t ) ，r无风险利率， σ 为年波动率。对风险资产在T时刻的模拟常常借助下面的代码： 该模拟显然是要借助正态随机数代码的，但是正态随机数的产生过程容易产生不均衡现象，也就是说： 通过上面的代码可以验证，在产生的10，000个标准正态随机数中，n1，n2的个数很难保证是1:1，有时差距还很大。 基于此，统计学家们构造出一种新的随机数产生方案，将密度函数分成有限个区间，然后在每个区间上按照各自的概率比重产生一定数量的随机数，最后用这些数进行Monte-Carlo模拟，该方法被称为分层Monte-Carlo方法(以下称为分层MC方法)，最终通过在每个小区间上的抽样来实现提高模拟精度这一目的。 2.2. 分层抽样原理 分层抽样的原理是对所输入数据进行处理，使其更加规律，从而使经验概率与理论概率更加匹配(Glasserman和Gaussia [ 11 ] )。下面通过一个简单的推导来说明分层抽样的性质。 假设拟估计： Y = E [ f ( U ) ] = ∫ 0 1 f ( y ) d y = μ 其中 U ∼ U i n f [ 0 ， 1 ] ，假设对样本分为两层，即： Y S = 1 2 f ( V 1 ) + 1 2 f (V2) 其中 V 1 ∼ U n i f [ 0 ， 1 2 ]                                                           V 2 ∼ U n i f [ 1 2 ， 1 ] 现在比较变量Y的粗糙估计量 Y c r u d e 和分层抽样估计量 Y S 之间的关系： Y ^ c r u d e = 1 n ∑ i = 1 n f ( μ i )                                                       Y ^ S = 1 n / 2 ∑ i = 1 n / 2 Y S i = 1 n [ ∑ i = 1 n / 2 ( f ( V 1 i ) + f ( V 2 i ) ) ] (1) 接下来(1)式期望和方差为： E [ Y ^ S ] = 1 n [ ∑ i = 1 n / 2 f ( V 1 i ) + f ( V 2 i ) ] = ∫ 0 1 f ( y ) d y = μ (2) V a r [ f ( V 1 i ) ] = ∫ 0 1 2 2 f 2 ( x ) d x − [ ∫ 0 1 2 2 f ( x ) d x ] 2 ≡ 2 ν 1 − 4 m 1 2 (3) V a r [ f ( V 2 i ) ] = ∫ 1 2 1 2 f 2 ( x ) d x − [ ∫ 1 2 1 2 f ( x ) d x ] 2 ≡ 2 ν 2 − 4 m 2 2 (4) 其中 ν 1 = ∫ 0 1 2 f 2 ( x ) d x                                                                   m 1 = ∫ 0 1 2 f ( x ) d x ν 2 = ∫ 1 2 1 f 2 ( x ) d x                                                                     m 1 = ∫ 1 2 1 f ( x ) d x 根据公式(3)、(4)得： V a r [ f ( V 1 i ) ] + V a r [ f ( V 2 i ) ] = 2 ( ν 1 + ν 2 ) − 4 ( m 1 2 + m 2 2 ) = 2 ( ν 1 + ν 2 ) − 2 ( m 1 + m 2 ) − 2 ( m 1 − m 2 ) 2 = 2 V a r ( f ) − 2 ( m 1 − m 2 ) 2 (5) 根据公式(1)、(5)，推导出公式(6)： V a r [ Y ^ S ] = 1 n 2 [ ∑ i = 1 n / 2 V a r ( f ( V 1 i ) ) + V a r ( f ( V 2 i ) ) ] = 1 n [ V a r ( f ) − ( m 1 − m 2 ) 2 ] (6) V a r [ Y ^ c r u d e ] = V a r ( f ) / n                                 V a r [ Y ^ c r u d e ] = V a r [ Y ^ S ] − ( m 1 + m 2 ) 2 / n (7) 整理公式(7)，得： V a r [ Y ^ S ] − V a r [ Y ^ c r u d e ] = ( m 1 + m 2 ) 2 / n (8) 假设 f ( U ) = e U                                                                                               n = 2 E [ f ( U ) ] = 1.718                                     V a r [ f ( U ) ] = 0.240 根据公式(6)、(7)计算，得： Y ^ c r u d e = ( e u 1 + e u 2 ) / 2                                                             V a r [ Y c r u d e ] = 0.121 m 1 = 0.6487                                                                                                                             m 2 = 1.069 V a r [ Y ^ S ] = V a r [ Y ^ c r u d e ] − ( m 1 + m 2 ) 2 / n = 0.0375 则 V R E R S = 3.227 从方差缩减比率 V R E R S 的值可看出分层抽样方差缩减技术的效果比较明显。
system: [SYSTEM]

user: 3. 均匀分布与非均匀的分层 根据第1节我们了解当分层层数为 n = 2 的情形，方差缩减的大小为(8)式。接下来我们将研究当分层层数 n > 2 的情形时均匀分布和非均匀分布分层抽样的估计量。方差缩减的大小在后面的章节给出证明。 3.1. 均匀分布的分层 假设把[0，1]分为n层，记第一层为 A 1 = [ 0 ， 1 / n ) ，第二层为 A 2 = [ 1 / n ， 2 / n ) ，依次类推， A n = [ ( n − 1 ) / n ， 1 ) ， 在均匀分布的条件下，落在每一个小区间的概率为 1 n ，假设样本空间为n，定义如下： V i = i − 1 + U i n ∼ U n i f [ i − 1 n ， i n ] 其中 U i ∼ U n i f ( 0 ， 1 ) ， i = 1 ， 2 ， ⋯ ， n 。因此 V i 为给定 U ∈ A i 时U的条件分布， V 1 ， V 2 ， ⋯ ， V n 构成了(0，1)均匀分布的一个分层样本。如果 Y = f ( U ) ，则有 E [ f ( U ) ] 的分层抽样估计为： Y S = 1 n ∑ i = 1 n f (Vi) 3.2. 非均匀分布的分层 给定概率 p 1 ， p 2 ， ⋯ ， p K ，       ∑ i = 1 K p i = 1 ，并定义K个分层分别为： A 1 = ( a 0 ， a 1 ) ， A 2 = ( a 2 ， a 3 ) ， ⋯ ， A K = ( a K − 1 ， a K ) 其中 a 0 = − ∞ ， a 1 = F − 1 ( p 1 ) ， a 2 = F − 1 ( p 1 + p 2 ) ， ⋯ ， a K = F − 1 ( p 1 + p 2 + ⋯ + p K ) ，F为概率分布函数。假设随机变量Y的概率分布函数为F，则有 P r ( Y ∈ A i ) = F ( a i ) − F ( a i − 1 ) = p i 。为了利用所产生的分层抽样 A 1 ， A 2 ， ⋯ ， A K ，需要生成随机变量 Y | Y ∈ A i 的样本，利用逆变换方法可实现。假设 V = a i − 1 + U ( a i − a i − 1 ) ，则 F − 1 ( V ) 具有 Y | Y ∈ A i 分布。
system: [SYSTEM]

user: 5. 总结 本文详细地分析了利用分层MC抽样技术模拟提高其计算精度与运行效率，及推导了分层抽样方差缩减的效果(也就是方差缩减比率 V R E R S 的值)，而且我们根据推导发现分层概率加权方法、分层匹配样本方法与最优分层方法这三种方法中，而且我们从(30)式看出最优分层方法方差缩减效果比分层匹配样本方法的方差缩减效果好。并把所研究的成果应用于欧式看涨期权的定价之中，其对于期权理论的研究和完善具有重要的基础作用。作用分析各种理财产品现值，为投资者和相关人员机构提供一种较好的与预测工具和手段。
system: [SYSTEM]

user: 摘 要 当飞行器以高超声速飞行时，其前方的气体剧烈压缩和粘性阻滞引起热障问题。由于高超声速流中非平衡效应对气动热的影响很大，化学非平衡效应成为高超声速气动热领域研究的关键问题之一。本文简要介绍了非平衡效应和真实气体效应，在此基础上，综述了化学非平衡模型下气动热数值模拟的研究进展，包括真实气体效应、化学非平衡效应、化学反应模型、近壁面第一层网格高度对气动热数值模拟结果的影响，以及化学非平衡流条件下圆球高超声速绕流的激波脱体距离随马赫数的变化规律，该规律有待相关风洞试验的验证。文章最后指出湍流对化学非平衡效应的影响及辐射效应对气动热数值模拟的影响是值得关注的研究方向。 关键词 :高超声速，气动热，化学非平衡流，网格划分，激波脱体距离 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. 引言 高超声速飞行器一般指飞行马赫数大于5的飞行器 [ 1 ] [ 2 ] [ 3 ] [ 4 ] 。当飞行器以如此高速度飞行时，强激波会包裹在飞行器周围并产生剧烈的气动加热。气动热结果的准确预测对高超声速飞行器的热防护系统(TPS)的设计至关重要 [ 5 ] 。气动热的研究方法分为三种，即飞行试验、地面试验和数值模拟。飞行试验能够提供充分可靠的信息，例如：研究地球大气层的物理特性、高超声速钝头体再入时的热环境和电子分布的RAM-CII飞行试验 [ 6 ] ，研究热流密度和壁面温度分布的Fire II飞行试验 [ 7 ] ，X-43飞行试验 [ 8 ] ；但是，飞行试验需要花费大量的人力、财力和物力。地面试验可模拟部分空中飞行，美国的CALSPAN的LENS系列激波风洞和膨胀管风洞代表最高水平 [ 9 ] [ 10 ] ，但不同的地面设备满足不同的要求，无法提供飞行全程的信息。上世纪八十年代以来，计算流体力学方面的长足进步和计算机技术的迅猛发展，使得数值模拟方法弥补了以上不足，采用CFD数值模拟方法成为预测高超声速飞行器气动特性的有效手段。随着研究的不断深入，高超声速飞行器的气动数值模拟经历了边界层方程(BL)、粘性激波层方程(VSL)、抛物化方程(PNS)和全NS方程的求解过程，其关键的气体状态方程从量热完全气体模型、化学平衡流模型、化学非平衡流模型、热力学非平衡流模型等不断深化，气体状态方程越来越接近实际情况。 飞行器表面的材料特性、飞行器的形状和飞行的弹道特性是影响气动热的主要因素。高温下，飞行器表面的防热材料对边界层内气体的化学反应的催化特性不同，不同的壁面催化条件使得气动热的变化十分巨大。飞行器的形状会影响流动特征和激波形状，特别是复杂外形高超声速飞行器的舵、翼、防热瓦缝隙等局部区域的流动状态非常复杂，对气动热有显著的影响。飞行的弹道特性包括飞行的速度、飞行的高度和飞行的姿态等。在高马赫数下，高温真实气体效应影响显著；在不同的飞行高度上，高超声速飞行器的热环境是完全不同的，由高空到低空气体处于自由分子流、过渡流和连续流状态，这就需要对不同高度下的高超声速飞行器热环境进行预测和计算；来流雷诺数还会使得气体存在层流到湍流的转捩现象，而湍流边界层的气动热更加严重；另外，当飞行高度很高时气体还会受到热力学非平衡的影响，低密度会降低气体粒子之间的碰撞频率，而高速度会缩短空气分子运动的时间，这些过程使得气体分子的平动、转动、振动和电子能激发等模式的热力学能产生一个滞后于平衡状态的效应。随着飞行高度的降低，例如70 km以下，热力学非平衡影响并不明显，也就是说，气体组分的所有内能模式是与重粒子的平动能模式处于平衡的。对于各类高超声速飞行器的绕流流动，仅需要模拟化学非平衡和粘性的影响 [ 11 ] - [ 16 ] 。如上所述的因素都与高超声速气动热及化学非平衡效应密切相关，本文综述了化学非平衡模型下气动热数值模拟的研究进展，包括真实气体效应、化学非平衡效应、化学反应模型、近壁面第一层网格高度以及激波脱体距离。
system: [SYSTEM]

user: 4. 高超声速环境下激波脱体距离的数值模拟 在高超声速流动中，考虑真实气体效应后最明显的就是激波脱体距离的变化 [ 43 ] [ 44 ] [ 45 ] [ 46 ] [ 47 ] 。高超声速飞行器头部的激波脱体距离与激波层内气体的热化学状态密切相关。高温气体的化学反应使激波后气体的密度增加，而激波脱体距离与驻点线上的平均密度成反比，因此激波脱体距离能够受到明显的非平衡现象的影响，其大小可以一定程度上反映出化学反应的激烈程度 [ 48 ] [ 49 ] [ 50 ] [ 51 ] 。图21为激波脱体距离受化学反应影响而发生变化的示意图，其中ΔP、ΔCR分别表示量热完全气体模型和化学反应模型的激波脱体距离 [ 52 ] 。对于高超声速流动的CFD数值模拟来说，激波脱体距离是常用于验证结果准确性的参数之一 [ 53 ] - [ 61 ] 。 图21. 化学反应对激波脱体距离的影响 钝体绕流是高超声速流动中的基础问题。在实验中，脱体激波是较直观且易测量的，图22为Nonaka实验中脱体激波的照片 [ 62 ] 。脱体激波和物面之间的区域为激波层，二者之间的距离即激波脱体距离。激波脱体距离是高超声速流动中一个很重要的参数，这是因为激波位置和激波形状直接决定了钝体表面的压强分布，从而决定了飞行器的气动力特性，且深入研究激波脱体距离有助于高超声速飞行器热防护系统的设计 [ 63 ] [ 64 ] [ 65 ] 。 图22. Nonaka实验脱体激波照片 [ 62 ] 胡雨濛以5组分化学非平衡流模型为基础，来流条件选择30 km~45 km处的大气条件，具体数值见表4，马赫数从24变化到34，来流中O 2 和N 2 的质量分数分别为0.21和0.79。 Table 4 Parameters of the free stream at different altitude H km ρ ∞ kg/m 3 T K μ ∞ ×10 − 5 kg/m∙s c m/s 30 0.01841000 226.5 1.4753 301.7 35 0.00846340 236.5 1.5287 308.3 40 0.00399570 250.3 1.6009 317.2 45 0.00198650 264.2 1.6712 325.8 表4. 不同高度下的来流条件 图23给出了当高度从30 km增加到45 km时，化学非平衡流与化学平衡流模型模拟的激波脱体距离随着马赫数增加的对比结果。在化学非平衡流动中，当大气层高度、钝体半径和来流密度相同时，激波脱体距离依然遵循随着马赫数的增加先减小后增大的变化规律，最小值用大图标标记。其中，在45 km的情况中，28与29 Ma计算出的激波脱体距离值大致相同，可以认为最小激波脱体距离对应的马赫数(记为“拐点马赫数”)在这二者之间。当马赫数低于拐点马赫数时，激波脱体距离随着马赫数的增加很快减小，高于拐点马赫数时，激波脱体距离会慢慢变大。在相同的大气高度处，非平衡流的最小激波脱体距离对应的马赫数均比平衡流的大。这是因为非平衡效应使得化学反应以有限速率进行，气体经过激波后的碰撞次数减少，气流需要经过一段距离才达到平衡状态，这种不稳定的状态使得气体仍然能够被继续压缩，即在较高马赫数下才能达到激波脱体距离的极限。在相同的马赫数下，平衡流中的部分热能用于化学反应和气体能量模式的激发，密度受化学反应的影响较大，平衡流激波后的密度大于非平衡流的，因此化学平衡流下的激波脱体距离都比化学非平衡流的小，这是高温非平衡效应下有限速率使化学反应进行的不充分造成的。由图可知，当马赫数相同时，随着大气高度的增加，平衡流与非平衡流计算的激波脱体距离的差值在变大，因为随着大气高度的增加气体的密度降低，气体分子之间的碰撞频率减少，非平衡效应会更加明显，在相同来流条件下达到平衡态所用的时间有所增加。在未来近空间高超声速飞行器的发展过程中，随着飞行器设计速度的增加，引起的激波脱体距离的新变化会影响到飞行器的整体设计，这种与以往认知相悖的现象是非常值得关注的，目前的风洞试验无法模拟如此高马赫数的来流条件，该规律有待后续马赫数达20以上的风洞试验验证。 图23. 不同高度处，随着马赫数增加，化学非平衡流和化学平衡流模型计算的激波脱体距离的对比
system: [SYSTEM]

user: 5. 结束语 “真实气体效应”是高超声速飞行器设计过程中必须考虑的问题，由于国外对我国的技术限制，自主研究和发展化学非平衡流数值模拟程序是至关重要的。本文综述了化学非平衡效应以及化学非平衡效应对计算结果的影响，简要如下： 1) 考虑真实气体效应后，激波层内温度降低，密度增大，而速度和压强与量热完全气体区别不大。 2) 化学非平衡激波后的温度高于化学平衡流的，且随着气体向近壁面流动，温度逐渐下降，密度逐渐增加，呈现一个明显的非平衡变化过程，而压强和速度与化学平衡流的区别不大。 3) 对于化学反应模型的选取应依照实际研究的问题，以避免浪费计算时间。 4) 第一层网格高度对壁面热流结果影响巨大，在计算条件允许下，应尽可能采用小的第一层网格高度。 5) 化学平衡下的与化学非平衡下的激波脱体距离均是随着马赫数的增加先减小后增大，存在激波脱体距离变化的拐点，但该现象有待后续试验验证。 6) 化学非平衡下的最小激波脱体距离出现的马赫数比平衡流下的大，说明非平衡条件下气体的可压缩性增加。 从本文综述可见，目前的高超声速气动热的数值模拟主要集中在钝体绕流的层流模型，没有考虑湍流对化学非平衡产生的影响，同时，数值模拟结果显示流场气体温度达到9000 K，需要考虑辐射对气动热的影响，这些方面的理论、实验研究是值得关注的方向。
system: [SYSTEM]

user: 摘 要 本文给出了向量值Bergman空间上块Toeplitz算子T F T G = 0的一个充分必要条件，其中F为一般本性有界向量值函数，G为调和多项式。 关键词 :向量值Bergman空间，块Toeplitz算子，零积，Mellin变换 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. 引言 记IN为自然数集，Z为整数集。定义D是复平面C上的开单位圆盘。dA表示D上正规化的面积测度，即 d A ( z ) = r π d r d θ = 1 π d x d y . L 2 ( D ， d A ) 是D上关于dA平方可积的函数全体构成的Hilbert空间，可表示为： L 2 ( D ， d A ) = { f : ∫ D | f ( z ) | 2 d A ( z ) < + ∞ } ， 对任意的 f ， g ∈ L 2 ( D ， d A ) ， L 2 ( D ， d A ) 空间中内积为： 〈 f ， g 〉 = ∫ D f ( z ) g ( z ) ¯   d A ( z ) . Bergman空间 L a 2 ( D ， d A ) 为 L 2 ( D ， d A ) 上由全体解析函数构成的闭子空间。定义D上由关于dA的本性有界可测函数全体构成的Banach空间为 L ∞ ( D ， d A ) ，P是从 L 2 ( D ， d A ) 到 L a 2 ( D ， d A ) 上的正交射影，则有 ( P f ) ( z ) = 〈 f ， K z 〉 = ∫ D f ( w ) K z ( w ) ¯   d A ( w ) ， K z ( w ) = 1 ( 1 − z ¯ w ) 2 ( z ， w ∈ D ) 为 L a 2 ( D ， d A ) 的再生核，令 k z ( w ) = K z ( w ) ‖ K z ‖ 为正规化的再生核。 定义1.1：设 λ ∈ L ∞ ( D ， d A ) ，以 λ 为符号的Toeplitz算子定义如下： T λ f = P ( λ f ) ，     f ∈ L a 2 ( d ， d A ) . 说明：设 λ ， μ ∈ L ∞ ( D ， d A ) ，则Toeplitz算子具有以下性质： 1) ∀ a ， b ∈ C ， T a λ + b μ = a T λ + b T μ ， 2) T f ∗ = T f ¯ 。 记 M n × n 为C上的 n × n 矩阵的全体。定义 L 2 ( D ， C n ) = L 2 ( D ， d A ) ⊗ C n 为D上关于dA平方可积的向量值函数空间。 L a 2 ( D ， C n ) = L a 2 ( D ， d A ) ⊗ C n 为D上的向量值Bergman空间，其中“ ⊗ ”为Hilbert张量积。以 Ψ ( z ) ∈ L ∞ ( D ， d A ) ⊗ M n × n 为符号的Toeplitz算子定义为 T Ψ h = P L a 2 ( D ， C n ) ( Ψ h ) ，   h ∈ L a 2 ( D ， C n ) ， 其中 P L a 2 ( D ， C n ) 是从 L 2 ( D ， C n ) 到 L a 2 ( D ， C n ) 的正交射影。 “ ⊕ ”表示直和 L ∞ ( D ， C n ) = L ∞ ( D ， d A ) ⊕ ⋯ ⊕ L ∞ ( D ， d A ) ， L a 2 ( D ， C n ) = L a 2 ( D ， d A ) ⊕ ⋯ ⊕ L a 2 ( D ， d A ) ，那么块Toeplitz T Ψ 有下列矩阵表示： T Ψ = [ T Ψ 11 ⋯ T Ψ 1 n ⋮ ⋱ ⋮ T Ψ n 1 ⋯ T Ψ n n ] ， 其中 Ψ = [ ψ 11 ⋯ ψ 1 n ⋮ ⋱ ⋮ ψ n 1 ⋯ ψ n n ] . 块Toeplitz算子和块Toeplitz算子的截断行列式存在于在数学和物理学的各个分支中，并有着重要的应用，尤其是在量子力学的理论研究中。例如，经典的二聚体模型 [ 1 ] 和不相交行走模型 [ 2 ] 的研究应用了块Toeplitz算子的截断行列式理论；盖尔芬德迪基层次的研究 [ 3 ] 应用了块算子理论等。 Brown与Halmos证明在单位圆周上的Hardy空间上，若有 T f T g = 0 ，则符号f或g必定有一个为零 [ 4 ] 。1984年，美国数学家Axler访问四川大学时提出一个问题，若Hardy空间上n个Toeplitz算子的乘积为零，是否必定有其中一个为零? 1994年，郭坤宇 [ 5 ] 证明了 n = 5 时结论是对的。Gu证明结论对于 n = 6 也对 [ 6 ] 。 n > 6 的情况后来也被解决 [ 7 ] 。Bergman空间上的函数论与Hardy空间不一样，因此Bergman空间上的算子理论与Hardy空间上的算子理论也不一样。即使在D上的Bergman空间上，有些看来非常简单的问题解决起来也很困难。例如具有有界调和函数符号的两个Toeplitz算子的乘积为零是否必有一个为零?这个问题就悬置了很长时间。直到2001年才由Ahern和Čučković解决并发表在文 [ 8 ] 中。文 [ 8 ] 证明了若f和g是有界调和函数，使得 T f T g = 0 ，则f或g恒等于零。近年来，许多学者开始研究多圆盘和单位球上的Bergman空间。例如，卢玉峰等 [ 9 ] [ 10 ] 分别研究了多圆盘和单位球Bergman空间上的有界Toeplitz和Hankel乘积。于涛和朱若卫研究了单位球向量值Bergman空间上的Toeplitz乘积有界的充要条件 [ 11 ] 。 本文通过Mellin变换研究Bergman空间上两个Toeplitz算子乘积的有限和为零问题，并以此为基础研究向量值Bergman空间上两个块Toeplitz算子的零积问题。 定义1.2：对 φ 是 [ 0 ， 1 ] 上的可积函数，可定义Mellin变换： φ ^ ( z ) = ∫ 0 1 φ ( r ) r z − 1 d r ， 此时，它是半平面 { z : Re z > 1 } 上的有界解析函数。其中 L 2 ( D ， d A ) 可分解为： L 2 ( D ， d A ) = ⊕ k ∈ Z e i k θ ℜ ，   ℜ = { u : D → C : ∫ 0 1 r | u ( r ) | 2 d r < ∞ } ， 因此对任意 f ∈ L 2 ( D ， d A ) ，都有 f ( r e i θ ) = ∑ k = − ∞ ∞ f k ( r ) e i k θ ，     f k ∈ ℜ . 引理1.1 [ 8 ] ：若 f ∈ L 2 ( D ， d A ) ，对 z = R e i θ ，则可定义f的Berezin变换 B f ： ( B f ) ( z ) = 〈 f k z ， k z 〉 = 2 ( 1 − R 2 ) 2 ∑ k = − ∞ ∞ R | k | [ ∑ n = 1 ∞ n ( n + | k | ) f ^ k ( 2 n + | k | ) R 2 ( n − 1 ) ] e . i k θ 引理1.2 [ 12 ] ：若 φ 为 { z : Re z > 1 } 上的有界解析函数，若存在自然数序列 { n k } k = 0 ∞ 使得 φ ^ ( n k ) = 0 且 ∑ k = 0 ∞ 1 n k = ∞ 则 ϕ = 0.
system: [SYSTEM]

user: 3. 向量值Bergman空间上两个块Toeplitz算子的零积 定理3.1：假设 F = ( f p q ) N × N ∈ L ∞ ( D ， d A ) ⊗ M N × N ， G = ( δ p q ) N × N = A z j + B z ¯ l ，其中 A = ( α p q ) N × N ∈ M N × N ， B = ( β p q ) N × N ∈ M N × N ， α p q ， β p q ∈ C ， j ， l ∈ I N ，且满足 ( A ′ + B ′ ) F ′ = 0. 那么则有 T F T G = 0 当且仅当 A ′ F ′ = 0 ( A ′ 表示矩阵A的转置)。 证明：必要性。由 F = [ f 11 f 12 ⋯ f 1 N f 21 f 22 ⋯ f 2 N ⋮ ⋮ ⋱ ⋮ f N 1 f N 2 ⋯ f N N ] ，     G = [ δ 11 δ 12 ⋯ δ 1 N δ 21 δ 22 ⋯ δ 2 N ⋮ ⋮ ⋱ ⋮ δ N 1 δ N 2 ⋯ δ N N ] ， 得 T F = [ T f 11 T f 12 ⋯ T f 1 N T f 21 T f 22 ⋯ T 2 N ⋮ ⋮ ⋱ ⋮ T f N 1 T f N 2 ⋯ T f N N ] ，     T G = [ T δ 11 T δ 12 ⋯ T δ 1 N T δ 21 T δ 22 ⋯ T δ 2 N ⋮ ⋮ ⋱ ⋮ T δ N 1 T δ N 2 ⋯ T δ N N ] ， 于是有 T F T G = [ T f 11 T f 12 ⋯ T f 1 N T f 21 T f 22 ⋯ T 2 N ⋮ ⋮ ⋱ ⋮ T f N 1 T f N 2 ⋯ T f N N ] [ T δ 11 T δ 12 ⋯ T δ 1 N T δ 21 T δ 22 ⋯ T δ 2 N ⋮ ⋮ ⋱ ⋮ T δ N 1 T δ N 2 ⋯ T δ N N ] = [ T f 11 T δ 11 + T f 12 T δ 21 + ⋯ + T f 1 N T δ N 1 ⋯ T f 11 T δ 1 N + T f 12 T δ 2 N + ⋯ + T f 1 N T δ N N ⋮ ⋱ ⋮ T f N 1 T δ 11 + T f N 2 T δ 21 + ⋯ + T f N N T δ N 1 ⋯ T f N 1 T δ 1 N + T f N 2 T δ 2 N + ⋯ + T f N N T δ N N ] = 0. 进而，对 1 ≤ p ， q ≤ N ，有 ∑ m = 1 N T f p m T δ m q = 0. 由 ( A ′ + B ′ ) F ′ = 0 ，得 ∑ m = 1 N ( α q m + β q m ) f m p = 0   ( 1 ≤ p ， q ≤ N ) ，运用定理1.1的结论得 ∑ m = 1 N α q m f m p = 0 对 1 ≤ p ， q ≤ N 成立。即 A ′ F ′ = 0. 充分性。由 ( A ′ + B ′ ) F ′ = 0 ，得 ∑ m = 1 N ( α q m + β q m ) f m p = 0   ( 1 ≤ p ， q ≤ N ) ，又 A ′ F ′ = 0 ，得到 ∑ m = 1 N α q m f m p = 0 ，     ∑ m = 1 N β q m f m p = 0 . 于是 ∑ m = 1 N T f m p T δ q m = ∑ m = 1 N T f m p ( T α q m z j + β q m z ¯ l ) = ∑ m = 1 N T α m q f m p T z j + T β m q f m p T z ¯ l = T ∑ m = 1 N α m q f m p T z j + T ∑ m = 1 N β m q f m p T z ¯ l = 0. 所以 ( T F T G ) ′ = 0 ， 故 T F T G = 0 ，定理证毕。
system: [SYSTEM]

user: 摘 要 笔者在洛南县多年的竹子生产实践中，基本摸清了竹子的生物学特性，从中总结出了竹子的栽植及管理技术要点，对竹子生产有一定的指导作用。 关键词 :竹子，栽植技术，管理 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. 引言 洛南县位于陕西省商洛市东北部，地处秦岭东段南麓，南洛河上游，属黄河水系，南北高，中间低，以山地为主，川原皆具，是一个“八山半水分半田”的土石山区。东西海拔670~2646米，平均海拔为1500米左右。属暖温带季风性湿润气候，气候差异较大四季分明，春季干旱多风，夏季炎热多雨；秋季多雨湿润，气温下降较快；冬季寒冷少雪。年平均气温为11.1℃，1月极端最低气温−18℃，7月份极端最高气温37.1℃。全县日平均气温 ≥ 0℃的活动积温4152.3℃；日平均气温 ≥ 10℃的活动积温3453.6℃。年平均日照时数为2045.0小时。竹子作为洛南的乡土林种得到了当地政府和人民的高度重视，2015年庵巡路栽植31.2公里102万株，2016至今已栽植180万株。因缺乏栽植经验及技术，导致栽植成活率、保存率、繁殖率低。对此，笔者通过4年时间的调查研究分析摸清了洛南县竹子的生物学特性，掌握了气候条件、土壤因子对竹子生长发育的影响，通过严把竹子的种苗选择、栽植方法、栽植时间、栽后管理4个关口，促其成活率、保存率、繁殖率大幅提高。
system: [SYSTEM]

user: 2. 种苗的选择 选择洛南当地1~2年生生长健壮、颜色深绿，无病虫害感染和机械损伤，分支较低，枝叶繁茂，鞭色鲜黄，鞭芽饱满，鞭根健全的全株，地径大于2 cm，高度3米以上的竹子作为种苗 [ 1 ] 。
system: [SYSTEM]

user: 3. 栽植 3.1. 起苗 起苗时，首先要视种源墒情提前2~3天向种源林地灌足透水，确保渗水深度达30 cm以上。其次要确定竹鞭的方向，一般竹鞭走向多数和第一枝盘方向一致。起苗时，在距母竹40 cm处轻轻挖开土层，找到竹鞭，再沿竹鞭方向呈椭圆形挖好，来鞭留30 cm，去鞭留40 cm。也可不予断鞭，保留长鞭连株起运(无竹鞭或鞭不足20 cm的竹子不能调运)。第三，减少机械损伤，要保护好竹鞭，尽量做到少伤竹鞭、鞭根和笋芽，更不能强行用力摇动种竹，以免损伤连接点。第四，起苗时要尽量多带“土球”，挖出的多株连体竹不能分开，确保土球完整。第五，及时遮阴，用遮阴网进行遮阴减少水分散失，保证正常生理代谢 [ 2 ] 。 3.2. 种苗运输 种苗采挖做到随挖随运，不得滞留到第二天。装运时要保护好竹鞭鞭芽与“螺钉”，要用如稻草、塑料薄膜等具有蓄水保湿作用的材料包扎好竹蔸和竹鞭。装卸时要轻拿轻放，使其枝叶舒展以防枝叶受损。长途运输时根系必须包扎，种苗要遮阴，中途进行洒水，做好通风降温等护理工作。 3.3. 栽前技术处理 种苗运到栽植地点后依据规划要求按照种苗优劣、高度进行分级栽植。如不受当时栽植效果影响均进行定干截枝栽植，其定干点位为基枝以上3~5层枝，且距下竹节15 cm左右处，然后灌足水，用黄泥封口，再用塑料袋或胶带包扎封护；枝上留3~5个叶进行截枝。 3.4. 整地 栽植竹子宜采用开沟整地，单行栽植时沟深、宽为60 × 60 cm，双行栽植可挖深宽为60 × 100 cm的条形沟，单株栽植可挖60 × 60 cm鱼鳞坑。对土壤瘠薄、石砾较多的胶泥土(粘土)要进行因地更换表层壤土或沙壤土，尽量避免更换犁底层土壤。山坡或立地条件差的可采用鱼鳞坑及“抽槽换土”的方法进行整地挖坑 [ 2 ] 。 3.5. 栽植方法 竹子要做到随起随栽，深挖沟浅栽竹，下拥紧(土)上松盖(土)，以栽植后鞭根舒展为原则。一般要求栽植前先垫15 cm种植土(壤土或黄棕壤土的表皮层土壤)，有条件的可以在栽植穴底部铺施一层10~15 cm有机肥，覆土部位应比地表稍高3~6 cm，踩踏后再浇足定根水，待水渗入后加盖一层松土，以减少水分蒸发。
system: [SYSTEM]

user: 4. 栽植时间 因受气候、海拔差距的影响，巡检镇、洛源镇当年5~8月、11月至翌年土壤解冻前不能栽植，其余区域当年9月至翌年4月底前均可栽植。
system: [SYSTEM]

user: 5. 栽后管理 竹子栽植后要做到“三分栽七分管”，这是提高成活率、保存率、繁殖率的关键。 5.1. 支护 栽植后要及时利用不合格或损坏的竹杆横向在竹子主干距下部1/3处用扎丝逐株扎绑，并且将竹子主干拨摆匀称。 5.2. 浇水 栽植后要及时浇水，做到“首浇透，三补水，七补墒”，即：首次要浇透水，确保土壤与根系充分接触不留缝隙及时供应养(水)分，三天、七天后要进行两次补水保障生理代谢。天气过于干燥时要向枝叶喷水，直到确定成活之后，再逐步减少或停止。 5.3. 抚育 防止杂草与竹子争夺养(水)分，要视杂草生长情况及时进行松土除草，土壤贫瘠的要适当追肥，补充养分。 5.4. 管护 对人畜活动频繁的区域要做好防护措施，建立护竹制度，严禁人畜入林践踏；固定专人落实责任，做好管护巡查。
system: [SYSTEM]

user: 6. 结语 本文通过对竹子的生物学特性的研究分析，结合洛南当地的气候特征，因地制宜，并对起苗、运输、栽植、管理等环节进行可行性分析、评估、研究，最终总结出指导洛南县竹子全程栽植的科学方法，该方法不仅使洛南县竹子栽植的成活率、保存率、繁殖率大幅提升，同时该理念、方法或措施对周边区县乃至西北同纬度地区的竹子栽植、园林绿化技术也有一定的指导和推广作用。
system: [SYSTEM]

user: 摘 要 随着我国综合国力的提升，航空运输以其自身优越性越来越受到重视。目前航班运行主要是依赖于航班计划，制定的航班计划是否合理直接关系到航空公司的整个的运行效率和经济效益，因此要改进整个航空运输体系至关重要的是优化航班信息。通过分析影响航空公司的经济效益的影响因素，改善不正常航班调度管理状况。目前，有些航空公司为了追求利益最大化，忽略了航班过站的时间裕度，导致系统的鲁棒性不是很理想。这就使得研究如何对航班计划进行鲁棒性优化很有现实意义。本篇论文主要对航空公司航班计划的合理编排问题，即对航班频率，航班时刻以及机型指派等问题进行了研究与探讨。主要研究内容是通过主成分分析找出航班计划编排中影响经济效益的主要因素，通过对主要因素的提取完善航班计划的制定；改变航班时刻，通过0~1整数规划制定航班计划使利益满足最大化；为了应对不正常航班出现，在考虑航空公司自身拥有的飞机数量和一天合理的停场维修时间的前提下，文中在航班计划中引入鲁棒性。 关键词 :航班计划，航班编排，主成分分析，0~1整数规划，鲁棒性 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. 引言 本文建立航线网络经济性分析、航班频率、时刻的顺序优化和综合优化的数学模型；确定航班网络经济性分析模型，对航线和航线网络的经济性进行探索性分析。难点在于航班计划的“鲁棒性”中如何设定合理的时间裕。 对待处理的数据进行航线收益分析，找出影响收益的主要因素，并根据分析结果提出针对亏损航线的整改措施。假定现有飞行航线和航空公司的营销能力是稳定的（航线、平均客座率、平均折扣率不变）并且假定各航线的航班时刻可以根据需要变动，制定一份航空公司的航班计划，使航空公司的收益最大化。航班计划的“鲁棒性”：重新制定一个带有“鲁棒性”约束的最优航班计划。 本文以2015年湖南省研究生数学建模竞赛为背景 [ 1 ] ，通过主成分分析，0~1整数规划，在航班计划中引入鲁棒性，对航空公司航班计划的合理编排问题，即对航班频率，航班时刻以及机型指派等问题进行了研究与探讨。