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

user: 4. 实验结果和分析 为了估计算法的性能，使用的是标准的512 × 512的Lena图像，混沌系统的初始值 x 0 = 0.456 ， y 0 = 0.242 ，参数 u 0 = 5.432 ， r 0 = 20 ， N 0 = 1000 ，加密结果和解密结果如图4所示。 图4. 图像的加解密图 4.1. 密钥空间 为了图像加密的安全性，图像密码系统的密钥空间应该足够大，从而可以有效地抵抗穷举攻击。在我们提出的密码系统中，包括的密钥有256位的哈希值，混沌系统的初始值 x 0 ， y 0 ， u 0 ， r 0 ，迭代次数 N 0 ，如果计算机的精度是10 −14 ，提出的混沌系统的密钥空间是10 108 ，远大于要求的20 100 。 4.2. 直方图分析 图像的直方图反映了一幅图像中的像素值的分布情况。一个安全的加密方案应使加密后的图像的直方图是均匀的才可以抵制任何统计攻击。在图5中给出了明文图像和密文图像的直方图，可以清晰的看到密文图像的像素值分布均匀，而明文图像的直方图跌宕起伏，所以加密后很好的隐藏了明文图像信息。 图5. 直方图分析 4.3. 相关性分析 明文图像相邻像素间的高相关性是加密后的密文图像应该克服的缺点，所以要求加密后的密文图像的相关性应小于0.1，为了测试提出的加密方案中相邻像素的相关性，我们从明文图像和加密后的密文图像中随机选取了1000对相邻的像素值，利用公式(16)~(18)进行了相关性的计算。从图6中可看出，明文图像在水平，垂直，对角上的相邻像素点均具有较强的相关性，而密文图像中的相邻像素点没有相关性。表1是加密前后相邻像素值的相关系数。 r x y = cov ( x ， y ) D ( x ) × D ( y ) (16) cov ( x ， y ) = 1 N ∑ i = 0 N ( x i − E ( x ) ) ( y i − E ( y ) ) (17) E ( x ) = 1 N ∑ i = 0 N x i ， D ( x ) = 1 N ∑ i = 0 N ( x i − E ( x ) ) 2 (18) 图6. 明文图像和密文图像相邻像素相关性 Table 1 Correlation coefficient of adjacent pixel 方向 水平 垂直 对角 明文图像 0.9878 0.9777 0.9710 密文图像 0.0040 −0.0170 0.0120 文献 [ 4 ] −0.0140 −0.0241 0.0509 表1. 相邻像素相关系数 4.4. 敏感性分析 4.4.1. 密钥敏感性分析 为了保证密码系统的安全性，一个好的密码系统应该对密钥有较强的敏感性。所谓的密钥敏感性是指当密钥发生微小变化后，让其加密相同的明文图像后得到的密文图像与原始密文图像之间的差别情况。如果这两个密文图像有明显的差别，则表明使用的密码系统的密钥敏感性强。若两个密文图像的差别很小，则认为密钥敏感性差。我们以明文图像Lena为例，图7(a)是正确密钥解密出的图像。在这里，我们将原始密钥 x 0 改为 x 0 = 0.45600000000001 ，其它密钥保持不变，采用新密钥解密明文得到的解密图像，如图7(b)所示。同样的，让 y 0 = 0.24200000000001 ，其它密钥不变，解密出的图像为图7(c)所示，让 u 0 = 5.43210000000001 ，解密出的图像如图7(d)所示。对比发现微小改变密钥不能解密出图像，因此，它具有很好的密钥敏感性。 图7. 密钥敏感性分析 4.4.2. 明文敏感性分析 加密后的图像对明文越敏感，表明算法的性能越好，越能抵抗各种攻击 [ 4 ]。通过对明文图像做微小改变，如只将明文的一个像素值改变，然后比较改变后得到的密文图像与原来的密文图像间的差别。比较的方法是使用NPCR (像素数改变率)和UACI (平均改变强度)这两个定性变量。它们的表示方法如式(19)~(21)。 NPCR = 1 M × N ∑ i = 1 M ∑ j = 1 N D ( i ， j ) × 100 % (19) UACI = 1 M × N ∑ i = 1 M ∑ j = 1 N | C 1 ( i ， j ) − C 2 ( i ， j ) | 255 × 100 % (20) 其中 D ( i ， j ) 的定义如下 D ( i ， j ) = { 0 ， C 1 ( i ， j ) = C 2 ( i ， j ) 1 ， C 1 ( i ， j ) ≠ C 2 ( i ， j ) (21) M，N是图像的大小， C 1 和 C 2 是加密后的两个密文图像，它们只在一个像素点处的灰度值是不同的。通过随机选取明文图像中的像素值使其像素值加1，进行测试后，得出的NPCR和UACI的值分别是99.6091%和33.4579%，从中可以看出，提出的算法具有较好的明文敏感性。 4.5. 信息熵 信息熵是另一种反映图像信息不确定性的一种测试标准，灰度值的分布越均匀，信息熵越大，不确定性也就越大.计算公式如(22)所示： H = − ∑ i = 0 L p ( i ) log 2 p ( i ) (22) 其中L = 256，P(i)表示灰度值i出现的概率。一幅L = 256个灰度级的灰度图像，其信息熵最大值8在图像是完全均匀分布的随机图像时达到。一幅图像的信息熵越接近最大值，像素灰度值的分布越均匀。根据式(22)计算得到的明文图像Lena的信息熵仅为7.4434，而密文图像的信息熵结果达到7.9993，比已有的文献 [ 5 ] [ 13 ] 加密后得到的密文图像的信息熵结果更好，更接近最大值8，如表2所示。实验结果表明加密后的图像的灰度值分布均匀，有很好的信息熵。 Table 2 Information entrop 明文 密文 文献 [ 5 ] 文献 [ 13 ] 7.4434 7.9993 7.9973 7.9985 表2. 信息熵
system: [SYSTEM]

user: 5. 总结 在这篇文章中，我们利用改进后的混沌系统来进行图像加密。在密码系统中，我们首先将明文图像进行比特位分解，接着进行比特级置乱和扩散后，经过比特级合并，得到新的图像有效的打乱了原始图像的相关性，使其具有随机性。最后再进行了一次像素级的块置乱算法得到最终的密文图像。整个密码系统在进行一轮加密后就能获得良好的加密效果。实验结果表明，提出的加密算法能够抵抗各种攻击。
system: [SYSTEM]

user: 摘 要 植物转录因子可以调控众多下游基因的表达，在植物的生长发育、响应外界环境刺激等方面起着重要的调控作用。随着转录因子在模式植物和作物中基因转录调控的研究不断深入，转录因子研究方法也在持续创新。我们主要从蛋白质与蛋白质互作和靶基因与蛋白质互作这2个方面阐述植物转录因子研究方法，为植物转录因子的相关研究提供参考方法。蛋白质与蛋白质互作的方法主要有：酵母双杂交系统、荧光共振能量转移技术、双分子荧光互补技术、下拉实验和免疫共沉淀技术；靶基因与蛋白质互作的方法主要有：酵母单杂交系统、染色质免疫沉淀技术、凝胶电泳迁移技术和荧光素酶报告系统。 关键词 :植物转录因子，转录调控，蛋白质与蛋白质互作，靶基因与蛋白质互作 Copyright © 2020 by author(s) 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. 引言 转录因子(Transcription Factor， TF)也称反式作用因子，其可直接或间接与基因启动子区域中顺式作用元件特异性结合、对基因转录进行调控的蛋白质 [ 1 ]。其主要功能是激活或抑制基因的转录效应。典型的转录因子一般具有4个功能结构域，即DNA结合结构域(DNA-Binding Domain)、转录调控结构域(Transcription Regulation Domain)、核定位信号区(Nuclear Location Signal)和寡聚化位点区(Oligomerization Site)。DNA结合结构域特异性结合顺式作用元件；转录调控结构域激活或抑制靶基因表达；核定位信号区是转录因子中富含精氨酸和赖氨酸残基的核定位区域，转录因子进人细胞核的过程受该区段控制；寡聚化位点区是不同转录因子借以发生相互作用的功能域。这些结构域共同调节靶基因的转录起始 [ 1 ] [ 2 ]。 转录因子在植物的生长发育及对外界环境刺激等方面发挥着重要的调控作用。从1987年Paz-AerS首次报道玉米转录因子基因的克隆以来，相继从高等植物中分离出的调控干旱、高盐、低温、激素、病原反应及生长发育等相关基因表达的转录因子已达数百种 [ 3 ]。植物中约有58个转录因子家族 [ 4 ]，其中有6个转录因子家族据研究报道参与了生物和非生物胁迫反应，如AP2/ERF (APETALA2/乙烯响应因子)，bHLH (基础螺旋–环–螺旋)，MYB (与成髓细胞相关)，NAC (无顶端分生组织(NAM))，拟南芥转录激活因子(ATAF1/2)，杯状子叶(CUC2)，WRKY和bZIP (碱性亮氨酸拉链) [ 5 ] [ 6 ]。 随着转录因子在模式植物和作物中基因转录调控的研究不断深入，转录因子研究方法也在持续创新。我们主要从蛋白质与蛋白质互作和靶基因与蛋白质互作这2个方面阐述植物转录因子研究方法。
system: [SYSTEM]

user: 2. 蛋白质与蛋白质互作 2.1. 酵母双杂交系统 酵母双杂交系统(Yeast two-hybrid system， Y2H)于1989年首次开发，彻底改变了寻找和鉴定互作蛋白的过程 [ 7 ]。该系统根据酵母GAL4转录因子的特性，由DNA结合结构域(Binding Domain， BD)和转录激活结构域(Activation Domain， AD)组成。DNA结合结构域由GAL4 N端1-147位氨基酸残基组成，转录激活结构域由C端768-881位氨基酸残基组成 [ 8 ]。这两个结构域之间由大的结构区域隔开，独立且功能稳定。一般情况下，单独的BD结合DNA上游激活序列(Upstream Activating Sequence， UAS)但不引起转录，单独的AD不能与UAS结合，只有当两个相互作用的蛋白质分别与BD和AD融合时，它们的相互作用产生嵌合型转录因子从而引起启动子下游基因的转录。 依据此特性，将编码已知蛋白(诱饵蛋白，Bait)的基因构建到含有BD序列的载体上，在酵母中表达产生融合蛋白BD-Bait，将编码未知蛋白(猎物蛋白，Prey)的基因构建到含有AD序列的载体上，在酵母中表达产生融合蛋白AD-Prey，当BD-Bait与AD-Prey在酵母核中相互作用时，DNA-BD和AD接近以恢复功能性GAL4转录激活因子，前者与报告基因(如lacZ、ADE2、HIS3和MEL1)的UAS结合，从而启动报告基因的表达 [ 8 ]。因此可以通过判断酵母克隆能否在与报告基因对应的营养缺陷型培养基上生长来证明两个蛋白是否存在相互作用。 与其它的研究蛋白质相互作用的方法相比，假阴性、假阳性以及互作蛋白的强制性核定位是Y2H的局限，当然Y2H有其独特的优点。首先，Y2H不仅可以检测稳定的相互作用蛋白，还可以检测弱而短暂的蛋白相互作用 [ 9 ]。第二，由于Y2H是在酵母细胞核内进行的，因此检测的蛋白质可能以天然构象存在 [ 8 ]。第三，Y2H是对生化方法(如免疫共沉淀和质谱分析)的补充 [ 9 ] [ 10 ]。近年来Y2H得到了很大的改进，不仅可以应用于蛋白质-DNA相互作用 [ 11 ] 和酵母三杂交 [ 12 ]，还可以用于膜蛋白，DNA结合蛋白和RNA结合蛋白的相互作用研究 [ 13 ]。 2.2. 荧光共振能量转移技术 高端共聚焦显微镜和荧光蛋白的发展促进了亚细胞尺度的蛋白质分析。荧光共振能量转移(Fluorescence Resonance Energy Transfer， FRET)是指处于激发态的供体荧光蛋白将能量以非辐射的方式转移到附近的受体分子 [ 14 ]。发生FRET必须有三个先决条件：1) 供体发射和受体吸收的光谱重叠区域必须有有效的能量转移；2) 荧光基团之间的距离必须小于10 nm；3) 荧光基团的双极取向必须以理想的平行方式对齐。虽然FRET常用于蛋白质相互作用，但其实际上是通过理论计算来检测 [ 15 ]。 有很多方法可以检测到FRET，每种方法都各有好处和缺陷。目前有三种最常用的技术：1) 敏化发射测量(FRET-SE)；2) 受体光漂白(FRET-AB)；3) 荧光寿命成像(FRET-FLIM)。 敏化发射测量(FRET-SE)主要用于动态检测活细胞中的FRET效率 [ 16 ]。FRET-SE的限制在于每次实验需要利用至少3个对照样本对光学检测系统以及荧光基团的光学性质进行事先校正，一旦校正完成，后续实验不能再改变任何系统参数，且每次实验都要重新校正，这极大地增加了实验的工作量和难度。 受体光漂白(FRET-AB)通过检测光漂白受体漂白前后供体的荧光强度来获得FRET效率 [ 17 ]。这是一种直接的方法，不需要高端显微镜，也不受光谱影响，可以自己控制样品。 荧光寿命成像显微术(FRET-FLIM)主要通过测量受体存在和不存在时的供体分子的荧光寿命来确定FRET效率 [ 18 ]。这是一种不受光谱，供体或受体浓度差异，激发强度变化和光漂白影响的技术 [ 19 ]。其限制因素主要是系统价格昂贵，而且需要熟练的操作人员进行操作。 2.3. 双分子荧光互补技术 双分子荧光互补技术(Bimolecular Fluorescence Complementation， BiFC)最初从大肠杆菌发展而来，后来广泛适用于哺乳动物和植物系统 [ 20 ]。其原理是将荧光蛋白在合适的位点切开形成不发荧光的N端和C端片段，这2个片段借助融合于其上的目标蛋白的相互作用，彼此靠近重新形成具有活性的荧光蛋白，并在该荧光蛋白的激发光激发下发射荧光。尽管BiFC被广泛应用，但它受到两个主要限制：1) 估计半衰期为10年的分裂荧光蛋白片段的不可逆互补排除了相互作用动力学的分析；2) 大多数分裂片段可以不受阻碍地自发重新组装。这些缺陷增加了假阳性的可能。 2.4. Pull down和免疫共沉淀技术 Pull down技术主要在体外验证两种蛋白之间是否存在直接的相互作用。主要有两种常用的pull down方法：一种是用谷胱甘肽巯基转移酶(Glutathione s-Transferase， GST)标签(GST-tag)的亲和纯化柱进行的，称为GST-pull down；一种是用带组氨酸(Histidine， His)标签(His-tag)的亲和纯化柱进行的，称为His-pull down [ 21 ]。以GST-pull down为例，首先是获得带有GST标签的融合蛋白，然后将带有GST标签的融合蛋白亲和到谷胱甘肽层析柱上，接着将与融合蛋白发生相互作用的蛋白被亲和到谷胱甘肽层析柱上，最后切掉标签，洗脱得到相互作用的蛋白。 免疫共沉淀技术(Co-Immunoprecipitation， Co-IP)以抗体和抗原之间的专一性结合作用为基础，用于检测和确定生理条件下蛋白质之间的相互作用。其基本原理是将细胞在非变性条件下裂解，在细胞裂解液中加入特异性的抗体，这样抗体可以和已知的抗原形成抗原–抗体复合物，如果在系统中存在与已知抗原相互作用的蛋白质，该蛋白质也会以复合物的形式沉淀下来，经过洗脱就可以得到与已知抗原相互作用的蛋白质，然后进行SDS-PAGE及Western blotting分析 [ 22 ]。该技术研究的相互作用的蛋白质都是经翻译后修饰的处于天然状态的，因此可避免人为的影响，同时可分离得到天然状态的相互作用蛋白复合物。但Co-IP难以检测到弱或瞬时相互作用。此外，多个相互作用蛋白质复合物的存在也可引起Co-IP中的相互作用。
system: [SYSTEM]

