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

user: 2. 电商大数据时代的现状 2.1. 电商时代发展简述 1996年前后，电子商务开始在中国起步。在21世纪初期，中国的电子商务迎来了第一次发展高潮，以当当网、卓越网为第一梯队，京东商城亦在中关村发迹。2010~2011年，团购网站也是如雨后春笋搬发展起来。自2013年开始，在中国互联网巨头中，以百度和阿里为代表开启了并购的风潮。2015年，根据全球十大电商市场数据盘点，中国的电商市场遥遥领先，阿里巴巴进军曼谷，口碑网投100亿推“全面开店”计划等等都说明了目前我国的电商发展还处于黄金时期。 2.2. 大数据时代电商面临的挑战 大数据时代的来临，对电子商务既是极大机遇，也是一系列挑战。 在数据集非常庞大的情况下，如何处理这些数据是基础。首先必须要对数据进行筛选和清洗，以便得到有价值和相关联的数据。然后才能进行客观实际的全面分析，摸清消费者的喜好。最后立足于这些喜好，设计并生产出迎合市场的产品。只有以消费者需求为中心，顺应市场，就能够真正实现销量与利润的双增长，最终促使电商得到更好的发展。目前一些企业已开始重视并开发大数据。
system: [SYSTEM]

user: 3. 网络指数的选取与概念说明 3.1. 网络指数的选取 网络指数有很多种，由于百度和阿里巴巴作为我国互联网巨头，而且电子商务的运营方式基本相同，数据采集和指数定义也相似，因此笔者只选择其中最具有代表性的百度指数和阿里指数作为研究对象。 百度指数，即百度搜索次数。如某词汇的百度指数为10万，那么就意味着网民通过百度搜索该词汇10万次。 阿里指数比较真实的反映了电子商务对数据的挖掘与分析，通过这些数据可以预测顾客的需求以及预测某一物品未来的需求程度。再大一点来说，可以反映电子商务的发展状况。与此相同，网络指数是对生活中各种现象的分析，是对社会民生的观察分析，可以预测人们生活状况和社会发展状况，对国家的各种政策的制定有重大的影响。 3.2. 两种网络指数的概念说明 3.2.1. 百度指数 百度指数是以全球最权威的中文检索数据为基础，根据搜索量和媒体检索量进行过滤和加权，通过科学、标准的运算，并且以直观的图形界面展现，帮助用户最大化的获取有价值信息。 百度指数能综合反映关键词在过去1天用户对它的关注和媒体对它的关注的一个参考值。任意关键词的百度指数都是该关键词在比较期的数值除以该关键词在基期的数值。比较期的数值和基期的数值是通过当天的用户搜索量和百度新闻中过去30天相关的新闻数量相比得来。 3.2.2. 阿里指数 阿里指数是阿里巴巴为了解电子商务平台市场动向而研发的数据分析平台，2012年11月26日，阿里指数正式上线。它是由阿里巴巴网站每日运营的5项基本数据统计计算得出，包括每天网站浏览量、每天浏览的人次、每天新增供求产品数、新增公司数和产品数这5项。 阿里指数主要包括如下三项： 1) 淘宝采购指数：根据在淘宝市场(淘宝集市 + 天猫)里所在行业的成交量计算而成的一个综合数值，指数越高表示在淘宝市场的采购量越多。 2) 1688采购指数：根据在1688市场里所在行业的搜索频繁程度计算而成的一个综合数值，指数越高表示在1688市场的采购量越多。 3) 1688供应指数：根据在1688市场里所在行业已上网供应产品数计算而成的一个综合数值，指数越高表示1688市场的供应产品越多。
system: [SYSTEM]

user: 4. 实证分析 4.1. 数据来源 在阿里指数官网上搜集并整理了智能手机行业2016年1月1日至2016年3月3日阿里指数里的淘宝采购指数、1688采购指数和1688供应指数 [ 7 ] 。 4.2. 统计建模分析 三组指数的数据是随时间变化的，因此笔者利用SAS软件编程 [ 8 ] 对数据进行时间序列分析。 4.2.1. 绘制时序图 图1中最上面红色表示的是淘宝采购指数的时间序列，绿色表示的是1688采购指数的时间序列，蓝色表示的是1688供应指数的时间序列。 4.2.2. 平稳性检验 为了判断序列是否平稳，除了需要考察时序图的性质，还必须对自相关图进行检验。下面笔者着重对淘宝采购指数的时间序列的平稳性进行分析。 图2横轴表示自相关系数，纵轴表示延迟阶数。从图2中发现序列的自相关系数衰减到零的速度非常缓慢，在比较长的延迟时期里，自相关系数一直为正，而后才开始变为负数。在自相关图中显出了一定的三角对称性，这是具有单调趋势的非平稳序列的一种典型的自相关图形式。这和图1中红色的曲线显示的显著的单调性是一致的，说明该指数的序列是非平稳序列。用同样的方法可以得出1688采购指数和1688采购指数的序列是平稳时间序列，从图1中也可以直观的看出，序列基本上在一个数值上随机波动。现将数据的简单描述汇总至表1。 4.2.3. 纯随机性检验 纯随机检验又称白噪声检验，此处显著性水平 取0.05。检验结果显示，无论是6阶，12阶还是24阶，检验统计量的 值都非常小(<0.0001)，所以可以以相当大的把握(置信水平 > 99.99%)断定1688采购指数和1688供应指数的时间序列为非白噪声序列(图3、图4)。 图1. 时序图 Table 1 Data analysi 指数 均值 方差 平稳性 淘宝采购指数 8892.3 2560.4 非平稳 1688采购指数 2226.0 378.1 平稳 1688供应指数 2101.7 55.8 平稳 表1. 数据分析 4.3. 模型拟合和预测 4.3.1. 模型拟合 1688采购指数模型定阶输出结果见表2，所以该模型为AR(1)模型时拟合最好。同时1688供应指数模型也同样为AR(1)模型时拟合最好。 经过检验，所有的P值均小于0.05，可以认为所有参数均显著。 图2. 淘宝采购指数自相关图 图3. 1688采购指数白噪声检验 图4. 1688供应指数白噪声检验 Table 2 Model fitting of 1688 procure inde Lags MA 0 MA 1 MA 2 MA 3 MA 4 MA 5 AR 0 11.80232 11.76387 11.6846 11.56017 11.4251 11.27005 AR 1 9.708679 9.774373 9.825117 9.852718 9.912487 9.955944 AR 2 0.768793 9.795009 9.84508 9.894549 9.901933 9.940308 AR 3 9.799704 9.833512 9.899254 9.950976 9.962964 10.00535 AR 4 9.808772 9.873049 9.934122 9.989545 10.01521 10.05808 AR 5 9.857694 9.918336 9.977967 9.995878 10.05778 10.07904 表2. 1688采购指数模型拟合 Table 3 Model error analysi 1 日期 3.4 3.5 3.6 3.7 3.8 2 真实值 2340 2342 2176 2317 2333 3 预测值 2366 2369 2272 2375 2378 4 相对误差 0.011 0.012 0.044 0.025 0.019 5 真实值 2044 2049 2054 2074 2043 6 预测值 2024 2053 2025 2083 2058 7 相对误差 0.009 0.002 0.014 0.004 0.007 表3. 模型误差分析 由此得到1688采购指数AR(1)模型为 (1) 1688供应指数AR(1)模型为 (2) 4.3.2. 模型预测 通过软件计算得到两个指数AR(1)模型2016年3月3日之后5天的指数预测值和真实值比较如表3所示。2~4行为1688采购指数模型，5~7行为1688供应指数模型。 模型拟合预测的相对误差 = |预测值 − 真实值|/真实值，从表3可知，相对误差均较小，说明模型拟合不错，预测效果比较好。
system: [SYSTEM]

