Regression-Based Correction and I-PSO-Based Optimization of HMCVT’s Speed Regulating Characteristics for Agricultural Machinery

: To improve the speed regulating characteristics of continuously variable transmission for agricultural machinery, in order to meet the engineering and technical requirements of precision agriculture and intelligent agriculture, the paper researches and proposes a method combining the analysis of speed regulating characteristics, regression-based correction, and the improved particle swarm optimization (I-PSO) algorithm. First, the paper analyzes the degree of deviation between the linearization degree and the theoretical value of the speed regulating characteristics of the variable-pump constant-motor system of agricultural machinery according to the measurement results of the bench test. Next, the paper corrects the speed regulating characteristics and compares the regression results based on four models. Finally, the paper proposes a design method for the expected speed regulating characteristics of agricultural machinery and it completes the optimization of speed regulating characteristics and the matching of transmission parameters with the I-PSO algorithm. Results indicate that the speed regulating characteristics of the variable-pump constant-motor system show high linearization (with a coefﬁcient of determination of 0.9775). The theoretical and measured values of the speed regulating characteristics have a certain deviation (with a coefﬁcient of determination of 0.8934). Therefore, correcting the speed regulating characteristics of the variable-pimp constant-motor system is highly necessary. In addition, the second reciprocal function model proposed has the highest correction precision (with a coefﬁcient of determination of 0.9978). The I-PSO algorithm is applicable to the design and application of hydro-mechanical continuously variable transmission (HMCVT) for agricultural machinery. The new method proposed can improve the HMCVT’s speed regulating characteristics efﬁciently and quickly. It also ensures that the speed regulating characteristics are highly consistent with the expected design characteristics (with a mean error of 1.73%). Thus, the research offers a theoretical direction and design basis for the research and development of continuously variable transmission units in agricultural machinery.