user: 3. 靶基因与蛋白质互作 3.1. 酵母单杂交系统 酵母单杂交系统(Yeast one-hybrid system， Y1H)由酵母双杂交系统发展而来，是一种在酵母细胞内分析与鉴定转录因子和DNA顺式作用元件相互结合的有效方法 [ 12 ]。依据此原理构建目的基因与转录激活结构域融合的载体，将顺式作用元件和报告基因分别连在启动子的上游和下游，转录激活域融合蛋白与特异的DNA序列结合激活启动子下游报告基因表达。因此可以通过判断酵母克隆能否在与报告基因对应的营养缺陷型培养基上生长来证明目的蛋白能否与DNA序列结合。由于酵母单杂交系统常常受所用报告基因的影响，较高的假阳性一直是该技术的“瓶颈”问题。此外，在酵母单杂交系统中，假阴性一直被忽略。如某些外源基因的蛋白表达产物可能对酵母本身有毒性作用，影响酵母细胞的生长和表型，而这其中可能就包含着与靶DNA序列相互作用的蛋白。 3.2. 染色质免疫沉淀技术 染色质免疫沉淀技术(Chromatin Immunoprecipitation， ChIP)是一种利用抗原抗体反应的特异性，在染色质水平上真实反映体内基因组DNA与蛋白质结合情况的转录调控技术。其原理是在生理条件下把细胞内的DNA与蛋白质交联在一起，通过超声波将染色质随机切成小片段，利用抗原抗体的特异性识别反应，将与目的蛋白结合的DNA沉淀下来，检测目的片段得到蛋白质与DNA相互作用信息 [ 23 ]。Chip实验涉及众多的实验步骤，包括细胞固定、染色质断裂、染色质免疫沉淀、解交联、DNA的纯化以及DNA的鉴定，而且结果的重复性较低，因此对Chip实验过程的每一步都应设计相应的对照。 近年来发展起来的Chip-seq技术将染色质免疫沉淀技术(Chip)与芯片技术和第二代测序技术(seq)相结合，用来进行靶基因的高通量筛选，为研究目的蛋白与整个基因组相互作用提供了可能。 3.3. 凝胶电泳迁移技术 凝胶电泳迁移技术(Electrophoretic Mobility Shift Assay， EMSA)又称为电泳迁移率改变分析技术和凝胶阻滞实验。这是一种利用探针体外分析DNA与蛋白质相互作用的凝胶电泳技术 [ 24 ]。其原理是蛋白质和DNA结合后增加了相对分子质量，在非变性聚丙烯酰胺凝胶电泳中，DNA-蛋白质复合物与单一的DNA片段相比，前者迁移速度明显较慢。具体方法是将含有假定蛋白质结合位点的已标记DNA片段与纯化蛋白进行孵育，然后将反应产物进行非变性聚丙烯酰胺凝胶电泳分析。传统的32P标记探针的凝胶阻滞实验，具有很高的灵敏度，但会接触到危险性的放射性同位素且不容易定量分析。即使现在有很多研究报道利用非放射性标记的凝胶阻滞实验，但非放射性标记探针的凝胶阻滞试剂盒的费用很高。 3.4. 荧光素酶报告系统 荧光素酶报告系统(Luciferase reporter assay)指的是以荧光素为底物来检测荧光素酶活性的检测系统 [ 25 ]。荧光素酶是一种广泛使用的报告基因，其可以催化荧光素氧化成氧化荧光素，在荧光素氧化的过程中，会发出生物荧光(bioluminescence)。荧光素和荧光素酶这一发光系统具有极高的灵敏性。荧光素酶不仅可以作为一种报告基因进行检测分析，还可以用于化学污染物的监测与分析，此外还可以用于RNA干扰的研究、蛋白质之间作用的研究及微生物检测等，应用前景十分广阔 [ 26 ]。
system: [SYSTEM]