user: 5. 结语 首先，淘宝采购指数的序列是非平稳时间序列，而1688采购指数和1688供应指数是平稳时间序列，进一步分析发现这两个还是非白噪声序列，说明还存在可提取有价值的信息，能够进行模型拟合和预测。 其次，笔者进行了模型拟合和未来五天的预测，得到了较好的预测效果。 最后，根据本文研究所得的结论和预测的结果，对1688采购商和1688供应商具有一定的参考价值。比如说淘宝采购商可以在淘宝采购指数较低时进行商品采购，这时采购人数一般较少。 当然，本文的研究还有很多不完善的地方。首先，本文所选取的样本量不是太多；其次虽研究了阿里指数和百度指数，但对智能手机只单纯考虑了阿里指数的数据，未能将两者有机地结合起来。这些不足对未来该领域的研究也许提供了一定的方向。
system: [SYSTEM]

user: 摘 要 本文在某型通用飞机的基础上设计了一个油电混合动力系统，并用MATLAB软件模拟验证该系统的实用性及优越性。根据运动学公式，使用MATLAB软件编程分别对原型机和该混合动力系统飞机的起飞距离、陆上行驶能量消耗、陆上行驶路程等性能进行模拟分析。结果表明，使用该动力系统的油电混合动力飞机，在上述各项性能中均优于原型机。 关键词 :油电动力系统，MATLAB，通用航空，混合动力
system: [SYSTEM]

user: 1. 引言 目前，活塞式航空发动机技术已趋于成熟，在环境和能源问题上很在再取得较大的改进 [ 1 ] 。而电动飞机却恰恰相反，具有环保的油电，但是由于电池技术发展缓慢，能量密度不高，很慢满足飞机队动力系统功重比的要求，因此还不能大规模应用。由于在汽车行业 [ 2 ] ，混合动力汽车在一方面融合了传统汽车成熟的技术和现有的工业基础，另一方面融合了电动汽车低油耗和低排放的特点，有效降低污染，并能节省燃油消耗，达到改善环境和节约能源的目的。所以，以电机作为辅助动力，活塞式发动机作为主动力源的混合动力飞机出现了。
system: [SYSTEM]

user: 2. 油电混合动力系统设计 该飞机动力装置为一台100 hp的Rotax公司生产的912ULS四冲程活塞发动机。本设计就是研究一个适合的油电混合动力系统，取代912ULS发动机，从而达到节能减排的目的。本次设计的混合动力系统框架图如图1所示，主要部件为发动机、蓄电池、电动机、离合器以及螺旋桨。电动机和车轮相连，发动机和螺旋桨相连。电机可作为单独动力驱动车轮行进。
system: [SYSTEM]

user: 3. 发动机、电机和电池的选择 3.1. 发动机的选择 在本文设计的混合动力系统中，活塞发动机依然是主动力源，发动机的选择需要考虑多方面的因素。其中最重要的是能够单独提供飞机巡航所需要的功率，且所选择的活塞发动机质量要求尽可能小。飞机平飞时的发动机功率计算公式 [ 3 ] ： (1) 式中， 为螺旋桨效率，K为最大升阻比， 为最小平飞速度。计算得所需最小功率为31.79 Hp，则本文所选的发动机的功率需大于31.79 Hp。对比各种现有发动机性能、质量与应用环境，最终选择81 Hp的ROTAX912-A发动机。 3.2. 电机的选择 电机在滑跑起飞阶段时提供辅助动力，功率则需满足： (2) 图1. 油电混合动力框架图 式中， 为活塞发动机功率， 为电动机功率， 为滑跑时阻力所做的功。 (3) (4) 将数值代入计算得， 。综合考虑电机外形尺寸、功率等因素，本次设计选用电机为斯洛文尼亚的EMRAX207电机，该电机最大持续功率为30 kW。 3.3. 电池的选择 根据此次动力系统的设计要求，电机在起飞滑跑时，提供辅助动力。且在陆上行驶时，提供全部或者部分动力。在飞机爬升和平飞时，电机不参与工作。所以电池储存的能量只需满足在起飞滑跑和陆上行驶的消耗即可。起飞滑跑阶段，为了缩短滑跑距离，动力系统必须以最大功率滑跑，此时电机和发动机均在其各自的最大功率状态下，电池容量可按照公式计算： (5) 式中， 为电机最大持续功率， 和 分别为起飞滑跑时间和陆上行驶时间， 为电。效率。为了不影响电池寿命，电池内不许保证有20%的剩余电量，因此在公式计算中加入了除以80%。设定陆上由电机驱动行驶的时间为2小时，由该公式计算得出的电量为460 Ah，电池可以从生产厂家订购，参数为60 V，460 Ah，重量约为40 KG。
system: [SYSTEM]

user: 4. 混合动力飞机滑跑起飞距离模拟 本文以飞机最大起飞的重量 进行计算。飞机所受的空气阻力为 ，地面摩擦力为 ，升力为 [ 3 ] [ 4 ] 。 (6) (7) (8) 为阻力系数， 为升力系数。本文所得升、阻系数可由红嘴鸥飞机阻力极曲线查得。 螺旋桨可用 拉力为： (9) (10) 式中， 为动力装置功率， 为螺旋桨直径。 当海拔过高引气空气稀薄，活塞发动机缸内燃烧过程因此而受到抑制，根据经验公式(11)，可以求出在不同海拔和不同温度下较为接近的发动机实际功率。 (11) 式中， 为发动机功率， 、 分别为实际压强和气温，下标“0”和“H”分别表示为发动机处于海平面和实际工作高度H处。 在海拔6000米以下，海拔高度与大气密度有密切关系，可以按照如下公式求得实际飞行高度下的大气密度： (12) 飞机在跑道滑跑起飞时的平均阻力为： (13) 式中 ， 为飞机离地瞬间的升阻比。根据动力学公式，飞机在跑道滑跑起飞时的平均加速度为： (14) 所以，可得飞机起飞滑跑距离为： (15) 式中， 为飞机滑跑起飞时瞬时速度。根据以上公式运用MATLAB软件，对原型机和混合动力飞机的起飞滑跑距离进行编程计算可以得到不同机场海拔高度时的滑跑距离，如图2所示。
system: [SYSTEM]

