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

user: 3. 县区逐日可更新水资源量计算 3.1. 我国县区可更新水资源量数据库 本文介绍县区可更新水资源量评估方法，并计算1961~2009年中国大陆各县区逐日可更新水资源量数据集，随机选择其中一个北方某县区为代表，进行实例分析该方法的优势。 本文采用的shp格式文件(2004版)从DIVA-GIS网站获取，网址如下：http://www.diva-gis.org/gdata。假设第j县区的可更新水资源量用 表示，那么，各县区逐日可更新水资源量计算方法如下： (1) 式中： 表示第j个县第t天的降水量， 是第j个县区总面积， 是这个县区的耕地面积。 和 分别是耕地和非耕地面积对降水的利用率。 和 可能随着不同地区气候条件和下垫面参数而有所差异，这两个系数是集成了裸土蒸发、土壤水动力和产流过程等水文过程可以通过模型模拟而得到的系数 [ 10 ] [ 11 ] 。 通常可将耕地面积 上的有效降水利用系数 取为0.7 [ 12 ] ，非耕地面积 上有效降水利用系数 取0.3 [ 13 ] 。各县区耕地面积和非耕地面积数据可采用EOS数据库统计的数据，也可以根据各级政府发布的统计年鉴进行统计 [ 14 ] 。 县区逐日面雨量的计算，采用从中国气象数据共享网获取的全国730个气象站点的逐日降雨量数据进行空间插值，经分析，部分站点资料序列年份较短或者缺省值较多，研究中选择了其中588个数据系列较长、数据质量较高的站点数据开展研究，少量缺省值用两个最邻近站点值取其平均值代替。 本研究采用泰森多边形方法 [ 15 ] 来进行空间插值计算，泰森多边形法利用离散分布的气象站的降水量来计算行政区面平均降水量。 将所有相邻气象站连成三角形，作这些三角形各边的垂直平分线，于是每个气象站周围的若干垂直平分线便围成一个多边形。然后，以某个县级单位作为目标，假设其总面积为 ，那么这个县级单位可以被划分为多个多边形，每个多边形内包含的唯一一个气象站，这个气象站的降水强度可以用来表示这个多边形区域内的降水强度。最后，将县区边界与多边形边界叠加在一起，计算出某个县区第i个多边形的面积 ，第i个多边形对应的唯一一个气象站的面积权重系数为 ，那么这个县区对应的面雨量 可以用县区内多边形对应的降水强度加权平均求得，公式如下： (2) (3) 泰森多边形的优点是充分考虑了站点分布不均，计算出了每个站点可以代表的面积，计算精度有所提高，因此被广泛应用于面雨量的计算当中。应用中，可以利用ARCGIS Toolbox，依次找到Analysis Tools→Proximity→ Create Thiessen Polygon，双击打开生成泰森多边形对话框。以588个站点为基础，做出泰森多边形。然后再用依次找到Analysis Tools→Overlay→Intersect得出每个县级单位包含的对应的泰森多边形面积所占这个县级单位总面积的比例权重。每个泰森多边形对应一个气象站点，也即对应这个站点的降水数据序列。由此，可以根据公式(3)得到各县区逐日降水量数据库。 3.2. 某县区逐日可更新水资源量的应用及优势分析 以北方某县区为例，分析年、月、日尺度水资源供需平衡分析中，日尺度供需水平衡分析的优势。 图1展示了该县区年、月尺度的供需平衡分析，在年尺度的供需平衡分析仅可以看出，该县区2009年的日平均需水量为70.75万m 3 ，而日平均可更新水资源量为54.55万m 3 ，依靠本区降雨作为内生可更新水资源量可以满足当地水资源需求的76%，剩下24%的用水需求需要依靠河流取水、抽取地下水来满足。然而，并不清楚该县区缺水程度的年内分布，无法精准的对缺水时段进行补给。在月尺度的供需水平衡分析中可以看出，该县区2009年8、11月份日平均可更新水资源量满足当地水资源需求后，仍分别有33.3万m 3 、40万m 3 节余。而其余月份日平均可更新水资源量不能满足当地水资源需求，其中以1、5、7月份缺水情况尤为严重，月尺度供需平衡分析比年尺度可以看出供需水年内分配情况，可以针对可更新水资源量有所节余的月份进行雨水资源化利用或者建立水库蓄水以备缺水月份使用。 图2展示了该县区日尺度的供需平衡分析及缺水累积效应，可以看出由于该县区2009年逐日需水量在40~100万m 3 之间波动，逐日可更新水资源量主要受降雨量影响在0~1300万m 3 之间波动，可以统计出旱季或者连续无雨日持续时间水资源量短缺的累积情况，展示出被月尺度分析均化的4月18日、4月19日和4月23 图1. 年、月尺度水资源供需平衡分析 图2. 日尺度水资源供需平衡分析及缺水累积效应 日，以及9月29日和9月30日逐日可更新水资源量远超当日需水量的节点，更可以反映该县区缺水累积效应的动态发展。
system: [SYSTEM]

user: 4. 结论 以行政单元作为水资源配置和考核的情境下，简单有效的评估行政单元可更新水资源量尤其必要。在水资源管理决策过程中，仅用月平均或者年平均水资源供需平衡分析指标进行水资源指标分配等决策没有考虑到内生可更新水资源量逐日波动变化特征。本文为逐日内生可更新水资源量的评估以及日尺度供需水平平衡提供了新方法、新思路，为区域水资源动态配置提供更为科学的依据。
system: [SYSTEM]

user: 摘 要 为给2016年汛期堰塞体整治、溢洪洞、发电厂房和库区复建交通隧道等项目施工的防洪度汛工作提供决策依据，最大限度减轻施工期间的洪灾风险，根据堰塞湖坝址各设计标准的入库洪水过程、原红石岩电站引水隧洞和溢洪洞联合泄流曲线及堰塞湖库容曲线，分析了当堰塞湖发生200、100、50、30、20、10、5年一遇洪水时，堰塞湖的洪水下泄和库水位变化情况，并选取了2015年汛期发生的三场洪水进行水量平衡分析。分析表明：2015年堰塞湖出、入库水量平衡性较好，实际下泄流量与泄流曲线及天生桥站的水位、流量关系符合性较好，洪水下泄通畅，安全可控。原红石岩电站引水隧洞和右岸溢洪洞(上半洞)能够合并下泄200a一遇以下标准的洪水，坝前水位低于堰塞体泄流槽底部高程。 关键词 :泄洪设施现状，施工期洪水分析，红石岩堰塞湖，牛栏江，鲁甸“8.03”地震
system: [SYSTEM]