user: 4. 结论 转录调控在拟南芥、水稻等模式植物和作物中的研究越来越多。转录调控技术的应用有助于更好地阐明植物内部的转录调控机制。在这篇综述中，我们从蛋白质与蛋白质互作和靶基因与蛋白质互作这2个方面总结了植物转录因子主要分析方法，以期为植物转录因子的相关研究提供参考方法。
system: [SYSTEM]

user: 摘 要 守恒律是非线性方程可积的一个重要指标，在偏微分方程系统的研究中扮演着重要的角色。本文基于Lax对，研究了一维散焦Hirota方程在非零背景下的无穷守恒律及其高阶守恒量。 关键词 :Hirota方程，Lax对，守恒律 Copyright © 2020 by author(s) 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. 问题背景及主要结论 随着现代科学的发展，线性微分方程已经无法确切的描述现实社会中各种现象出现的个体内在关系，随之而来的，非线性科学得以快速的发展，并在数学、物理、流体等领域中广泛应用。也正源于此，越来越多科学家投身到非线性偏微分方程的研究之中，并取得巨大的成就。而在这之间，数学家发现许多具有广泛应用背景的方程如：KdV方程、Schrödinger方程、sine-Gordon方程等都存在孤立子解，于此，孤立子理论得以产生并迅速发展。同时，因为光孤子在光纤通信中的重要作用，也使得众多的数学家开始致力于描述光孤子运动的非线性发展方程的研究，而这其中最具有代表性的方程就是Hirota方程。 广义的Hirota方程如下： i ( ψ t + μ 1 ψ + ( ν 1 + μ 1 x ) ψ x ) + ( ν 2 + μ 2 x ) ( ψ x x + 2 | ψ | 2 ψ )   + 2 μ 2 ( ψ x + ψ ∫ − ∞ x | ψ | 2 d x ) + i ν ( ψ x x x + 6 | ψ | 2 ψ x ) = 0 。 (1.1) 针对广义的Hirota方程，Poresezian [ 1 ] 、田播 [ 2 ] 等人分别利用Bäcklund、Hirota双线性方法得出了方程的孤子解，并分析了该孤子解的性质。而当方程(1.1)中系数 μ 1 = ν = ν 1 = 0 时候，方程就退化为广义的非线性Schrödinger方程： i ψ t + ( ν 2 + μ 2 x ) ( ψ x x + 2 | ψ | 2 ψ ) + 2 μ 2 ( ψ x + ψ ∫ − ∞ x | ψ | 2 d x ) = 0 。 (1.2) 该方程对于描述连续极限条件下不同磁性相互作用的动力学特性具有明显的效果，该Schrödinger方程的Lax对在 [ 3 ] 文献中也已经给出。当方程(1.1)中系数 μ 1 = μ 2 = ν 1 = 0 ， ν 2 = 1 ， ν = ϵ ，方程(1.1)就变成了变系数Hirota方程： i ψ x + ψ x x + 2 | ψ | 2 ψ + i ϵ ( ψ x x x + 6 | ψ | 2 ψ x ) = 0 。 (1.3) 方程(1.3)的可积性在文献 [ 4 ] 中给出，同时，Radhakrishnan等人 [ 5 ] 利用Hirota双线性方法得出了该方 程的孤子解。当方程(1.1)中系数 μ 1 = μ 2 = ν 1 = 0 ， ν 2 = 1 2 ， ν = 2 4 ，得到方程： i ψ x + 1 2 ( ψ x x + 2 | ψ | 2 ψ ) + 2 4 i ( ψ x x x + 6 | ψ | 2 ψ x ) = 0 。 (1.4) 对式子(1.4)作坐标变换： x = 1 2 x ，即可得到方程： i ψ x + ψ x x + | ψ | 2 ψ + i ( ψ x x x + 3 | ψ | 2 ψ x ) = 0 。 (1.5) 相似于方程(1.5)，本文深入研究了下列一维Hirota方程： i ψ x + ψ x x − | ψ | 2 ψ + i ( ψ x x x − 3 | ψ | 2 ψ x ) = 0 。 (1.6) 其中： ψ = ψ ( x ， t ) ∈ ℂ ， ( x ， t ) ∈ ℝ × ℝ 。这是一个散焦立方非线性方程，包括三阶色散项和非线性色散项。守恒律反映了某物理量不随时间改变的一种现象。而在孤子理论中，人们发现这样一个事实，如果非线性发展方程具有孤立子解，那么几乎也具有无穷守恒律。因此，就一个孤子系统而言，寻找其无穷守恒律，对于证明此系统的可积性，甚至在证明其适定性方面有重大意义。下面我们将基于Lax对 [ 6 ]，来研究散焦Hirota方程(1.6)在非零背景下的无穷守恒律及其高阶守恒量。 定义1.1 [ 7 ] 给定一般的非线性偏微分方程： F ( x ， t ， u ) = 0 ， (1.7) 其中： u = u ( x ， t ) 是关于x和t的函数，而 F ( x ， t ， u ) 是 x ， t ， u 及其导数的函数。若存在一对连续可微函数 ω ( x ， t ， u ) 和 J ( x ， t ， u ) ，使得u按照方程(1.7)发展时候满足关系 ∂ t ω ( x ， t ， u ) = ∂ x J ( x ， t ， u ) ， (1.8) 则称此关系式为式(1.7)的守恒律，而 ω ( x ， t ， u ) 和 J ( x ， t ， u ) 分别称为守恒密度与连带流。 如果当 | x | 趋于无穷时，密度与流充分快地趋于零，则将守恒律在整个数轴上对x积分得： d d t ∫ − ∞ ∞ ω ( x ， t ， u ) d x = 0 ，(1.9) 则积分 ∫ − ∞ ∞ ω ( x ， t ， u ) d x 与时间t无关，是方程(1.7)的守恒量。 注：设二阶矩阵 M = ( m 1 m 2 m 3 m 4 ) ， 我们定义如下符号意义： M o f f = ( 0 m 2 m 3 0 ) ，     M d i a g = ( m 1 0 0 m 4 ) 。 定理1.1 HIROTA方程(1.6)具有无穷守恒律。 推论1.1 已知方程(1.6)，则有 S ( ψ ) = ∫ ℝ [ | ψ x x | 2 + 3 | ψ | 2 | ψ x | 2 + 1 2 ( ψ ¯ ψ x + ψ ψ ¯ x ) 2 + ( 1 − | ψ | 2 ) 2 ( 1 + 1 2 | ψ | 2 ) ] d x     。 (1.10) 是其高阶守恒量。
system: [SYSTEM]