user: 5. 混合动力飞机陆上行驶性能模拟 假定飞机在 的水泥公路，且天气环境对飞机的行驶状态无影响，行驶时由于速度的不同，所需功率也会不同，电池消耗电量和发动机的耗油量也随之不同。 在陆上行驶时，功率计算公式为， (16) (17) 式中， 为陆上行驶时的阻力， 为陆上行驶时的速度。设定发动机功率转化为螺旋桨功率的效率70%。 由螺旋桨单独驱动飞机行进和由螺旋桨和电机共同提供动力驱动飞机行进时所能提供的速度模拟结果如图3所示。 由于行驶时由于速度的不同，所需功率也会不同，电池消耗电量和发动机的耗油量也随之不同。根据燃油消耗量计算公式 [ 5 ] ， (18) 式中，B为燃油消耗量，b为燃油消耗率，P为输出功率，d为燃油比重。在起飞过程中，功率P为最大 图2. 不同海拔高度下的滑跑距离 图3. 不同速度下的功率 功率。当飞机由螺旋桨驱动，在陆上行驶时，功率为， (16) (17) 式中， 为陆上行驶时的阻力， 为陆上行驶时的速度。 由螺旋桨牵引飞机在陆上行驶总路程计算公式： (19) (20) 式中，Q为油箱总量。本文只是针对动力系统进行改装，并未变油箱大小。 由螺旋桨和电机共同驱动飞机在陆上行驶总路程计算公式： (21) (22) (23) 在陆上行驶的过程中，由ROTAX 912ULS活塞发动机作为动力源，由螺旋桨牵引飞机行进；与由ROTAX 912-A活塞发动机与电机同时作为动力源，共同驱动飞机行进的路程计算如图4所示。 当在陆上行驶过程中，由于油量消耗而释放出来的能量为： (24) 式中， 为质量， 为单位热值。 为燃油完全燃烧释放的热量。 当在陆上行驶过程中，电机所消耗的能量为： 图4. 不同速度下的总路程 (25) 式中， 为发动机实际功率， 为发动机实际功率的工作时间。 当单独由电机提供动力，飞机在陆上行驶时的功率为： (16) (17) 单独由螺旋桨提供动力在陆上行驶，与单独由电机提供动力在陆上行驶时，不同速度状态下能量消耗与路程的关系模拟结果如图5所示。 理想状态下，飞机陆上性能模拟结果如表1所示。
system: [SYSTEM]

user: 6. 结论 从上述模拟结果得出，在飞机载重不变的情况下，飞机的起飞滑跑距离缩短了12.5%。陆上行驶最 图5. 不同速度下能量消耗与路程的关系 Table 1 Simulation results contras 原型机 混合动力机 起飞滑跑距离/m 128 112 陆上行驶最大功率/kw 51.45 53.28 空重/kg 280 321.6 有效载荷/kg 270 220 最大起飞重量/kg 550 550 陆上行驶总路程/km 110 196 表1. 模拟结果对比 大功率虽然相差不大，但是在相同的燃油消耗量的前提下，陆上行驶总路程增大了78.1%。虽然有效载荷略微减小，但是对机场跑道长度要求下降了。并且可由电机单独驱动在陆上行驶，消耗较少的能量，行驶更大的路程。通过模拟结果，本次设计的混合动力飞机与原型机相比，起飞滑跑距离短，能源消耗少。适用于飞行训练、旅游观光，且节能减排。
system: [SYSTEM]

user: 摘 要 铜作为一种重要的有色金属，对其他有色金属的价格具有较强的影响力，因此期货铜市场的稳定发展不仅对现货铜市场具有重要作用，同时对其他有色金属市场的稳定，促进国民经济的又好又快发展，都具有重大意义。本文运用两状态的Markov状态变换模型对期货铜的收益率和大幅波动情况进行了预测，取得了良好的效果。 关键词 :铜期货，Markov态变换模型，复合似然函数
system: [SYSTEM]

user: 1. 引言 铜作为一种重要的工业重金属，被广泛地应用于电气、轻工、机械制造、建筑工业等领域，在国民经济的发展中有着重要的地位。但是，自2008年金融危机以来，铜的价格波动十分频繁，给铜的生产和销售环节带来了巨大的风险，对我国经济的平稳运行产生了巨大的冲击。因此，对于铜的价格和波动预测具有十分重要的理论和现实意义。 Markov状态换模型由Hamilton (1989) [ 1 ] 首先提出，他在研究美国的经济周期时，运用三状态两阶滞后的Markov机制转换模型研究了美国1953~1984年间季度实际产出增长的波动，模型回归效果良好。Cai (1994) [ 2 ] 、Hamilton、Susmel (1994) [ 3 ] 几乎同时提出了将Markov状态换模型引入ARCH模型中的方法，这一创举使得对变结构金融资产序列的波动性建模取得了显著的进步和发展。孙金丽、张世英(2003) [ 4 ] 将MS-GARCH模型运用到中国股票市场的分析中，证明具有结构转换功能的MS-GARCH模型的分析结果由于单一状态条件下的GARCH模型。朱钧钧 [ 5 ] (2011)利用马尔科夫状态转换–门限GARCH模型的MCMC估计对我国股市波动率的双重不对称性进行了解释，得到了稳定并且可信的结果。 本文共分为以下几个部分：第一章、前言，介绍本文的背景、框架和内容；第二章、理论方法综述，介绍本文所采用方法的相关理论介绍；第三章、模型分析，介绍模型的分析结果，重点分析模型的收益率、涨跌方向预测准确率、大幅波动的预测情况；第四章、模型的缺陷和不足，介绍模型中仍存在的缺陷和未来改进的方向；第五章、总结，对分型进行概括和总结。
system: [SYSTEM]

user: 2. 理论方法综述 2.1. Markov链 设随机过程 只取整数值 ，假定， 的取值只与 有关，而与t − 1期前S的取值无关，那么此随机过程被称为N状态下的马尔科夫链，用公式可表示为： 转移概率 代表状态 条件下 的概率，若该过程服从一阶Markov链过程，可用转移矩阵P来表示转移概率： 其中， ， 。 2.2. Markov模型 下面以正态分布为例，构建期货黄金指数的Markov模型，假定t时期的期货黄金指数为 ， 。假定存在一个不可观测变量 代表t时刻期货黄金指数所处状态，当 ， ， 的边缘分布为 ， 服从N状态下的一阶马尔科夫链。 Markov模型中需要估计的参数有转移概率 ， ， 。我们采用EM算法去进行参数估计，关于似然函数的选择有两种方式：完全似然函数Zucchini (2009) [ 6 ] 、复合似然函数Zou (2013) [ 7 ] 。不过在Zou (2013) [ 7 ] 已经说明复合似然函数效果要好于完全似然函数。 在正态分布条件下，复合似然函数可表示为： 其中， 表示 时刻为 状态而 状态。 表示 属于状态μ状态的概率。 每时刻的状态都是未知的，我们可以将其作为缺失变量。因此可以利用EM算法估计参数 ， ， 。利用所预测出来的参数，通过平滑概率，可以计算出当t = T时， 属于每个状态的概率。我们也可以通过参数P、 、 与第t时刻的观测值 可以预测出当 时， 属于每个状态的概率，均值和方差。 本文运用惩罚下复合似然函数Chen (2012) [ 7 ] 作为似然函数 在这里 为序列的样本方差。由Zou (2013) [ 7 ] 可知在这种方法下迭代步骤分为E-STEP和M-STEP两个步骤： 在E-STEP，我们运用贝叶斯公式计算条件期望值 其中 。在EM算法中，我们假定这个变量为隐藏变量。 在M-STEP，我们最大化似然函数，则可得： ， ， 其中 和 ，其中 。
system: [SYSTEM]