user: 1. 概述 2014年8月3日16时30分，云南省昭通市鲁甸县龙头山镇发生6.5级地震，震区鲁甸县火德红镇李家山村和巧家县苞谷垴乡红石岩村交界处，牛栏江干流红石岩段右岸山体大面积塌方，形成最高危险等级(I级)的大型堰塞湖。堰塞体位于原红石岩水电站取水坝下游1.2 km处，顶部左岸高，右岸低，迎水面最低高程点1222.0 m，左岸最高点为1270.0 m，下游最低点高程为1091.7 m，堰塞体总方量约1000万m 3 ，是唐家山堰塞体的4倍。堆积体以白云岩为主，夹带中至薄层状的砂、页岩及灰岩。 在应急和后续处置阶段采取了在堰塞体开挖泄流槽、拆除原红石岩水电站调压井施工支洞堵头检修闸门、利用引水隧洞调压井井筒自由泄流、采取措施防止原电站进水口被漂浮物堵塞及新建应急泄洪洞等五种工程排险措施，至2014年10月4日晚8时堰塞湖内蓄水已全部排空，入流、出流基本平衡，恢复原河道正常过流，红石岩堰塞湖由此转入后期整治阶段。 牛栏江属长江流域金沙江下段右岸一级支流，发源于昆明市嵩明县杨林海，流域呈南北向狭长形，地理坐标系统位于东经102˚53'~104˚05'、北纬25˚02'~27˚24'之间。河流由南向北流经嵩明、沾益、会泽、巧家、鲁甸等县后，在昭阳区麻耗村附近汇入金沙江，河长469 km，流域面积13，672 km 2 。牛栏江位于昭通市境内的流域面积2012 km 2 ，占流域总面积的14.7%。拟建红石岩堰塞湖水库正常蓄水位1200.0 m，相应库容1.41 × 10 8 m 3 ；校核洪水位1208.06 m，总库容1.85 × 10 8 m 3 。堰塞湖上游有牛栏江–滇池补水工程的德泽大(二)型水库、黄梨树电站、小岩头电站、象鼻岭电站(在建中)和入库控制站大沙店水文站；下游有红石岩电站(本次地震受损)、天生桥水文站(“8.03”地震抢险期间新建)、天花板和黄角树电站。堰塞湖下游两岸分布有鲁甸县龙头山等10个乡镇，涉及3万余人和3.6万亩耕地。
system: [SYSTEM]

user: 3. 结语 本文采用堰塞湖水位、导流隧洞下泄流量和天生桥站实测水位对选取的洪水场次进行了分析，通过以上三场入库洪水的变化分析看出：堰塞湖导流隧洞的泄洪比较正常，下游天生桥站实际监测到的水位所查出的流量与下泄流量偏差均在5%以内，与泄流曲线及天生桥站的水位、流量关系符合性较好。还考虑到天生桥站下游约300 m处的沙坝河汇入牛栏江干流时对天生桥水位产生的顶托影响，认为下游天生桥站的出流量是合理的。2015年堰塞湖入库洪水除10月10日的量级大于2年一遇外，其余的均小于2年一遇。8月30日和10月10日的两场洪水均已漫过原红石岩电站取水坝，在堰塞体前形成大量积水，但水位较低，溢洪洞未过流，导流隧洞来水量下泄通畅。尽管2015年堰塞湖的入库流量小于5年一遇，但经过入库洪水预测、现场实际监测和一系列的资料分析等工作，较好的验证了堰塞湖的水量平衡关系，可为来年汛期堰塞湖出、入库洪水的分析工作提供借鉴和参考，也为2016~2017年堰塞体防渗墙、溢洪洞和库区复建交通隧道等项目施工期间的防洪度汛工作提供可靠的技术支撑。 2015年施工单位的作业面位置较高，超过导流隧洞下泄能力的洪水虽然形成新的堰塞湖，但采取了拦污设施对漂浮物进行了打捞，过流通畅，积水在1~2天内即可全部下泄完毕。2016~2017年汛期，堰塞湖整治工程防洪风险的压力仍然很大，溢洪洞洞身已全面转入衬砌作业阶段，堰塞体整治、库区复建交通隧道等项目施工的作业面位置越来越低，施工人员增多，机械设备多，临水低洼处施工组织和防洪避险难度也将增大，堰塞湖施工范围内的防洪形势将更加严峻。建议在加强施工洪水预报工作的同时，充分利用堰塞湖上下游即将建设的水情测报系统，实际监测上游来水量、库水位和天生桥站的水位、流量数据，及时进行堰塞湖入流量和出流量的平衡分析工作。在编制堰塞湖整治工程项目防洪度汛应急预案的基础上，继续在导流隧洞上游采取拦污漂拦截和打捞漂浮物等措施，保障导流隧洞和溢洪洞过流通畅。本方案也有待在来年的防洪度汛工作中作进一步优化和必要的调整，以便更好的将其应用于堰塞湖整治工程防洪度汛工作的实践当中去。
system: [SYSTEM]

user: 摘 要 利用Apriori关联规则挖掘出消费者购房意向中重点考虑的因素，并且通过分析频繁项集，以消费者个人因素为前项，分析不同群体的多元化购房意向，通过了解用户个人特征推测其购房意向，为开发商定位产品提供建议，为销售商精准营销提供方法。 关键词 :消费者，购房意向，关联规则，个人因素，频繁项集
system: [SYSTEM]

user: 1. 研究背景及现状 近年来，开发商数量不断快速增长，同时一些问题也开始显现出来。一方面房价过高让不少国人望“房”兴叹，一方面住宅滞销也屡屡出现。越来越多的业内人士意识到不能盲目兴建楼盘，而是要将开发与消费需求联系起来。消费者需求的研究越来越被重视。研究者们也指出房地产开放商要想使目标顾客的需要得到满足，就必须从顾客的角度出发来思考问题。 向长江等(2002)通过对星汇园潜在购房者的问卷调查及其统计分析，探讨了品牌形象、本体素质、户型设计和地段等因素对潜在购房者感知价值的影响，以及户型设计、室外环境、感知价值等因素对购买意向的影响 [ 1 ] ；吴殿廷等(2007)以北京市和瓦房店市为例，对比了大都市和中小城市居民购房决策影响因素的异同点，指出大都市居民更重视物业和生态环境，而中小城市居民更重视商业环境 [ 2 ] ；于丹等(2007)以大连市为例，利用方差分析和因子分析的方法讨论了影响用户购买意向的因素 [ 3 ] ；蒋绪军(2014)结合客户购买决策模型和商品房客户购房实际情况，构建商品房客户购买决策一般模型 [ 4 ] ；李莹等(2015)以广东省的“85后”消费者作为研究对象，调查他们的购房意愿以及影响因素，在购房决策所考虑的房屋属性方面，“85后”消费者最重视房屋自身的地理位置、价格以及交通这三方面的价值，表现出较强的实用性导向 [ 5 ] 。
system: [SYSTEM]

user: 2. 问题提出 对于消费者购房影响的研究中，多是基于描述给出了单个因素对结果的影响，而并没有分析出影响因素与结果之间的关联关系，也没有刻画出组合因素对消费者购房意向的影响。为更加全面的反映不同消费者群体的购买意向，本文利用Apriori关联规则挖掘的方法，以消费者个人因素为前项，以消费者购房需求为后项。通过住房自身特征、住房所在小区特征和小区周边环境特征三方面表征消费者的购房需求。进而挖掘出具有不同个人特征的消费者群体与不同购房需求之间的关系。为开发商的产品定位和消费者需求探索提供合理可靠的建议。
system: [SYSTEM]

user: 3 理论基础 1) Apriori算法基本思想 首先找出所有的频繁项集，这些项集出现的频繁性至少和预定义的最小支持度一样，然后由频集产生强关联规则，这些规则必须满足最小支持度和最小可信度，再使用第1步找到的频集产生期望的规则，产生只包含集合项的所有规则，一旦这些规则被生成，那么只有那些大于用户给定的最小可信度的规则才被留下来。 2) 算法分析步骤 第一步：通过迭代的方法，找出事物数据库中的全部频繁项集，即支持度大于或等于用户设定的阈值的项集； 第二步：利用频繁项集构造出满足最小置信度的规则 [ 6 ] 。
system: [SYSTEM]

