Parameter Optimization of the Harvest Method in the Standardized Hedge Cultivation Mode of Lycium barbarum Using Response Surface Methodology

: In a previously published study, to optimize the vibrating and comb-brushing harvesting, the main factors and their parameter values were obtained based on the FEM and RSM. However, the study was based on the extensive cultivation mode which need to be improve. To realize the mechanization of the harvesting of Lycium barbarum L., as well as to face the standardized hedge cultivation mode, a vibrating and comb-brushing harvester machine was designed, which was primarily composed of an execution system, a motion system, and a control system. The mathematical model between the harvest effect index and the operation parameters was established based on response surface methodology (RSM). The effects of various parameters on the harvest index were analyzed, and the best parameter combination was determined: a vibration frequency of 38.73 Hz, a brush speed of 14.21 mm/s, and an insertion depth of 26.07 mm. The ﬁeld experiment showed that the harvesting rate of ripe fruit was 83.65%, the harvesting rate of unripe fruit was 7.22%, the damage rate of the ripe fruit was 11.49%, and the comprehensive picking index was 87.85. The ﬁndings provided a reference for the development of L. barbarum harvesting mechanization in a standardized hedge cultivation mode. models of the harvesting rate of the ripe fruit, the harvesting rate of the unripe fruit, the damage rate of the ripe fruit, and the comprehensive picking index were established. Furthermore, an optimized parameter combination was determined based on RSM with a vibration frequency of 38.73 Hz, a brush speed of 14.21 mm/s, and an insertion depth of 26.07 mm. The ﬁeld experiment showed that the harvesting rate of the ripe fruit was 83.65%, the harvesting rate of the unripe fruit was 7.22%, the damage rate of the ripe fruit was 11.49%, and the comprehensive picking index was 87.85. These results provide a design basis for future research on the vibrating and comb-brushing harvest method of L. barbarum . However, the harvesting effect still needs to be improved. In further studies, according to the growth characteristics and cultivation mode, the mechanism of vibrating and comb-brushing harvesting should be studied. The harvesting time, area, and location could be optimized to increase the harvesting rate of the ripe fruit, reduce the harvesting rate of the unripe fruit, and the damage rate of the ripe fruit to achieve a better harvesting effect.


