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

user: 3. 讨论 卡氏肺孢子菌肺炎(Pneumocystis carinii pneumonia， PCP)，是由寄生于人体肺泡上皮细胞表面的卡氏肺囊虫大量繁殖造成上皮细胞损害，阻碍气体交换的一种机会性感染 [ 1 ]。该病是由T细胞和体液免疫反应(在较小程度上)介导的。在免疫抑制的患者中，感染会导致由单核细胞、CD8淋巴细胞和激活的巨噬细胞介导的免疫反应功能失调，从而导致弥漫性肺损伤 [ 2 ]。此外，吸入的肺孢子虫滋养体抑制肺泡内的上皮修复过程，导致严重的肺损伤 [ 3 ]。PCP感染后可以表现为轻微的炎症至重症肺炎，严重者导致呼吸衰竭甚至死亡 [ 4 ]，临床特征是干咳、发热和呼吸困难。胸部X线片上，最常见弥漫性双肺间质浸润，高分辨率计算机断层扫描可能显示出广泛的磨玻璃影及囊性病变 [ 5 ] [ 6 ]，是一种临床病死率极高的疾病。 本病多发生在免疫功能缺陷及长期应用免疫抑制剂的患者中。HIV阴性的患者PCP发病率较低，约为0.01%~0.10% [ 1 ]。但值得注意的是，随着免疫抑制剂的逐步使用，非HIV感染的免疫抑制人群的发病率也随之增加 [ 5 ] [ 7 ] [ 8 ]。在HIV阴性患者中，炎症性疾病占PCP的20% [ 9 ]。对于接受免疫抑制药物治疗的自身免疫性疾病患者，合并PCP目前还没有具体的共识指南 [ 10 ] [ 11 ] [ 12 ]。已经明确的是，炎症性疾病及其并发症和免疫抑制治疗方案增加了这类患者患PCP的风险 [ 13 ]。 肠系膜脂膜炎(Mesenteric panniculitis，MP)是一种临床上十分罕见的慢性肠系膜炎症性疾病。MP的病因可能是腹部手术、外伤、缺血、药物、过敏或自身免疫性疾病 [ 14 ] [ 15 ] [ 16 ] [ 17 ]。临床上根据MP患者的具体情况，通常选择保守治疗，如抗感染和/或免疫抑制治疗、糖皮质激素、非甾体类抗炎药物、秋水仙碱、黄体酮、环磷酰胺药物治疗等，但疗效报道不一致 [ 15 ] [ 16 ] [ 17 ] [ 18 ]。该患者炎症性疾病本身及糖皮质激素、环磷酰胺、吗替麦考酚酯等免疫抑制药物的使用，可能都是促进PCP的因素。 非HIV感染患者合并PCP的死亡率(27%~62%)高于HIV感染者(4%~15%) [ 6 ] [ 10 ] [ 19 ] - [ 24 ]。成人数据显示，与HIV感染者可能缓慢且隐匿的表现不同，非HIV感染的免疫功能受损者可能会出现快速进展性的病程，通常在几天内出现呼吸衰竭 [ 1 ] [ 5 ] [ 21 ] [ 25 ]。此外，非HIV感染患者合并PCP的临床表现也取决于所使用的免疫抑制剂，其中最常提到的有：长时间应用中到高剂量糖皮质激素、环磷酰胺、利妥昔单抗、阿林单抗和肿瘤坏死因子拮抗剂 [ 26 ] [ 27 ]。糖皮质激素可以减少CD4淋巴细胞，进而可能导致PCP更易发生。在非HIV感染免疫抑制人群中，使用糖皮质激素已成为PCP的主要致病因素，当强的松龙达到或超过16 mg时，风险尤其增加 [ 28 ]。有趣的是，在一些研究中，相当比例的PCP只在皮质类固醇逐渐减量过程中出现 [ 29 ]。本例病人在肠系膜脂膜炎治疗过程中先后应用环磷酰胺、吗替麦考酚酯等免疫抑制剂及糖皮质激素，且正是在糖皮质激素减量过程中出现PCP感染。值得注意的是，并非所有免疫抑制患者都会出现淋巴细胞减少症 [ 29 ]。与HIV感染者一样，监测CD4+ T淋巴细胞似乎不足以预测PCP发生风险 [ 30 ]。支气管镜、支气管灌洗或经支气管活检诊断率高于痰液诊断的诊断率，但不易获得 [ 31 ]。非HIV感染的PCP患者少见，缺乏特异有效且易行的诊断手段，是导致该病诊断困难乃至误诊、漏诊的重要原因。该病例CD4 + T淋巴细胞未见明显下降，支气管镜检测及肺泡灌洗液均未检见卡氏肺孢子球菌，但不能据此排除PCP感染。PCP患者临床病程凶险，病死率极高，尤其非HIV感染的PCP患者病情进展迅速，提示我们临床诊治过程中，辅助检查不十分支持的情况下仍不应忽视该病可能。另外，肠系膜脂膜炎罕见，对肠系膜脂膜炎的病因及治疗认识不足，导致未对患者免疫抑制药物及糖皮质激素用药史引起足够重视也是PCP诊断延误的原因之一。 PCP的一线治疗是甲氧苄啶–磺胺异甲噁唑静脉注射，因其具有良好的组织渗透性、快速反应和低成本等优点 [ 5 ]。不良反应包括白细胞减少、血小板减少和包括史蒂文斯–约翰逊综合征在内的皮疹 [ 32 ]。然而，这些反应在儿童中似乎不太常见 [ 33 ]。也有研究支持卡泊芬净用于PCP的治疗。目前已有TMP-SMZ治疗PCP失败换用卡泊芬净治疗成功的报道 [ 34 ]，亦有卡泊芬净联合TMP-SMZ成功治疗PCP且减少TMP-SMZ不良反应发生率的报道 [ 35 ]。PCP早期确诊治疗反应较好，多数可以得到恢复，更严重的病程往往认为继发于更具播散性的肺部炎症反应和诊断延迟 [ 5 ] [ 36 ]。所以，临床工作中，应充分强调早期诊断的重要性，提高对非HIV感染合并PCP患者，乃至肠系膜脂膜炎合并PCP患者的认识。同时，提倡对于临床上疑诊PCP的患者，在第一时间，而不必须等待病原学证据确诊，果断、及时地给予经验性治疗。
system: [SYSTEM]

user: 摘要 以轨道交通常用的6005A铝合金焊接接头为研究对象，采用200℃、230℃和250℃对6005A铝合金焊接接头进行火焰调修，研究不同调修温度对铝合金焊接接头微观组织及力学性能的影响。研究结果表明：随着调修温度的升高，6005A铝合金焊接接头的抗拉强度和断后伸长率均呈一定程度的下降。不同调修温度下，焊接接头的微观组织变化不大，焊接接头的硬度分布规律相同，相较于其他调修温度接头，250℃调修接头软化区宽度略有增加。
system: [SYSTEM]