user: 5. 结论和建议 1) 首先，从整体看来，从住房自身特征来看，现房、房价小于等于10，000元/m 2 、房屋面积在51~120 m 2 的房屋更受消费者欢迎；从住房所在小区特征来看，小区配套公众设施有游泳池、健身设施和商场超 图7. 小区周边环境频繁项 Table 8 The association rules table of the personal factors and the characteristics of the surrounding environment of communit 前项 方向 后项 支持度 置信度 提升度 {B1 = 0， B3 = 1， B4 = 0， B6 = 1} ≥ {A12 = 3} 0.010 1 18.18 表8. 个人因素与小区周边环境特征的关联规则表 注：对于置信度为1的情况是由于样本量较小。 市、小区物业管理包括家政服务、水电通讯服务和保安服务的小区更被关注，这些体现出现代消费者对健康和便捷的需求，开发商可以依此为卖点来吸引消费者；从小区周边环境特征来看，除传统的医院和银行外，消费者对学校的需求也十分高，那么学区房的市场可以保证。 2) 其次，从开发商的角度来说，顾客群的划分也十分重要，从分析可知中老年人群相较于青年选择大面积、高房价的房子的趋势更明显。年龄、收入和家庭结构会在一定程度上影响消费者选择，且购房者均多考虑未来生活安排，所以需要挖掘其可能的潜在需求。若开发商计划针对某一群体进行楼盘售卖，那么便可以推出各方面条件适合该群体的楼盘；若开发商计划开发多元化的楼盘，那么便需要各种类型取交集或互补开发以吸引客户。如针对青年、低收入、家庭成员小于等于两人(即单身、同居或结婚暂时无子女)的群体可以开发小户型、安保设施齐全、有家政服务并且距离学校较近的楼盘。从销售商的角度来说，通过了解有购房需求的消费者的个人和家庭特征，通过关联规则挖掘的结果将其购房意向具体化，有针对性的进行推荐，有助于提高房屋售卖率。如消费者特征为中老年、家庭成员小于等于2人，那么便可以为其推荐91~120 m 2 、三室两厅两卫的房屋。当然，这并不是绝对的规律，还要根据购房者具体意愿进行分析。 3) 最后，本研究发现，仅通过某个单独的个人特征来对消费者购房意向进行定位和推荐并不全面，多种因素组合才能还原出消费者的真实面貌，进而为开发商提供决策建议。通过关联规则挖掘可以通过划分出特征的细微不同，进而实现较为购房意向精准的预测和推荐，从而提高房屋出售率。需要说明的是，本次调查问卷样本量较小，可能出现个别指标与实际有一定差距，但并不影响总体特征的反映。
system: [SYSTEM]

user: 摘 要 本文考察一类含参变量的函数，揭示函数在参数不同取值下的性态，并利用matlab进行了仿真，结果验证结论的正确性。 关键词 :稳定点，零点，最终单调递增(减)
system: [SYSTEM]

user: 1. 引言 从现有的文献来看，对含参变量函数性态的研究甚少，而多数都是讨论含参变量函数的积分、Riemann边值问题或Hilbert边值问题 [ 1 ] - [ 3 ] ，梅宏则在文献 [ 4 ] 中给出了一类含参变量积分的常差分方程计算方法。 本文对一类特殊的带参变量函数的进行分析，揭示该类函数在参数不同取值下的性态。 本文考察如下函数： ， 其中参数 。
system: [SYSTEM]