Introduction
There is a wide variety of agriculture and forestry machinery. The tractor is one of the most important pieces of agricultural working machinery applied in the fields of agriculture and forestry [1][2][3]. Moreover, the agriculture working machinery also includes the grain harvester combine, the cotton picker, and so on. The forestry working machinery includes the skidder, the forest fire truck, the harvesting-cultivating combination machine, and so on [4]. The agriculture and forestry machinery generally works in the conditions of severe environments and variable loads [5][6][7][8]. Running reliably at the required speed is one of the most basic performance requirements of agriculture and forestry machinery. It requires a wider range of speed changes. the theoretical analysis formulas of kinematics and mechanics combined with the Newton-Raphson method for parameter matching. Cheng et al. [35] performed a mechanism analysis of the HMCVT's speed regulating characteristics, composed of the multi-planetaryrow compound transmission, and designed a non-equal-ratio-transmission HMCVT using the improved genetic algorithm. The research matched the required range of tractor working speed with the change range of each HMCVT section transmission ratio. This helped to improve the flexibility of the HMCVT's speed regulating characteristics design.
To sum up, using the theoretical relational expression of speed regulating characteristics to design the HMCVT is common. However, the method depends on the precision of the model of the speed regulating characteristics. The HMCVT is composed of the mechanical system (generally composed of the fixed-shaft gear pair and the planetary gear train) and the hydraulic system (generally composed of the pump and the motor). In the process of transmission, it can be considered that the actual transmission characteristics are consistent with the theoretical transmission characteristics of the mechanical system. However, the hydraulic system is greatly affected by the environment, thus causing a certain deviation between the actual and theoretical transmission characteristics. Therefore, only when the transmission characteristics of the mechanical and hydraulic systems are both correct, the HMCVT's speed regulating characteristics have certain precision. In particular, it is difficult to improve the precision of the HMCVT's speed regulating characteristics model in the research, development, and design stage. Generally, the comparative verification of the speed regulating characteristics can only be done after the HMCVT is made, thus resulting in increasing research and development costs and duration. In addition, current studies have a serious deficiency in this respect.
To solve the problem described above, the paper proposes a regression-based correction method and an I-PSO-based optimization method for the HMCVT's speed regulating characteristics. The paper mainly researches the following three aspects. Firstly, the paper describes a test of the speed regulating characteristics of the hydraulic system using the test bench of a variable-pump constant-motor system. Then, based on the measurement results of the test, the paper analyzes the linearization degree of the measured output speed and the deviation degree between the measured and theoretical values. The method of least squares is used for the regression analysis of the speed regulating characteristics of the variable-pump constant-motor system. The regression model is considered as the correction model of the HMCVT hydraulic system's speed regulating characteristics. Secondly, according to agricultural machinery's working speed requirement, engine characteristics, and the corrected model of the HMCVT's speed regulating characteristics, the paper uses the I-PSO algorithm for the optimization design of the speed regulating characteristics and matching of transmission parameters. Thirdly, using the coefficient of determination R 2 and the mean absolute percentage error MAPE as evaluation indexes, the paper compares before and after optimization and before and after correction of the speed regulating characteristics. Comparison results show that the corrected model of the speed regulating characteristics of the hydraulic system has high precision; the speed regulating characteristics based on the I-PSO algorithm are highly consistent with the expected characteristics and the matching of the transmission design parameters has a good result. According to the method proposed, it is only necessary to test the output speed of the pump-motor system chosen to complete the optimization design of the speed regulating characteristics and the matching of the transmission parameters of the HMCVT effectively in the research, development, and design stage. Figure 1 shows the transmission scheme of the tractor HMCVT researched in the paper. The continuously variable transmission unit has eight gear pairs (suppose that their transmission ratios are i 1 , i 2 , i 3 , i 4 , i 5 , i 6 , i 7 , and i 8 , respectively), two planetary rows P 1 and P 2 (suppose that the parameters of planetary rows are k 1 and k 2 , respectively), three wet clutches C 0 , C 1 , and C 2 , and one variable-pump constant-motor system.  Figure 1 shows the transmission scheme of the tractor HMCVT researched in the paper. The continuously variable transmission unit has eight gear pairs (suppose that their transmission ratios are i1, i2, i3, i4, i5, i6, i7, and i8, respectively), two planetary rows P1 and P2 (suppose that the parameters of planetary rows are k1 and k2, respectively), three wet clutches C0, C1, and C2, and one variable-pump constant-motor system. When a tractor runs normally and works (corresponding to the HMCVT hydraulicmechanical-power mixed working mode), the overall power of the transmission system is output from the engine and then input into the variable-pump constant-motor system (the hydraulic power flow) and the planetary gear train (the mechanical power flow). The divided power converges through the planetary gear train and then is output through P1's gear ring or P2's planetary carrier. When a tractor starts up (corresponding to the purely hydraulic working mode of HMCVT), the overall power of the transmission system is output from the engine and then input into the variable-pump constant-motor system directly, and drives gear pair i8 directly to offer the power to the output end. The core part of the continuously variable unit (i.e., the power confluence mechanism) of the tractor is the Simpson planetary gear transmission mechanism and has three working sections (including the purely hydraulic working section H0 and the hydraulic-mechanical mixed working sections HM1 and HM2). Figure 1 also shows the power output routes of working sections, in which sections H0, HM1, and HM2 are shown in blue, green, and red, respectively.

Working Principle of Tractor HMCVT Researched
The calculation formulas for the speeds of the confluence mechanism planetary carrier, gear ring, and sun gear are as follows [36,37]: in which r n is the output speed of the gear ring; c n is the input speed of the planetary carrier, and s n is the input speed of the sun gear.
Theoretically, the transmission relationship between the variable pump and constant motor is as follows [38]: When a tractor runs normally and works (corresponding to the HMCVT hydraulicmechanical-power mixed working mode), the overall power of the transmission system is output from the engine and then input into the variable-pump constant-motor system (the hydraulic power flow) and the planetary gear train (the mechanical power flow). The divided power converges through the planetary gear train and then is output through P 1 's gear ring or P 2 's planetary carrier. When a tractor starts up (corresponding to the purely hydraulic working mode of HMCVT), the overall power of the transmission system is output from the engine and then input into the variable-pump constant-motor system directly, and drives gear pair i 8 directly to offer the power to the output end. The core part of the continuously variable unit (i.e., the power confluence mechanism) of the tractor is the Simpson planetary gear transmission mechanism and has three working sections (including the purely hydraulic working section H 0 and the hydraulic-mechanical mixed working sections HM 1 and HM 2 ). Figure 1 also shows the power output routes of working sections, in which sections H 0 , HM 1 , and HM 2 are shown in blue, green, and red, respectively.
The calculation formulas for the speeds of the confluence mechanism planetary carrier, gear ring, and sun gear are as follows [36,37]: in which n r is the output speed of the gear ring; n c is the input speed of the planetary carrier, and n s is the input speed of the sun gear. Theoretically, the transmission relationship between the variable pump and constant motor is as follows [38]: in which n m is the speed of the constant motor; n p is the working speed of the variable pump, and ε is the displacement ratio of the variable-pump constant-motor system. Combining Equation (1) with Equation (2), we then obtain the relational expression of the 3-section HMCVT's speed regulating characteristics (i.e., the relational expression between the displacement ratio of the pump-motor system and the transmission ratio of HMCVT) as follows [39]:  Figure 2 shows the variable-pump constant-motor test bench built for the HMCVT for a tractor.  Table 1 gives the ranges of speed and torque. The test bench uses the variable pump of the Linde HPV-02 model with 55 cm 3 /rev displacement and 75 kW continuous working power, and the constant motor of the Linde HMF-02 model with 55 cm 3 /rev displacement and 93 kW continuous working power. Test bench of 'variable-pump and constant-motor' system. Note: 1 Engine (DEUTZ TCD2013L062V); 2 Gear box combined by the 'variable-pump and constant-motor' system; 3 Speed torque sensor of ZJ-5000A model; 4 Auxiliary gearbox; 5 Electrical eddy current dynamometer of DW250 model; 6 Hydraulic system (realizing lubrication, cooling, and other functions).