Introduction
Lycium barbarum L. (L. barbarum) is an infinite inflorescence deciduous shrub plant belonging to Lycium barbarum in the eggplant family [1][2][3][4]. Its ripe fruit are rich in polysaccharides, carotene, and other important physiological active ingredients [5][6][7][8]. It has the functions of protecting the liver, kidney, eyesight, and nerves as well as fighting off cancer [9,10]. At present, the manual harvesting method of L. barbarum has the problems of low efficiency and a high cost [11,12]. With the increase in planting area and a reduction in the labor force, the development of mechanized harvesting technology is necessary [11,[13][14][15][16].
L. barbarum is a species that is distributed worldwide. However, only China has cultivated and exploited it to a large scale, and, as such, its harvest research is mainly concentrated in China [15,[17][18][19][20][21][22]. The current harvest research of L. barbarum mainly focuses on designing different mechanical structures and obtaining the best harvest parameters [20,23,24]. However, there are still several problems such as diverse and inconsistent evaluation indexes, which make it difficult to quantitatively evaluate the comprehensive harvest effect [1,10,24]. Among the various harvesting methods of L. barbarum, the vibrating harvester has a high efficiency but inflicts great damage. By contrast, the comb-brushing

Standardized Hedge Cultivation Mode
In this study, we investigated the cultivation mode of L. barbarum in Ningxia Zhengqihong Industry Development Co., Ltd., Guyuan City, Ningxia Hui Autonomous Region (36 • 17 32.9" N, 106 • 6 41.5" E). The main variety planted was Keqi No. 2, which exhibits good growth conditions, more fruit, no diseases or insect pests, and no obvious defects. For the convenience of harvesting, the shrubs were artificially pruned into a standardized hedge cultivation mode (Figure 1) when the shrubs were 2-3 years old. The data of the cultivation mode are shown in Table 1.  The structure of the vibrating and comb-brushing harvester is shown in Figure 2. structure consisted of an executive system, a motion system, and a control system. execution system included a vibrating comb-brush bar, a vibration source, and a sh absorber connector. The motion system consisted of two cross ball screw sliding tab where the sliding tables used 42BYG stepper motors as a power source to drive the ex utive system movement within two degrees of freedom. The control system include battery, a CNC SHIELD V3 engraving machine expansion board, an A4988 driver boa and an Arduino UNO R3 MCU (manufactured by Shenzhen Two Trees Technology C Ltd., Shenzhen, China).  The structure of the vibrating and comb-brushing harvester is shown in Figure 2. The structure consisted of an executive system, a motion system, and a control system. The execution system included a vibrating comb-brush bar, a vibration source, and a shock absorber connector. The motion system consisted of two cross ball screw sliding tables where the sliding tables used 42BYG stepper motors as a power source to drive the executive system movement within two degrees of freedom. The control system included a battery, a CNC SHIELD V3 engraving machine expansion board, an A4988 driver board, and an Arduino UNO R3 MCU (manufactured by Shenzhen Two Trees Technology Co., Ltd., Shenzhen, China).

Operating Principle
When the harvester worked, the control system managed the insertion depth and brush speed of the execution system through the motion system. The executive system moved horizontally to insert the vibrating comb-brush bar into the branches. The executive system then drove the branches and fruit to vibrate, and then moved vertically downward at a certain speed to achieve the combined vibrating and comb-brushing effect.

Execution System
The structure of the execution system is shown in Figure 3. The vibration was transmitted to the branches and fruit through the flexible vibrating comb-brush bar. With the downward movement of the execution system, the fruit could be harvested under the action of vibration and comb-brushing.
where the sliding tables used 42BYG stepper motors as a power source to d utive system movement within two degrees of freedom. The control system battery, a CNC SHIELD V3 engraving machine expansion board, an A4988 and an Arduino UNO R3 MCU (manufactured by Shenzhen Two Trees Tec Ltd., Shenzhen, China).

Operating Principle
When the harvester worked, the control system managed the insertion depth a brush speed of the execution system through the motion system. The executive syste moved horizontally to insert the vibrating comb-brush bar into the branches. The exec tive system then drove the branches and fruit to vibrate, and then moved vertically dow ward at a certain speed to achieve the combined vibrating and comb-brushing effect.

Execution System
The structure of the execution system is shown in Figure 3. The vibration was tran mitted to the branches and fruit through the flexible vibrating comb-brush bar. With t downward movement of the execution system, the fruit could be harvested under the a tion of vibration and comb-brushing. The execution system provided vibration through the centrifugal force generated the rotating eccentric block. The exciting force F, vibration amplitude S, and excitati acceleration a could be described by the following: The execution system provided vibration through the centrifugal force generated by the rotating eccentric block. The exciting force F, vibration amplitude S, and excitation acceleration a could be described by the following: where F is the exciting force of the vibrating comb-brush bar (N), m is the mass of the eccentric block (kg), r is the distance between the centroid of the eccentric block and the center of rotation (m), ω is the angular frequency of the motor rotation (rad/s), S is the vibration amplitude of the vibrating substance (mm), G is the mass of the vibrating substance (kg), and a is the excitation acceleration of the vibrating substance (N/kg).

Motion System
The motion system consisted of two ball screw slides mounted vertically. The screw slides were driven by a stepper motor. The speed of the stepper motor could be adjusted by the pulse of the MCU. The movement length L and rotational speed V r could be obtained as follows: where L is the movement length of the screw slides (mm), h is the lead of the screw (mm), R is the rotational angle of the motor ( • ), V r is the rotational speed of the stepper motor (r/s), p is the step angle after subdivision ( • ), and f is the pulse frequency (Hz).
As the lead of the screw was 8 mm, the speed of the sliding table V l was obtained as follows: where V l is the speed of the screw slide table (mm/s).

Control System
The circuit diagram of the control system is shown in Figure 4. The rated speed of the vibration motor in the execution system was 3200 rpm and the speed was controlled by a PWM DC governor (manufactured by Shenzhen Quanqiuyi Technology Co., Ltd., Shenzhen, China). The screw slide platform was driven by a 42BYG stepper motor (manufactured by Shenzhen Global Yi Technology Co., Ltd., Shenzhen, China) and powered by a DC-24120 large-capacity lithium polymer battery (manufactured by Shenzhen Dipuwei Technology Co., Ltd., Shenzhen, China). The motor was controlled by an Arduino MCU, a CNC SHIELD V3, and an A4988 (manufactured by Shenzhen Two Trees Technology Co., Ltd., Shenzhen, China). The development environment was Arduino IDE. The Arduino MCU adopted an acceleration and deceleration calculation method. The moving length or speed of the sliding table could be controlled by changing the rotational angle or pulse frequency of the stepping motor.
Horticulturae 2022, 7, x FOR PEER REVIEW 5 of 14 The motion system consisted of two ball screw slides mounted vertically. The screw slides were driven by a stepper motor. The speed of the stepper motor could be adjusted by the pulse of the MCU. The movement length L and rotational speed Vr could be obtained as follows: where L is the movement length of the screw slides (mm), h is the lead of the screw (mm), R is the rotational angle of the motor (°), Vr is the rotational speed of the stepper motor (r/s), p is the step angle after subdivision (°), and f is the pulse frequency (Hz).
As the lead of the screw was 8 mm, the speed of the sliding table Vl was obtained as follows: where Vl is the speed of the screw slide table (mm/s).

Control System
The circuit diagram of the control system is shown in Figure 4. The rated speed of the vibration motor in the execution system was 3200 rpm and the speed was controlled by a PWM DC governor (manufactured by Shenzhen Quanqiuyi Technology Co., Ltd., Shenzhen, China). The screw slide platform was driven by a 42BYG stepper motor (manufactured by Shenzhen Global Yi Technology Co., Ltd., Shenzhen, China) and powered by a DC-24120 large-capacity lithium polymer battery (manufactured by Shenzhen Dipuwei Technology Co., Ltd., Shenzhen, China). The motor was controlled by an Arduino MCU, a CNC SHIELD V3, and an A4988 (manufactured by Shenzhen Two Trees Technology Co., Ltd., Shenzhen, China). The development environment was Arduino IDE. The Arduino MCU adopted an acceleration and deceleration calculation method. The moving length or speed of the sliding table could be controlled by changing the rotational angle or pulse frequency of the stepping motor.

Performance Experiment
The main goal of the vibrating and comb-brushing harvest method is to achieve ripe fruit harvesting, reduce ripe fruit damage, and reduce unripe fruit harvesting as well as

Performance Experiment
The main goal of the vibrating and comb-brushing harvest method is to achieve ripe fruit harvesting, reduce ripe fruit damage, and reduce unripe fruit harvesting as well as flower and leaf shedding. Therefore, based on the research of [1,[10][11][12][23][24][25], the harvesting rate of ripe fruit (I 1 ), the harvesting rate of unripe fruit (I 2 ), and the damage rate of ripe fruit (I 3 ) were selected as the evaluation indexes for the experiments and calculated according to the following: where n 1 is the amount of ripe fruit harvested, n 2 is the amount of ripe fruit unharvested, n 3 is the amount of unripe fruit harvested, n 4 is the amount of unripe fruit unharvested, and n 5 is the amount of damaged ripe fruit harvested.
To comprehensively evaluate the harvesting effect, weight distribution was carried out on the measurement results to obtain a single comprehensive picking index (I), which could be described by the following: where I 2s is the score of unripe fruit harvested and I 3s is the score of damaged ripe fruit harvested. Through a principal analysis and literature review, it could be determined that the main factors affecting the performance of the harvester were vibration frequency, brush speed, and insertion depth of the comb-brush bar into the branch. The value range of each factor was determined by a preliminary test: the vibration frequency (X 1 ) was 20~40 Hz, the brush speed (X 2 ) was 10~40 mm/s, and the insertion depth (X 3 ) was 10~50 mm. The insertion depth and brush speed were adjusted by the motion system and the vibration frequency was measured by a VICTOR 6236P rotation speed/line speed meter (manufactured by Shenyang Zizun Science and Technology Co., Ltd., Shenyang, China).
A three-factor and three-level quadratic orthogonal rotation combination test was used in the experiment conducted for the present study. The factor codes are shown in Table 2 and the test schemes and results are presented in Table 3. A total of 17 groups of tests were carried out; a group of 5 tests was carried out and the average value of the 5 test results was taken as the test result of the group [27][28][29]. Design-Expert 12 software was used to design the experimental scheme and analyze the results.

Regression Analysis
After an analysis from the implementation of Design-Expert 12 software, a regression mathematical model was obtained that took the harvesting rate of the ripe fruit as the response function and the coding value of each factor as the independent variable.
A variance analysis (ANOVA) was performed on the harvesting rate of the ripe fruit, as shown in Table 4. The results showed that the regression model was p = 0.0003 (<0.05), indicating that the model had a statistical significance. The factors X 1 , X 3 , X 1 X 2 , and X 3 2 were significant (p < 0.05) whereas the other factors were not. The lack of fit was p = 0.5526 (>0.05), indicating that none of the factors were irrelevant. As above, the regression mathematical model of the harvesting rate of the unripe fruit using the codes of the factor as the independent variable was as follows: Horticulturae 2022, 8, 308 8 of 14 The ANOVA was performed on the harvesting rate of the unripe fruit, as shown in Table 5. The results showed that the regression model was p = 0.0115 (<0.05), indicating that the model had a statistical significance. The factors X 1 , X 2 , X 3 , and X 1 X 2 were significant (p < 0.05) whereas the other factors were not. The lack of fit was p = 0.1848 (>0.05), indicating that none of the factors were irrelevant. As above, the regression mathematical model of the damage rate of the ripe fruit using the codes of the factor as the independent variable was as follows: The ANOVA was performed on the damage rate of the ripe fruit, as shown in Table 6. The results showed that the regression model was p = 0.0147 (<0.05), indicating that the model had a statistical significance. The factors X 1 , X 2 , and X 3 were significant (p < 0.05) whereas the other factors were not. The lack of fit was p = 0.1799 (>0.05), indicating that none of the factors were irrelevant. As above, the regression mathematical model of the comprehensive picking index using the codes of the factor as the independent variable was as follows: The ANOVA was performed on the comprehensive picking index, as shown in Table 7. The results showed that the regression model was p = 0.0155 (<0.05), indicating that the model had a statistical significance. The factors X 1 , X 2 , and X 3 were significant (p < 0.05) whereas the other factors were not. The lack of fit was p = 0.3131 (>0.05), indicating that none of the factors were irrelevant.

Response Surface Analysis
The response surface methodology was used to analyze the influence of various factors on the harvesting rate of the ripe fruit and the response surface of the regression equation is presented in Figure 5. It could be seen from Equation (13) and Table 3 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and brush speed. The interaction between the vibration frequency and the brush speed was significant. As shown in Figure 5a, with an increase in the vibration frequency, the impact times and impact force on the fruit increased and the harvesting rate of the ripe fruit gradually increased. However, the vibration frequency could not continue to increase. It was found in the actual test that when the vibration frequency exceeded 50 Hz, the amplitude noticeably decreased and tended to disappear with an increase in the frequency. Due to inertia, the vibrating comb-brush bar was unable to keep up with the rapidly changing excitation force so the system no longer vibrated. As shown in Figure 5b, with a decrease in the insertion depth, the harvesting rate of the ripe fruit gradually increased. This was because the vibrating comb-brush bar could be regarded as a simple supported beam. Therefore, the further away from the fixed end of the force, the greater the amplitude. As shown in Figure 5c, with a decrease in the brush speed, the harvesting rate of the ripe fruit rapidly increased at first but then slowly decreased.
The response surface methodology was used to analyze the influence of various factors on the harvesting rate of the unripe fruit and the response surface of the regression equation is presented in Figure 6. It could be seen from Equation (14) and Table 4 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and brush speed. The interaction between the vibration frequency and the brush speed was significant. As shown in Figure 6a, with an increase in the vibration frequency, the harvesting rate of the unripe fruit gradually increased. As shown in Figure 6b, with a decrease in the insertion depth, the harvesting rate of the unripe fruit gradually increased. As shown in Figure 6c, with a decrease in the brush speed, the harvesting rate of the unripe fruit gradually increased. This was because the unripe fruit were mostly concentrated at the end of branches; an increase in the vibration frequency and comb-brushing time could lead to the fracture of the end branches.
quency. Due to inertia, the vibrating comb-brush bar was unable to keep up with the rapidly changing excitation force so the system no longer vibrated. As shown in Figure 5b, with a decrease in the insertion depth, the harvesting rate of the ripe fruit gradually increased. This was because the vibrating comb-brush bar could be regarded as a simple supported beam. Therefore, the further away from the fixed end of the force, the greater the amplitude. As shown in Figure 5c, with a decrease in the brush speed, the harvesting rate of the ripe fruit rapidly increased at first but then slowly decreased.  Response surface and contour plots of the harvesting rate of ripe fruit affected by the vibration frequency X 1 and brush speed X 2 (a), the vibration frequency X 1 and insertion depth X 3 (b), and the brush speed X 2 and insertion depth X 3 (c).
Horticulturae 2022, 7, x FOR PEER REVIEW 10 of 14 The response surface methodology was used to analyze the influence of various factors on the harvesting rate of the unripe fruit and the response surface of the regression equation is presented in Figure 6. It could be seen from Equation (14) and Table 4 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and brush speed. The interaction between the vibration frequency and the brush speed was significant. As shown in Figure 6a, with an increase in the vibration frequency, the harvesting rate of the unripe fruit gradually increased. As shown in Figure  6b, with a decrease in the insertion depth, the harvesting rate of the unripe fruit gradually increased. As shown in Figure 6c, with a decrease in the brush speed, the harvesting rate of the unripe fruit gradually increased. This was because the unripe fruit were mostly concentrated at the end of branches; an increase in the vibration frequency and combbrushing time could lead to the fracture of the end branches. The response surface methodology was used to analyze the influence of various factors on the damage rate of the ripe fruit and the response surface of the regression equation is presented in Figure 7. It could be seen from Equation (15) and Table 5 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and the brush speed. The interaction effect of each factor was not significant. As shown in Figure 7a, with an increase in the vibration frequency, the damage rate of the ripe fruit gradually increased. As shown in Figure 7b, with a decrease in the insertion depth, the damage rate of the ripe fruit gradually increased. As shown in Figure 7c, with a decrease in the brush speed, the damage rate of the ripe fruit gradually increased. This was because the number of hits to the ripe fruit increased.  Response surface and contour plots of the harvesting rate of unripe fruit affected by the vibration frequency X 1 and brush speed X 2 (a), the vibration frequency X 1 and insertion depth X 3 (b), and the brush speed X 2 and insertion depth X 3 (c).
The response surface methodology was used to analyze the influence of various factors on the damage rate of the ripe fruit and the response surface of the regression equation is presented in Figure 7. It could be seen from Equation (15) and Table 5 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and the brush speed. The interaction effect of each factor was not significant. As shown in Figure 7a, with an increase in the vibration frequency, the damage rate of the ripe fruit gradually increased. As shown in Figure 7b, with a decrease in the insertion depth, the damage rate of the ripe fruit gradually increased. As shown in Figure 7c, with a decrease in the brush speed, the damage rate of the ripe fruit gradually increased. This was because the number of hits to the ripe fruit increased.
The response surface methodology was used to analyze the influence of various factors on the comprehensive picking index and the response surface of the regression equation is shown in Figure 8. It could be seen from Equation (16) and Table 6 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and the brush speed. The interaction effect of each factor was not significant. As shown in Figure 8a, with an increase in the vibration frequency, the comprehensive picking index gradually increased. As shown in Figure 8b, with a decrease in the insertion depth, the comprehensive picking index rapidly increased at first but then slowly decreased. As shown in Figure 8c, the brush speed had little effect on the comprehensive picking index; with a decrease of the brush speed, the comprehensive picking index slowly increased. all factors, the vibration frequency had the greatest influence, followed by the insertion depth and the brush speed. The interaction effect of each factor was not significant. As shown in Figure 7a, with an increase in the vibration frequency, the damage rate of the ripe fruit gradually increased. As shown in Figure 7b, with a decrease in the insertion depth, the damage rate of the ripe fruit gradually increased. As shown in Figure 7c, with a decrease in the brush speed, the damage rate of the ripe fruit gradually increased. This was because the number of hits to the ripe fruit increased.  Response surface and contour plots of the damage rate of ripe fruit affected by the vibration frequency X 1 and brush speed X 2 (a), the vibration frequency X 1 and insertion depth X 3 (b), and the brush speed X 2 and insertion depth X 3 (c).
Horticulturae 2022, 7, x FOR PEER REVIEW 11 of 14 The response surface methodology was used to analyze the influence of various factors on the comprehensive picking index and the response surface of the regression equation is shown in Figure 8. It could be seen from Equation (16) and Table 6 that, among all factors, the vibration frequency had the greatest influence, followed by the insertion depth and the brush speed. The interaction effect of each factor was not significant. As shown in Figure 8a, with an increase in the vibration frequency, the comprehensive picking index gradually increased. As shown in Figure 8b, with a decrease in the insertion depth, the comprehensive picking index rapidly increased at first but then slowly decreased. As shown in Figure 8c, the brush speed had little effect on the comprehensive picking index; with a decrease of the brush speed, the comprehensive picking index slowly increased.

Field Experiment Verification
Based on the comprehensive picking index, the optimal parameter combination of each factor was obtained: a vibration frequency of 38.73 Hz, a brush speed of 14.21 mm/s, and an insertion depth of 26.07 mm. The field experiment ( Figure 9) was repeated 10 times to reduce the impact of random errors. The result showed that the harvesting rate of the ripe fruit was 83.65%, the harvesting rate of the unripe fruit was 7.22%, the damage rate of the ripe fruit was 11.49%, and the comprehensive picking index was 87.85.  Response surface and contour plots of the comprehensive picking index affected by the vibration frequency X 1 and brush speed X 2 (a), the vibration frequency X 1 and insertion depth X 3 (b), and the brush speed X 2 and insertion depth X 3 (c).

Field Experiment Verification
Based on the comprehensive picking index, the optimal parameter combination of each factor was obtained: a vibration frequency of 38.73 Hz, a brush speed of 14.21 mm/s, and an insertion depth of 26.07 mm. The field experiment ( Figure 9) was repeated 10 times to reduce the impact of random errors. The result showed that the harvesting rate of the ripe fruit was 83.65%, the harvesting rate of the unripe fruit was 7.22%, the damage rate of the ripe fruit was 11.49%, and the comprehensive picking index was 87.85.
Compared with the relevant experiments [1,12,25], it was found that the harvesting method by vibrating and comb-brushing was effective. In all cases, the harvesting rate of the ripe fruit was more than 80% and both the harvesting rate of the unripe fruit and the damage rate of the ripe fruit were less than 15%. Therefore, the vibrating and combbrushing harvesting method could be used for harvesting L. barbarum. We also found a few interesting phenomena in the field experiment. First, the higher the maturity of the ripe fruit, the easier it was to fall off. This could be because the connecting force of the stalk gradually decreased with the ripening of the fruit. Therefore, a reasonable selection of harvesting times may effectively improve the harvesting rate of ripe fruit. Second, because the ripe fruit and unripe fruit on the branches grew in turn, the harvesting could be stopped after harvesting in the ripe fruit area, thus reducing the harvesting rate of the unripe fruit. Finally, the vibrating and comb-brushing harvesting depended on the vibrating comb-brush bar to hit and brush the fruit to make it fall off, which led to a high damage rate of the ripe fruit. If the fruit could be harvested without touching, the rate of ripe fruit damage could be reduced.

Field Experiment Verification
Based on the comprehensive picking index, the optimal parameter combination of each factor was obtained: a vibration frequency of 38.73 Hz, a brush speed of 14.21 mm/s, and an insertion depth of 26.07 mm. The field experiment (Figure 9) was repeated 10 times to reduce the impact of random errors. The result showed that the harvesting rate of the ripe fruit was 83.65%, the harvesting rate of the unripe fruit was 7.22%, the damage rate of the ripe fruit was 11.49%, and the comprehensive picking index was 87.85. Figure 9. Field experiment. Figure 9. Field experiment.

Conclusions
In this study, faced with the standardized hedge cultivation mode of L. barbarum, a vibrating and comb-brushing harvester was designed. To optimize the vibrating and combbrushing harvest method, its main factors and parameter value ranges were obtained by the analysis of each system and a preliminary test of the harvester. By using the quadratic orthogonal rotation combination test, the mathematical models of the harvesting rate of the ripe fruit, the harvesting rate of the unripe fruit, the damage rate of the ripe fruit, and the comprehensive picking index were established. Furthermore, an optimized parameter combination was determined based on RSM with a vibration frequency of 38.73 Hz, a brush speed of 14.21 mm/s, and an insertion depth of 26.07 mm. The field experiment showed that the harvesting rate of the ripe fruit was 83.65%, the harvesting rate of the unripe fruit was 7.22%, the damage rate of the ripe fruit was 11.49%, and the comprehensive picking index was 87.85. These results provide a design basis for future research on the vibrating and comb-brushing harvest method of L. barbarum. However, the harvesting effect still needs to be improved. In further studies, according to the growth characteristics and cultivation mode, the mechanism of vibrating and comb-brushing harvesting should be studied. The harvesting time, area, and location could be optimized to increase the harvesting rate of the ripe fruit, reduce the harvesting rate of the unripe fruit, and the damage rate of the ripe fruit to achieve a better harvesting effect.

Conflicts of Interest:
The authors declare no conflict of interest.