密度、氮肥及硼肥对紫色辣椒凤紫运1号的产量效应分析
隋益虎, 胡能兵, 唐金宝
安徽科技学院 农学院,安徽 凤阳 233100

作者简介:隋益虎(1970—),男,安徽定远人,博士,教授,主要从事园艺遗传育种学、生物统计学的教学与研究。E-mail: suiyh@ahstu.edu.cn

摘要

采用三因素五水平二次回归正交旋转组合设计的方法,研究不同密度( x1)、氮肥( x2)及硼肥( x3)对凤紫运1号紫色辣椒产量的影响,结果表明,三个因素线性项对辣椒产量影响为正向极显著,效应大小顺序是 x2> x3> x1,二次项影响为负向极显著,顺序是x12>x32>x22, x1 x3互作效应正向极显著, x1的边际效应最大。当定植密度4 123株·667m-2,氮肥40.020 kg·667m-2,硼肥1.000 kg·667m-2时产量最高,为3 206.678 kg·667m-2。综合因子频数分析表明,95%置信度下凤紫运1号产量高于2 500 kg·667m-2的高产栽培的因子优化组合为:密度3 903~4 172株·667m-2,氮肥施量31.186~36.081 kg·667m-2,硼肥施量0.724~0.880 kg·667m-2

关键词: 辣椒; 密度; 氮肥; 硼肥; 产量
中图分类号:S641.3 文献标志码:A 文章编号:1004-1524(2018)03-0413-07
Effect of planting density, fertilizer level of nitrogen and boron on yield of purple pepper Fengziyun No.1
SUI Yihu, HU Nengbing, TANG Jinbao
College of Agriculture, Anhui Science and Technology University, Fengyang 233100, China
Abstract

In order to optimize the main cultivation factors that affect the production of pepper crops, the method of three-factor five-level quadratic orthogonal rotation design was used to study the effects of different density ( x1), nitrogen fertilizer level ( x2) and boron fertilizer level ( x3) on the yield of purple pepper Fengziyun No. 1. The results showed that under the experimental conditions, the linear effect of the three factors on the yield of pepper was positive and very significant, and x2 had the greatest impact on production, followed by x3, and x1 had the least impact. The effect of the quadratic term was negative, andx12 had the greatest impact on production, followed byx32, andx22 had the least impact. The interaction effect of x1 x3 was very significantly positive, the marginal effect of x1 was the largest. When the planting density was 4 123 plants per 667m2, the nitrogen fertilizer was 40.020 kg per 667m2 and the boron fertilizer was 1.000 kg per 667m2, the fresh pepper fruit weight was the highest of 3 206.678 kg per 667m2. The results of comprehensive factor analysis showed the optimized factor combination of high yield of Fengziyun No. 1 with a yield of more than 2 500 kg per 667m2 at 95% confidence was that density 3 903-4 172 plants per 667m2, nitrogen fertilizer 31.186-36.081 kg per 667m2, boron fertilizer 0.724-0.880 kg per 667m2.

Keyword: pepper; density; nitrogen fertilizer; boron fertilizer; yield

辣椒(Capsicum annuum L.)是一种世界性蔬菜作物, 在我国辣椒种植面积位居蔬菜作物之首[1]。目前广泛种植的辣椒为常规绿果(未成熟果)类型, 鲜有紫果类型。紫色辣椒相对绿色辣椒有较高的抗逆性, 且未成熟果富含的花青素具有保健等特点。安徽科技学院辣椒课题组经过十多年努力, 选育出了三系配套的紫色线椒新品种— — 凤紫运1号, 2013年已通过安徽省级审定/登记。该品种母本是通过种间远缘杂交创制的绿色辣椒雄性核质互作不育系(CMS) 1110A[2], 父本是含有高辣椒素的紫色辣椒纯系9007-2。F1代植株的子叶、真叶、茎、花和未成熟果均呈紫色, 尤其是未成熟果中含有较高的花青素; 同时, 具有抗病毒病和疫病, 耐高低温和强光等特点, 因此, 可有效减少农药的使用, 提高商品果实的安全性[3]