Test of Speed Regulating Characteristics for 'Variable-Pump and Constant-Motor' System
The test bench uses the speed torque sensors of the ZJ-2000A model and ZJ-5000A model of Lanling Jiangsu, China. Table 1 gives the ranges of speed and torque. The test bench uses the variable pump of the Linde HPV-02 model with 55 cm 3 /rev displacement and 75 kW continuous working power, and the constant motor of the Linde HMF-02 model with 55 cm 3 /rev displacement and 93 kW continuous working power.

Test of Speed Regulating Characteristics for 'Variable-Pump and Constant-Motor' System
The test aims to measure the speed regulating characteristics of the variable-pump constant-motor system in the real environment, and offers measured data to the studies on the linearization degree of the motor output speed, the deviation degree between the measured value and theoretical value of output speed, and the correction of the speed regulating characteristics based on regression.
The test requires us to fix the variable pump's input speed and record the constant motor's output speed by changing displacement ratio ε.
According to Equation (2), the change range of the absolute value of the displacement ratio is 0~1. Therefore, to cover the change range of displacement ratio ε completely, the test chooses displacement ratios of 0.2, 0.25, 0.3, 0.375, 0.5, 0.625, 0.75, and 1 (in total, 8 groups of tests), and considers the average of stable constant motor output speeds in each group of tests as the output speed in the current working condition. In addition, the research records the input speed data of the variable pump to analyze and determine whether the volatility of the input speed is reasonable. The analysis aims to alleviate the influence of the input speed volatility on the system's output speed.

The Analysis Method for Speed Regulating Characteristics of 'Variable-Pump and Constant-Motor' System
The paper calculates and obtains n p (the average input speed of variable pump) and σ p (the standard deviation) to determine whether the volatility of the input speed is reasonable. The calculation formulas of average input speed n p and standard deviation σ p are as follows: in which N p is the total number of measured data and n pi is the value of the ith measured datum. The higher the linearization degree of the HMCVT's speed regulating characteristics, the better the adjustability of the tractor's driving speed [24]. According to the expression of the speed regulating characteristics of the HMCVT used, i.e., Equations (3)-(5), the HMCVT's output speed and displacement ratio ε have a linear variation relationship theoretically. The gear system transmission is reliable, and its theoretical and actual transmission characteristics are basically consistent. Therefore, the linearization degree of the variable-pump constant-motor system's speed regulating characteristics plays a decisive role in the linearization degree of the HMCVT's speed regulating characteristics.
The paper proposes to perform linear fitting to the motor output speed based on the least squares method and using R 2 , the coefficient of determination of linear fitting, as the determination basis of the linearization degree. The following is the calculation formula of R 2 : (8) in which y ideal is the theoretical value, y measured is the measured value of the bench test, and N is the total number of data to be fitted.
The paper uses the coefficient of determination R 2 to evaluate the degree of consistency between the measured and theoretical output speeds of the pump-motor system.