user: 3. 模型分析 本章我们将利用Markov状态转换模型对2001年1月2日至2015年12月21日的上海期货铜指数日收益率进行分析。我们利用过去交易日的数据结合Markov状态转换模型对下一交易日上海期货铜指数日收益率进行预测，分别预测出下一交易日上海期货铜指数日收益率的均值、方差、上涨(下跌)概率，并结合前一交易日的持仓量和交易量对下一交易日大幅波动的情况进行预测 3.1. 日收益率数据分析 我们使用上海期货铜指数日收益率作为研究对象，在这里我们用指数益率作为收益率的代表，令 ，其中， 代表第t天上海期货铜指数的收盘价，为此我们可以得到关于日收益的时间序列图，见图1所示。 可以看出，2001年~2015年上海期货铜日收益率大都在0附近徘徊，同时，收益率存在着一定的波动聚集性，这主要是指在金融时间序列中，一个较大的波动后往往会跟着一个较大的波动，而一个较小的波动后往往跟着一个较小的波动，这反映的是金融波动的正相关性和正反馈效应。因此，通过过去一个 图1. 2001~2015年上海期货铜日指数收益率 时间段上海期货铜的日交易数据对下一交易日上海期货铜的日收益率进行预测具有一定的合理性。 3.2. 模型构建 我们利用Markov状态转换模型对上海期货铜的日指数收益率进行分析，构建模型。我们假定每一个交易日的上海期货铜的指数收益率服从正态分布，且在训练期间内，上海期货铜的指数收益率存在两种波动状态，即上海期货铜的指数收益率服从两个不同的正态分布，且两种状态之间的转换过程服从一阶马尔可夫过程。 此时，复合似然函数可表示为： 其中， 。 通过复合似然函数，利用EM算法，可以估算出转移概率 ， ， 。利用预测出来的参数，通过平滑概率，可以计算出当t = T时， 属于每个状态的概率，再结合转移概率 ，可以预测出当 时， 属于每个状态的概率，通过状态参数 、 可以预测数 的均值和方差，通过均值和方差，我们可以计算出 大于或小于0的概率(即上涨或下跌的概率)以及 大于0.02或小于−0.02大概率(即出现大幅波动的概率)。 3.3. 模型效果 我们通过以下3个方面来衡量Markov状态转换模型的分析效果，一、模型的收益率情况，二、模型预测的准确率情况，三、模型对大幅波动的预判情况。 关于模型训练区间的选取并没有特定的标准，由于模型当中需要估计6个参数，训练区间不宜选取太短；同时，由于假定上海期货铜的指数收益率存在两种波动状态，训练区间不宜选取太长。为此，我们分别考虑100个交易日、150个交易日、200个交易日的情况。 3.3.1. 模型的收益率 为了考察模型的收益率，我们按照模型预测出上涨或下跌的概率进行投资，当上涨概率大于50%时买入，上涨概率小于50%时卖出，模型2001年1月21日~2015年10月21日的收益率如表1所示，在3种不同训练期间情况下，模型的平均收益率均在10%以上，在训练期间等于150个交易日时，平均收益最高为18.16%。 3.3.2. 模型的准确率 为了考察模型的收益率，我们依然考查的是通过Markov状态转换模型预测出的下一交易日上涨或下跌的概率，我们设定一个阈值Ɵ (Ɵ ≥ 0.5)，当上涨(下跌)概率 ≥ Ɵ时，我们认为下一交易日存在上涨(下跌)，然后，我们将实际情况与我们的预测情况进行比较，分析模型的预测准确率。图2~图4所示为Ɵ在0.5~0.6时，模型预测的准确率的情况。 Table 1 2001~2015 three model yield 2001~2005 2006~2010 2011~2015 平均收益率 100交易日 21.39% 33.51% 0.28% 17.57% 150交易日 21.78% 26.28% 7.28% 18.16% 200交易日 25.98% 12.69% 4.60% 14.08% 表1. 2001~2015三种模型收益率情况 图2. 训练期间为100个交易日时模型的预测准确率 图3. 训练期间为150个交易日时模型的预测准确率 图4. 训练期间为200个交易日时模型的预测准确率 可以发现，模型对下一交易日上涨行情的预测要优于下跌行情的预测，但随着阈值的不断提高，无论是总体预测准确率还是上涨行情和下跌行情的准确率都不断提高，最终模型预测的总体准确率都在60%以上。 3.3.3. 模型对大幅波动的预判 为了考察模型对大幅波动的预判，我们运用预测出的下一交易日收益均值、方差、以及由此计算出的下一收益率的绝对值大于2%的概率、前一交易日的交易量和持仓量，通过随机森林模型对下一交易日收益率的绝对值大于2%的情况进行预测。 我们以训练期间为150交易日的模型为例，训练集选用2014年前150个交易日，测试集为2014年8月13日~2015年10月21日。表2所示为模型对测试集大幅波动预测的情况。 可以发现，模型对大幅波动的预测准确率为65.22%，整体准确率为85.2%，表明利用预测出的下一交易日收益均值、方差、以及由此计算出的下一收益率的绝对值大于2%的概率、前一交易日的交易量和持仓量等信息可以在一定程度上对下一交易日的大幅波动进行预警。
system: [SYSTEM]

user: 4. 模型的缺陷和不足 虽然模型在收益率、准确率和大幅波动的预警等方面表现良好，同时值得注意的是，模型当中依然存在着一些不足，例如，我们模型中假设上海期货铜的指数收益率服从正态分布，但金融数据中往往存在着“尖峰厚尾”的情况，我们的模型会在一定程度上低估模型的风险。如图5所示，我们模型预测出的收益率较实景情况中的收益率更多的集中在0附近，对大幅波动的预测存在一定的偏差。为此，模型可以通过将正态分布调整为t分布来更好地拟合这一情况。
system: [SYSTEM]

user: 5. 总结 本论文建立了两状态的马尔可夫转换模型，同时利用复合似然函数方法以及EM算法来对铜期货的数据进行模拟预测，模型的分析结果良好，同时模型具有以下优势： 第一，模型的收益率良好。我们用模型对过去100天，150天和200天的交易日的数据进行训练，对下一交易日收益率大于0的概率进行预测，利用这个概率进行投资，从2001~2015年15年的交易数据来看，3个模型的年平均收益率均在10%以上，以150个交易日为训练期间的模型表现最为良好，模型的年平均收益率为18.16%。 第二，模型的准确率比较高。随着阈值的不断提高，模型的准确率也不断上升，从最初的百分之五十左右一直上升到了百分之六十之上，也就是说通过这个马尔可夫模型预测出来的第二天涨跌的概率具有一定的实际和理论分析意义。 第三，同时我们利用随机森林的方法，结合通过马尔可夫模型预测出的下一交易日收益均值、方差、以及由此计算出的下一收益率的绝对值大于2%的概率、前一交易日的交易量和持仓量，对第二天波动幅度超过2%的情况的预测。我们利用2014年前150个交易日，预测2014年8月13日~2015年10月21日的交易 Table 2 The predict situation of model to large fluctuation 波动情况 交易天数 预测准确率 总体情况 288 85.42% 波动浮动超过2% 23 65.22% 波动浮动未超过2% 265 87.17% 表2. 模型对大幅波动的预测情况 图5. 实际收益率与预测收益率比较 情况，预测结果表示，模型可以准确的预测出其中65.22%的大幅波动情况，因此，我们可以认为这个模型对出现大幅波动的预测情况良好。 但值得注意的是，模型依然存在着一定的缺陷，由于我们的模型假设是正态分布，通过图5可以看出，通过模型所预测出来的大部分收益率的波动情况跟实际情况而言较窄，也就是说我们的模型在一定程度上低估实际情况中期货市场的各种风险。由于低估了这些风险，会使得我们的模型对第二天收益率的估计造成一定的风险。我们可以通过将模型假设改为t分布或其他的厚尾分布从而改进这个模型。 综上所述，虽然这个模型有着一定的缺陷，但是该模型的收益率良好，准确率也达到了百分之六十以上，同时对期货市场出现大幅波动的预测结果良好，我们可以认为这个模型在预测铜期货的收益率有着一定的实际用处，同时可以根据所预测出来的大幅波动来对市场的投资进行一定的调整。
system: [SYSTEM]