user: 2. 主要结果 首先，函数 对 求导数 ， 令 ，则有 (1) 同样，函数 对 求导数 ， 再令 ，则有 (2) 当 时，记 ，显然 ，因而 是函数 的稳定点。 引理1 对于函数 和 ，下述成立： (i) 当 时，对于任意 ，有 ，从而 关于 在 上严格单调递增。 (ii) 当 时，对于任意 ，有 ，从而 关于 在 上严格单调递增。 (iii) 当 时，有表1成立。 Table 1 The monotonicity of the function at the cas 表1. 函数 在 下的单调性 且 。 证明 (i) 当 时，对于任意 ，有 ，由(2)知， 关于 在 上严格单调递增。 (ii) 当 时，对于任意 ，有 ，由(2)知， 关于 在 上严格单调递增。 (iii) 当 时，注意到对任意 ，有 ，所以 在 上关于 严格单调递减。而 ，所以存在唯一的 使得 ，再结合(2)知表1成立。而 显然成立。注意到 ，故 。因而 。 定理1 对于函数 和 ，下述成立： (i) 当 时， 存在零点 ，有表2成立。 且 。 (ii) 当 时，有表3成立。 (iii) 存在 满足 ，且 (a) 当 时，存在 及 使得表4成立。 (b) 当 时，有表5成立。 (c) 当 时，存在 及 使得表6成立。 证明 (i) 当 时，注意到 ， ，再结合(1)、(2)及引理1 (i)可知，存在 使得 ，从而表2成立。显然 ， 。 (ii) 当 时，注意到 ， 及 ，再结合(1)、(2)及引理1 (ii)可知表3成立。 (iii) 当 时，记 ，由引理1 (iii)知 。由 ，解得 ，而 Table 2 The monotonicity of the function at the cas 表2. 函数 在 下的单调性 Table 3 The monotonicity of the function at the cas 表3. 函数 在 下的单调性 Table 4 The monotonicity of the function at the cas 表4. 函数 在 下的单调性 Table 5 The monotonicity of the function at the cas 表5. 函数 在 下的单调性 Table 6 The monotonicity of the function at the cas 表6. 函数 在 下的单调性 ，从而有表7成立。 令 ，可得 满足 ，且有表8成立。 (a) 当 时，由表8可知 ，再结合引理1 (iii)及(1)知，存在 及 使得表4成立。 (b) 当 时，由表8可知 ，再结合引理1 (iii)及(1)知表5成立。 Table 7 The monotonicity of the function I(α 表7. 函数I(α)的单调性 Table 8 The value of the functio 表8. 函数 的取值 (c) 当 时，由表8可知 ，再结合引理1(iii)及(1)知，存在 及 使得表6成立。 定理2 对于函数 ，下述成立： (i) 当 时， 。 (ii) 当 时， 。 证明 (i) 当 时，由定理1知， 是关于变量 最终单调递增的，所以不妨设 。显然 ，下证 。反证，假设 ，则 ， 与假设矛盾，因而 。 (ii) 当 时，由定理1知， 是关于变量 最终单调递减的，注意到对任意 有 ，所以不妨设 。显然 ，下证 。 ， 即 。 定理3 对于 有 且 ， 。 证明 当 时，方便起见，记 ， 。则 。 将 对 分别求一阶导数和二阶导数 ， ， 令 ，得 。注意到 ，所以有表9成立。 即 在 上严格单调递减。而 ， ，所以存在唯一一点 ，使得 ，亦即 。而 ， ， ，从而有表10成立。 因此存在唯一一点 使得 ，且当 时， ；当 时， 。由 ， ，知 ， 。 Table 9 The monotonicity of the functio 表9. 函数 的单调性 Table 1 0 The monotonicity of the functio 表10. 函数 的单调性
system: [SYSTEM]

user: 3. 数值仿真 这一部分我们用matlab描绘在参数 不同取值下函数 的图像。 当 时，函数 在参数 不同取值下的图像如图1所示。 当 时，函数 在参数 不同取值下的图像如图2所示。 图1. 在参数 下的图像 图2. 在参数 下的图像
system: [SYSTEM]

user: 4. 小结 本文引入函数 和 ，通过考察其性质(引理1)进而考察函数 的性质(定理1、定理2)。接着考察了函数 的稳定点 (定理3)，为分类参数 的区间进而准确描绘函数 的图像做了准备。
system: [SYSTEM]

user: 摘 要 本文研究了微分方程对称方法在非线性偏微分方程边值问题中的应用。首先，基于微分特征列集算法确定了给定非线性偏微分方程边值问题的多参数对称；其次，利用对称将非线性偏微分方程边值问题化为常微分方程初值问题；最后，利用龙格-库塔法求解了常微分方程初值问题的数值解。 关键词 :非线性偏微分方程边值问题，微分特征列集算法，对称方法，龙格–库塔法
system: [SYSTEM]

user: 1. 引言 Lie对称是公认的普适性方法之一，它以诸多传统方法为其特例，如：分离变量法、行波变换、相似变换等 [ 1 ] 。对称理论在现代数学、物理和力学等科学中有重要的理论和实际意义，并且已有了广泛的应用 [ 2 ] [ 3 ] 。但是除了在文 [ 1 ] [ 4 ] [ 5 ] 中研究者做了一些对称群在边值问题的应用外，这方面的研究还很少，所以利用对称群研究偏微分方程(简记为PDEs)边值问题是对称理论应用新的研究领域。在诸多求解非线性的问题中，常常利用相似变换对PDEs进行约化或降阶。我们知道PDEs的Lie变换群可以产生更一般形式的相似变换，并且这些变换有更多的数学和物理意义，所以利用对称方法研究非线性PDEs边值问题比直接用相似变换更有其优越性。 最近我们基于微分特征列集算法 [ 6 ] [ 7 ] 研究了对称方法在非线性PDEs边值问题中的应用，并且将对称方法、对称分类与数值方法、同伦分析方法进行结合解决了非线性PDEs边值问题 [ 8 ] - [ 11 ] 。也有研究者基于微分特征列集算法，将对称方法与变分迭代法、同伦摄动法结合，求解了PDEs边值问题 [ 12 ] [ 13 ] 。本文中，我们将基于微分特征列集算法，有效结合对称方法和数值方法求解一个非线性PDEs边值问题。
system: [SYSTEM]

user: 2. 一个非线性PDEs边值问题的对称约化及其数值解 考虑一个如下非线性PDEs边值问题 (1) 其中 为常数。相应的边界条件为 ， ，若 (2) ，若 (3) 下面引入流函数 (4) 将上面的关系代入(1)~(3)，得到 (5) 满足边界条件 ， ，若 (6) ，若 (7) 其中函数 ， 在后面根据边界条件在对称作用下的不变性来确定。 2.1. 边值问题的对称约化 下面利用对称方法约化边值问题(5)~(7)。 第一步：确定PDE边值问题(5)~(7)的多参数对称。 假设非线性PDE(5)的对称的无穷小向量为 其中 ， ， 为该对称的无穷小生成函数。 1) 产生关于无穷小生成函数的确定方程组。 根据Lie算法，利用产生确定方程组的算法 [ 7 ] 可得到方程(5)的对称的确定方程组，但是难于人工求解。我们可以利用微分特征列集程序包进行计算得到与确定方程组等价的特征列集对应的方程组，即 (8) 2) 确定多参数对称。 通过求解方程组(8)，得到无穷小生成函数 即可得到对称的无穷小向量为 (9) 其中 为任意常数， 为任意函数。 第二步：约化边值问题。利用对称(9)将边值问题(5)~(7)约化为常微分方程初值问题。 对称(9)的特征方程为 由 ，得到不变量 同理可从特征方程得到 (10) 将(10)代入方程(5)，可得到常微分方程 (11) 根据PDEs边值问题的边界条件在对称群作用下的不变性 [ 1 ] 知，边界条件(6)，(7)在对称(9)的延拓作用下不变，故有 ，当 (12) ，当 (13) ，当 (14) 其中 为对称 的一阶延拓的无穷小向量， 为对称 的二阶延拓的无穷小向量，即 由关系式(12)，(13)，(14)可以确定函数 ，即 ， 其中 为任意常数。 为了对应边界条件(6)和(7)，取 ，故有 当 时， 。 根据边界条件(6)，(7)和关系(10)，得到初值条件 (15) 2.2. 边值问题的数值解 下面我们利用四阶龙格–库塔法求解常微分方程初值问题(11)，(15)的数值解。 为了使用龙格–库塔法，首先将(11)，(15)化为一阶常微分方程初值问题。令 则(11)，(15)可以化为 (16) 相应的有初值条件 (17) 令 其中 都是 的函数，那么可以建立龙格–库塔公式 其中 ， 为步长。 常微分方程(16)中的常数 和初始条件(17)中的 取其不同值时，我们也可以得到相应的数值解。下面我们取 ，并且 时，借助Mathematica符号系统可以得到函数 在区间 上的数值解，如图1所示。 图1. 在 上的数值解
system: [SYSTEM]

user: 3. 结论 本文研究了微分方程对称方法在非线性PDE边值问题中的应用，也探索了利用对称方法与数值方法的有效结合求解PDE边值问题的新途径。我们首先基于微分特征列集算法确定了给定非线性PDE边值问题的多参数对称，并利用对称将非线性PDEs边值问题约化为常微分方程初值问题。其次我们充分考虑到对称方法与数值方法的彼此互补性，借助于Mathematica符号系统利用四阶龙格－库塔法求解了上一步得到的常微分方程初值问题的数值解。得到的结果充分体现了对称方法在偏微分方程应用中的优越性，并且我们在本文中有效结合对称方法和数值方法求解了边值问题，这项研究也推广了对称方法的应用范围。
system: [SYSTEM]

user: 摘 要 Catalan数是指通项公式为 的序列中的 这些数，其最早是由我国清代数学家明安图开始研究的。本文运用Catalan数与生成函数法来解决一类特殊数学结构的计数问题，构造出该数学结构解个数的显性表达式。最后还给出了该数学结构计数问题的另一种解决方案，从另一个角度也利用了广义的Catalan数来解决问题。 关键词 :Catalan数，组合计数，生成函数法
system: [SYSTEM]