user: 1. 引言 6005A铝合金属Al-Mg-Si系铝合金，因其具有中等强度，挤压性和耐蚀性良好等特点，非常适于制造轨道车辆车体主体结构用的复杂截面多孔中空型材，从而被广泛应用于车体的生产制造 [ 1 ] [ 2 ]。铝合金由于热导率和热膨胀系数高，焊接后容易出现焊接变形问题，影响车体结构的稳定性。在现有高速铝合金车体生产制造过程中，虽然采取了很多焊前和焊接过程控制工艺措施，但焊后仍会存在一定的焊接变形，当残余变形量超过生产技术要求时，往往需要在焊后进行机械矫正和热调修矫正 [ 3 ] [ 4 ]。热调修矫正作为一种焊接变形的矫形方法，是在焊接构件选定的区域按一定的方向进行火焰加热，使该区域的金属产生压缩塑性变形，冷却时产生收缩，进而调整焊件最终形状。这种方法操作简单，在高速列车车体生产中有广泛的应用。焊后火焰调修可以控制焊接变形，同时也会引起材料组织性能和残余应力的变化，热调修不当甚至会造成焊接构件的使用性能降低 [ 3 ] [ 4 ] [ 5 ] [ 6 ]。高速列车用6005A铝合金属于可热处理强化铝合金，热调修方法必须严格控制额外的热输入，避免产生焊接接头严重软化，影响车体的服役安全 [ 5 ] [ 6 ] [ 7 ]。不同温度下火焰调修的铝合金焊接接头组织及性能如何变化是亟待解决的问题，因此研究调修温度对6005A铝合金焊接接头组织与性能的影响，对指导铝合金车体生产制造具有重要的工程实际意义。
system: [SYSTEM]

user: 2. 材料及试验方法 试验材料采用6005A-T6铝合金中空型材(见图1)，板厚3 mm，单边自带有35˚ V型坡口，型材挤压成型时已有1 mm钝边，组对间隙为0~1 mm，且自带永久垫板。焊接方法为MIG焊，所用焊丝牌号为ER5356，焊丝直径为1.2 mm，焊接保护气为高纯氩。焊接工艺参数为：焊接电压21~22 V，焊接电流220~230 A，焊接速度800 mm/min，单道焊接成型。6005A-T6铝合金和焊丝的化学成分见表1。焊后采用氧乙炔火焰加热焊接接头区域，加热温度分别为200℃、230℃和250℃，为了减小热量对工件的影响及降低热量在工件上的传递范围，冷却方式采用水冷降温。热调修后，利用线切割机对焊接接头进行金相、拉伸、弯曲、硬度试样的加工取样。采用型号为AX10 ZEISS的蔡司光学显微镜观察焊接接头金相组织。拉伸试验参照GBT 2651-2008焊接接头拉伸试验方法和GB/T 228-2010金属材料室温拉伸试验方法制成矩形横截面拉伸试样。采用DNS300微机控制电子万能试验拉伸机进行焊接接头的拉伸试验和弯曲试验，拉伸速度为3 mm/min，按标准GB/T 2653-2008焊接接头弯曲试验方法根据试板的厚度选取试验跨距和压头直径，弯曲角度为180˚。在距离焊接接头上表面1 mm处布置测试线，用HVS-30D维氏硬度计测量焊接接头各区域的硬度值，测量载荷为1 kg，载荷持续时间为15秒，硬度点之间的间隔距离为1 mm，测试方向为从焊缝中心向母材方向。 图1. 6005A-T6挤压中空型材结构示意图 Table 1 Chemical composition of 6005A aluminum alloy and ER5356 welding wire (wt% 牌号 Si Mg Fe Cu Mn Cr Zn Al 6005A-T6 0.6-0.9 0.4-0.6 <0.35 <0.30 <0.10 <0.10 <0.10 余量 ER5356 ≤0.25 ≤0.10 ≤0.10 0.05~0.2 4.5~5.5 0.05~0.2 ≤0.10 余量 表1. 6005A铝合金与ER5356焊丝的化学成分(wt%)
system: [SYSTEM]

user: 3. 试验结果及分析 3.1. 调修温度对接头组织的影响 图2和图3分别为不同调修温度下焊接接头焊缝及熔合区的微观金相组织。从图2可以看出，未调修和不同调修温度下的接头焊缝均为树枝状铸态结晶组织，在光学显微镜下各温度加热调修后的焊缝组织与未调修焊缝组织相比没有出现明显差别。这主要是因为在焊接过程中，焊缝中心经历了极高的热循环，热调修温度在没有超过焊缝中心之前的温度情况下，不会对焊缝金相组织造成明显的影响。从图3可以看出，各调修温度下熔合区靠近焊缝侧的组织均呈铸态枝晶状，靠近HAZ侧晶粒较大，但从所观察到的区域来看，与未调修接头相比，3种调修温度下靠HAZ侧的晶粒尺寸未发生明显的长大。 图2. 不同调修温度条件下焊缝金相组织 图3. 不同调修温度条件下熔合区金相组织 3.2. 调修温度对接头拉伸性能的影响 表2为未调修和不同调修温度条件下的6005A-T6铝合金焊接接头的拉伸试验结果(表中抗拉强度平均值和断后伸长率平均值均为3个平行样的平均值)，图4为未调修和不同调修温度条件下接头拉伸断口的扫描电镜(SEM)形貌。表2可以看出，经过不同温度调修的接头抗拉强度和断后伸长率均低于未调修接头，且随着调修温度的升高，接头抗拉强度和断后伸长率程下降趋势，但仍然高于标准155 MPa的要求。所有的焊接接头均断于热影响区(HAZ)，距离焊缝中心约10 mm处，虽然焊缝的熔敷金属强度略低于母材，但是MIG焊和焊后热调修均会使热影响区产生多次热循环，导致双重的软化现象致使焊接接头其强度降低较大，故而热影响区成为整个6005A-T6铝合金焊接接头最为薄弱的环节。图4可以发现，所观察的4种条件下的接头断口均呈现为典型的韧窝断口，未调修、200℃和230℃调修温度下的接头断口形貌并无特别明显的区别，250℃调修温度的接头断口韧窝直径要小于未调修接头，一般韧窝尺寸越大，说明在形成韧窝的过程中发生的塑性变形更严重，变形更充分，那么断裂过程中会吸收更多的能量，材料的韧性就越好。 Table 2 Tensile test result 调修温度 抗拉强度平均值/MPa 断后伸长率平均值/% 断裂部位 未调修 203 7.01 HAZ 200℃ 195 6.52 HAZ 230℃ 190 6.23 HAZ 250℃ 181 5.98 HAZ 表2. 拉伸试验结果 图4. 不同调修温度条件下拉伸断口SEM形貌 3.3. 调修温度对接头硬度的影响 图5为未调修和不同调修温度接头的硬度分布曲线。从图中可以发现，4种接头的硬度曲线分布规律基本相同，即母材区硬度 > 热影响区淬火区硬度 > 焊缝区硬度 > 过时效软化区硬度。随着远离焊缝中心，接头的硬度由焊缝区过渡到热影响区淬火区硬度值升高，由热影响区淬火区过渡到过时效软化区硬度值降低，再由过时效软化区过渡到母材区硬度值升高。其中过时效软化区的硬度最低值约为60 HV，低于焊缝区的硬度值，成为整个焊接接头最薄弱的区域 [ 8 ]，接头拉伸试验的断裂部位正好也位于该区域，从而验证了这一点。这是由于6005A-T6铝合金为可热处理强化铝合金，在焊接过程中会产生软化现象，主要原因是热影响区在焊接高温作用下强化相会脱溶析出并聚集长大，使强化效果减弱，产生“过时效”现象，并形成软化区 [ 9 ]。因为随着距焊缝中心的距离增加，距离焊缝中心较远的热影响区的温度会逐渐降低，热循环温度虽然比原有时效处理温度高，但又未达到固溶温度，强化相高温溶解不充分且析出数量少，并且原有未溶解的过渡强化相β′(Mg 2 Si)也因为焊接热循环的影响而聚集长大，转变为稳定的β(Mg 2 Si)强化相，由于析出的和长大的强化相粒子粒径不均匀，强化效果弱，从而形成过时效软化区，就使6005A-T6铝合金焊接接头的强度和硬度明显降低 [ 10 ]，所以该区域的显微硬度最低，拉伸时也因强度降低而最先断裂。从图5中还可以发现250℃调修接头的软化区宽度要略大于未调修和其他调修温度接头的软化区宽度，说明调修温度越高，受到热循环的影响，6005A-T6铝合金接头的软化区会变宽。这是因为调修温度越高，高温停留时间越长，强化相越容易脱溶析出并聚集长大，同时受到调修温度影响的区域越宽，使得软化区变宽。 图5. 不同调修温度接头硬度分布曲线 3.4. 调修温度对接头弯曲性能的影响 通过对未调修和不同调修温度的焊接接头进行面弯和背弯试验，发现所有焊缝均能达到弯曲试验标准，未发现明显的弯曲裂纹，说明接头的弯曲性能良好。这说明在本文中的调修温度下，热调修不会影响焊接接头的弯曲性能。
system: [SYSTEM]

user: 4. 结论 通过采用200℃、230℃和250℃对6005A铝合金焊接接头进行火焰调修，并与未调修接头的组织和基本力学性能进行对比，得出以下结论： 1) 随着调修温度的升高，与未调修接头相比，焊缝和熔合区的微观金相组织未发生大的变化。 2) 经过不同温度调修的接头抗拉强度和断后伸长率均低于未调修接头，且随着调修温度的升高，接头抗拉强度和断后伸长率呈一定程度的下降趋势，但仍然高于标准155 MPa的要求。 3) 在所采用的调修温度中，硬度曲线分布规律基本相同。但250℃调修接头的软化区要略大于未调修和其他调修温度接头的软化区宽度。 4) 从实验的结果来看，200℃、230℃和250℃调修温度均可用于生产实际的变形调修中，不过调修温度采用200℃最佳。
system: [SYSTEM]