user: 摘 要 本文得到了带有弱化表面的颗粒增强复合材料的本构模型，其可以描述颗粒的脱黏破坏以及基体的塑性和颗粒对复合材料损伤的影响。此损伤模型是基于Eshelby等效夹杂方法和自洽方法，假设基体、完整颗粒和带有弱化表面的颗粒都属于复合材料的夹杂体，根据Eshelby细观力学的方法，通过数值计算估算出三相颗粒增强复合材料的有效弹性模量。得出在复合材料损伤过程中，颗粒体积分数变化对颗粒复合材料有效刚度的影响，并且通过分析计算说明颗粒的体积分数与尺寸对复合材料应力–应变关系的影响。 关键词 :颗粒复合材料，弱化表面，体积分数，等效刚度，表面脱黏
system: [SYSTEM]

user: 1. 引言 复合材料的有效弹性模量在分析和设计复合材料中起着很重要的作用，并且在复合材料细观力学分析中得到了广泛的应用 [ 1 ] ，因此很多的研究人员尝试着应用数值方法和理论方法来估算复合材料的有效弹性模量，解释了表面脱黏损伤是控制局部应力场和复合材料总的力学特性的重要因素 [ 2 ] 。复合材料中一个最基本的问题是预测单相或多相材料有效或平均宏观性能的特性和关系，为了研究表面脱黏对复合材料的力学性能的影响，一些学者对颗粒增强复合材料的三相模型进行了研究。本篇论文的三相模型假设是由基体、完整颗粒以及带有弱化表面的颗粒组成。 Togho等人 [ 3 ] 、Togho和Weng [ 4 ] 、Togho和Chou [ 5 ] 基于Eshelby等效夹杂原理和Mori-Tanaka平均场的概念建立了增量模型，解释颗粒的脱黏破坏过程，基体的塑性和颗粒的尺寸对复合材料的变形和破坏的影响。Bian和Wang [ 6 ] 、Wang等人 [ 7 ] 基于Eshelby等效夹杂原理和自洽理论研究了具有两相夹杂的混凝土复合材料中颗粒的体积分数对等效刚度的作用。Niordson和Tvergaard [ 8 ] 应用有限单元方法基于应变梯度塑性计算了单个颗粒的非连续性增强复合材料的颗粒尺寸对材料损伤的影响。Nan和Clarke [ 9 ] 推出了颗粒的尺寸和断裂对金属基复合材料的弹塑性变形的影响，并且延伸有效介质方法解释了颗粒的尺寸对复合材料基体的应力–应变关系的影响以及颗粒的破坏准则。 由于颗粒的断裂或者颗粒与基体之间的脱黏都会产生材料的变形，Prangnell等人 [ 10 ] 分析了金属基复合材料的弹性相应，他是通过非连续的金属基复合材料的变形来进行描述。Zhao和Weng [ 11 ] 、Ju和Lee [ 12 ] 通过考虑颗粒的表面脱黏相来研究复合材料细观力学的破坏模型。Zhong等人 [ 13 ] 研究了三相尺寸具有非完整界面颗粒的颗粒增强复合材料的微观力学模型，得出的结论是当位移通过表面牵引的比例增加时，若牵引保持连续，则非完整表面就以弹簧模型为特征。 微观力学模型的提出能够有效的预测颗粒增强复合材料的有效弹性性能和弱化表面的演变过程，假设材料在未受到外界力的作用时，复合材料中的颗粒都是完整的，当有外界应力时，材料发生变形进而产生带有弱化表面的颗粒。为了在复合材料中存在具有弱化表面的球形颗粒，本文采用了模型 [ 14 ] [ 15 ] ，其模型具有轻微的弱化表面椭圆夹杂相的Eshelby’s张量，基于Eshelby’s细观力学，可以获得三相颗粒增强复合材料的有效弹性模量。Bian和Cheng [ 16 ] 延伸到三相复合材料(颗粒、基体和界面)，分析得出复合材料的有效弹性模量。Lee [ 17 ] 描述了三相颗粒复合材料(基体、完整球形颗粒以及带有弱化表面的球形颗粒)的微观本构模型，利用Weibull [ 18 ] 概率函数计算出非完整颗粒的体积分数，画出单轴、双轴和三轴载荷作用下颗粒1增强复合材料应力–应变关系。 在本篇论文的模型中，利用Eshelby [ 19 ] 等效夹杂原理和自洽原理建立了颗粒强化复合材料的有效体积和剪切模量。基于上述两种方法说明了完整颗粒到非完整颗粒转变的过程(脱黏破坏)，得出颗粒的体积分数对等效刚度的影响。通过估算优化的等式来描述颗粒的脱黏破坏对等效刚度改变的规律，同时也说明了颗粒的体积分数和尺寸对复合材料应力–应变关系的影响。
system: [SYSTEM]

user: 2. 复合材料的模型分析 2.1. 模型的建立 图1(a)表示了颗粒增强复合材料在单向拉伸时的起始状态，图1(b)则表示复合材料在单向拉伸受力时的破坏情况。假设复合材料颗粒的原始体积分数为 ，当外界应力作用后完整颗粒的体积分数为 ，体积和剪切模量分别为K1、G1，带有弱化性的体积分数为 ，体积剪切模量分别为K2、G2，脱黏的体积分数为 。 2.2. 等效模型的增量理论 2.2.1. 等效模型 图2(a)作为最初的原始模型是一个均匀各项同性无限大的基体中包含了单个颗粒，基体和颗粒的弹性模量分别是 和 。图2(b)则基于Eshelby等效夹杂的方法，利用 来代表颗粒的等效应变并且基体的弹性模量还是 。在图2(c)模型中，我们把基体和颗粒都当做复合材料的夹杂体，经过修改过的模型就变成了等效基体的弹性模量是 (基体颗粒弹性模量的平均值)，等效应变是 。本文中复合材料模型包含的三相分别是完整颗粒、带有弱化表面的颗粒和基体，首先本文得到 图1. 单向拉伸作用下颗粒增强复合材料：(a) 原始状态；(b) 破坏状态 图2. 复合材料不同单元的细观力学模型 (1) 式中：0——基体，1——完整颗粒，2——非完整颗粒， ——等效应变。 (2) (3) (4) (5) 式中： 分别是基体、完整颗粒的初始体积分数，K、G分别是复合材料的体积和剪切模量。 因为有弱化表面相的存在，所以组成颗粒强化复合材料的Eshelby张量有两个，即 (6) (7) (8) (9) (10) (11) (12) (13) 式中： ——颗粒的直径， ——基体的拉梅常数， ——基体的剪切模量， ——泊松比， ， ——代表了切向和法相的界面的参数。 把式(7)~(13)带入到公式(6)即可获得 的表达式。 由公式(1)得： (14) (15) (16) (17) 式中： ， ， ——应变及应力球张量， ——应变及应力偏张量。 2.2.2. 在外载荷作用下的增量本构关系 基于前面模型的假设，颗粒的有效弹性模量可以写成： (18) 由式(14)和式(18)可得： (19) 根据胡克定律得： (20) 式中： ， 是颗粒等效刚度的体积和剪切模量。 根据公式(18)得： (21) 由式(19)及式(20)可得： (22) (23) 当完整颗粒和带有弱化表面颗粒完全脱黏时，就会产生脱黏引起的空洞，由公式(1)得： (24) (25) 即 。 把公式(1)、(24)代入公式(21)中简化为： (26) 同理 。
system: [SYSTEM]