user: 1. 引言 Catalan数有着非常丰富的历史背景，是组合数学中一类使用广泛，拥有众多组合意义的组合数，常用来解决在各类组合计数问题。一般记为 。其名字来源于18世纪比利时的数学家Eugene Charlie Catalan。其在1838年研究Catalan问题：不满足乘法结合律的n个因子的乘积，在保持因子顺序的结合数？最终获得了Catalan数的显性表达式： 。 C n 是该数列的第n + 1项。这是贾宪三角中垂线上的数依次除以其相对应的自然数之后所得一个数列。其与Fibonacci sequence，Stirling数一起作为组合数学之中的典型例子。但实际上卡Catalan并不是第一个研究Catalan数的数学家。在十七世纪三十年代时，我国清代数学家明安图在研究无穷级数时得到了Catalan数的三种不同算法 [ 1 ] ，并作为计算特定级数的方法在之后的研究中使用，建立了一个暗含了Catalan数的几何模型 [ 2 ] 。其中有两个算法公式在当代的组合数学书籍以及论文之中都没有涉及到。有关于明安图计算无穷级数时所得的一系列结果都详细记载在《割圆密率捷法》中，在这一本书中给出了Catalan数的卷积型递推公式为 。 其中另一个重要的公式是获得Catalan数的组合型递推公式 。 在这之后不久的1758年，大数学家欧拉与赛格纳讨论这么一个问题：用一个凸多边形彼此之间不相交的对角线把凸多边形划分成若干个三角形区域的方法数。针对于此问题，划分n + 1条边的多边形的方法数就是Catalan数。而之后的董祐诚等清代数学家受明安图的影响下，对Catalan数都做出了一定的研究成果，其研究时间都早于Catalan研究Catalan数。之后比较值得注意的是1839年Bi Noto，其通过一系列证明，最终推出了Catalan数的生成函数的表达式： 。 从上文可以知道，Catalan数最早并不是Catalan提出的，但由于一种习惯或约定，现在已经变成数学名词。 在19世纪，数学家们得到两个Catalan数的递推公式： 随着越来越多的数学家研究Catalan数，在越来越多的领域里发展和拓广Catalan数的应用和新的组合解释，如：非结合代数问题、随机游动问题、图形剖分问题、特定非负整数解问题、出栈问题、几何、代数问题，网格路径等组合解释。初步估计下来应该在200种以上，他们散落在各种文献之中。初文昌在《关于Catalan序列》一文中给出了25种不同的组合解释 [ 3 ] 。 近年来，对于Catalan数的研究主要是得到了Catalan数各种拓广形式和其相对应的枚举结果。主要是以下几个方面的：高阶Catalan数、广义Catalan数 [ 4 ] 、连续q-模拟、Catalan数的同余(或整除)性质，二项式系数问题 [ 5 ] ，Motzki数 [ 6 ] ，Riordan群 [ 7 ] 等方面等。本文的主要工作就是利用Catalan数与生成函数法来解决一类特殊数学结构的计数问题，构造出该数学结构解的显性表达式。最后还给出了该数学结构计数问题的另一种解决方案，从另一个角度也利用了广义的Catalan数来解决问题。
system: [SYSTEM]

user: 2. 一个广义Catalan数问题 问题1： 关于该数学结构计数问题，是一类广义Catalan数问题，涉及到线性不等 式组的整数解以及整数规划问题。是一类特定非负整数解问题。现在能在网上找到的解答都是在已知结果的情况下，运用数学归纳法来解决的。其前提都是在已知 解结构的情况下来完成的。但并不知道这个解是怎么来的。接下来将用两种不同方法来构造出该数学结构具体解的形式。它们分别涉及网格路径法、生成函数法两种不同的数学方法以及广义Catalan数。在求解该数学机构计数问题过程之中，可以领悟到Catalan数运用的广泛性及其深刻性。
system: [SYSTEM]

user: 3. 预备知识 先给出在之后的证明过程中需要的一般概念，定义和相关结论。 生成函数法的思想其实就是让幂级数与离散的数列一一对应，把数列中的第n项与幂级数中指数为n的项相对应起来，再利用幂级数的运算性质去确定离散数列的结构。 定义1. 对于一个有限或无限的序列 [ 8 ] 称 为 的生成函数，称 为 的指数型生成函数。 定义2.不满足乘法结合律的n个因子的乘积，在保持因子顺序的情况下，有 种结合方法。这里的 就等于Catalan数序列中的第n项。 以Catalan数为系数的形式幂级数的和为 。 引理1. 网格图中从点 到点 ( )的最短路径数为 。 证：从网格的任一结点处向上走一格记为y；向右走一格记为x；那么，从点O走到点M要向上走b格，向右走a格。从点O到点M的每一条最短路径都对应一个由a个x和b个y组成的长为 的字，且两者是一一对应的。因此，只需要在 个字母中指出哪a个是x (或哪b个是y)就对应一条最短路径。所以， 引理2. 网格图中从点 到点 ( )的最短路径满足 的路径数为 。 证： 等于从点 到点 的最短路径减去从点A到点M穿过直线 的最短路径数。后者又与从点 到点M的最短路径数相等。即 。 有关的组合恒等式： 。 。
system: [SYSTEM]

user: 4. 生成函数法解决 这里的话，要先试着去找出 它的边界条件和递推关系式来。 先就 的取值进行一个讨论 1) 若 ，则问题转化为求 ， 2) 若 ，则原问题转化为 。 从而问题可转化为 。 综合1)，2)可以得到这样一个关于 的递推关系式。 。 但这时边界条件还没有得到。接下来，来找到 的边界条件，就 ， 分开讨论，得出其所对应的边界条件。 3) 当 时，原问题可转化为 ；此时解的个数相当于从 这 个数中选出一个数的种数 。故 。 4) 当 时，原问题可转化为 又等价于 来思考其几何意义，其表示从网格图上点 出发达到直线 上的各点最短路径数(满足 )的和。又根据引理2：网格图中从点 到点 ( )的最短路径满足 的路径为 。 故根据定义我们可以得到 的计算公式 。 5) 当m = 1时，原问题可转化为 这与4相似，类比4得到： 。 6) 现知 的递推关系式 边界条件 ； ； 。 下面利用生成函数法来解出该递推关系式。 构造一个二元幂级数 (1) 根据边界条件， 时的 的表达式是已知的，故把其分成 和 两种情况来分开计算 (2) 同样是根据边界条件， 时的 的表达式是已知的，故把 与 两种情况也分开计算 (3) 其中为了方便计算把(3)式中的第三项单独拿出来，进行变形化简之后再放回(3)式中进行计算 (4) 将 的递推关系式 代入(4)中可得 (5) 先将(5)中第一项的式子里提出一个 ，第二项里提出一个 ，使得 中n，m可以与x，y的指数相互对应起来得到 (6) 通过改变指标的取值范围，可以将 与x，y的指数相互对应起来得到了 (7) 由(1)式的关系，把(1)进行适当的拆分，移项后，可以得到： 将这两个式子带入到(7)中代替第一项和第二项化简后可以 (8) 适当改变n，m的取值范围，变形(8)式，可以得到下面这样一个式子 (9) 将(9)式带入到(3)式之中去，已知 又 ； 。 我们将其代入后计算可以得到： 化简后的(3)式为 (10) 将 代入(10)式 (11) 又(1)式等于(11)式，我们将其相等后化简可得下式： (12) 在这里 为Catalan数，根据已知的结论， 。 将其带入到(12)式之中，可得到 (13) 运用泰勒展开公式 ； 则(13)式展开后变形化简之后，可得到下式 (14) 根据生成函数的理论可知 至此该数学结构解的个数就被我们所构造出来了。在这之中运用了广义Catalan数来帮助进行计算，体现了Catalan数应用的广泛性，在这里使用生成函数法，给出了 的一个显式的表达式。在这里的关键是能够把 这一序列与生成函数 建立一个一一对应，然后再利用生成函数的性质和运算，得到另一个序列，从而两个序列相比较得到了 的表达式。现在已经用生成函数法得到了该数学结构的计数的显性表达式，但这个过程中也发现了，计算比较的繁琐，接下来将通过解析该数学结构的几何意义后，通过网格路径法来给出一种完全不同的解决方案。
system: [SYSTEM]