user: 摘要 课程思政是课程建设的一个趋势，本文针对数控设备应用与维护专业的特点，以《数控机床电气控制与PLC》为例，探索思政元素的挖掘，进行有效的单元设计，通过课程思政的实践使学生在接受教育的同时树立正确的专业科学观、价值观，达到全方位育人的目的。
system: [SYSTEM]

user: 1. 前言 2016年12月习近平总书记在出席全国高校思想政治工作会议并发表重要讲话，指出“做好高校思想政治工作，要用好课堂教学这个主渠道思想政治理论课，要坚持在改进中加强，提升思想政治教育亲和力和针对性，满足学生成长发展需求和期待，其他各门课都要守好一段渠，种好责任田使各类课程与思想政治理论课同向同行，形成协同效应” [ 1 ]。这一纲领性的文件对专业课任课教师、高校课程建设提出了更高的要求，对于做好新时代高校思政工作具有重要的意义，为高校践行办好中国特色社会主义大学的专业课程建设指明了方向，是实现“立德树人”这一根本目标的新思路。要推进高校课程思政建设工作，工科专业课程是最核心和最难的一部分：首先，工科专业课程在高校课程具有重要的地位，具有为国家工业新时代发展培育接班人的历史使命，这是工科专业课程思政实践的必要性；其次，目前工科专业课教师思政教育意识、能力匮乏，生搬硬套出现“思政”与“课程”“两张皮”的现象，受教育者在受教育过程中无法形成核心科学的价值观，甚至出现反作用，损害原本的知识体系及思政的育人功能，这是工科专业课程思政实践的紧迫性；最后，工科专业课程中蕴含了大量的思政元素，思政载体丰富，实践空间大，平台、形式、方法灵活多元，已经具备了实践课程思政的条件，这是工科专业课程思政实践的可行性 [ 2 ] [ 3 ] [ 4 ]。因此，工科专业课程思政实践研究是新时代专业建设的重要内容，对全面推进高校课程思政工作，形成思政教育合力具有重要的意义。
system: [SYSTEM]

user: 2. 数控机床电气控制与PLC课程的思政意义 数控机床电气控制与PLC课程是数控设备应用与维护专业非常重要的一门专业课，涵盖了基础的机床电气原理、连接调试、PLC等内容，是是高职工科专业中具有实践意义且较为典型的专业课程，对知识及技能相结合的要求较高，知识点较多，技能水平要求高，也是企业实际生产需求较为迫切的一项专业技能 [ 5 ] [ 6 ]。因此在教学过程中需要引导学生对课程具有高度的认同感，需要学习课程时要有严谨的态度、精益求精的追求、与其他同学团结协作的意识以及高度的专业认同感，才能较好地完成课程的学习，养成良好的职业素养。课程思政可以充分发挥全员、全课、全过程的协同效应，形成价值体系，内容体系，评价体系，将思想政治教育的显性育人功能和其它类课程的隐性育人功能充分有机渗透与融合，真正实现知识传授、能力培养和价值引领相统一，着力培养学生精益求精的大国工匠精神，激发学生科技报国的家国情怀和使命担当。
system: [SYSTEM]

user: 3. 整体教学设计的创新 课程教学设计包括制定教学目标、教学内容的选取、教学环节的设计、教学方法的创新、教学环境的创建等等 [ 7 ]。通过“课程思政”深入挖掘各门学科专业课程的价值观教育资源，以课程知识体系为本体，价值体系为指引，以点、线、面、拓展的形式多方面挖掘资源，并形成全面的内容资源体系，思路如图1所示。将思政元素融入课堂就是需要从企业需求、国家职业标准出发，制定与学生匹配相应的教学目标，其次是教学内容的选取，选取与企业实际岗位需求、科技最新发展息息相关的内容，同时系统地将课程内容串起来，完整地实现课程的教学知识及技能、素养目标。针对某一教学内容的课堂教学环节的设计，需考虑如何将思政元素很好地融入其中 [ 8 ] [ 9 ]。教学方法上灵活运用，以学生为中心充分运用学生讨论、探索、实践、反思等方式充分调动其学习的积极性。教学环境将实训室、多媒体课堂、企业实训基地三地以及线上、线下二维空间来进行教学。 图1. 课程思政内容体系的构建
system: [SYSTEM]