user: 3. 复合材料的有效弹性模量 复合材料总的本构关系可以写成： (27) 表示复合材料总的有效弹性模量， 表示空洞的等效弹性模量。 基于Reuss理论，总的有效模量 可以写成： (28) 式中： ， ——分别是基体和颗粒的等效弹性模量。 根据公式(22)、(23)和(26)得： (29) 把公式(29)代入(28)中可得： (30) 由公式(19)得出： (31) (32) (33) (34) 根据本模型可以有效的预测复合材料的有效弹性模量。把公式(2)~(3)、(14)~(18)、(21)~(23)、(26)、(31)~(34)代入公式(30)可得： (35) (36)
system: [SYSTEM]

user: 4. 理论分析结果 为了说明我们提出的颗粒增强复合材料细观力学弹性破坏模型，我们考虑了单轴拉伸的SiC颗粒与Al合金的应用计算分析。我们所采用的Al合金基体是6061-T6，具体的材料性能参数见表1。 表1指出了复合材料两种夹杂体的杨氏模量和泊松比。本文的基本模型虽然是三相夹杂复合材料，但由于完整颗粒和带有弱化项颗粒本身的力学性能没有发生变化，因此只显示两种夹杂。假设原始颗粒 的体积分数为 ， ， ，已知基体的剪切模量 。 Table 1 Material properties of particle composite Parameter (参数) Matrix (6061-T6合金) Particles (SiC) Young’s modulus (杨氏模量)/GPa 68.3 490 Possion’s ratio (泊松比) 0.33 0.17 表1. 颗粒复合材料的力学特性 根据公式(35)可以得到有效体积模量与体积分数之间的变化关系，具体如图3所示。 图3(a)表示的是有效体积模量与空洞的体积分数和带有弱化表面颗粒的体积分数之间的关系。当我们只考虑带有弱化表面颗粒和空洞的体积分数对有效体积模量的影响时，结果如图3(b)和图3(c)所示。图3(b)表示的是当给定空洞的体积分数时等效弹性模量 与 之间的关系，很明显当空洞的体积分数 是定值时，等效弹性模量 随着 的增大而减小。而且带有弱化表面颗粒的分数越大， 减小的就越明显。图3(c)表示的是当给定带有弱化表面颗粒的体积分数为常数时，等效弹性模量 与 的关系。当 给定时， 随着 的增大而减小，且 越大 减小的越迅速。 根据公式(36)可得有效剪切模量与体积分数的关系，具体如图4所示： 图4(a)表示的是复合材料等效剪切模量与空洞的体积分数和带有弱化表面颗粒的体积分数之间的关系。如果我们只考虑单个体积分数(带有弱化表面颗粒或空洞)对等效剪切模量的影响时，就如图4(b)和图4(c)所示。图4(b)表示的是当给定空洞的体积分数时复合材料的等效剪切模量 与 之间的关系。显然当空洞的体积分数 是定值时，等效剪切模量 随着 的增大而减小。而且带有弱化表面颗粒的分数越大， 减小的就越明显。图4(c)表示的是当给定带有弱化表面颗粒的体积分数为常数时，等效剪切模量 与 的关系走向。当 给定时， 随着 的增大而减小，并且 越大， 减小的越迅速。
system: [SYSTEM]

user: 5. 颗粒尺寸与体积分数的影响 根据Ramburg-Osgood等式，复合材料基体的应力–应变关系如下 [ 9 ] ： (37) 式中： 分别表示等效应力、应变， 表示基体的杨氏模量、等效材料的屈服应力和工作强度指数， 。其中， 表示纯基体的初始屈服应力(86 MPa)及基体的强度。 (38) 上述公式中 是基体的剪切模量，b (=0.283 nm)表示柏氏矢量，d为颗粒的直径尺寸， 是一些常数，1.0，0.4，1.0，2.0， 表示温度的改变量， ( )热膨胀的不匹配度 [ 20 ] 。 基体的等效塑性应变 [ 4 ] ， 表示材料常数3/7。如果考虑基体为Al ( ， )，颗粒SiC ( ， )，Ramburg-Osgood参数A356-T4的 ，n = 0.212。颗粒的原始体积分数为 。 图5和图6结果显示了Togho和Itoh [ 3 ] 、Nan和Clarke [ 9 ] 等的结果与本文计算方法的比较。图5说明了颗粒的尺寸对复合材料的影响情况，本文的模型是基于Eshelby等效夹杂和自洽方法的理论基础 图3.有效体积模量与空洞体积分数和带有弱化表面颗粒的体积分数之间的关系 图4. 等效剪切模量与空洞体积分数和带有弱化表面颗粒的体积分数之间的关系 图5. 颗粒尺寸对材料的应力–应变关系的影响 图6. 不同的计算方法的本构关系与实验结果的对比 建立的。此模型中颗粒的尺寸越大对材料的应力、应变关系的影响越小，尺寸越小影响就显得明显许多。其中ID model是Togho等人 [ 3 ] 所描述没有尺寸影响的模型。 图6显示了本文的计算方法相较于其他的结果与实验结果更吻合。首先从建模的角度来说：Togho和Nan等人所采用的模型都是基体、颗粒和最后脱黏而形成的空洞，空洞则是颗粒的脱黏形成。本研究论文的模型则采用的是三相模型：基体、完整颗粒、弱化表面相颗粒以及脱黏形成的空洞。此文中的空洞是有两部分组成，完整颗粒的脱黏形成空洞和弱化表面相颗粒脱黏形成的空洞。再者从运用的方法不同来说：Togho所运用的主要是Mori-Tanaka’s法与等效夹杂的相结合，而Nan等人所运用的是断裂破坏理论。本篇论文所采用的思想是等效夹杂和自洽方法的结合来计算复合材料的本构关系，根据自洽的基本定义，对此模型的修正即把颗粒与基体都作为复合材料的夹杂相。Togho和Nan等人考虑了颗粒脱黏直接到空洞的过程，而本文考虑了由完整颗粒到弱化表面颗粒最后到空洞脱黏的整个过程。所以，基于此模型建立的基本理念，理论分析显示本文计算结果更接近于实验数据。
system: [SYSTEM]