The Regression-Based Correction Method of HMCVT's Speed Regulating Characteristics for Tractor
According to Equation (2), the relational expression of i pm , the transmission ratio of the variable-pump constant-motor system (i.e., n p /n m , the ratio of the input speed of pump and the output speed of motor), and the displacement ratio is as follows: According to the variation trend of measured data, the research uses the linear and nonlinear least squares methods for the regression analysis of the speed regulating characteristics of the variable-pimp constant-motor system. The polynomial regression [40,41] has good generalization ability and strong applicability in the field of engineering technology. In addition, according to Equations (2) and (9), the research adds items to the numerator and the denominator based on the original theoretical model to correct the model. To sum up, the paper uses four models in the regression analysis, including the 2-order polynomial model, the 3-order polynomial model, and two ε reciprocal function models. The four models' forms are as follows.
The 2-order polynomial model of speed regulating characteristics: Agriculture 2022, 12, 580 The 3-order polynomial model of speed regulating characteristics: The first reciprocal function model of speed regulating characteristics: The second reciprocal function model of speed regulating characteristics: The paper uses the coefficient of determination R 2 to evaluate and determine the final variation expression of the speed regulating characteristics of the variable-pump constant-motor system, and substitutes the final form of i pm into Equations (3)-(5) to obtain the corrected change law of the speed regulating characteristics of the tractor HMCVT. Figure 3 shows the technical route of the correction method proposed.

The I-PSO-Based Optimization Design Method of HMCVT's Speed Regulating
Characteristics for Tractor Figure 4 shows the original speed regulating characteristics of the HMCVT researched. The HMCVT of this model has three sections (H0, HM1, and HM2), of which the speed regulating characteristics intersect mutually. Moreover, there are two synchronous section-shift transmission ratios. The design aims to ensure the continuity of the transmission ratio of the HMCVT.

The I-PSO-Based Optimization Design Method of HMCVT's Speed Regulating
Characteristics for Tractor Figure 4 shows the original speed regulating characteristics of the HMCVT researched. The HMCVT of this model has three sections (H0, HM1, and HM2), of which the speed regulating characteristics intersect mutually. Moreover, there are two synchronous section-shift transmission ratios. The design aims to ensure the continuity of the transmission ratio of the HMCVT. To make the HMCVT system easier to control, the two synchronous section-shift transmission ratios should be distributed in the marginal area with the displacement ratio ε of −1~1 as much as possible because, in this case, the displacement ratio ε has a wider range of change. If the actually required displacement ratio ε has a small change range, the requirement to control precision is high. For instance, suppose that the designed HMCVT transmission ratio i cvt has a change range of 1~11; then, if displacement ratio ε changes in the range of −1~1, ε changing by 0.1 corresponds to i cvt changing by 0.5, on average; if displacement ratio ε changes in the range of −0.5~0.5, ε changing by 0.05 corresponds to i cvt changing by 0.5 on average. This indicates that the wider displacement ratio ε's available range is, the more applicable it is to practical engineering.