user: 4. 课程思政实践 专业认真组织“课程思政建设”专题教研活动，通过学习各项课程思政建设工作会议核心精神及指导意见，统一认识思想政治教育在工科专业课程融入的必要性，统一思政课堂革命的正确方向。其次，邀请学校思政课教师共同研讨思政元素的发掘及融入课程的方法，重点关注工科专业课程的特色内容“工程伦理教育、精益求精的大国工匠精神、培养和践行社会主义核心价值观、科技报国的家国情怀和使命担当等民族精神”，明确课程的思政目标，根据课程的知识体系充分挖掘思政元素，并以渗透的方式融入教学环节中 [ 10 ]。 工科专业课程思政的教学实践 对照课程思政方法，通过多方主体研讨，确定课程融入思政目标的整体教学目标，充分挖掘课程内容中的具体思政元素，从知识、技能、价值为一体的目标出发，重构教学环节、创新教学方法、多元化教学手段、拓展教学空间，以有效的融合为评价标准进行教学设计的创新，形成课程思政特色教学整体设计、单元设计、拓展模块等材料。以《数控机床电气控制与PLC》课程的第一章为例，课程思政的发掘及切入点如表1所示。 Table 1 The integral design of the first chapte 序号 项目/知识点 “课程思政”有机渗透切入点 1 我国数控机床的发展 1) 我国数控机床的发展史，从进口到自己生产，从低端设备到高端设备的研发，是国家工业水平的提升，也是综合国力的提升，逐渐向发达国家的技术水平看齐，激发学生名族自豪感，认可自身的专业。 2) 我国在机床行业依然存在很多技术瓶颈，需要很多代人的热情及努力去突破，明确发展中的差距，培养学生的忧患意识及使命感。 3) 机床行业的发展变化对数控人才的需求也产生变化，从机床操作到机床维护、机床研发的人才转变，让学生对这一课程及自身的发展有比较明确的方向。 2 数控机床XK630电气柜的认识与制图 1) 通过电气柜内规范的布局及连接，强化学生的职业素养，培养耐心细致的工作作风和严肃认真的工作态度。 2) 制图单规范化回执及签字强化学生的职业道德意识，恪守信用、尊重规则的职业道德与修养。 3 低压元器件的认识(一) ——保护元器件 1) 通过教学中对于低压元器件的认识、连接、选型的要求，契合数控装调高级工的电气基础要求，培养学生自我追求的意识，及严谨的工作态度；以国家职业标准来指导学习。 2) 保护元器件的选型及维护不当会造成严重的电路安全事故，引导学生在学习时务必以求真的态度，规范操作，养成良好的职业素养。 4 低压元器件的认识(二)——控制元器件 控制元器件是实现电路控制要求的关键元器件，且种类繁多，规格、型号复杂，在这一部分，“授人以鱼不如授人以渔”，以问题引导的方式，使学生掌握通过查阅资料、实践等方式来寻求解决方案，培养学生独立解决问题的能力，以及务实求真的态度。 5 电气连接的工具简介及使用 电气连接工艺是通过工具将元器件按要求进行连接，通过实际大量的训练，观摩企业实际操作视频，来磨炼自己的工艺技能，在每一个实践项目都进行严格的工艺检测，培养学生孜孜不倦、精益求精的工匠精神。 6 XK630低压元器件的 选型与连接 通过一章的学习，学生能独立进行元器件的选型及掌握连接工艺技能，建立学生的专业自信，实践创新的工匠精神。 表1. 数控机床电气结构的认识 课堂教学单元设计，以《数控机床电气控制与plc》课程为例，选取电气连接的工具简介及使用为例，将“工匠精神”中的精益求精、严谨、独立的态度充分地融入到教学单元设计中，其中充分地使用了线上讨论、讲授实践、岗位模拟、总结等教学方法，同时将教学环境拓展到校外的企业实际工位上，实现“工匠精神”的内化于心的教学过程。单元设计如表2所示。 Table 2 Tool introduction and use of electrical connection 项目名称 电气连接的工具简介及使用 课时 1 教学目标 知识目标：1.掌握电气连接所需的各种工具的作用；2.掌握电气连接的标准及方法 技能目标：1.能熟练使用各种工具；2.能根据要求完成连接工艺 思政目标：1.精益求精的工匠精神；2.严谨、独立的学科文化 教学环节的设计 内化于心 线上导入 观看企业电气连接车间的规范化操作视频 认知与目标的树立 讨论(引导) 结合电气装调职业标准总结电气连接的标准 1) 根据视频总结电气连接标准 2) 根据电气装调职业标准要求进行补充完整 严谨的科学态度 讲授及实践 教师工具使用演示及作用讲解， 指导学生独立操作 1) 每位学生独立操作 2) 独立完成知识工单 独立学习、自信 岗位模拟 1.制作U形片连接头 1) 对于工艺不合格的必须重新连接； 2) 小组比赛制体现追求更好的工艺作品； 3) 对于能力较弱的同学，将工匠精神的要求延伸至“第二课堂” 孜孜不倦、精益求精的工匠精神 2.小组赛，制作6个U形片连接头，并将保护元器件进行连接 3.教师进行工艺检查及小组互评，反馈工艺问题及解决方法 4.重新审视，根据反馈意见重新制作并连接 5.对于工艺未达标的学生通过课后时间进一步练习，提交标准的工艺连接作品 总结 电气连接的工具名称及作用； 电气连接的标准及方法 引导学生多方面总结，做出全面的回答； 对于有疑问处进行纠正及补充 严谨的科学态度 校外延伸 参观高端设备电气连接车间、 建立工艺精湛的电气连接作品展示台等 自我追求及专业认可 岗位实践 在实训基地进行岗位实践 践行价值观、工匠精神 表2. 电气连接的工具简介及使用
system: [SYSTEM]

user: 5. 总结 以《数控机床电气控制与plc》课程为例的课程思政教学实践实现了教学目标的重设，以课程内容为基础充分地发掘思政元素，并以课堂全部环节为载体融入思政元素，创新多主体参与、教学手段的多样化及教学空间的多方拓展，对于工科人才培养具有一定的参考意义。
system: [SYSTEM]

user: 摘要 基于服从负偏态分布数据的一种水平镜像算法，本文提出一种改进Box-Cox变换——镜像Box-Cox变换，并进行数值实验，实验结果显示，与传统的Box-Cox变换相比较，镜像Box-Cox变换在处理正偏态分布数据的效果上与传统Box-Cox变换处理效果相同的基础上，其处理负偏态分布数据的效果要优于传统Box-Cox变换。再进行模拟回归模型实验，实验结果表明，经过镜像Box-Cox变换的数据建立的回归模型的拟合和预测效果有所提高，且效果优于使用传统Box-Cox变换后的数据。
system: [SYSTEM]