user: 6. 结论 本文通过修改的复合材料细观力学模型，基于Eshelby等效夹杂和自洽方法，可以估算出三相复合材料的等效弹性模量。同时根据自洽方法提出的等效模型的有效等式，利用该等式可以预算出颗粒增强复合材料的体积和剪切模量。 当复合材料在没有受到外界应力或者受到的力不足以使颗粒发生损伤破坏时，颗粒的体积分数增加有助于复合材料整体的力学性能提高。当受的应力足以使颗粒发生损伤破坏，其破坏产生的带有弱化表面颗粒的体积分数和空洞的体积分数对复合材料的力学性能有很大的负效应。当带有弱化表面颗粒的体积分数增大时，复合材料的等效刚度会随完整颗粒的体积分数的减小而减小。当完整颗粒和带有弱化表面颗粒完全脱黏时，这时产生空洞的体积分数就会明显的增加，使得复合材料的力学性能减小。 在复合材料的损伤过程中，应变会随着损伤颗粒体积分数 的增大而减少，即随着破坏的加大，材料的承载能力也会逐渐下降，当复合材料完全破坏时( )直至材料完全损伤。并且，计算结果表明相较于其他的模型本文模型与实验的结果更吻合。
system: [SYSTEM]

user: 摘 要 [目的]：观察头盆环在极重度脊柱畸形患者的治疗过程中的安全性及有效性。[方法]：2011年9月~2012年6月住院患者7例，其中脊柱结核后凸3例(静止期)，青少年特发性脊柱侧凸4例。Cobb’s角在136˚~180˚，平均(168.8˚)。身高在128 cm~145 cm (平均135.4 cm)。截骨矫正前行头盆环牵引。牵引70~213天(平均128.4天)后，行截骨矫正术后30天，观察角度身高及临床症状。[结果]：7例患者全部在头盆环牵引下顺利完成了软组织松解，截骨矫形融合内固定术，脊柱后凸Cobb’s角由治疗前的136˚~180˚ (平均168.8)矫正至0˚~42˚ (平均33.5˚)，矫正率81%；身高由牵引前的128 cm~145 cm (平均135.4 cm)增高至138~160 cm (平均150.7 cm)患者身高差9~27 cm。患者外观、行走姿势均得到明显改善，无脊髓神经损伤等严重并发症发生，2例合并脊髓受压(不全瘫)，患者肌力有III~级增至IV级。[结论]：头盆环牵引、软组织松解在治疗极重度脊柱畸形过程中的一种安全有效的方法之一。它可通过缓慢牵引使极重度脊柱畸形变为中、轻度脊柱畸形，这样减小矫正的手术风险，使手术矫正更为安全。 关键词 :头盆环，极重度脊柱畸形，软组织松解牵引，手术
system: [SYSTEM]

user: 1. 引言 脊柱畸形在儿童和青少年中较为常见的一种病变，由于发病年龄小，心、肺功能严重受损，甚至瘫痪，给患者健康造成严重危害。一般来说，畸形角度大于135˚，就将其归为极重度脊柱畸形 [ 1 ] 。极重度脊柱畸形患者由于畸形严重，心肺功能受损，手术耐受性差，直接矫形风险极大，临床治疗上用单纯器械内固定矫形无法取得满意的疗效。认为Cobb’s大于120˚则不宜手术 [ 2 ] 。应用头盆环的缓慢牵引、相应部位的软组织松解，使极重度脊柱畸形变为中、轻度脊柱畸形，这样避免了脊髓的损伤，降低了手术风险，使后期脊柱的矫形更为安全。我科利用头盆环牵引、软组织松解后行截骨矫形融合内固定术治疗极重度脊柱畸形患者7例，取得了较好的效果，现报道如下。
system: [SYSTEM]

user: 2. 材料与方法 2.1. 一般资料 2011年9月~2012年6月收治7例患者，男性6例，女性1例，年龄：13岁~25岁，平均21岁；其中脊柱结核后凸3例(静止期)，青少年特发性脊柱后凸4例。Cobb’s角在136˚~180˚平均(168.8˚)。身高在128~145 cm平均(135.4 cm)。术前肺功能检查均有中~重度不同程度的受损，其中有2例合并脊髓受压(不全瘫)，肌力III - 级。在患者及家属深入了解病情及治疗方案后同意治疗方案并签署同意书后开始治疗。 2.2. 方法 2.2.1. 牵引前 如图1所示，患者来院畸形严重，除常规X片正位片、侧位片(如图2所示侧弯120˚、后凸165˚)检查，了解畸形严重程度外，还要行脊髓MRI检查，确定无禁忌症(脊髓空洞、脊髓分裂等) [ 3 ] 存在； 2.2.2. 治疗方法 常规局麻下行头盆环安置术。先后取左侧卧位、右侧卧位，在局部浸润麻醉后用直径4 mm斯氏针分别从两侧髂前上棘沿髂前、髂后上棘连线穿入。再改为仰卧位置头环：局部浸润麻醉下在双侧眉弓中外1/3处、耳廓后上方1.5 cm处钻入特制螺钉直接穿皮肤拧入颅骨外板，固定头环。如图3所示，在安装头盆环手术后观察1周左右，无头痛头晕、无腹痛恶心等症状后安置盆环和支撑杆，保持头部于中立位，颈部牵引后均匀受力，根据患者病情，牵引速度3 mm~5 mm/次，当出现颈部疼痛或四肢运动感觉异常等过牵症状时停止牵引，症状仍无缓解时减低牵引高度5~10 mm。平稳后复查X片、CT、MRI，根据影像结果了解影响牵引的软组织的张力位置。在全麻下行相应部位的软组织松解手术，同时使椎体达到相对的稳定性。术后1周左右在患者身体情况平稳下再次缓慢牵引，牵引速度1 mm~3 mm/次不等(如图4所示，牵引后复查X片后凸70˚，侧弯45˚)。同时让患者每天吹气球、下蹲、爬楼等动作行肺功能锻炼；口服钙剂增加骨量，减少头钉松动可能。在再次牵引中反复出现过度牵引症状后停止牵引，适当降低高度并维持牵引一段时间，待症状完全恢复后再次复查CT、MRI，牵引70~213天平均(128.4天)后行截骨矫形融合内固定术(如图5所示，矫形后在头盆环牵引保护下复查X片畸形后凸42˚，侧弯40˚)。 2.3. 评价方法 测量手术前手术后身高及身高差；按照Cobb’s角测量方法，根据手术前、手术后X片，测量截骨部位头侧椎体上终板与尾侧椎体下终板之间的Cobb’s角度数及Cobb’s角度差。术后30天复查X片测量角度，随访6月~15月。
system: [SYSTEM]

user: 3. 结果 7例患者全部在头盆环牵引下顺利完成了软组织松解，两例患者在牵引过程中肌力恢复至IV级。后行截骨矫形融合内固定术，患者脊柱后凸Cobb’s角由治疗前的136˚~180˚ (平均168.8)矫正至0˚~42˚ (平均33.5)，矫正率81%，Cobb’s角度差116˚~180˚；身高由牵引前的128~145 cm (平均135.4 cm)增高至138~160 cm (平均150 cm)，身高增高差9~27 cm，如表1。
system: [SYSTEM]