HMCVT transmission ratio
The paper sets the tire radius of the tractor as 0.9 m and the other transmission system's transmission ratio as 9. Figure 5 shows the external characteristics of the diesel engine used.  To make the HMCVT system easier to control, the two synchronous section-shift transmission ratios should be distributed in the marginal area with the displacement ratio ε of −1~1 as much as possible because, in this case, the displacement ratio ε has a wider range of change. If the actually required displacement ratio ε has a small change range, the requirement to control precision is high. For instance, suppose that the designed HMCVT transmission ratio cvt i has a change range of 1~11; then, if displacement ratio ε changes in the range of −1~1, ε changing by 0.1 corresponds to cvt i changing by 0.5, on average; if displacement ratio ε changes in the range of −0.5~0.5, ε changing by 0.05 corresponds to cvt i changing by 0.5 on average. This indicates that the wider displacement ratio ε 's available range is, the more applicable it is to practical engineering. The paper sets the tire radius of the tractor as 0.9 m and the other transmission system's transmission ratio as 9. Figure 5 shows the external characteristics of the diesel engine used. In Figure 5, the values of torque in the positions marked with the left and right red lines (i.e., the speeds of the diesel engine are 1085 and 1992 r/min, respectively) correspond to 90% of the maximum torque of the diesel engine. The design in the paper uses Wang's analysis of the research results of Resch and Renius for reference [42]. In the whole life cycle of the tractor, the proportion of time for which the tractor works in the speed section of 4~20 km/h is approximately 76~93% (in which the proportion of time in 4~12 km/h is 61~68% and the proportion of time in 12~20 km/h is approximately 15~25%).
Therefore, the paper designs section HM1 as the farmland working section of the tractor and section HM2 as the driving and transportation section of the tractor. The design can also avoid the frequent switch problem of the wet clutch. To sum up, the transmission ratios of the sections of the tractor HMCVT are shown in Table 2.
Engine torque/Nm In Figure 5, the values of torque in the positions marked with the left and right red lines (i.e., the speeds of the diesel engine are 1085 and 1992 r/min, respectively) correspond to 90% of the maximum torque of the diesel engine. The design in the paper uses Wang's analysis of the research results of Resch and Renius for reference [42]. In the whole life cycle of the tractor, the proportion of time for which the tractor works in the speed section of 4~20 km/h is approximately 76~93% (in which the proportion of time in 4~12 km/h is 61~68% and the proportion of time in 12~20 km/h is approximately 15~25%).
Agriculture 2022, 12, 580 9 of 18 Therefore, the paper designs section HM 1 as the farmland working section of the tractor and section HM 2 as the driving and transportation section of the tractor. The design can also avoid the frequent switch problem of the wet clutch. To sum up, the transmission ratios of the sections of the tractor HMCVT are shown in Table 2. Section H 0 is the starting section of the tractor, so, in the design, it should have a synchronous section-shift transmission ratio with section HM 1 , and the minimum vehicle speed corresponding to the synchronous section-shift transmission ratio should be equal to or greater than 4 km/h. Moreover, when the displacement ratio ε = 0, in section HM 1 , only the gear system is transmitting the power, so the power transmission efficiency is the maximum in this case. To ensure that section HM 1 fully plays its role in the case where the displacement ratio ε is 0, the paper considers the speed of the diesel engine corresponding to 90% of maximum torque as the critical speed for calculation, and calculates the range of vehicle speed corresponding to the critical speed. Meanwhile, we further calculate the proportion of the vehicle speed range in 4~12 km/h, and take the HMCVT transmission ratio corresponding to 80% of the maximum proportion as the designed value of the expected transmission ratio of section HM 1 in the case where displacement ratio ε is 0. The corresponding calculation formula (engine output speed through the HMCVT and other transmission systems to cause the tire to produce translation speed) of the tractor's speed and the HMCVT transmission ratio is as follows (0.377 is the coefficient used for unit conversion): u a = 0.377 r d n e i 0 i cvt (14) in which u a is the running speed of the tractor, r d is the radius of the wheel, n e is the working speed of the engine, and i 0 is the overall transmission ratio of the other transmission system. According to the calculation result, when displacement ratio ε = 0, the value range of i cvt , the HMCVT transmission ratio corresponding to section HM 1 , is 5.37~7.82.
According to Equations (3)-(5), the HMCVT's speed regulating characteristics have a nonlinearity characteristic and many characteristic parameters (i.e., many transmission parameters to be designed). We classify the transmission parameters of the HMCVT researched and obtain the following seven transmission parameters to be designed: i 1 and i 2 ; i 3 , i 4 , and i 5 ; i 6 ; i 7 ; i 8 ; k 1 ; and k 2 . If we use the enumeration method to match the transmission parameters and optimize the speed regulating characteristics (suppose that the precision of the transmission parameter of the mechanical part is 0.01), it is necessary to carry out 350 5 + 150 2 = 5.25 × 10 12 times of calculation and matching processes (the range of gear transmission ratio is 0.5~4 and the value range of planetary row's characteristic parameter is 2.5~4). In addition, there are constraint conditions for the transmission parameters in the matching process. Meanwhile, a great amount of work is required if replanning the HMCVT's speed regulating characteristics to match the working requirements of other agricultural machinery. The reasons above cause the difficulty in optimizing the speed regulating characteristics using the enumeration method. The heuristic intelligent optimization algorithm has obvious advantages in performance optimization, parameter matching, and identification, and has been used to accurately solve a series of complex engineering problems [43][44][45][46]. Therefore, the paper proposes using the I-PSO algorithm for the optimization of the speed regulating characteristics of the tractor HMCVT and the matching of the transmission param-eters. As for the I-PSO algorithm used, the paper uses the I-PSO algorithm process [47,48], proposed in previous research and verified for engineering applications, for reference.
The optimization of the speed regulating characteristics of the HMCVT is a multi-object optimization problem. The objective function f itness of the I-PSO algorithm proposed and used in the research is as follows: f itness = i HM 1 (ε = −1) − 14.14 + i HM 1 (ε = 1) − 4.34 + i HM 2 (ε = 1) − 4.34 (15) in which i HM 1 (ε = −1) and i HM 1 (ε = 1) are the transmission ratios of the HMCVT working in section HM 1 with the displacement ratios of −1 and 1, respectively, and i HM 2 (ε = 1) is the transmission ratio of the HMCVT working in section HM 2 with the displacement ratio of 1. According to Equation (15), the paper considers the sum of three sub-objective functions with equal weight as the overall objective function. Each term of a sub-objective function considers the absolute value of error between the design value and expected value as the calculation formula. Meanwhile, the constraint conditions of parameter matching are as follows: i HM 2 (ε = −1) ≤ 2.17 (16) 5.37 ≤ i HM 1 (ε = 0) ≤ 7.82 (17) i H 0 (ε = c) = i HM 1 (ε = c) ≥ 8.48 (18) in which i H 0 (ε = c) is the transmission ratio of the HMCVT working in stage H 0 when the displacement ratio is equal to c.
To sum up, Figure 6 shows the optimization process of the speed regulating characteristics based on the I-PSO algorithm proposed. To sum up, Figure 6 shows the optimization process of the speed regulating characteristics based on the I-PSO algorithm proposed.