影响辣椒生产的农艺因子很多, 为了使不同品种辣椒的栽培产量或经济效益最大化, 前人在单一因子或多种因子组合等方面做了大量研究[4, 5, 6, 7, 8, 9, 10]。而本课题组根据供试土壤环境条件及辣椒新品种的生长发育特点, 试验选择密度、氮肥和硼肥为主要优化因子, 筛选最佳因子水平组合, 旨在为特色辣椒新品种凤紫运1号的高产栽培与大面积推广提供参考。

1 材料与方法
1.1 材料与实验条件

实验于2015年1— 9月在安徽科技学院蔬菜基地中心进行。实验地前茬为小白菜, 黏性黄棕壤土。0~20 cm 耕层含有机质1.21%, 全氮0.09%, 碱解氮78 mg· kg-1, 速效磷27 mg· kg-1, 速效钾152.5 mg· kg-1, 热水溶性硼0.28 mg· kg-1。供试辣椒品种为凤紫运1号紫色辣椒新品种, 供试肥料的氮肥为尿素, 微肥硼肥为四硼酸钠。

1.2 试验设计

采用三因素五水平二次回归正交旋转组合设计, 选取密度(x1)、氮肥(x2)、硼肥(x3)为自变量, 五水平分别为1.682、1、0、-1、-1.682, 以辣椒鲜果实产量(y)为目标函数。试验共设置23个处理组合, 小区(畦)长度2.75 m, 宽度1 m(带沟), 面积为2.75 m2。每畦双行, 行距0.5 m, 每穴单株种植, 试验区四周设有保护畦。小区和折合667 m2的密度水平、氮肥和硼肥施用量及其编码列于表1

表1 试验因子及其水平编码 Table 1 Coding of experimental factors and levels
1.3 田间管理与数据记录

2015年1月在高效节能温室内育苗, 4月9日进行小拱棚覆盖栽培, 16日根据设计量, 按小区均匀施用硼肥; 5月10日按小区均匀施用尿素, 并撤除小拱棚。其后的水分、中耕除草等管理均按常规生产进行并严格控制作业的一致性, 减少试验误差。每个小区的植株都按株编号, 6月3日进行第一次采摘, 以后每1周采摘1次, 直至8月30日止, 每个被采果实长度要求大于10 cm。每次采摘后立即用电子天平称重并按株分别记录鲜果质量, 最后对整个生长期每株鲜果质量累加求和, 再计算小区产量, 折算667m2产量。

2 结果与分析
2.1 辣椒产量的二次回归模型

根据每株辣椒的多次产量记录结果算得的小区与667m2产量列于表2。试验的三因素二次回归正交旋转组合设计结构矩阵及其试验结果(小区产量)部分计算过程列于表3

表2 二十三个处理组合的辣椒果实产量 Table 2 Pepper fruit yields of 23 treatment combinations
表3 三因素二次回归正交旋转组合设计结构矩阵及其试验结果计算 Table 3 Structure matrix of three-factor quadratic regression orthogonal rotational combination designs and processing of the experimental data

按二次回归正交旋转组合设计的分析方法[11], 利用Excel软件计算, 得小区产量(y, kg· 2.75m-2)依密度(x1)、氮肥(x2)和硼肥(x3)变化的函数模型为:

y=10.310 3+0.655 1x1+1.108 2x2+0.792 2x3+0.163 1x1x2+0.649 9x1x3+0.270 2x2x3-0.765 4 x12-0.405 5 x22-0.410 1 x32

2.2 回归模型的显著性检验

根据变异来源将全试验总平方和分解为回归、失拟和误差项, F测验表明:F1=失拟均方/误差均方=0.641 9, Pr=0.675 5, 失拟不显著, 说明回归方程的拟合度好, 主要影响因子已被考查。F2=回归均方/剩余均方=41.57, Pr=9.104× 10-6, 回归极显著, 说明试验指标产量依考查因子(x1x2x3)回归在总体上达到极显著水平。进一步将模型中回归项平方和再分解为各因子一次项、交互项及二次项, F测验表明:除交互项x1x2x2x3回归没有达到显著, 其余项均达极显著水平(表4)。

表4 辣椒产量三因素二次回归正交旋转组合试验方差分析 Table 4 Variance analysis of pepper yield based on three-factor quadratic orthogonal rotation combination design
2.3 单因素效应