user: 5. 网格路径法解决问题 我们现在来思考一下问题1的几何意义。先来再看一次这个不等式组 。 从上到下来开，从原点出发，左式每加上一个x i ，就向右移动一个单位，向上移动 个单位。那么最终来看就是向右移动了 个单位，向上移动了 个单位。在做了这样的转换之后，得到了下面这样一个等价问题： 网格图中从点 到 上各点 ( )的最短路径经过任一节点 时，满足 (即路径与直线 )的路径数的和为 。 定理1. 记网格图中从 到点 ( )的最短路径经过任一节点 时，满足 (即路径不穿过直线 )的路径数 [ 9 ] 证明：从点 到点 的不穿过直线 的最短路径数等于从点 到点 的最短路径数减去从点O到点M穿过直线 的最短路径数。后者又等于从点 到点 的最短路径数减去从点O到点M碰到或穿过直线 的最短路径数。 下面建立从点 到点 的最短路径数减去从点O到点M碰到或穿过直线 的最短路径(记做 )与从点 到点M的最短路径(记做 )之间的一一对应关系。 设 与 的第一个交点为P，从点O到点P的这段路径OAP关于直线 的对称路径记为 ，以 代替OAP得到从点 到点M的一条路径 。反之，从点 到点M的任一条路径与 必有交点。记第一个交点为P，取 关于直线 的对称路径 ，以 代替 ，得到从点O到点M的一条路径 ，即 与 是一一对应的。于是 从 的几何意义我们可以得到 将 ， 代入可以得到 。 由此得到了满足该数学结构的非负整数序列 的个数 。
system: [SYSTEM]

user: 摘 要 柑橘类果汁在加工过程，在风味上受柠檬苦素影响甚巨，本研究针对台湾柚类柠檬苦素葡萄糖转移酶基因(limonin glucosyltransferase gene， lgt)，选殖出一1600 bp片段，并对其作用活性位置预测出一段转录的44胺基酸序列，得到一个完整的开放框架(open reading frame)，此段功能性基因有助于未来进行构筑与选殖，应用于柑橘类果汁在风味改进的研究。本研究亦探讨台湾的文旦柚和蜜柚lgt在开花至成熟时期的表现情况。结果得知，蜜柚中仅有第二型存在，且果实中含有大量lgt表现，而文旦同时具有两型lgt存在，但表现量相较于蜜柚低，此结果可应证前人之结论，即文旦柚于果实成长后期采收，因 lgt表现量不足，而导致其果实中仍含有大量苦味化合物仍有存在，导至于加工过程中苦味化合物的生成。本研究利用內果皮以快速蛋白纯化系统(FPLC) HiTrap Q管柱进行纯化，得分子量大约40 kDa之LGTase，并以2，2'-二辛可宁酸(bicinchoninic acid)定量蛋白质为0.6，1 mg/ml，此方法有助提高柚类加工过程中果皮废弃物的利用性。 关键词 :柚类，文旦，蜜柚，柠檬苦素葡萄醣转移酶，加工废弃物
system: [SYSTEM]

user: 1. 引言 柑橘(Citrus spp.)属芸香科常绿果树，为全球产量和销售最大宗之柑橘果汁，因颜色鲜艳和香气浓厚之特性而广受喜爱，然而柑橘果汁中常因具苦味而造成其质量之下降，而在加工过程中加以去除以提高果汁之风味；近年来研究报告指出，苦味化合物包括柠檬苦素类(limonoid)与黄烷酮配糖体(flavonoid glycoside)两大类 [ 1 ] [ 2 ] 。到目前为止，已发现有36种柠檬苦素化合物 [ 3 ] 。柠檬苦素类三帖类化合物，主要存在于芸香科和楝科的植物中，其中柠檬苦素(limonin)为果汁苦味主要来源之一，也是造成果汁制造质量降低的重要原因。90年代，大部份的学者针对如何将引起苦味之化合物以分离的方法去除，或者使用微生物将此些化合物进行代谢分解 [ 4 ] [ 5 ] ；晚近，由于相继找到此类化合物在抑制癌症的效果 [ 6 ] - [ 9 ] ，因而研究方向转而希望利用酵素对其进行作用，以达到降低苦味及保留此类化合物之目的，然而纯化酵素之成本高，若将其用于工业生产上可能不符合经济效益。故2000年左右，学者开始将其相关的酵素基因进行转殖，期望可生产出成本低之酵素，以供工业上使用 [ 10 ] 。过去学者研究之柑桔果汁去苦味方法主要有几类，如以微生物将苦味物质代谢与分解、以各式管柱对其进行吸附，或以酵素将其作用为不具苦味之化合物 [ 11 ] 。柠檬苦素之含量若高于6 ppm则会造成些微苦味的感觉，若是高于18 ppm则苦味会较重 [ 12 ] 。此外，如脐橙果实在鲜食时并不具苦味，但若进行果汁加工之后会产生苦味物质，主要是鲜食时多数的柠檬苦素是以不具苦味的前驱物形态存在，如limonoic acid A-ring lactone，但在加工过程中其在酸性环境下被limonin D-ring lactone hydrolase作用，则会形成具苦味之化合物 [ 13 ] ，由于柠檬苦素多存在于内果皮，因此在榨汁时亦会进入果汁中造成苦味，在酸性环境下经过酵素作用，则使柑桔果汁呈现具苦味之情况，故又被称为后生苦味。本年度研究重点于建立柠檬苦素葡萄糖转移酶蛋白质纯化及基因选殖系统，并将纯化后之柠檬苦素葡萄糖转移酶之基础物化性狀研究。研究中除寻求最大酵素活性來源外，也将藉由mRNA差異表现探讨其可能之表现情形。
system: [SYSTEM]