The Comparative Analysis Method of Speed Regulating Characteristics
To compare the optimization results of the speed regulating characteristics better, the paper uses three comparative analysis methods.
First, we calculate the errors of i HM 1 (ε = −1), i HM 1 (ε = 1), and i HM 2 (ε = 1) of the speed regulation of the tractor HMCVT after optimization with their expected values to verify the optimization result of the speed regulating characteristics. Meanwhile, we check whether the transmission ratios of i HM 2 (ε = −1), i HM 1 (ε = 0) and the intersection position of section H 0 and section HM 1 are in the ranges specified by Equations (16)- (18).
Next, we compare and analyze the difference in speed regulating characteristics before and after optimization by calculating and generating figures.
Finally, based on the optimization results, we compare and analyze the difference in the speed regulating characteristics before and after correction by calculating and generating figures.
The research uses the coefficient of determination R 2 and the mean absolute percentage error MAPE for the analysis and evaluation of different types of speed regulating characteristics. See Equation (8) for the calculation formula of the coefficient of determination R 2 .
The following is the calculation formula of the mean absolute percentage error MAPE. (19) in which y ideal is the ideal expected value of the speed regulating characteristics of the HMCVT of the tractor, y design is the designed result of the speed regulating characteristics of the HMCVT of the tractor, and n is the total number of data of the speed regulating characteristics. Figures 7 and 8 show the test results of the speed regulating characteristics of the variable-pump constant-motor system. In the test, n p (the mean working speed of variable pump) is 718.57 r/min and σ p (standard deviation) is 0.88 r/min. The results indicate that in the whole test process of the speed regulating characteristics, the variable pump's working speed remains essentially stable, with low volatility. The variable pump's working speed ultimately will not affect the constant motor's output speed (i.e., it will not affect the speed regulating characteristics of the pump-motor system).   In Figure 8, the output speed of the constant motor (showing the law of step rising) corresponds to the cases where the displacement ratio ε is 0.2, 0.25, 0.3, 0.375, 0.5, 0.625, 0.75, and 1, respectively, from left to right. We use the least squares method for the linear fitting of the variation characteristics of displacement ratio ε , and obtain the following result:

The Results and Analysis of Speed Regulating Characteristics for 'Variable-Pump and Constant-Motor' System
2 R , the coefficient of determination of linear fitting, is 0.9775. The result shows that the variable-pump constant-motor system researched has a high degree of linearization, consistent with the original transmission law of the variable-pump constant-motor system. Meanwhile, it indicates that the variation characteristics between the tractor's speed and displacement ratio ε also show a linear relationship in the adjustment of the system. Figure 9 shows the comparison between the measured and theoretical values of the output speed of the variable-pump constant-motor system. In Figure 8, the output speed of the constant motor (showing the law of step rising) corresponds to the cases where the displacement ratio ε is 0.2, 0.25, 0.3, 0.375, 0.5, 0.625, 0.75, and 1, respectively, from left to right. We use the least squares method for the linear fitting of the variation characteristics of displacement ratio ε, and obtain the following result: R 2 , the coefficient of determination of linear fitting, is 0.9775. The result shows that the variable-pump constant-motor system researched has a high degree of linearization, consistent with the original transmission law of the variable-pump constant-motor system. Meanwhile, it indicates that the variation characteristics between the tractor's speed and displacement ratio ε also show a linear relationship in the adjustment of the system. Figure 9 shows the comparison between the measured and theoretical values of the output speed of the variable-pump constant-motor system. From Figure 9, we can see that the coefficient of determination 2 R between the measured and theoretical values is 0.8934. The result shows that the variation characteristics of the measured and theoretical values are basically consistent, but there is still a certain deviation. Therefore, correcting the model of the speed regulating characteristics of the variable-pump constant-motor system is particularly important. From Figure 9, we can see that the coefficient of determination R 2 between the measured and theoretical values is 0.8934. The result shows that the variation characteristics of the measured and theoretical values are basically consistent, but there is still a certain deviation. Therefore, correcting the model of the speed regulating characteristics of the variable-pump constant-motor system is particularly important.