由回归方程中已无量纲化处理和消除相关性的偏回归系数可判断, 三个因素线性项对辣椒产量影响为正向, 影响大小的顺序是x2> x3> x1, 二次项影响为负向, 顺序是 x12> x32> x22。应用降维法, 将上述回归模型中任意两个因素固定在零水平, 得到仅有其余某一个因素(xi)对产量(yi)的一元偏回归子模型:y1=10.310 3+0.655 1x1-0.765 4 x12; y2=10.310 3+1.108 2x2-0.405 5 x22; y3=10.310 3+0.792 2x3-0.410 1 x32。当x1x2x3的水平分别在0.427 9、1.366 5和0.965 9时, 可获得2 534.71、2 684.36和2 593.51 kg· 667m-2的最高鲜果实产量(图1)。

图1 三因素对辣椒产量的效应Fig.1 Effects of three factors on the yield of pepper

2.4 两因素互作效应

将其中一个因素固定在零水平, 求出其余两个因素与试验指标产量之间的二元回归子模型。通过此模型的有关计算值寻找出该两因素相互作用而引起辣椒果实产量变化的规律。本实验因x1x2x2x3均不显著, 因此只分析x1x3的交互效应。将回归方程中x2取零水平, 得到x1x3与辣椒果实产量y之间的二元回归子模型为:y1, 3=10.310 3+0.655 1x1+0.792 2x3+0.649 9x1x3-0.765 4 x12-0.410 1 x32

x1x3的各水平编码值代入方程中计算相应互作的辣椒果实鲜质量值y1, 3, 比较表明, x1x3存在正向互作效应, 不论种植密度在何水平下, 增施硼肥都能带来辣椒产量的增加; 随着密度增大, 硼肥施用水平亦需相应提高。但在不同密度水平下, 各硼肥水平的增产作用大小不同。低密度水平(-1~-1.682)配合零水平硼(0.5 kg· 667m-2), 高密度水平(1~1.682)配合高水平硼(1.0 kg· 667m-2)均能获得最佳互作效应(表5)。

表5 密度(x1)与硼肥(x3)对辣椒果实鲜质量的互作效应 Table 5 Interaction effect of planting density (x1) and boron fertilizer (x3) on the fresh pepper fruit weight
2.5 边际效应

为了研究辣椒果实产量随各因素水平增减的变化率, 对产量数学模型求一阶偏导数, 同时将其他2个变量分别固定在“ 0” 水平, 得到dy1/dx1=0.655 1-1.530 8x1, dy2/dx2=1.108 2-0.811 0x2, dy3/dx3=0.792 2-0.820 2x3, 可见密度(x1)、氮肥(x2)和硼肥(x3)的边际产量均有明显变化, 其中密度(x1)的边际效应最大; 当各变量分别取较低水平时, 增产效应较大, 反之则较小; 当三因素超出适宜范围即越过产量最佳点后, 边际产量的效应降低, 并最终为负值。

2.6 模型的综合统计选优

在本研究约束范围x∈ [-1.682, 1.682]内, 辣椒果实的最高产量为3 206.678 kg· 667m-2, 其处理因子水平组合为x1=1, x2=1.682, x3=1.682(即密度4 123株· 667m-2, 氮肥40.020 kg· 667m-2, 硼肥1.000 kg· 667m-2)。同时对125个处理组合进行产量模拟得到超过2 500 kg· 667m-2产量的组合有24个, 说明凤紫运1号紫色辣椒的稳产性与高产性均较好。

进一步对综合因子进行频数统计分析, 获得高产栽培的优化区间(95%置信度), 表明在本实验条件下, 按常规技术标准进行管理, 凤紫运1号产量在2 500 kg· 667m-2以上的优化组合模式为:密度为3 903~4 172 株· 667m-2, 氮肥施量为31.186~36.081 kg· 667m-2, 硼肥施量为0.724~0.880 kg· 667m-2(表6)。

表6 凤紫运1号产量高于2 500 kg· 667m-2的综合因子水平频数分析 Table 6 Comprehensive frequency analysis of different factor levels combination on fresh pepper fruit weight of Fengziyun No.1 higher than 2 500 kg· 667m-2
3 结论与讨论

在本实验的土壤、管理和设计的因子水平条件下, 紫色辣椒凤紫运1号的果实产量依密度(x1)、氮肥(x2)和硼肥(x3)三因子的回归极显著, 回归方程可用于预测产量并指导高产栽培。本研究中, 当处理组合密度为4 123株· 667m-2 (x1=1), 氮肥40.020 kg· 667m-2 (x2=1.682), 硼肥1.000 kg· 667m-2 (x3=1.682)时, 辣椒果实的产量最高, 达3 206.678 kg· 667m-2