user: 摘 要 青光眼滤过术后合并硬核白内障，同时角膜内皮细胞密度低，前房条件差，经典白内障手术风险大，并发症发生率高，我们收治1例，经后路超声粉碎摘除白内障，一期植入后房型人工晶体，术后效果良好，患者不仅脱盲而且脱残，现与同仁分享。患者，女，76岁，主诉：右眼青光眼滤过术后10余年，视力下降5年加重2年。眼科检查：视力：右0.01，矫正无助，左无光感，眼压：右10.3 mmHg，左10.8 mmHg。右眼上方球结膜见扁平滤过泡，前房中央深度约2.5 CT，周边约1/4 CT，前房上方周切口附近见半透明胶样组织堆积，且局限性虹膜前粘连，瞳孔缘见机化膜，4:00~5:00后粘连，晶状体皮质及核棕黑色混浊，眼底窥不入。右眼角膜内皮细胞计数751个/mm 2 ，B超示玻璃体轻度混浊伴后脱离。诊断：1) 白内障OD；2) 青光眼滤过术后OU。常规术前准备，于2019年3月28在神经阻滞麻醉下行右眼瞳孔成形联合经后路晶状体超声粉碎吸除及后房型人工晶体植入术，手术顺利。术后抗生素及糖皮质激素滴眼液点眼预防感染、减轻炎症反应。术后第四天右眼视力0.3，眼压13 mmHg，角膜透明，人工晶体位正，治愈出院。术后6月随诊，右眼视力0.25，矫正视力0.4，眼压12 mmHg，角膜透明，人工晶体位正，复查角膜内皮细胞计数740个/mm 2 。结论：青光眼滤过术后合并V级核白内障，同时角膜内皮少，前房条件差，采用经后路超声粉碎摘除白内障，一期植入人工晶体是一种行之有效的方法，值得推广。 关键词 :硬核，白内障，经后路晶状体超声粉碎，青光眼，角膜内皮 Copyright © 2019 by author(s) 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. 引言 青光眼滤过术后常并发白内障或白内障加重，当白内障发展到一定程度影响视力时，需要手术治疗。青光眼滤过术后合并白内障的手术属于复杂性白内障手术，因为青光眼滤过术后常常伴有角膜内皮数不足，前房浅，虹膜弹性差，术中又需注意青光眼滤过泡保护等等。因此年轻的眼科医生或经验不足的眼科医生遇到青光眼滤过术后白内障常常望而怯步，同时患者因手术风险大难以接受手术治疗，当患眼视力严重障碍时，白内障的发展由软核变为硬核，手术难度和风险进一步增加，甚至失去经典白内障超声乳化联合人工晶体植入术的手术时机。今年3月我院收治一例单眼青光眼滤过术后V级核白内障患者，术前视力0.01，且眼前段条件较差，采用经睫状体扁平部超声粉碎摘除白内障联合后房型人工晶体植入术，取得了良好的效果，兹报告如下。
system: [SYSTEM]