user: 2. 材料与方法 2.1. 材料 样品取得于新竹县柑橘果园，于3月开花期至6月成熟果每隔30天取样一次。取得样品将果办膜分离储藏于−80℃，抽取DNA、RNA与纯化LGTase。 2.2. 方法 2.2.1. DNA萃取 新鲜叶子于液态氮下磨碎，取0.1 g粉末分别装于微量离心管中，以Wizard ® Genomic DNA Purification Kit (Promega， USA)进行DNA萃取。 2.2.2. 限制酶之DNA图谱分析 将DNA溶液(0.1~1.0 μg)，加入反应缓冲溶液(10×) 1 μL，加入限制酶1 μl，及去离子水依序混合，并调整反应之总体积为10 μl，将混合后之溶液置于37℃下作用两小时，反应完成后以65℃水浴10分钟来终止反应，即可进行电泳分析。 2.2.3. RNA萃取 将样品于液态氮下磨成粉末，取0.1 g装于1.5 ml离心管中，以RNeasy Plant Mini Kit (QIAGEN， Germany)进行RNA萃取。 2.2.4. cDNA合成 将萃取出RNA以Reverse Transcriptases and RT-PCR Kits (Promega， USA)进行cDNA合成。 2.2.5. 实时定量PCR 设计LGTase基因专一性引子(表1)片段大约为208 bp。于罗氏 ® Light-Cycler系统上机操作，PCR反应条件如下95℃ 3 min；95℃ 3 min，50℃ 30 sec，72℃ 1.5 min，35 cycles；72℃ 7 min。 2.2.6. LGTase纯化 取果实内果皮部位与溶剂(pre-cooled 0.1M saline，含0.5% PVP，1 mM PMSF在0.1M stock solution DMSO) 1:4，于约600 rpm转速均质3 min，将均质后的泥浆物再搅拌2 h，以滤网过滤。于4℃下离心10，000 rpm，30 min，75%(NH 4 ) 2 SO 4 使之饱和状态，将蛋白质沉淀物进行离心收集，沉淀物进行溶解使用pre-cooled 10 mM Tris-HC1 buffer，pH 7.8透析12 h，同时置换四次。透析完经过离心去除不溶性物质，上清液为粗酵素萃取物粗萃液进行硫胺划分35%~75%，将硫酸胺饱和蛋白质沉淀物进行透析10 mM Tris–HC1 buffer，pH 7.8，以FPLC HiTrap Q管柱进行纯化。 2.2.7. LGTase活性测定 反应溶液1 mM limonin (withan open D-ring)，2 mM UDPG-2Na，5 mM MnCl 2 ，20mM Tris-HC1 buffer混合后于37℃下反应1 hr。反应终止：沸腾加热100℃ 5 min，取上清液并离心10，000 rpm，15 min。HPLC Nucleosil 7 C 18 reversephase column (4.6 × 250 mm)进行分析，洗提15%~30% methanol线性梯度，35 min，流速0.4 ml/min。比对limonin glucoside标准品，计算单位活性pkat。 2.2.8. 蛋白质定量 以BCA (bicinchoninic acid)方法进行测定，当BCA以水溶性钠盐形式存在时，对cuprous ion (Cu 2+ )是1种靈敏、稳定性高且具高特異性之试剂。在碱性环境下，蛋白质可还原Cu 2+ 成为Cu + ，接着Cu + 再和BCA试剂反应形成紫色的生成物，于波长562 nm下测定其吸光値，以已知浓度BSA (bovine serum albumin) Table 1 Real-time RT-PCR amplification temperatures on type I and II LGTas Primer Sequence(5’-3’) Amplicon length (bp) Optimum Tm (˚C) LGTase1 5’ CGAGAATAATGGGAACTG 3’ 1600 bp 46 LGTase2 5’TGAAAGGGAATAGATAAGTAC 3’ LGTF2 5’ACTCGTAGCATCGGGGGTGCCGGT 3’ 208 bp 52 LGTA1 5’ TCAGCGCGTTCTCCTCCAGCGCCA 3’ LGTF2 5’ACTCGTAGCATCGGGGGTGCCGGT 3’ 208 bp 52 LGTA2 5’ TCAGCGCGTCCTTCA GCTCCG 3’ 表1. 实时定量RT-PCR扩增一型与二型产物引子及其最适链合温度 为蛋白质标准品，绘制标准曲线，即可定量蛋白质。 2.2.9. SDS PAGE测定 SDS-PAGE胶体的制备參考Laemmli (1970)之方法，以10% polyacrylamide gel及Mini-Protein II迷你电泳槽进行电泳，以2比1的比例混合样本和处理液(2.5 ml 1.5 M Tris-HCl， pH 8.8， 4 ml 10% SDS， 2 ml glycerol， 1 ml β-mercaptoethanol， 292 ml 1% bromophenol blue， 2 mL 1 M DTT)，进行电泳条件为100 V 2小时。 2.2.10. LGTase生化特性测定 1) 最适pH、温度 将酵素分别加入pH 3.5~10.5反应液(Tris-HCl)作用24小时，利用10 mM Tris-HCl buffer，pH 7.0透析12小时去除盐类，离心后取上清液进行蛋白质活性测定，最适温度测定以10℃~80℃做测定。 2) 最稳定pH、温度 将酵素分别在pH 3.5~pH 10.5反应液中于4℃作用12 h，去除盐类，离心取上清液进行蛋白质活性测定。
system: [SYSTEM]

user: 4. 结论 本研究选殖出台湾文旦柚和蜜柚之柠檬苦素葡萄糖转移酶基因，此1600 bp片段可转录出一含44胺 图5. 以HiTrap Q管柱纯化文旦LGTase FPLC系统纯画图SDS-PAGE Table 4 Purification of wentan PUMMELO lgt enzym Purification steps Total protein (mg) Total activity (unit) Specific activity (unit/mg) Yield% Purification fold Crude enzyme 261.76 259 0.99 100 1.0 35%-75%(NH4) 2 SO 4 126 250 1.98 96 2 FPLC 0.6 107 178 41 180 表4. 文旦柚lgt酵素纯化 Table 5 Purification of miyu pummelo lgt enzym Purification steps Total protein (mg) Total activity (unit) Specific activity (unit/mg) Yield% Purification fold Crude enzyme 264 367 1.4 100 1.0 35%-75%(NH4) 2 SO 4 135 281 2.08 76 1.5 FPLC 1 188 188 51 134 表5. 蜜柚lgt酵素纯化 图6. HPLC测定酵素活性图谱 图7. 文旦柚LGTase最适pH 基酸序列，此段功能性基因有助于未来进行构筑与选殖之用。探讨台湾的文旦柚和蜜柚lgt在开花至成熟时期的表现情况时，蜜柚中仅发现有第二型存在，而文旦柚虽同时具有两型lgt存在，但表现量相较于蜜 图8. 文旦柚LGTase最适温度 柚低，导致其果实中仍含有大量苦味化合物仍有存在，导至于加工过程中苦味化合物的生成。本研究亦以管柱进行进一步纯化，得分子量大约40 kDa之LGTase，经酵素活性换算得文旦柚LGTase纯化产率达41%、51%，电泳分析得分子量大约40 kDa。文旦内果皮纯化出LGTase其最适与最稳定pH为7.8与5.5~8.5，而蜜柚内果皮纯化出LGTase其最适与最稳定pH为7.8与4.5~8.5。
system: [SYSTEM]

user: 摘 要 研究了黄绿绿僵菌对桃蚜GSTs和AChE 2种酶系的影响，结果表明：24 h以后，桃蚜体内AChE的活力急骤下降，直到40 h时达到最低，抑制率为52.81%，其后一直维持在该水平；桃蚜体内GSTs的活力在染菌40 h后，其活力快速下降，直到56 h降到最低，抑制率为38.91%，并保持在该水平。 关键词 :黄绿绿僵菌，桃蚜，酶活性
system: [SYSTEM]

user: 1. 引言 昆虫病原真菌穿透寄主体壁进入血腔中，必然会遇到昆虫防御系统的抵抗，即昆虫受到刺激后会产生应激反应，引起的一系列代谢变化 [ 1 ] 。外源物进入昆虫体内后会与多种代谢酶系发生作用，其中GSTs等解毒酶是在昆虫体内代谢过程中起着重要作用的酶类，AChE是有机磷杀虫剂和氨基甲酸酯类杀虫剂的作用靶标，这类酶在对外源化合物的解毒代谢方面起着重要作用，药剂对上述酶系的影响均会造成昆虫代谢紊乱而影响正常的生理生化过程。 本研究以桃蚜为试虫，研究了黄绿绿僵菌对GSTs和AChE 2种酶系的影响，以期从生理生化角度探索提取物对试虫作用的靶标，及昆虫的解毒机制，阐述黄绿绿僵菌的杀虫作用机理，为其开发应用提供理论依据。
system: [SYSTEM]