高产原因分析表明, 除三因子一次项回归正向极显著外, x1x3互作亦极显著, x2x3存在弱正向互作(P=0.069 9), 这可能与充足的单株营养状况以及硼肥促进授粉受精有关。有研究[5, 12, 13]表明, 种植密度通过影响通风、透光从而影响辣椒株高、分枝数、结果数量和单果质量等, 合理密植是辣椒高产和高效栽培的重要农艺措施之一。氮是辣椒必需的大量元素, 能够改善其多项生理功能[14], 氮的作用大于磷、钾[15]。李子双等[16]研究表明, 氮肥对产量的影响大于磷肥与硅钙肥, 能改善辣椒吸收营养元素等功能, 但应适量施用, 过多施用反而会造成产量降低、品质下降等现象。硼是植物的必需微量元素, 一些学者研究证实硼营养对油菜[17]、大豆[18]及辣椒[8]等植物的生长发育及产量均有较大的影响。硼参与细胞生长许多代谢过程[19, 20], 尤其在花粉母细胞减数分裂敏感阶段, 缺硼引起雄性不育[21, 22], 最终导致结实率下降、落果、果实畸形, 严重降低了产量[23]。凤紫运1号是三系配套的杂种一代, 增施硼肥可能有助于恢复基因在F1中发挥作用, 改善其花粉粒活性从而提高产量。此外, 前人还证实硼能促进某些作物对大量元素氮、钾、磷的吸收[24, 25]以及对中、微量元素Ca2+、Mg2+、Zn2+、Cu2+的吸收[26], 从而平衡这些矿质营养使植株生长良好。本研究中由于氮肥(x2)和硼肥(x3)均取了最高水平(1.682), 所以实际生产中随着两者施用水平的提高, 还有更高产的可能。

三因子综合试验模拟得到了高产(≥ 2 500 kg· 667m-2)处理组合24个, 占总组合数的比例达19.2%, 说明凤紫运1号在生产上具有较好的稳产、高产性。在95%置信度下, 凤紫运1号高产优化农艺措施模式为:密度3 903~4 172 株· 667m-2, 氮肥施量31.186~36.081 kg· 667m-2, 硼肥施量0.724~0.880 kg· 667m-2。此模式也为类似的果数与果质量并重型辣椒杂交品种的高产栽培提供了借鉴。

The authors have declared that no competing interests exist.