user: 2. 病例介绍 患者，女，76岁，以“右眼青光眼滤过术后10余年，视力下降5年加重2年”为主诉于2019年3月27日至我院就诊。既往高血压病病史5年。眼科检查：视力：右0.01，矫正无助，左无光感，眼压：右10.3 mmHg，左10.8 mmHg。双眼上方球结膜见滤过泡，右眼扁平，左眼囊状，双眼角膜透明。右眼中央前房深度约2.5 CT，左眼中央前房深度约3.5 CT，双眼周边前房约1/4 CT；右眼前房上方周切口附近见半透明胶样组织堆积，且局限性虹膜前粘连。右眼瞳孔直径约3 mm，对光反应迟钝，瞳孔缘见机化膜，4点至5点处后粘连；左眼瞳孔大小约4 mm × 3 mm，欠圆，向上移位，对光反应迟钝。右眼晶状体皮质及核棕黑色浑浊(Emery核分级V级)，左眼人工晶体位正。眼底：右眼窥不入；左眼视盘界清色苍白，C/D = 1.0，视网膜红润、在位。辅助检查：角膜内皮细胞计数：右751个/mm 2 ，左654个/mm 2 ；B超：双眼玻璃体轻度混浊伴后脱离。诊断为：1) 白内障OD；2) 青光眼滤过术后OU；3) IOL植入术后OS；4) 高血压病。 入院后术前常规准备，3月28在球后神经阻滞麻醉下行右眼瞳孔成形联合经后路晶状体超声粉碎吸除联合人工晶体植入术。手术方法：11点角膜缘隧道切口，前房注入粘弹剂，分离粘连虹膜，瞳孔成形。建立标准的玻切三通道(3点、5点、10点)，用两把MRV刀自10点、3点巩膜通道口从两侧晶状体赤道后部刺入晶状体核内，保留前囊膜，机械性碎核，然后超声粉碎晶状体核碎块并吸除，清除坠入玻璃体腔的晶状体核块，玻切头吸刮前囊膜下上皮细胞。检查视盘界清，色粉红，C/D < 0.3，视网膜红润。经角膜缘隧道切口，前房注入粘弹剂，将AR40e折叠人工晶体植入后房，襻固定于睫状沟内，灌洗前房，水密前房，术毕。术后第一天右眼视力0.03，眼压T + 1，角膜轻度水肿，上皮下小水泡，房水光斑(+)，人工晶体位正。处理：经角膜隧道口前房放液至眼压17 mmHg，抗生素及糖皮质激素滴眼液滴眼预防感染、减轻炎症反应。术后第二天，右眼视力0.2，眼压12 mmHg，角膜透明。术后第四天右眼视力提高为0.3，眼压13 mmHg，角膜透明，治愈出院。术后6月复诊，右眼视力0.25，矫正 + 1.50 DS − 3.50 DC × 95 ° = 0.4，眼压12 mmHg，角膜透明，人工晶体位正。复查角膜内皮细胞计数740个/mm 2 ，右眼30˚中心视野检查大致正常。
system: [SYSTEM]

user: 3. 讨论 目前，青光眼滤过术后硬核(III级以上硬度)白内障的治疗方法有三种，即经典白内障超声乳化吸除联合囊袋内人工晶体植入术、小切口白内障囊外摘除联合囊袋内人工晶体植入术及经睫状体扁平部晶状体超声粉碎摘除联合后房型人工晶体植入术 [ 1 ]。经典白内障超声乳化联合囊袋内人工晶体植入术的适应症为白内障核硬度III~IV级，角膜内皮细胞数在1300个/mm 2 以上，瞳孔散大≥4 mm以上，前房中等深度等。小切口白内障囊外摘除联合人工晶体植入术的适应症为核硬度V级，角膜内皮细胞数大于1000个/mm 2 ，而小于1300个/mm 2 ，切口不影响功能性滤过泡，虹膜无广泛的前粘连等。经睫状体扁平部超声乳化摘除白内障联合后房型人工晶体植入术适应症为白内障核硬度在IV级或V级，角膜内皮数≤1000个/mm 2 ，前房偏浅、瞳孔不能散大或前后粘连。本文病例术前经过详细询问病史，发现青光眼滤过术后10余年，曾数次到其他医院就诊，因眼部条件差，未给予手术治疗，随着病情发展，生活不能自理，患者十分痛苦，来我院就诊，仔细眼部检查，白内障核硬度V级，角膜内皮细胞751个/mm 2 (<1000个/mm 2 )，前房浅，瞳孔不能散大且后粘连等。若采用前两种手术方法摘除白内障手术风险大，效果可能差，甚至出现严重并发症，故选择非常规摘除白内障的方法即经睫状体扁平部超声乳化摘除白内障联合后房型人工晶体植入，术中保留前囊膜，穿刺刀机械碎核，在玻璃体腔内超声粉碎核块并吸除，减少超声乳化对角膜内皮的损伤及对前房的影响，然后经角膜切口植入折叠式人工晶体，襻固定于睫状沟，术程顺利。本文病例术后效果良好，术眼视力由术前0.01恢复到0.4，青光眼滤过泡没有受到影响，眼压正常，角膜透明。患者生活质量得以提高，患者非常满意。该术式的特点为术中能最大限度地减少对角膜内皮细胞的影响 [ 2 ] [ 3 ] [ 4 ]，本病例仅丢失11个/mm 2 ，与Eom、Diddie等学者的研究结果一致 [ 2 ] [ 5 ]。可是术中也增加了视网膜损伤的风险，手术费用也高，增加了患者的经济负担。我们建议该手术应由经验丰富的前后节医生完成，减少术中创伤，降低手术风险。此外，青光眼滤过术后当合并白内障影响视力时，应尽早行经典白内障超声乳化摘除联合囊袋内人工晶体植入，这时晶状体核硬度低，手术风险小，效果好，否则，手术风险增大，甚至影响手术效果。
system: [SYSTEM]

user: 4. 同意书 该病例报道已经获得患者的知情同意。
system: [SYSTEM]