user: 1. 引言 1.1. 研究背景 现实中，我们会遇到的数据纷繁复杂，不同的数据根据我们所做的假设的不同，需要进行不同的变换，以便我们能够在已有理论上对其进行分析。例如：股票收益率等数据的特殊性，不可观测的误差可能是和预测变量相关的，但其不服从正态分布，于是给线性回归的最小二乘估计系数的结果带来误差，为了满足线性回归的四个假设条件而又不丢失信息，有时需要对数据进行处理变换；又例如方差分析需要试验误差具有独立性、无偏性、方差齐性和正态性的条件，若不满足这些条件就需要对数据进行处理 [ 1 ]。 Box-Cox变换是George Box和David Cox在1964年提出的一种参数化广义幂变换方法 [ 2 ]，其主要特点是引入一个参数 λ ，通过数据本身估计该参数 λ ，从而确定应采取数据变换形式 [ 3 ]。常用于稳定方差、减少数据在统计建模中的非正态性和增强关联性度量的有效性。 基于正偏和负偏的数据互为镜像关系，本文提出的镜像Box-Cox变换，提高偏态分布数据正态化效果，且便于运算，并进行数值实验以验证该结论。 1.2. 正态性检验及回归模型评价指标说明 1) Shapiro-Wilk检验 [ 4 ] (W检验) W检验是用来检验数据是否符合正态分布的。可计算得到一个相关系数，它越接近1就越表明数据和正态分布拟合得越好。且W检验还会给出一个P值，若P值大于0.05，就无法拒绝其符合正态分布。若统计量W值接近1，但P值小于0.05，我们仍然拒绝其符合正态分布。W检验计算公式为： W = ( ∑ i = 1 n a i y i ) 2 ∑ i = 1 n ( y i − y ¯ ) 2 其中， y i 为样本的次序统计量， a i 为一个待估常量。 2) MAPE [ 5 ] (Mean Absolute Percentage Error，平均绝对百分比误差) MAPE常用于描述准确度，它是一个百分比值，因此比其他统计量更容易理解。MAPE的值越小，说明预测模型拥有更好的精确度。其数学表达式为： MAPE = 1 n ∑ i = 1 n | y i − y ⌢ i y i | × 100 %
system: [SYSTEM]