user: 4. 讨论 头盆环最初由Perry和Nickel [ 4 ] [ 5 ] 在1959年详细报道利用牵引原理固定颈椎不稳，后研究发现利用头盆环纵向撑开力，可对中重度脊柱后凸和后侧凸的矫正有显著疗效 [ 6 ] 。我们收治的患者在发病后多年就诊，已出现严重的畸形，两例出现神经症状，下肢肌力III - 级，利用头盆环牵引的持续、慢性撑开力对抗体重、肌肉收缩以及椎体、椎间盘、韧带等软组织拉力，阻止畸形的加重，使脊柱周围软组织蠕变松解延长，挛缩的肌肉、韧带得到松解，各椎骨间的韧带、小关节松动，弯曲的脊柱逐渐伸直 [ 5 ] ，避免了对脊髓的直接牵拉、损伤。通过头盆环牵引、调节，患者身高均有不同程度的增高，使畸形逐渐得到控制、纠正，避免了直接手术导致截瘫的风险 [ 4 ] ，再结合截骨矫形、器械固定，畸形得到纠正，两例患者下肢肌力III − 级恢复至IV级。头盆环牵引是治疗严重脊柱畸形的有效手段 [ 7 ] 。从表1所示，在青少年 Table 1 Contrast before and after tractio 术前身高 (cm) 术后身高 (cm) 身高差 (cm) 术前Cobb’s角 (度) 术后Cobb’s角 (度) Cobb’s角度差 (度) 年龄 (岁) 发病至就诊时间 (年) 病例1 128 138 10 180 0 180 13 5 病例2 140 149 9 136 (侧弯135) 20 (侧弯75) 116 (侧弯60) 23 11 病例3 133 160 27 165 (侧弯120) 42 (侧弯40) 123 (侧弯80) 16 15 病例4 129 147 18 180 37 143 21 17 病例5 145 155 10 180 40 140 24 22 病例6 132 150 18 161 33 128 26 10 病例7 141 156 15 180 29 151 25 20 表1. 牵引前后信息对比表 图1. 来院外观 图2. X片：后凸165˚，侧弯120˚ 至成年的患者，畸形都得到了纠正，身高均有不同程度增长，身高差10~27 cm，如图1~6治疗前后比较所示患者身高差27 cm。牵引后软组织蠕变松解后给我们纠正其畸形成为可能。其中，成年患者牵引的身高增长相对青少年、青春期患者的身高增长小，畸形纠正角度相对较小，而在青少年和青春期的患者(2例)，畸形矫正角度相对较大，由180˚矫正至0˚，后凸165˚、侧弯120˚矫正至后凸42˚、侧弯40˚。头盆环牵引最佳时期是在青少年期还是青春期？后期矫正的畸形是否会随着生长发育增长而丢失？由于本组病例数较少，以后在治疗过程中需密切注意此类患者的随访以寻求一个最佳的牵引时期。 在安装头环、盆环时，需密切观察患者头痛、腹痛腹胀等情况。安装头环的特制螺钉主要在颅骨的 图3. 安装头盆环 图4. 牵引后X片：后凸70˚，侧弯45˚ 图5. 矫形术后头盆环保护下X片：后凸42˚，侧弯40˚ 外板内达到稳定，防止外板在牵引中出现切割及螺钉进入颅骨内板。一旦穿过颅骨内板，进入蛛网膜下腔，则会出现头痛、颅内出血甚至颅内高压；安装盆环的斯氏针，主要是从髂前上棘穿过髂骨至髂后上棘，稍有不慎，进入腹腔，则会损伤肠道。我们经验是：安装头环、盆环后禁食观察7天左右，防止腹膜炎、颅内出血等并发症的出现 [ 8 ] 。待患者无头痛、恶心、腹痛、心慌等症状，才进行头环、盆环的安装。7例患者中有1例在安装头环钉后感局部胀痛，观察无恶心呕吐，行头颅CT检查未见颅内出血，对症治疗后症状缓解，后顺利行头盆环的牵引。 图6. 矫形术后并去除头盆环外观 极重度脊柱畸形患者在病变过程中，黄韧带、后纵韧带出现钙化、骨化，脊柱出现旋转，关节突关 节的融合，以增强病变脊柱的稳定性，同时也会使脊髓受压增加，导致患者出现神经症状。在牵引过程中，椎管在后伸时容积减少，会进一步增加脊髓受压的几率。所以在牵引过程中，在纠正患者脊柱冠状面、矢状面、横断面不同受力的同时需注意患者有无颈部不适，吞咽困难、呛咳、视物模糊、下肢乏力、感觉过敏等过牵症状 [ 9 ] ，若出现症状，则达到牵引的最大限度，需停止牵引，适当降低高度并维持牵引一段时间，待症状完全恢复后行软组织松解术并保持脊柱的相对稳定。后再根据病情再次牵引或行内固定手术。本组患者中有1例出现呛咳，2例出现下肢肌力下降表现，经降低高度维持牵引一段时间后症状消失。 在头盆环牵引过程中，让患者进行下蹲、扩胸、爬楼、辅助行走等运动，改善患者肺功能，减轻畸形引起的腹部受压，增加消化功能 [ 10 ] ，增强患者体质。手术过程中发现可降低软组织张力，减小Cobb’s角，降低手术风险。 牵引过程中，常规注意钉道有无感染、患者有无过度牵引表现 [ 9 ] 外，特别需注意头环钉有无松动，一旦头环钉松动，将直接影响头盆环牵引。由于极重度脊柱畸形患者骨质较正常人明显下降，我们加强补钙，增加骨强度，增加稳定性。但有无更好的增加该类患者的骨强度，还需进一步研究。 田慧中等报道 [ 2 ] 当脊柱畸形Cobb’s角大于120˚以上者，无论用否头盆环牵引，其矫正效果和矫正率总是不满意的。我们临床通过在头盆环牵引详细观察、处理，结果表明：7例患者，后凸角度均大于135˚，为极重度畸形，在正确应用头盆环牵引，适当松解软组织，松解术后再次头盆环牵引或行内固定手术，均得到了良好的效果，身高差9~27 cm，有明显改变。4例从180˚后凸畸形矫正为0˚至40˚不等，1例从161˚后凸畸形矫正至33˚，2例后凸、侧凸畸形的患者也得到了良好的矫正，达到良好的效果。后期患者是否会随着年龄增长而出现矫正的角度丢失、身高是否会受限等还需进一步随访研究。 极重度脊柱畸形患者在矫形治疗中头盆环牵引、软组织松解是有效、安全的方法之一，它可通过缓慢牵引改变畸形的严重程度，降低患者后期矫正手术的风险，使手术矫正更为安全。虽然病例数较少，但为青少年或青春期的极重度脊柱畸形患者行头盆环牵引、软组织松解、后期矫形达到治愈提供一种有效的方法。
system: [SYSTEM]

user: 摘 要 “京津冀一体化”的提法在政策上由来已久，最早可追溯到1982年的《北京城市建设总体规划方案》，该方案明确提出要建设双重“首都圈”的设想。2014年2月26日，习近平总书记在听取京津冀协同发展工作汇报时又强调了京津冀协同发展的重要性。如今两年过去，本文在两年前原有文献的基础上，加入新增文献，梳理了有关京津冀市场一体化的研究。根据这些研究，可以发现随着时间的发展，京津冀内部以及之间的商品市场和要素市场的一体化都在不断完善。但在京津冀内部的不同地区，同一市场的一体化进程是不均匀的；而且在相同的地域范围内各行业、各要素种类的市场一体化进程也是不均匀的。 关键词 :京津冀，市场一体化，综述