user: 摘 要 首次采用结构方程模型(SEM)就城市化与工业化对房价的不同作用路径与模式进行了分析。探究结果表明：城市化与工业化对省际房价的作用路径与模式存在着明显的差别，虽然二者对房价均具有正向的直接影响，但二者的间接作用路径差异较大。其中，工业化通过对空气污染状况的作用来影响房价的变化，城市化通过对居民生活状况(收入、人均GDP)与废气排放的角度来影响房价。除此之外，工业化还通过调控产业结构对房价产生影响，而城市化中人口比重对房价的作用是“倒U型”，市区人口比重对房价作用存在阈值(转折点为0.615)，即在一般情况下，当市区人口比重超出该阈值，将会导致该区域房价下降。因此，从调控房价的角度来看，应该将重心放在城市化的作用路径上。 关键词 :房价，城市化，结构方程模型，工业化 Copyright © 2019 by author(s) 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. 引言 随着我国经济的快速稳定发展以及城市化的不断加深，使得我国房价呈现出快速的上升势头以及人们对于住房的需求的增强，且该现象在经济发展程度较高的城市中尤为严重。随着我国城市化进程的不断加深，使得住房价格持续攀高，这也进一步加深了地区房价极端分化的现象。区域房价的分化以及虚高难已对经济发展产生较为重要的推力，甚至可能会导致更为严重的金融风险。因此，探究影响房价的因素及其作用路径对于稳定房价与制定相关政策具有重要的意义。在本文中，从城市化与工业化的视角出发，探究的重点在于城市化与工业化是否对房价具有显著作用，二者是通过哪种路径来影响房价。而从诸多研究的结果中可以发现，二者是影响城市房价的主要原因。但是，在现有的探究中并没有深入分析二者对我国房价的作用路径。一般情况下，某一要素对房价的作用路径主要有4种：1) 只有直接影响，即该影响因素直接作用于房价；2) 既有直接作用又有间接作用，即该影响因素不仅对房价具有直接作用，还通过其他媒介(中间因素)对房价产生影响；3) 只有间接影响，即该因素仅仅通过其他媒介对房价产生影响；4) 既没有直接影响也没有间接影响。而在现有的研究中多采用参数或非参数回归的方式来探究城市化、工业化对房价的影响，但该类处理方法只是探究因素间的直接联系，而无法反映出各影响因素的间接作用以及作用路径，更加无法反映出各因素影响途径间的差异性。 因此，本文首次尝试将结构方程模型(SEM)应用到房价探究领域中，为更进一步探究我国房价提供新思路。该方法是一种既包含显变量又包含潜变量的多方程模型，其不仅能够揭示变量间的直接影响路径还能够发现其间接影响路径，从而更加系统化地探究各变量对房价的作用机制。
system: [SYSTEM]

user: 2. 文献综述 国内外研究学者研究了城市房价的影响因素，Gabriel (1999)、Vandana (2012)、骆永民(2011)、邓翔(2013)和刘广平(2016) [ 1 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ] 通过实证分析研究影响城市房价的主要原因，结果表明城市化和工业化是影响城市房价的主要原因。关于城市化对房价的影响存在两种看法，如：陈石清(2008)、况伟大(2013)、王立平与李艳萍(2014) [ 6 ] [ 7 ] [ 8 ] 等通过实证分析提出城市化水平的提高对房价具有良好的促进作用；姜松和王钊(2014) [ 9 ] 通过建立省际面板模型却得出相反的结果，即城市化会在一定程度上抑制房价的上涨。关于工业化对房价的作用，从张协奎(2009)、黄淳(2009)以及于守华(2013) [ 10 ] [ 11 ] [ 12 ] 的研究中可以看出工业化对房地产业具有显著的促进作用；而楚尔鸣(2018) [ 13 ] 从全国的层次角度出发，提出工业化的中心是产业结构的调整，且产业结构对城市化既有直接作用也存在着通过房价的中介效用；另一方面，从王拉娣(2016)、刘刚(2016)、鞠方(2017)以及陈晓与张文杰(2017) [ 14 ] [ 15 ] [ 16 ] [ 17 ] 的研究中可以看出，城市化与工业化进程的推进下，会引起居民收入的增加、人口的流动以及就业机会的增加，使得城市人口的比重不断上升，进而导致对住房需大于供，使得房价不断的上涨。关于污染状况对房价的作用，陈永伟与陈立中(2012) [ 18 ] 提出地区空气质量对该区域房价具有较高的关联性，且人们愿意提供降低空气质量而支付更多的房租；刘晓罡(2017) [ 19 ] 也提出相类似的结论，并且指出该现象在中高端购房者中更为明显。而从杜雯翠(2013) [ 20 ] 的研究中可以发现城市化对空气状况具有双向效应；从Dinda (2004) [ 21 ] 的分析结果中则可以看出提工业化对空气状况既有直接作用。 整合上述文献不难看出，工业化与城市化对房价的作用路径并不仅仅是直接作用，还有可能通过其他中介途径(如：居民收入与就业、空气污染状况、产业结构的变化等等)产生间接作用。更甚是可能仅存在对房价的间接作用，而直接作用会在引入间接因素后消失，即所得到的直接作用仅仅是一种虚假作用。但现有的研究中多是分析对房价的直接影响，并没有考虑到间接作用。除此之外，大多文献在模型中常常会遗漏一些起到“中间桥梁”作用的因素，使得所得到结果出现虚假作用。
system: [SYSTEM]

user: 3. 研究方法与数据处理 3.1. 研究方法 为避免上述现象的出现，本文采用SEM分析进行探究城市化、工业化对房价的影响路径，具体理论模型见图1。 图1. 影响因素路径分析理论模型 SEM模型：结构方程模型简记为SEM，与传统计量方法的差异主要在于其可以容纳可观测到的显性变量，也可以容纳无法直接观测获取的潜在变量。该模型主要由两部分组成：测量方程与结构方程，分别反映的是潜变量与指标、潜变量间的联系。具体模型如下： 测量方程： y p × 1 = Λ y p × m η m × 1 + ε p × 1 x q × 1 = Λ x q × n ξ n × 1 + δ q × 1 结构方程： η m × 1 = Β m × m η m × 1 + Γ m × n ξ n × 1 + ε m × 1 其中，x、y分别表示外生指标与内生指标； 、 Λ y 为因子载荷； η 、 ξ 分别为内生与外生潜变量； 、 表示路径系数值； δ 、 ε 表示误差项。 而在文中选择使用标准化法对测量模型进行设置，并对根据模型的性质判断测量方程的合理性，当测量方程具有较高的合理性时，采用CFI、TLI、AIC以及RMSEA等多项指标评估结构模型。 3.2. 变量选择与来源 在本文中，根据理论模型(图1)以及参考相关研究文献，本文选取各潜变量对应的指标进行分析，具体见表1。 Table 1 Index of housing price and its influencing factor 房价及其影响因素 指标 房价 住宅商品房平均价格 城市化 市辖区人口比重与人口密度 居民生活状况 城镇居民人居收入、就业人员数以及人均GDP 工业化 第二产业增加值所占的比重 废气排放状况 SO 2 、NO x 以及粉尘排放密度 产业结构 第三产业增加值与第二产业增加值的比例 空气质量 AQI、PM 2.5 、PM 10 、SO 2 、CO、NO 2 浓度年均值 自然状况 年平均气温、日照时数、降水量 表1. 房价及其影响因素指标 城市化的重点在于人口的流动。而人口流动主要表现在2点：城镇人口密度与比重的变化。因此，在本文中选取市区人口密度与比重作为度量标准。工业化的中心是产业结构，尤其是第二产业，因此，本文选取第二产业增加值所占比例作为度量指标。 居民收入水平高低以及就业人员数量的多少会影响人民购房的意愿。一般情况下，当人们整体具有较高的收入以及稳定的工作时，该区域的房价基本上都会偏高，反之亦反。因此，在文中选取城镇居民人均收入以及就业人员数作为居民收入与就业的度量指标。而对于房价则选择使用住宅商品房平均价格作为替代变量。 空气污染状况与废气排放状况均为潜变量，在本文中选取我国所使用的6种污染物浓度(AQI、PM 2.5 、PM 10 、SO 2 、CO、NO 2 以及O 3 )作为评判空气污染状况的指标；选取SO 2 、NO x 以及粉尘排放密度作为废气排放状况的指标，该指标采用各地区SO 2 、NO x 以及粉尘排放量与该地区土地面积的比值计算得出。除此之外，本文还加入了一些控制变量，例如：人均GDP与自然状况(平均气温、全年日照时数以及全年降水量)。 其中，除了空气污染状况数据来源于《中国空气质量在线监测分析平台》，其他数据均来自于2017年《中国统计年鉴》与《中国环境统计年鉴》。 为避免由于变量间存在量纲所带来的误差，本文首先对所收集的数据标准化处理；其次，对潜变量(空气污染状况与废气排放)进行单因子检验，使得测量方程中的因子载荷较大(具体见表2)。从表2中可以看出，SO 2 、O 3 的因子载荷值较小，因此剔除该两项指标。 Table 2 Single factor test of latent variable 指标 废气排放状况 空气污染状况 二氧化硫浓度 氮氧化物浓度 粉尘浓度 AQI PM 2.5 PM 10 SO 2 CO NO 2 O 3 因子载荷值 0.857 0.966 0.898 0.989 0.958 0.926 0.421 0.809 0.802 0.063 表2. 潜变量单因子检验
system: [SYSTEM]