The Results and Analysis of Regression-Based Correction of HMCVT's Speed Regulating Characteristics
We use the least squares method for the regression analysis of Equations (10)- (13). Figure 10 shows the results.
From Figure 9, we can see that the coefficient of determination 2 R between the measured and theoretical values is 0.8934. The result shows that the variation characteristics of the measured and theoretical values are basically consistent, but there is still a certain deviation. Therefore, correcting the model of the speed regulating characteristics of the variable-pump constant-motor system is particularly important.

The Results and Analysis of Regression-Based Correction of HMCVT's Speed Regulating Characteristics
We use the least squares method for the regression analysis of Equations (10)- (13). Figure 10 shows the results.  Table 3 shows the values and coefficients of determination 2 R of the parameters of the four correction models. According to Figure 10 and Table 3, the second reciprocal correction model has the proper number of parameters and the highest coefficient of determination. Therefore, we choose the second reciprocal correction model as the model of the speed regulating characteristics of the variable-pump constant-motor system. The research also shows that the correction based on the original law model produces a good result.
Therefore, the modified models of the speed regulating characteristics of the HMCVT of the tractor are shown as follows and the variation law of characteristics is shown in Figure 11.  Table 3 shows the values and coefficients of determination R 2 of the parameters of the four correction models. According to Figure 10 and Table 3, the second reciprocal correction model has the proper number of parameters and the highest coefficient of determination. Therefore, we choose the second reciprocal correction model as the model of the speed regulating characteristics of the variable-pump constant-motor system. The research also shows that the correction based on the original law model produces a good result. Therefore, the modified models of the speed regulating characteristics of the HMCVT of the tractor are shown as follows and the variation law of characteristics is shown in Figure 11.
Agriculture 2022, 12, x FOR PEER REVIEW 16 of 20 Figure 11. The comparison of speed regulating characteristics of transmission parameters before and after optimization. Figure 12 shows the iteration evolution curve of the optimization of the speed regulating characteristics based on the I-PSO algorithm.

The Results and Analysis of I-PSO-Based Optimization for Speed Regulating Characteristics
According to Figure 12, for the strongly nonlinear problem, the I-PSO algorithm produces a good optimization result. From the perspective of the overall iteration process, the I-PSO algorithm's optimization result can evolve continuously. It indicates that the I-PSO algorithm can avoid the problem of prematurity well.
After the optimization,    Figure 12 shows the iteration evolution curve of the optimization of the speed regulating characteristics based on the I-PSO algorithm.