user: 2. Box-Cox变换 假设样本里一共有n个数据点，分别是： y = ( y 1 ， y 2 ， ⋯ ， y n ) ′ ，我们把变换后新的数据点记为： y ( λ ) = ( y 1 ( λ ) ， y 2 ( λ ) ， ⋯ ， y n ( λ ) ) ′ 当 y ≥ 0 时，Box-Cox变换是对原始数据做如下变换： y ( λ ) = { y λ − 1 λ ，       λ ≠ 0 log y ，           λ = 0 (1) 当存在 y i < 0 ，   i = 1 ， 2 ， ⋯ ， n 时，Box-Cox变换是对原始数据进行如下变换： y ( λ ) = { ( y + β ) λ − 1 λ ，           λ ≠ 0 log ( y + β ) ，           λ = 0 (2) λ 是一个待定变换参数，对不同的 λ ，所做的变换自然就不同，所以这是一个变换族。我们将(1)式称为Box-Cox变换的基本公式；将(2)式称为Box-Cox变换的扩展公式。
system: [SYSTEM]

user: 3. 镜像Box-Cox变换 众所周知，在面对偏态分布数据时，我们需要使其变换为正态分布，要在保持数据大小次序不变的同时，减小数据之间的距离。对于处理正偏态分布的数据就需要一类增长率单调递减的函数，例如对数函数、平方根函数等等。而根据Box-Cox变换公式： y ( λ ) = { y λ − 1 λ ，         λ ≠ 0 log y ，             λ = 0 要使变换函数增长率单调递减，就需要变换的二阶导数小于0，即： ( y ( λ ) ) ′ ​ ′ < 0 ⇒ { λ ( λ − 1 ) y λ < 0 ，       λ ≠ 0 − 1 y 2 < 0 ，                         λ = 0 ⇒ λ ∈ [ 0 ， 1 ) 故Box-Cox变换在处理正偏态分布数据时，待定变换参数 λ ∈ [ 0 ， 1 ) ，同理可得，其在处理负偏态分布数据时，待定变换参数 λ ∈ [ 1 ， + ∞ ) 。 所以，Box-Cox变换在处理正偏态分布数据时，相较于其处理负偏态分布数据时，除了幂函数，还多了对数变换，且变换参数的取值范围较小，易于求解。 理论介绍 ∀ y i ∈ ℝ ，镜像Box-Cox变换是对原始数据做如下变换： y ( λ ) = { α ⋅ ( α ⋅ y + β ) λ − 1 λ ，           λ ≠ 0 α ⋅ log ( α ⋅ y + β ) ，           λ = 0 (3) 若原始数据服从正偏态分布，且 ∀ y i ≥ 0 ，则： α = 1 ， β = 0 ；若原始数据服从正偏态分布，且 ∃ y i < 0 ，则： α = 1 ， β = max ( | y | ) + 1 ；若原始数据服从负偏态分布，则： α = − 1 ， β = max ( | y | ) + 1 其中， max ( | y | ) 表示原始数据取绝对值后，数据中的最大数。 无论选择传统的Box-Cox变换还是镜像Box-Cox变换，最关键的问题在于怎样选定一个最优的 λ ，使得变换后的样本(及总体)正态性最好 [ 6 ]。求解最优参数 λ ，我们可以采用极大似然估计法和Bayes方法 [ 6 ] [ 7 ]。
system: [SYSTEM]

user: 4. 数值实验 我们将分三种情况进行数值实验：1) 数据全为正数，2) 数据全为负数，3) 数据有一部分正数和一部分负数。 本节使用Python软件随机生成上述三种情况的负偏态分布 [ 8 ] 数据，之后将此类数据分别使用传统的Box-Cox变换和本文提出的镜像Box-Cox变换进行处理，对变换后的数据进行偏度、峰度 [ 9 ] 和Shapiro-Wilk [ 10 ] 检验，并画出数据的频率直方图和P-P图，据此可以比较两种变换的效果。 4.1. 符号说明 本节将实验中处理的不同数据使用不同的符号表示，便于之后的实验结果的描述，符号说明如表1所示。 Table 1 Symbol descriptio 符号 说明 І − 型数据 表示全为正数的服从负偏态分布的原始模拟数据 П − 型数据 表示全为负数的服从负偏态分布的原始模拟数据 Ш − 型数据 表示部分负数、部分正数的服从负偏态分布的原始模拟数据 Mx数据 表示x型数据经过镜像Box-Cox变换后的数据，x取І − 、П − 、Ш − Tx数据 表示x型数据经过传统Box-Cox变换后的数据，x取І − 、П − 、Ш − 表1. 符号说明 4.2. 实验结果 4.2.1. 数据正态性检验图示结果 图1从左到右分别表示的是І − 型、П − 型、Ш − 型、TІ − 型、TП − 型、TШ − 型、MІ − 型、MП − 型和MШ − 型数据的直方图和P-P图。从图中可以看出经过镜像Box-Cox变换后的数据更加接近正态分布。 4.2.2. 数据正态性假设检验结果 从表2中可以看出无论是І − 型、П − 型还是Ш − 型数据，经过镜像Box-Cox变换后的数据在偏度、峰度和W值的表现均比传统Box-Cox变换后的数据更加接近正态分布。
system: [SYSTEM]

user: 5. 回归模型模拟 5.1. 实验分析 生成一组服从负偏态分布的回归因变量，数据中即包含正数、也包含负数，计算原始数据的偏度、峰度，并进行Shapiro-Wilk检验，之后对原始数据建立回归模型，进行预测并计算MAPE；再使用传统Box-Cox变换及镜像Box-Cox变换对原始数据进行处理，计算经过变换后的数据的偏度、峰度，并进行Shapiro-Wilk检验，之后使用处理后的数据建立回归模型，进行预测并计算MAPE。 图1. 三类数据的频率直方图和P-P图 Table 2 Data normality test results and optimal parameters ( λ 统计量 数据类型 偏度(skewness) 峰度(kurtosis) W值 P值 最优参数( λ ) І − 型数据 −0.9919 0.8031 0.9327 3.213e−14 TІ − 型数据 −0.2097 −0.7730 0.9805 3.0208e−06 2.5746 MІ − 型数据 0.0256 −0.2940 0.9974 0.6136 0.3642 П − 型数据 −1.1063 1.5032 0.9277 8.4380e−15 TП − 型数据 −0.2745 −0.7039 0.9809 3.8808e−06 2.3344 MП − 型数据 −0.0098 0.0173 0.9967 0.4014 0.3023 Ш − 型数据 −1.0621 1.7657 0.9346 5.4640e−14 TШ − 型数据 −0.2344 −0.6995 0.9827 1.1579e−05 2.5655 MШ − 型数据 0.0170 −0.1521 0.9977 0.7196 0.3909 表2. 数据正态性检验结果及最优参数( λ ) 5.2. 实验结果 从表3可以得到经过镜像Box-Cox变换后的数据更接近正态分布，并且其预测值与实际值的平均绝对百分比误差(RMSE)为2.802%，远小于传统Box-Cox变换后的数据的70.77%。 Table 3 Data normality test results and model fitting effect evaluatio 统计量 数据类型 偏度 (skewness) 峰度(kurtosis) W值 P值 RMSE 原始数据 −1.1638 2.4959 0.9379 5.342e−20 11.15% 传统Box-Cox变换数据 −0.0543 −0.1824 0.9983 0.4239 70.77% 镜像Box-Cox变换数据 0.0062 0.1481 0.9987 0.7126 2.802% 表3. 数据正态性检验结果及模型拟合效果评价
system: [SYSTEM]

user: 6. 结语 通过数值实验，我们发现无论何种类型(全为正数、全为负数和部分正数、部分负数)的负偏态分布数据经过镜像Box-Cox变换后，数据基本服从正态分布，且效果要优于使用传统Box-Cox变换。 进行模拟回归模型实验结果表明，使用传统Box-Cox变换后的数据建立的回归模型，进行预测后其RMSE为70.77%；使用镜像Box-Cox变换后的数据建立的回归模型，进行预测后其RMSE为2.802%，经过镜像Box-Cox变换的数据建立的回归模型的拟合和预测效果有所提高，且效果优于使用传统Box-Cox变换后的数据。 本文提出的镜像Box-Cox变换通过计算机易于实现，效果相较传统Box-Cox变换有所提高，可以作为处理非正态数据的一种可靠的方法。
system: [SYSTEM]

user: 摘要 该文提出一种基于stm32为主芯片的智能牛体刷控制器。该控制器采用stm32为主控制芯片和mpu6050为角度检测芯片，两者通过IIC总线进行通讯，采用可控硅方式控制单相交流异步电机。该控制器具备控制电机正反转功能，过电压保护功能，过电流保护功能。通过实际的应用反馈，该控制器可有效应用在牛场等环境中。
system: [SYSTEM]

user: 1. 引言 在奶牛养殖过程中，产奶量是其重要指标，直接影响养殖效益。为奶牛提供舒适安全的生活环境，是提高产奶量的重要手段。一般具体措施有适当提高生活温度、提高饲料和饮水标准、保持环境通风换气、刷拭牛体 [ 1 ]、保持环境卫生等。其中，擦拭牛体方式，即保持牛体干净清洁的手段，一般采用人工方式擦拭或者牛体自动设备擦拭 [ 2 ]。目前市场存在的智能自动刷设备，普遍存在控制不灵敏，保护不到位，功能不完善的问题。据此，本文提出一种基于stm32为主芯片的智能牛体刷控制器，具备正反转自动切换，过压欠压保护，过流保护，角度检测灵敏等实用功能。刷体机械部分和牛体刷控制器刚性连接固定。当牛体接触刷体，带动控制器离开基准角度一定幅度，控制器自动检测触发工作，控制电机转动带动刷体旋转。转动过程中，当发生牛体毛发过长卷入刷毛，导致电机堵转，控制检测回路电流超负荷，控制电机反向旋转，保证牛体安全。当错误持续发生时，具备及时切断控制回路的功能，安全有效。
system: [SYSTEM]

user: 2. 控制器总体设计 该牛体刷控制器包含软、硬件两部分。硬件部分包含电源转换单元、角度检测单元、电流监测单元、电压监测单元、控制输出单元和控制核心单元等部分。软件部分，主要实现通过对电压、电流、以及当前角度信息的采集，实现电机的控制逻辑。整个控制器的工作流程为，当控制器检测到本身离开基准偏移一定角度后，自动触发控制电机输出控制回路开关打开，使电机转动。转动过程，通过AD转换电路，实时采集电机的电流和电压信号，按照设定好的阈值参数，进行停机报警等操作。 该控制器的实物图，如图1所示。 图1. 牛体刷控制器实物图
system: [SYSTEM]

user: 3. 控制器的软硬件设计方案 3.1. 控制器的硬件设计 该控制器的硬件布局，如图2所示。 图2. 牛体刷控制器硬件布局图 电源转换单元：负责进行电源转换，见图3，为控制器其他控制器电路提供稳定的低压直流12 V、5 V、3.3 V。同时需要对电机端提供标准市电交流220 V。整个电源回路，首先对市电220 V交流电，通过EI30变压器转换成12 V交流，然后通过LLD06整流电路输出12 V直流，进而通过L78L05和SPX1117M3-3.3器件，输出稳定的5 V和3.3 V电压。 图3. 电源转换单元原理图 角度检测单元(见图4)：负责实时输出稳定的当前控制器板的姿态。该控制器采用MPU6050，全球首例9轴运动处理器，内部集成3轴MEMS陀螺仪，3轴MEMS加速计和以及1个可扩展的数字运动处理器DMP。考虑到产品的应用场景和精度问题，利用内部DMP [ 3 ] [ 4 ] [ 5 ] [ 6 ]，DMP (Digital Motion Processor)，数字运动处理器，通过融合陀螺仪和加速度计的数据，直接输出四元数，通过解算得到当前姿态欧拉角，减轻主芯片处理器的计算量。 图4. 角度检测单元原理图 电流监测单元、电压检测单元：电流监测电路，通过电流互感器器件，利用AD采样电路实时获取，电机运行时的电流信息。电压检测电路，通过对变压器二次侧的转换电压，利用AD采样电路，实时获取当前市电运行的电压信息。根据预先设定的电压、电流参数，判断当前是否正常运行状态。 控制输出单元：作为电机的驱动控制回路。该控制器采用可控硅的控制方式，控制电机的转动和停止。可控硅(Silicon Controlled Rectifier)简称SCR，是一种大功率电器元件，也称晶闸管。具有体积小、效率高、寿命长等优点。在自动控制系统中，可作为大功率驱动器件，实现用小功率控件控制大功率设备。它在交直流电机调速系统、调功系统及随动系统中得到了广泛的应用。 控制核心单元：控制核心主要负责对各路IO输入信号的计算处理，以及对控制输出功能。同时，具备参数存储的功能。核心控制器，采用高性能的ARM Cortex-M3核的32位处理器STM32F103C8T6，工作频率可高达72 MHz，内置高速存储器(20 KB SRAM、64 KBFLASH)，工作电压范围在2~3.6 V。 3.2. 控制器的软件设计 该控制器的软件实现，包括：初始化、读取当前角度信息、判断当前工作状态、控制电机开关等部分。 如下表1所示，软件预先的内部设定参数值。 如下图5所示，控制器软件运行流程图。 初始化包含控制器时钟初始化，参数初始化，通讯总线初始化、外设MPU6050的参数初始化，以及LED状态灯和数码管显示的初始化等。设置时钟单元，使能内部时钟单元，使主频配置在72 Mhz。然后，读取预先设定的电压、电流防护参数值，单次触发控制电机转动的周期参数值和触发控制器工作的偏移角度阈值参数。然后，初始化IIC总线与MPU6050的交互流程，对MPU6050的DMP数字运动处理器进行相应的初始化，包含设置DMP需要的陀螺仪和加速度计数据来源，设置数据输出的FIFO (First Input First Output)存储器，设置数据采样频率，设置陀螺仪的方向，加载DMP固件信息，设置数据输出频率，最后使能DMP计算。 Table 1 Pre-set internal parameter 参数名称 参数值 电压下限值 187 V 电压上限值 253 V 电流上限值 10 A 单次控制电机转动时间 90 S 最大偏移角度 15 表1. 预先设定的内部参数 图5. 控制器软件流程图 循环工作流程，通过IIC总线实时获取DMP当前的计算得出的四元数，进而解算得到欧拉角信息。主循环逻辑，判断当前的欧拉角与初始化水平姿态比较，当超过最大偏移角度阈值参数后，控制对应的管脚输出，控制对应正向可控硅导通，进而控制电机转动。同时，通过AD采样实时采集电机运行电流和当前电压值。当检测到电机电流超过对应参数时，认为可能存在“堵转”现象，主循环逻辑控制可控硅关闭，保护电机。等待计时2秒，控制反向可控硅导通，循环往复。当持续5次检测到电流超限时，控制器计时4分钟停止响应任何外部触发，进入到休眠模式。当检测到电压超过上下限位时，控制电机停止，进入到休眠模式，等待电压恢复正常后，恢复到正常工作模式。
system: [SYSTEM]

user: 4. 应用注意事项 根据原理图实物搭建如图6所示。 图6. 控制器实物连接图 应用中应着重安全问题，一是使用者安全的问题、因应用市电220 V，使用人员应着重注意安全问题；二是线路安全问题，所有电机的等线路的操作，必须在断电时操作。防止电路短路等意外发生。 根据原理搭建测试，系统运行符合预期。
system: [SYSTEM]

user: 5. 结束语 本文介绍了牛体刷控制器的软硬件构成，以及实际应用应当注重的安全问题。该控制器相较于市场已有的控制器，角度检测更加准确，功能更加全面。精确的角度检测，可准确识别是否是牛体触发，防止误触发，节约电能。自带正反转控制方式，外加可靠的电压电流监测模块，可有效保护牛体安全，避免对牛体造成伤害。随着智能养殖的进一步发展，机器取代人工，已经成为共识。通过长时间牛场的稳定测试运行，符合预期，已经取得不错的效果，具有一定的市场推广价值。
system: [SYSTEM]

user: 摘要 COVID-19在县级尺度内的病例时间变化与空间聚集特征及其影响因素在疫情防控工作中起到重要作用。本文尝试以县级尺度为单位，收集2020年1月26日至3月1日广东省确诊病例信息，采用空间自相关分析法、因子分析法、多元线性回归分析法等，探究广东省疫情在时空上的变化特征，并揭示其影响因素。结果表明，广东疫情特征主要表现为：时间上分为四个主要阶段即上升期、高发期、下降期、低水平波动期；空间上，病例呈现较强的空间聚集性且聚集强度逐渐增强，热点区县主要分布在东莞、广州、深圳市、佛山、中山、清远等地市；疫情情况与人口流动规模、客运量、人均GDP和限额以上餐饮业餐费收入显著相关。
system: [SYSTEM]

user: 1. 引言 新型冠状病毒(COVID-19)属于β属冠状病毒，主要传播途径为经呼吸道飞沫和密切接触传播，接触病毒污染的物品也可造成感染，在相对封闭的环境中长时间暴露于高浓度气溶胶情况下存在经气溶胶传播的可能 [ 1 ]。COVID-19的流行过程具有较强的聚集性特征，先在某地(如武汉华南海鲜批发市场)爆发成为一级传染点，感染者通过人口流动等形式携带病毒传播至下一级传染点，从而导致疫情在大范围内进行传播。因此，对于COVID-19病例的聚集区域(特别是县级尺度)的研究和原因分析，对城市开展疫情防控具有重要的借鉴意义。 目前，国内外已有大量关于COVID-19疫情的研究，主要集中在临床研究、流行病学调查研究、疫情对社会的影响等，如Backer J A [ 2 ] 等通过对早期疫情阶段在武汉市意外发现的88里确诊病例的旅行史和症状，分析COVID-19的平均潜伏期约为6.4天。中国疾病预防控制中心新型冠状病毒肺炎应急响应机制流行病学组 [ 3 ] 根据内地报告72，314例病例对COVID-19流行病学特征进行分析研究，新型冠状病毒肺炎传播流行迅速，随着人们返回工作岗位，需积极应对可能出现的疫情反弹。沈国兵 [ 4 ] 发现此次疫情已经给我国外贸和就业以及全球产业链和供应链带来严重不利冲击，加剧了我国巨量高校毕业生叠加农村城镇化的大量人口的就业压力。目前，针对病例聚集的研究尺度主要有国家、省、市、县等，如Kang D [ 5 ] 等从中国省级尺度研究疫情的空间相关性；曾永明 [ 6 ] 等从市级角度分析人口流动与城市早期新冠肺炎疫情的时空扩散和分布关系；张延炀 [ 7 ] 等从县级角度分析发现河南省COVID-19存在空间聚集性，传播风险经历弱–强–弱的变化过程。这些研究表明，COVID-19的时空聚集性在不同空间尺度上有不同特点，微观层面的COVID-19疫情空间聚集特点和因素分析对疫情防控工作的开展起到重要作用。 广东省截止2020年3月1日24时累计确诊病例已达1350例，累计出院1034例，累计死亡7例，其确诊病例仅次于湖北省，排全国第二。同时，广东省扼守中国南大门，是粤港澳大湾区重要组成省份，是我国重要交通枢纽和人流、物流中心，在全国疫情防控工作中具有战略性意义。因此，本文基于病例个体数据，研究分析广东省的疫情时间发展变化过程、病例性别和年龄构成特点，再基于县级尺度分析病例时间序列变化与空间聚集特征，并进一步分析空间聚集影响因素，为今后重大安全事件防控提供决策参考。
system: [SYSTEM]

user: 2. 数据与方法 2.1. 数据来源 病例数据全部来自广东省以及21个地级市卫健委官方网站和官方公众号新冠肺炎疫情情况通报的数据(2020年1月19日~2020年3月1日)。收集的病例信息包括所属县区、性别、年龄、确诊时间、是否有境外旅居史等。根据病例是否有境外旅居史以及感染原因，剔除境外输入病例，保留境内确诊病例。在2020年1月19日~3月1日，广东省共确诊1350例，经过数据处理与甄别，本文采用其中符合本次研究的1272例病例，约占总确诊病例的94%。 以全国地理信息资源目录服务系统提供的2019年版1:100万矢量地图为基础，根据需求利用ArcGIS 10.2提取出广东省县级行政区域矢量地图，包括64个市辖区、20个县级市、34个县、3个自治县，由于东莞市和中山市行政单位划分的特殊性，本文根据卫健委发布的病例归属信息，特将两个城市分别划分为34和24个镇。通过县级国标码将病例数据与广东省县区边界地理数据进行关联，利用EXCEL和ArcGIS 10.2建立广东省COVID-19的GIS数据库。 人口流动是疫情聚集性爆发的重要原因，由于数据的限制，本文采用广东省2019年常住人口和年末户籍人口的差来表示人口流动规模。2019年常住人口、2019年年末户籍人口、人均GDP、限额也上餐饮业餐饮费收入等数据均来自《广东省2020年统计年鉴》和各地级市统计年鉴。 2.2. 研究方法 本文主要从时间与空间两个维度对广东省疫情聚集性特征进行分析，在时间上基于logistic增长曲线拟合分析广东省确诊病例趋势变化过程，空间上分析广东省确诊病例聚集性变化趋势与聚集具体情况，最后结合因子分析与多元回归分析的方法，进一步探究疫情的空间分布影响因素。 2.2.1. 全局空间自相关 本文选取Moran’s I指数来识别新冠肺炎病例的全局空间自相关性，其公式为： I = n ∑ i = 1 n ∑ j = 1 n w i j ( x i − x ¯ ) ( x j − x ¯ ) n ∑ i = 1 n ∑ j = 1 n w i j ∑ i = 1 n ( x i − x ¯ ) 2 (1) 式(1)中，x i ，x j 为第i，j县区的新冠肺炎报告患病率，x为全省所有县区新冠肺炎平均报告患病率；w ij 为第i，j县区的空间矩阵；n为所有县区数。Moran’s I的取值范围为[−1， 1]，0表示空间不相关，[0， 1]表示空间正相关，[−1， 0]表示空间负相关。Z检验评价Moran’s I的统计学意义，计算公式如下： Z ( I ) = I − E ( I ) V A R ( I ) (2) 式(2)中：E(I)和VAR(I)分别是Moran’s I指数的期望和标准差。在给定的显著性水平下，若Moran’s I显著且为正值，表示确诊病例较多(或较少)的县区显著集聚，即存在正相关。若Moran’s I显著且为负值，表示县区间确诊病例数具有显著的空间差异，即存在负相关。当Moran’s I接近E(Ii)时，表示单元呈随机分布。 基于目前流行病学调查和研究结果，新型冠状病毒(COVID-19)潜伏期为1~14天，多为3到7天 [ 8 ]。广东省在2020年1月19日出现第一例新冠肺炎确诊病例，为深圳确诊，按照诊断报告日期，本文将新冠肺炎确诊病例在广东省的情况分为6个时段，以7天为一个周期，选取1月26日、2月2日、2月9日、2月16日、2月23日、3月1日这六天的广东省各县区累计确诊病例数进行空间自相关分析。 2.2.2. 热点分析 全局莫兰指数可以识别疫情在空时间上是否具有相关性，但是未能准确识别聚集区域和热点区域。故本文引入局域Getis-Ord G i * 指数来进一步识别广东省确诊病例在县级单位上的热点，主要利用ArcGIS 10.2进行计算Getis-Ord G i * 指数并输出图表，若指数得分为正值且较大，说明这个空间单元为高值聚类单元，即此单元为高病例单元，相反则为低值聚集单元。公式如下所示： G i * = ∑ j w i j x j − x ¯ ∑ j w i j ∑ j x j 2 n − x ¯ 2 n ∑ j w i j 2 − ( ∑ j w i j ) 2 n − 1 (3) (注：式(3)中所有符号的释义同式(1))。 2.2.3. 疫情空间分布影响因素分析 为了分析县级尺度下疫情的空间分布因素，从人口流动、交通客运、经济水平、餐饮零售中选取4个变量来进行实证研究(表1)。首先使用SPASS软件进行KOM值检验这4个因子的相关性，KOM值在0与1之间，当KOM值越接近1，意味着变量间的相关性就越强，越适合做因子分析；当KOM值越接近0时，意味着相关性越低，越不适合做因子分析。然后使用Bartlett球形度检验因子分析方法的适用性，sig值小于0.05则说明符合标准，数据呈球形分布，各变量相互独立，适合做因子分析。最后使用多元回归分析法，3个因子作为自变量，各县级累计确诊病例作为因变量，以此分析不同因子对疫情空间的影响力。 Table 1 Moran’s I statistics of the total number of confirmed cases in Guangdong Provinc 时间 Moran’s I 预期指数 P 方差 Z得分 2020年1月26日 0.175788 −0.005587 0.000016 0.001770 4.311434 2020年2月2日 0.231042 −0.005587 0.000000 0.001904 5.423624 2020年2月9日 0.244015 −0.005587 0.000000 0.001915 5.703119 2020年2月16日 0.249521 −0.005587 0.000000 0.001924 5.815945 2020年2月23日 0.253234 −0.005587 0.000000 0.001928 5.894653 2020年3月1日 0.253840 −0.005587 0.000000 0.001926 5.911127 表1. 广东省累计确诊病例数全局Moran’s I各统计量