user: 4. 实证分析 通过上述的分析，本文模型建立如下：首先，对房价具有直接影响的因素建立方程(模型1)；其次，在上述模型的基础上添加一阶间接影响，即对空气污染状况、居民收入与就业状况以及城市化进行解释(模型2)；然后，在模型2的中增添解释废气排放与产业结构的方程(模型3)；最后，对模型3进行修正优化，得到模型4 (具体见表2)。 Table 3 SEM model estimation result 因变量 自变量 观测指标 模型1 模型2 模型3 模型4 房价 工业化 第二产业增加值所占比重 0.845 *** 0.871 *** 0.88 *** 0.89 *** 城市化 市辖区人口比重 0.958 *** 0.964 *** 0.971 *** 0.975 *** 人口比重平方 1.032 *** 1.025 *** 1.019 *** 1.015 *** 市辖区人口密度 −0.113 −0.023 −0.034 人口密度平方 −0.058 0.008 −0.003 居民收入与就业状况 城镇居民人均收入 1.094 *** 0.957 *** 0.949 *** 0.935 *** 城镇就业人员 0.053 −0.025 −0.029 人均GDP 0.81 *** 0.925 *** 0.934 *** 0.948 *** 自然状况 年平均气温 0.916 *** 0.921 *** 0.922 *** 0.922 *** 全年日照时数 0.816 *** 0.823 *** 0.823 *** 0.823 *** 全年降水量 0.87 *** 0.859 *** 0.857 *** 0.856 *** 潜变量 空气污染状况 AQI −0.99 ** −0.99 *** −1.122 *** −0.99 *** PM 2.5 −0.958 *** −0.961 *** −0.99 *** −0.958 *** PM 10 −0.924 *** −0.93 *** −0.957 *** −0.926 *** CO −0.805 *** −0.819 *** −0.925 *** −0.81 *** NO 2 −0.803 *** −0.814 *** −0.808 *** −0.805 *** 空气污染状况 工业化 第二产业增加值所占的比重 0.942 * 0.651 * 0.535 * 城市化 市辖区人口比重 0.148 0.187 人口比重平方 0.158 0.197 人口密度 0.004 0.007 人口密度平方 0.001 −0.0006 产业结构 第三产业增加值与第二产业增加值的比例 1.213 0.825 潜变量 废气排放 二氧化硫浓度 0.291 *** 0.271 * 0.313 ** 氮氧化物浓度 1.304 *** 0.3191 * 0.372 ** 粉尘排放浓度 2.098 *** 0.285 * 0.328 ** 废气排放 工业化 第二产业增加值所占的比重 0.0959 城市化 市辖区人口比重 0.387 *** 0.457 *** 人口比重平方 −0.406 *** −0.476 *** 人口密度 0.0136 人口密度平方 0.0012 自然状况 年平均气温 −0.451 ** −0.483 ** 居民收入与就业状况 工业化 第二产业增加值所占的比重 −0.121 −0.074 城市化 市辖区人口比重 0.769 *** 0.833 *** 0.896 *** 人口比重平方 0.817 0.874 人口密度 −0.018 −0.0292 人口密度平方 0.006 −0.0026 产业结构 工业化 第二产业增加值所占的比重 −0.808 −0.797 * 城市化 产业结构 第三产业增加值与第二产业增加值的比例 0.578 *** 0.636 0.621 *** 优性判断指标 TLI 0.523 0.572 0.573 0.795 CFI 0.584 0.639 0.644 0.853 RMSEA 0.283 0.268 0.264 0.119 ACI 719.99 667.71 664.22 390.28 表3. SEM模型估计结果 从表3中的优性判断指标中可以看出模型4要优于其他3个模型，因此，本文针对模型4进行分析。 首先对影响房价的直接因素进行分析：1) 工业化对房价具有显著作用，且为正向的线性作用，该结果与诸多文献的结果保持一致。总体来说，在第二产业较为发达的省份，该区域的房价通常会偏高；而与之对应的是水平较低的省份，该区域的房价通常会偏低。2) 城市化中市区人口比重对省份房价具有正向显著作用，且其平方也显著，表明人口比重与房价间存在非线性关系，可能是由于随着城市人口的急剧增长，当超出某一阈值后，会导致人们愿意在市区购房的意愿降低，继而导致房价的下跌。但在控制其他变量不变后，市区人口密度及其平方项对房价均没有明显的作用。3) 居民状况(收入与人均GDP)对房价具有正向推动作用，这进一步说明随着人们生活水平的提高，物价(住房价格)也处于增长的现实。4) 地区所处自然环境(自然状况：降水、日照以及气温)对房价具有显著正向关系，表明人们自然条件优越的省份，房价也相对较高。5) 空气污染状况与房价负相关，表明人们在评估住房价格时已渐渐将当地空气状况列入考察范围内。 从间接影响的角度进行分析：1) 从空气污染状况方程、城市化方程以及产业结构方程来看，废气排放对空气污染状况具有直接作用，进而对房价产生影响；而工业化通过对产业结构产生影响，继而影响城市化，最终对房价产生作用。2) 从废气排放方程中可以看出，城市化中的人口比重对废气排放具有显著影响，表明城市化在对房价直接作用的同时，还通过废气排放对空气污染状况产生作用，进而对房价产生作用；而平均气温与废气排放间也存在正相关，主要是由于气温较低会增加煤炭的消耗，进而使得废气排放量增多。3) 从居民生活状况方程来看，工业化对居民收入的影响并不显著，相反城市化却与之正相关，表明城市化中的人口比重也通过提高居民生活状况，进而影响房价。 综上所述，城市化与工业化对房价的作用路径存在着差异，二者对房价的影响并不仅仅是直接作用也存在着一阶间接作用甚至是二阶间接作用。在城市化中，市区人口比重的影响力要远高于人口密度的作用力，具体的影响路径见图2。 图2. 城市化与工业化对房价的影响模式与路径
system: [SYSTEM]