参考文献
[1] 王立浩, 张正海, 曹亚从, . “十二五”我国辣椒遗传育种研究进展及其展望[J]. 中国蔬菜, 2016 (1): 1-7.
WANG L H, ZHANG Z H, CAO Y C, et al. Progress and prospect of genetics and breeding of pepper in China during the twelfth five-year plan period[J]. China Vegetables, 2016 (1): 1-7. (in Chinese) [本文引用:1]
[2] 隋益虎, 张子学, 胡能兵, . 辣椒雄性不育系1110A的选育[J]. 安徽科技学院学报, 2013, 27(2): 14-17.
SUI Y H, ZHANG Z X, HU N B, et al. Breeding of a new CMS capsicum line 1110A[J]. Journal of Anhui Science and Technology University, 2013, 27(2): 14-17. (in Chinese with English abstract) [本文引用:1]
[3] 胡能兵, 隋益虎, 张子学, . 紫色辣椒色素性状观察及其叶片色素提取体系优化[J]. 西北植物学报, 2012, 32(12): 2450-2456.
HU N B, SUI Y H, ZHANG Z X, et al. Observation of pigment in purple pepper and extraction system optimization of leaf purple pigment[J]. Acta Botanica Boreali-Occidentalia Sinica, 2012, 32(12): 2450-2456. (in Chinese with English abstract) [本文引用:1]
[4] 李德斌. 不同栽培密度对辣椒产量和经济性状的影响[J]. 现代农业科技, 2013 (16): 74.
LI D B. Effect of different planting density on fruit yield and economic characters of pepper[J]. Modern Agricultural Sciences and Technology, 2013(16): 74. (in Chinese) [本文引用:1]
[5] 肖莉, 郑元红, 王慧, . 不同栽培方式及种植密度对辣椒产量的影响[J]. 长江蔬菜, 2011 (4): 39-42.
XIAO L, ZHENG Y H, WANG H, et al. Effect of planting method and density on yield of hot pepper[J]. Journal of Changjiang Vegetables, 2011(4): 39-42. (in Chinese with English abstract) [本文引用:2]
[6] 吕长山, 王金玲, 于广建, . 氮肥对辣椒果实品质及产量的影响[J]. 东北农业大学学报, 2005, 36(4): 448-450.
LV C S, WANG J L, YU G J, et al. Effect of nitrogen fertilizer on quality and yield of pepper fruits[J]. Journal of Northeast Agricultural University, 2005, 36(4): 448-450. (in Chinese with English abstract) [本文引用:1]
[7] 隋益虎, 汪建飞, 张子学, . CO2气肥、N肥及B肥对辣椒的施用效应[J]. 安徽农业技术师范学院学报, 1999, 13(1): 32-37.
SUI Y H, WANG J F, ZHANG Z X, et al. Effect of fertilizer CO2, N and B on pepper[J]. Journal of Anhui Agrotechnical Teachers College, 1999, 13(1): 32-37. (in Chinese with English abstract) [本文引用:1]
[8] 左明玉. 硼肥用量及施用方式对辣椒产量、效益的影响[J]. 安徽农业科学, 2014, 42(9): 2619-2621.
ZUO M Y. Effects of boron fertilizer amount and application mode on pepper yield and benefits[J]. Journal of Anhui Agricultural Sciences, 2014, 42(9): 2619-2621. (in Chinese with English abstract) [本文引用:2]
[9] 曾长立, 康六生. 植物生长调节剂、氮肥与密度配伍对辣椒产量及品质的影响[J]. 江西农业大学学报, 2009, 31(4): 644-649.
ZENG C L, KANG L S. The effect of composite method of plant growth regulators, nitrogen fertilizer and planting density on the yield and quality of hot pepper[J]. Acta Agriculturae Universitatis Jiangxiensis, 2009, 31(4): 644-649. (in Chinese with English abstract) [本文引用:1]
[10] 张福建, 陈昱, 吴超群, . 外源脂肪酸对辣椒生长及根际土壤环境的影响[J]. 浙江农业学报, 2017, 29(5): 760-766.
ZHANG F J, CHEN Y, WU C Q, et al. Effects of exogenous fatty acids on growth and rhizosphere soil of pepper plants[J]. Acta Agriculturae Zhejiangensis, 2017, 29(5): 760-766. (in Chinese with English abstract) [本文引用:1]
[11] 张禄捷, 李荣, 姜子涛. 二次回归正交旋转组合设计优化茼蒿籽中总黄酮的微波提取工艺[J]. 浙江农业学报, 2015, 27(7): 1233-1238.
ZHANG L J, LI R, JIANG Z T. Optimization of microwave-assisted extraction technology for flavonoids from Chrysanthemum coronarium L. seeds by quadratic regression orthogonal rotary[J]. Acta Agriculturae Zhejiangensis, 2015, 27(7): 1233-1238. (in Chinese with English abstract) [本文引用:1]
[12] 赵贞祥, 杨永岗, 张二喜, . 旱地辣椒栽培中密度、氮、磷及钾肥因子的优化[J]. 土壤, 2013, 45(4): 628-632.
ZHAO Z X, YANG Y G, ZHANG E X, et al. Optimization of transplanting density, nitrogen, phosphorus and potassium fertilization on yield of pepper in drought area[J]. Soils, 2013, 45(4): 628-632. (in Chinese with English abstract) [本文引用:1]
[13] 岳振平. 不同种植密度对早春大棚辣椒产量和产值的影响[J]. 北方园艺, 2014 (5): 51-52.
YUE Z P. Effect of different planting density on fruit yield and output value of pepper in early spring[J]. Northern Horticulture, 2014 (5): 51-52. (in Chinese with English abstract) [本文引用:1]
[14] 付小松, 何建文, 韩世玉. 水分和氮素胁迫对辣椒产量形成及抗逆性评价[J]. 长江蔬菜, 2015 (4): 25-28.
FU X S, HE J W, HAN S Y. Effects of water and nitrogen stress on pepper yield and stress resistance evaluation[J]. Journal of Changjiang Vegetables, 2015(4): 25-28. (in Chinese with English abstract) [本文引用:1]
[15] 黄科, 刘明月, 吴秋云. 氮磷钾施用量与辣椒产量的相关性研究[J]. 江西农业大学学报, 2002, 24(6): 772-776.
HUANG K, LIU M Y, WU Q Y. A study on the correlation between combined application of N P K and the yield of hot pepper[J]. Acta Agriculturae Universitatis Jiangxiensis, 2002, 24(6): 772-776. (in Chinese with English abstract) [本文引用:1]
[16] 李子双, 王薇, 张世文, . 氮磷与硅钙肥配施对辣椒产量和品质的影响[J]. 植物营养与肥料学报, 2015, 21(2): 458-466.
LI Z S, WANG W, ZHANG S W, et al. Effect of nitrogen, phosphorus and silicon calcium fertilizer on yield and quality of pepper[J]. Journal of Plant Nutrition and Fertilizer, 2015, 21(2): 458-466. (in Chinese with English abstract) [本文引用:1]
[17] 张秀省, 沈振国, 沈康. 硼对油菜花器官发育和结实性的影响[J]. 土壤学报, 1994, 31(2): 146-152.
ZHANG X S, SHEN Z G, SHEN K. Effect of boron on floral organs development and seed-setting of rapeseed ( Brassica napus L. )[J]. Acta Pedologica Sinica, 1994, 31(2): 146-152. (in Chinese with English abstract) [本文引用:1]
[18] 刘鹏, 杨玉爱. 钼硼对大豆氮磷钾吸收及产量的影响[J]. 中国油料作物学报, 2000, 22(3): 57-60.
LIU P, YANG Y A. Effect of molybdenum or boron application on N, P, K absorption and yield in soybean[J]. Chinese Journal of Oil Crop Sciences, 2000, 22(3): 57-60. (in Chinese with English abstract) [本文引用:1]
[19] MIGUEL R N, LUIS P, REGUERA M, et al. Developmentally regulated membrane glycoprotein sharing antigenicity with rhamnogalacturonan Ⅱare not detected in nodulated boron deficient Pisum sativum[J]. Plant, Cell and Environment, 2007, 30(11): 1436-1443. [本文引用:1]
[20] CAMACHO-CRISTOBAL J J, REXACH J, GONZALEZ-FONTES A. Boron in plants: deficiency and toxicity[J]. Journal of Integrative Plant Biology, 2008, 50(10): 1247-1255. [本文引用:1]
[21] RAWSON H M, SUBEDI K D. Sterility in wheat in subtropical Asia: extent, causes and solutions[C]// Proceedings No 73 (Workshop held at Lumle Agricultural Research Centre, Nepal, 18-21 September 1995. ) Canberra, ACIAR, 1996: 1-154. [本文引用:1]
[22] HUANG L B, PANT J, DELL B, et al. Effects of boron deficiency on anther development and floret fertility in wheat ( Triticum aestivum L. ‘Wilgoyne’)[J]. Annals of Botany, 2000, 85(4): 493-500. [本文引用:1]
[23] 祖艳群, 林克惠. 硼在植物体中的作用及对作物产量和品质的影响[J]. 云南农业大学学报, 2000, 15(4): 359-364.
ZU Y Q, LIN K H. The role of boron in plants and its effect on the yield and quality of crops[J]. Journal of Yunnan Agricultural University, 2000, 15(4): 359-364. (in Chinese with English abstract) [本文引用:1]
[24] CAMACHO-CRISTOBAL J J, GONZALEZ-FONTES A. Boron deficiency causes a drastic decrease in nitrate content and nitrate reductase activity, and increases the content of carbohydrates in leaves from tobacco plants[J]. Planta, 1999, 209(4): 528-536. [本文引用:1]
[25] MARSCHER H, MARSCHER P. Mineral nutrition of higher plants[M]. 3rd edition. New York: Academic press, 2012: 158-483. [本文引用:1]
[26] SHAABAN M M, El-FOULY M M, ABDEL-MAGUID A A. Zinc-boron relationship in wheat plants grown under low or high levels of calcium carbonate in the soil[J]. Pakistan Journal of Biological Sciences, 2004, 7(4): 633-639. [本文引用:1]