user: 2. 材料与方法 2.1. 供试材料 黄绿绿僵菌Mf82菌株是从褐飞虱罹病虫体分离纯化后保存于昆虫病理实验室的菌种。在SDAY (SDAY配方含4%葡萄糖、1%蛋白胨、1%酵母粉及2%琼脂)斜面上接种培养10 d左右，培养条件：置于4℃冰箱中保存备用。用于饲养桃蚜的甘蓝叶片，采自烟草所温室大棚内播种的甘蓝。 2.2. 供试昆虫 桃蚜采自即墨试验农场，在烟草所的大棚温室内所栽的烟草上饲养繁殖3代以上，选取大小一致，活动正常的的蚜虫供试。 2.3. 试剂与仪器 甲醇(Methanol)，上海振企化学试剂有限公司生产；乙酸乙酯(Ethylacetate)，广东汕头市西陇化工厂生产；丙酮(Acetone)，上海凌峰化学试剂有限公司。所有溶剂均为分析纯。离心机(TGL8 Rrefrigerated Centrifuge)：Hema公司；分光光度计(DU730)：BECKMAN COULTER公司；恒温水浴锅(B-260)：上海亚荣生化仪器公司；冰箱(ULTRALOW)：SANYO公司；光照培养箱。移液枪，培养皿，广口瓶，毛笔等。 乙酰乙碱酯酶测定试剂盒、GSTs测定试剂盒、考马斯亮兰蛋白测定试剂盒。 2.4. 黄绿绿僵菌孢子悬浮液的配制 将保存在冰箱中的供试菌株转接马铃薯葡萄糖琼脂(PDA)平板培养基上(图1)，(25±1)℃、12 L：12 D的光照条件下培养7 d，待真菌充分产孢后，用灭菌的接种环轻轻将分生孢子粉刮到盛有5 mL 0.05%的 图1. PDA培养基上培养，菌落形状观察 Tween-80湿润剂的三角瓶中，在涡旋混合器上充分振荡10 min，待链珠状分生孢子断裂和分散后，以双层无菌纱布过滤菌液，获得孢子悬浮液母液。用血球计数板计数每毫升母液的分生孢子数，然后按比例稀释，配制成浓度为1.0 × 108个孢子/mL的孢子悬浮液供试。 2.5. 黄绿绿僵菌对桃蚜的接菌方法 用干净的小毛笔从大棚里饲养的桃蚜中采集大小一致的无翅若蚜，放入有甘蓝叶片的培养皿(d = 12 cm)中，用手动喷雾器将浓度为1.0 × 108个孢子/mL的孢子悬浮液均匀喷于供试虫体，每皿喷1.0 mL孢子液，以0.05%的Tween-80水溶液处理作为空白对照。用保鲜膜将培养皿口封住，用大头针在其上刺约20个小孔，既透气又可防止桃蚜逃逸。将处理的桃蚜置于(25±1)℃、相对湿度80%左右、光周期为L:D = 14:10 h的恒温光照培养箱内饲养。每次接菌的桃蚜数约为100头，重复3次。于处理后每间隔8 h取样一次，共取8次，以供制取酶液。 2.6. 桃蚜体内酶活性测试 2.6.1. 酶液的制备 试虫先用生理盐水漂洗2~3次。准确称取干净的待测试虫重量，按重量体积比(1:9)加生理盐水制成10%的组织匀浆，在4℃下以3500 r/min离心10 min，取上清液作为待测酶液。 以上酶液提取都是在冰浴中进行，所用生理盐水、匀浆器都需要进行遇冷处理，以防止酶失活。所提酶样均放入−20℃冰箱内储藏备用。 2.6.2. 酶活性的测定 AChE、GSTs和蛋白质含量的具体测定方法按照南京建成生物工程公司生产的试剂盒说明书进行。每处理重复3次。
system: [SYSTEM]

user: 3. 结果与分析 3.1. 黄绿绿僵菌对桃蚜AChE活力的影响 用108个孢子/mL的黄绿绿僵菌孢子悬浮液处理桃蚜，处理后每隔8 h取样一次，共取8次，测定桃蚜体内AChE的活性变化情况。由图2可以看出，在24 h内，处理组的蚜虫AChE的活力与对照无明显差异。24 h以后，其活力急骤下降，直到40 h时达到最低，抑制率为52.81%，其后一直维持在该水平。 3.2. 黄绿绿僵菌对桃蚜GSTs活力的影响 从图3可见，蚜虫体内GSTs的活力在染菌后的40 h内，与对照相比无明显变化，但从40 h后，其活力快速下降，直到56 h降到最低，抑制率为38.91%，并保持在该水平。 相比之下，黄绿绿僵菌对桃蚜体内AChE活性的抑制作用更大，抵制速度也更快。黄绿绿僵菌对桃蚜体内这两种酶活性的抑制是由于菌体本身侵染造成的，还是其代谢产物或者说绿僵菌毒素造成的？结合毒力测定现象和酶活测定结果推测，后一种的可能性更大，或者至少是二者共同作用的结果。在毒力测定中桃蚜的死亡速度较快，在染菌48 h甚至更早的时候后即出现死亡，对AChE活性的抑制速度较快，抑制程度较高，这些现象或结果说明，黄绿绿僵菌的毒素在对蚜虫的致死作用中起着重要作用。因此，该菌作为蚜虫的生防菌有着重大潜力。
system: [SYSTEM]

user: 4. 小结与讨论 桃蚜一般营全周期生活，也可以一直营孤雌生殖的不全周期生活，比如在北方地区的冬季，可在温室内的茄果类蔬菜上继续繁殖为害。桃蚜的繁殖很快，华北地区一年可发生10余代，长江流域一年发生20~30代 [ 2 ] 。本试验采用烟草上饲养繁殖3代以上的成蚜。 张仙红(2006)等 [ 3 ] 研究表明，被玫烟色拟青霉Paecilomyces fumosoroseus侵染的3、4龄菜青虫体内 图2. 桃蚜感染黄绿绿僵菌后体内AChE活性的变化 图3. 桃蚜感染黄绿绿僵菌后体内GSTs活性的变化 GSTs活力在感病前期显著高于同期未感染的菜青虫，而在感病后期明显低于同期未感染的菜青虫。苏筱雨等研究显示，感染球孢白僵菌后的桑天牛幼虫，其体内的GSTs活性则表现出先刺激后抑制的作用，而AChE活性影响则呈先抑制再激活再抑制再激活的趋势。王龙江(2010)等 [ 4 ] 研究发现，红火蚁感染白僵菌后，AChE的比活力呈先上升后下降的趋势，在接种后36 h时最高，之后不断下降。 本实验结果表明，黄绿绿僵菌处理对桃蚜AChE和GSTs活性的影响均表现为侵染初期抑制作用不明显，随着处理时间的延长均出现显著的抑制作用，并保持在较低水平，但最早出现抑制效应和达到最大抑制效应的时间不同，对于AChE分别出现在处理后的32 h后和40 h；而对于GSTs则分别出现在处理后的48 h和56 h。从毒力测定试验观察，黄绿绿僵菌处理后48 h，桃蚜开始逐渐死亡。这可能是因为处理后第24 h黄绿绿僵菌菌丝还没有侵染到桃蚜体内，48 h以后黄绿绿僵菌菌丝侵染桃蚜体内，并在桃蚜体内迅速生长增殖，黄绿绿僵菌菌丝造成虫体组织结构的破坏，显著抑制了AChE和GSTs的活力，最终导致桃蚜的死亡。 毒定测定显示，黄绿绿僵菌对桃蚜的入侵在48 h后即可引起桃蚜死亡，这既可能是由于菌丝的入侵抑制其体内重要代谢酶的活性，更可能是黄绿绿僵菌产生的毒素造成的。黄绿绿僵菌的入侵导致桃蚜体内AChE和GSTs的活力的活性受到抑制，是否与黄绿绿僵菌在入侵过程中伴随着毒素的产生有关？其具体机制还有待于进一步研究。这对于以后用黄绿绿僵菌对桃蚜进行生物防治有重要意义。
system: [SYSTEM]

user: 摘 要 本文中，基于单层多齿光栅结构，我们设计、优化了一种高性能的谐振光栅反射镜。严格耦合波分析算法表明，对于TE偏振光，此结构在1.47~2.1 μm宽谱范围内，其0级反射率大于98.5% (R 0 > 98.5%)，在中心波长1.8 μm处，角度带宽约为28˚。另外，我们还分析了此器件的工艺容差性能。 关键词 :衍射光栅，泄露模谐振，反射镜