user: 5. 结论 本文从结构性的视角出发，采用结构方程模型(SEM)分析了城市化与工业化对我国省份房价的作用路径与模式。结果表明： 1) 二者对省际房价的作用路径与模式存在着明显的差异。 2) 从直接作用路径来讲，城市化与工业化对房价均存在着直接作用，且二者的作用模式也表现的一致，均为正向关系。从一阶间接作用的角度来讲，工业化通过对空气污染状况的作用来影响房价的变化，城市化则是通过对居民生活状况(收入、人均GDP)的角度来影响房价。 3) 从二阶间接作用的角度来看，工业化不存在二阶作用，而城市化则会通过废气排放的途径对房价产生影响。除此之外，工业化还通过对产业结构调整的作用，来影响城市化，进而对房价产生作用；而平均气温也会通过对废气排放的影响来对房价产生二次作用。 4) 从城市化中人口比重对废气排放的作用是“倒U型”，进而对房价产生作用，市区人口比重对房价作用存在阈值(转折点为0.615)，也就是说，在一般情况下，当某省份的市区人口比重超出该阈值，将会引起其房价下降。因此，从调控房价的角度来看，应该将重心放在城市化的作用路径上。
system: [SYSTEM]

user: 摘 要 本文选取赣南地区18区县围屋数量、民间戏曲遗址数、民间小吃遗留种数等多项客家文化指标，采用赋值加权方法，得出各区县客家文化特征值，并结合莫兰指数对整个赣南地区区县尺度下客家文化空间相关性和空间分布进行评价。研究结果表明：赣南十八县客家文化西北地区分布较东南地区稀疏，整个靠近广东福建省界边缘的东南地区的客家文化指标分布较密集；赣南十八县客家文化空间相关性较高，西北部大部分地区客家文化值处在低水平集聚区；瑞金市处在显著HL区，形成文化极化效应；定南县在Lisa分布图的显著HH区，且周边地区客家文化水平较高。该研究结果可为保护和发展该地区客家文化提供相关针对性措施和决策建议。 关键词 :客家文化，空间分布，莫兰指数，文化值，赣南地区 Copyright © 2020 by author(s) 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. 引言 客家文化是历史遗留下来的文化瑰宝，表现在建筑、语言、饮食、生活等各个方面，部分文化堪称“活化石文化”，有较大的历史研究意义，能够带来一个广阔的旅游经济市场空间。客家文化具有时间分布长、空间分布广特点，发源于岭南地区，上古时期就已经有了一个雏形，秦始皇统一七国后，中国南北大门打开，客家文化得以大规模蔓延传播，到如今客家文化遍布中国各地，甚至越出国界，据不完全统计国外就有500万左右客家人定居海外，国内客家人密集分布在赣闽粤三省交界处。李凤英在对福建省非遗项目畲族和客家文化特征的综述中，以客家文化中的“上刀山、下火海”等习俗作为研究对象，为中国的文化多元化提出了发展建议 [ 1 ]。宋德剑以梅州客家文化产业为研究重点，认为在广东地区客家文化产业存在着巨大的发展空间，需从文化凝聚、政策维持、以及多元化筹资等多方面大力弘扬客家文化 [ 2 ]。赣南地区客家文化的传承方面，权亚楠以赣南客家人文精神为基础，提出了对乡村建设的建议 [ 3 ]。近年来，已有众多学者对客家文化进行了研究 [ 4 ]，但在客家文化的空间分布和相关性评价方面研究尚少。鉴于此，本文以赣南地区为例，采用赋值加权法 [ 5 ] 对当地客家文化进行评价，并结合莫兰指数和LISA散点图对赣南地区18个县区进行空间分布统筹性分析，得出赣南18县区尺度下的客家文化空间相关性和集聚性的关系，并给出不同层面的文化发展建议 [ 6 ] [ 7 ]。
system: [SYSTEM]

user: 2. 研究区概况和研究方法 2.1. 研究区概况 赣南地区位于江西省南部，与闽粤地区交界，赣江的发源地，是红色故都和革命摇篮，同时也是客家文化的传播区和遗存区。研究区赣州市面积39379.64 km 2 ，拥有18个区县 [ 8 ] [ 9 ] [ 10 ] [ 11 ] [ 12 ]。客家文化散布在赣州市各个地区，如客家围屋文化有龙南县的关西围屋、全南县雅溪围屋、安远县东生围屋。客家民间文学有南康市木根源，客家民间音乐有兴国县山歌、于都县唢呐《公婆吹》、会昌县踩鼓、崇义县竹洞畲族山歌等。客家民间舞蹈有定南县瑞狮、章贡区高跷、上犹县九狮拜象、大余县南安罗汉舞、赣县云灯等。客家传统戏剧有赣南采茶戏、信丰县手端木偶戏、宁都县中村摊戏、于都县江西提线木偶戏，兴国县端戏等。客家曲艺有南康市古文、于都县古文等。客家传统手艺有石城县砚制作技艺、瑞金市手工艺术模具雕刻铜模、崇义县豆腐制作技艺、全南县蓝巾帕制作、兴国县鱼香肉丝制作技艺等，客家民俗有龙南县围屋习俗、上犹县客家门匾习俗、兴国县三僚堪舆文化、宁都县石上曾坊桥帮灯、定南县定南客家哭嫁习俗、龙南县香火龙宁都县洛口南云竹篙火龙等。客家传统医药有寻乌县挑积、于都县敦本堂熊氏民间中医等。客家传统美术有章贡区客家竹雕、大余县核微雕技艺、兴国县篆刻等 [ 12 ] [ 13 ] [ 14 ] [ 15 ]。 2.2. 研究方法 本文将客家文化指标首先进行赋值加权分析，得出赣南18县的客家文化特征值。参照前人研究的基础上 [ 16 ] [ 17 ]，采用Moran’s I指数和Lisa分布图分析赣南地区的客家文化空间分布特征，该方法能够综合且较合理的揭示客家文化的空间聚散特征。 全局莫兰指数计算公式如下： I = ∑ i = 1 n ∑ i = 1 n ω i j ( x i − x ¯ ) ( x j − x ¯ ) S 2 ∑ i = 1 n ∑ j = 1 n ω i j (1) 其中I表示莫兰指数； x i 表示要素i处的值； x j 表示与i相邻要素j的值；n为要素的个数；S 2 表示x值和其均值 x ¯ 的方差； ω i j 为要素之间的空间权重。 局部莫兰指数的Z-score计算公式如下： Z ( I i ) = I i − E ( I i ) σ (2) 其中I i 表示局部莫兰指数； E ( I i ) 表示I i 的数学期望； σ 表示标均差。当 Z ( I i ) > 0时，表示要素值与邻近区域具有高相似的集聚特征，当 Z ( I i ) < 0时，代表要素值与邻近区域具有低相似的集聚特征。