The Results and Analysis of I-PSO-Based Optimization for Speed Regulating Characteristics
According to Figure 12, for the strongly nonlinear problem, the I-PSO algorithm produces a good optimization result. From the perspective of the overall iteration process, the I-PSO algorithm's optimization result can evolve continuously. It indicates that the I-PSO algorithm can avoid the problem of prematurity well.
After the optimization,  According to Figure 12, for the strongly nonlinear problem, the I-PSO algorithm produces a good optimization result. From the perspective of the overall iteration process, the I-PSO algorithm's optimization result can evolve continuously. It indicates that the I-PSO algorithm can avoid the problem of prematurity well.
After the optimization, i HM 1 (ε = −1) has an error of 1.08% compared with the expected value 14.14; i HM 1 (ε = 1) has an error of 2.32% compared with the expected value 4.3; i HM 2 (ε = 1) has an error of 1.79% compared with the expected value 4.34. This indicates that the optimization results of the speed regulating characteristics are consistent with the expected results.
After the optimization, i HM 1 (ε = 1) and i HM 2 (ε = 1) have an error of 0.51%; the transmission ratio characteristics of section H 0 and section HM1 have an intersecting point when displacement ratio ε is −0.39. This shows that after the optimization design, the tractor HM-CVT's speed regulating characteristics are continuous, without any point of discontinuity.
After the optimization, i HM 2 (ε = −1) = 0.95 < 2.17. The results show that the highest driving speed of the tractor after optimization meets the requirement of >40 km/h. Moreover, after the optimization design, i HM 1 (ε = 0) = 6.74, and the transmission ratio value in the intersection position of section H 0 and section HM 1 is 8.71. The optimization result completely meets the design requirement. Figure 11 shows the comparison results of the speed regulating characteristics of the HMCVT of the tractor before and after optimization.
According to Figure 11, after the optimization, the position of the intersection point of section HM 1 and section HM 2 moves to the position with a larger displacement ratio and finally stops in the position with the displacement ratio of 1. This improves the available range of the displacement ratio, completely consistent with the design expectation.
Based on the optimization result, Figure 13 shows the comparison results of the speed regulating characteristics before and after correction.  Figure 11 shows the comparison results of the speed regulating characteristics of the HMCVT of the tractor before and after optimization.
According to Figure 11, after the optimization, the position of the intersection point of section HM1 and section HM2 moves to the position with a larger displacement ratio and finally stops in the position with the displacement ratio of 1. This improves the available range of the displacement ratio, completely consistent with the design expectation.
Based on the optimization result, Figure 13 shows the comparison results of the speed regulating characteristics before and after correction. ) is 0.5139, 0.6257, and 0.6351, respectively, before and after correction; MAPE is 13.69%, 6.95%, and 10.24%, respectively. Meanwhile, section HM1 and section HM2 have no intersection point before correction. This shows that the speed regulating characteristics of the HMCVT are not even continuous. To sum up, if we use the speed regulating characteristics in ideal conditions for optimization design, the design results have certain errors compared with the actual transmission characteristics.

Conclusions
The paper proposes an analysis method for the linearization degree and deviation degree of the speed regulating characteristics of a variable-pump constant-motor system. The linear variation relationship between the actual output speed and displacement ratio of the tractor with the HMCVT is verified using the test bench results. Moreover, the research points out that the correction based on measured results is necessary for the speed regulating characteristics of the variable-pump constant-motor system. According to Figure 13, R 2 , the coefficient of determination, of the speed regulating characteristics of section H 0 (choosing the transmission data with the displacement ratio of −1 ≤ ε ≤ −0.4), section HM 1 (choosing the transmission data with the displacement ratio of −1 ≤ ε ≤ −0.5), and section HM 2 (choosing the transmission data with the displacement ratio of 0.5 ≤ ε ≤ 1) is 0.5139, 0.6257, and 0.6351, respectively, before and after correction; MAPE is 13.69%, 6.95%, and 10.24%, respectively. Meanwhile, section HM 1 and section HM 2 have no intersection point before correction. This shows that the speed regulating characteristics of the HMCVT are not even continuous. To sum up, if we use the speed regulating characteristics in ideal conditions for optimization design, the design results have certain errors compared with the actual transmission characteristics.

Conclusions
The paper proposes an analysis method for the linearization degree and deviation degree of the speed regulating characteristics of a variable-pump constant-motor system. The linear variation relationship between the actual output speed and displacement ratio of the tractor with the HMCVT is verified using the test bench results. Moreover, the research points out that the correction based on measured results is necessary for the speed regulating characteristics of the variable-pump constant-motor system.
The second reciprocal correction model proposed can be built easily (with only three parameters to be determined) and is highly consistent with the measured value of the speed regulating characteristics (the coefficient of determination is close to 1), so it should replace the original model of the speed regulating characteristics of the variable-pump constant-motor system. Moreover, the research shows that the correction of the original law model produces a good improvement result.
The method proposed for the optimization of the speed regulating characteristics and transmission parameter matching of agricultural machinery HMCVT is highly effective. Using the method combined with the I-PSO algorithm, the optimization design of the speed regulating characteristics and transmission parameter matching for the HMCVT can be completed effectively, merely by testing the output speed of the pump-motor system of the chosen model in the research, development, and design stage.

Data Availability Statement:
The data presented in this study are available on demand from the corresponding author or first author at (luzx@njau.edu.cn or chengzhun38@163.com).