Analyzing the Solution of Chemotaxis Equations with Logistic Source Term

When the parameters have the characteristics of sensitivity and specificity, while the traditional analysis method is the analysis of equations, under the influence of the characteristics of parameters, because the ability to consider the overall existence and limitations of the equation is weak, this leads to the deviation of the analysis results beyond the allowable range. With the chemotactic equations with logistic source terms, a new analysis method is proposed, on the basis of setting the optimal adjustment parameters of logistic source term, the method evaluates the well-posed behavior of the solution of the chemotactic equation system with logistic source term by quantitative analysis of the global existence and limitations of the solution. *e experimental results show that the proposed analysis method, after the analysis of the overall existence and limitations of the equation, has smaller deviation in the results of the appropriate qualitative evaluation and meets the analysis requirements of the solution problem of the chemotaxis equations. It can be seen that the analysis method is more practical.


Introduction
With the rapid development of science and technology, there are more and more cases that need to deal with high-dimensional data in human work and production. For these cases, we can use the method of variable selection to solve them. Whether it is regression or classification problem, the traditional analysis method is realized by using stepwise regression combined with information criterion. e practice proves that the traditional analysis method is practical, but there are many shortcomings at the same time. It is not stable to use stepwise regression method to select variables, and when there are many alternative variables of parameters, the calculation workload is huge, resulting in large calculation error, which cannot meet the requirements of the solution. erefore, overcoming the shortcomings of the traditional analytical methods has effectively become one of the hot spots in the study of the solution of chemotactic equations. e global existence of solutions is an important direction in the study of chemotaxis. e analysis method proposed in this paper is to analyze the solution of the system of chemotaxis equations with logistic source term and realize the accurate, systematic, and comprehensive problem analysis starting from the boundedness, global existence, and limitations [1].
is method provides an important analysis method for the future development of equation system solution and also provides a technical support for the study of equation system in other fields. In Brownian motion, a particle moves in response to the instantaneous difference in impact made by bombarding molecules on different sides of the particle. In chemotaxis, the motion of a cell is influenced by the molecules of the critical substrate through a chemical interaction rather than by the direct impact characteristic of Brownian motion [2]. In order to get advantages of advanced technologies and robotics, citizens have to take a move in the face of automation farming. Certain businesses take this step in provision of agricultural machinery for lease as well as for custom hiring [3]. e present-day FEM stands on three legs: mathematical models, matrix formulation of the discrete equation, and computing tools to do the numerical work.
e third is the one that has undergone the most dramatic changes [4,5]. Logistic regression analysis is a variant of traditional regression analysis for binary dependent variables which does not suffer from the weaknesses of the traditional approach to such data but which nevertheless shares many of its powerful features such as the ability to study quantitative as well as categorical independent variables and to include in the model interactions among the independent variables [6,7]. It is known that logistic source term regression is a kind of generalized linear regression model, which is used to classify 0/1 problem; that is to say, the prediction result belongs to "0" or "1" classification problem. Suppose there is a continuous explanatory variable Y′ in theory, which represents the possibility of the event, and the value range is W. When the value of the variable crosses a critical point j, for example, when j � 0, the event occurs. So, there are the following: when Y ′ > 0, Y � 1, and in other cases, Y � 0. Here, Y is the actual observed response variable, Y � 1 means the event occurred, and Y � 0 means the event did not occur. Suppose that there is a linear relationship between Y ′ and the explanatory variable X, namely,

Analytical
In the above relationship, α 0 indicates no event; α 1 indicates an event has occurred; y indicates an error term. en, the probability equation of the event can be described by the following equation: Generally, it is assumed that the error term y in (2) has logistic distribution or standard normal distribution. Since both the logistic distribution and standard normal distribution are symmetrical, we can get the following results: According to the above equation, K(α 0 + α 1 X) is the distribution function of error c, and the form of distribution function K(α 0 + α 1 X) depends on the hypothetical distribution of c [8]. When c is a logistic distribution, After simplification and treatment, the following equation is obtained: where e represents a natural number. Assuming the conditional probability P(Y � 1|X) � P ′ , a regression model with logistic source term can be obtained, as shown in the following equation: From the above calculation, it can be seen that P ′ is a nonlinear function of X, which can be transformed into a linear function. According to (6), the conditional probability that an event will not occur is defined as the following equation: Based on the above two equations, the ratio of the probability of occurrence to the probability of nonoccurrence is obtained as follows: is ratio is called the occurrence ratio of events and is recorded as O. In statistics, ln(P ′ /1 − P ′ ) is called logit transformation [9]. e logit transformation is carried out for the regression model of the logistic source term, as shown in Figure 1.
According to the curve trend in the figure above, the transformed result can be expressed in the following linear form: where ln( * ) represents the base function. When there are m explanatory variables in the above equation, the regression model of the corresponding logistic source term can be described by the following equation: where P ′ � P(Y � 1|X 1 , X 2 , . . . , X k ) is the conditional probability of Y � 1 under the action of m explanatory variables. On the basis of the above model settings, when estimating the optimal adjustment parameters of the logistic source term, it is also necessary to add a penalty term to the 2 Mathematical Problems in Engineering parameters to realize the estimation of parameters and the selection of variables [10]. Suppose the response variable Y is a binary variable, whose values are Y � 1 and Y � 0.
e m explanatory variables that affect the value of Y are X 1 , X 2 , . . . , X m . Suppose that there are n observations in the data set; that is, { }, and k � 1, 2, . . . , n. en, the following regression model of logistic source term is considered in the following equation: where P k ′ � P(y k � 1|x k ) is the conditional probability of the result y k � 1 under the given x k condition. Similarly, the conditional probability of y k � 0 is P(y k � 0|x k ) � 1 − P ′ . Combining (8), the probability of an observation value is obtained as follows: It is known that observations are independent of each other, so their joint distribution can be expressed as the product of each marginal distribution; then the likelihood function of n observations is as follows: In the above equation, F(α) represents the likelihood function for all events α. According to the above equation, the log likelihood function is obtained as follows: In the equation, P k ′ � e α 0 +α 1 x k1 +···+α p x kp /1 + e α 0 +α 1 x k1 +···+ α p x kp . k ′ represents the observation value adjacent to k observation value; α k′ represents the event corresponding to adjacent observation value; x kk′ represents the relationship parameter between two adjacent observation values. According to the related theory of Lasso method, the coefficient estimates α ′ in Lasso-logistic regression model are given by the minimum value of convex function of the following equation: In the equation, S(α) is the logarithmic likelihood function of the estimated value α′ . e coefficient estimates α′ can be written as follows: e above equation can be divided into two parts: in the first part, n k�1 (y k (α 0 + p k′�1 α k′ x kk′ ) − ln(1 + exp (α 0 + p k′�1 α k′ x kk′ ))) represents the fitting degree of the model; in the second part, λ p k�1 |α k | represents the punishment degree for the selected variables, where λ represents the adjustment parameter, and there is λ ≥ 0. Finally, the selected variable set should minimize the values of these two parts [11]. Using Lasso-logistic model for variable selection, the key is to adjust the parameter λ. e size of λ affects whether the coefficient of some variables can be compressed to 0 in the process of variable selection, which will directly affect the result of variable selection. K-fold cross validation, generalized cross validation, and BIC criterion are also suitable for Lasso-logistic model. In this analysis, K-fold cross validation is used to select the adjustment parameters.
Let L(y, g(x)) be the loss function. Firstly, the data set is randomly divided into K parts with approximately the same capacity. e first k(k � 1, 2, . . . , K) − th data is taken out as the verification set in turn, and the remaining K − 1 − th data is taken as the training set. Finally, the verification results on all K data are combined to obtain the CV value of K-fold cross validation; namely, the pretest error of K-fold cross validation is shown in the following equation: where n is the sample size, y k is the observation value of the k -th sample, k(k) is the number of the k -th sample, and g − k(k) is the model fitted on the remaining K − 1 − th data after removing the k(k) − th data. It is assumed that a group of models are fitted and expressed as g − k (x, λ), and the model contains the adjustment parameter λ [12]. en the prediction error obtained by K-fold cross validation is shown in the following equation: Finding λ, which makes CV(λ) the minimum, is the optimal adjustment parameter of Lasso-logistic model, so as to reset the optimal adjustment parameter of logistic source term.

Quantitative Analysis of the Global Existence and Solutions' Limitations of Chemotaxis Equations with Logistic
Source Terms. In order to analyze the solution of chemotaxis equations with logistic source term, it is necessary to obtain the global existence and limitation of the solution of chemotaxis equations by quantitative analysis. e known chemotaxis equations with logistic source term are as follows: In the above equations, Ω ⊂ * 3 is the bounded region with smooth boundary, n is the normal vector outside the unit on zΩ, μ ∈ * , σ > 0. x t represents the trend primary equation under the influence of time series parameters; z t represents the dependent variable affected by the result of the equation; Δx and Δz represent the difference value, respectively; T represents the periodic parameter; u represents the variable of the trend equation. e global existence and boundedness of initial boundary value solutions are discussed [13]. Assume that the initial value meets the conditions shown in the following equation: e classical solutions of homogeneous initial boundary value problems for chemotactic systems are global and bounded. If the above conditions are true, then when the adaptive parameter σ ≥ 3 + (23/3)‖z 0 ‖ 2 L ∞(Ω) is limited, the system of (19) is in the only bounded global classical solution (x, z), and there are constants at the same time. For all the time series parameters Q > 0, there are Among the above equations, L represents the constraint function of missing equation and W represents the maintenance function. When the constant Q satisfies the above equation, the assumption is true; that is, the initial value assumption is true, which can quantitatively analyze the global existence and limitations of the solution of the chemotactic equation system [14]. It is known that, for any conditional parameter t ∈ (0, T max ), there is a condition shown in the following equation: It is proved that in (19), both sides of the first equation are multiplied by x at the same time and both sides are integrated. For any t ∈ (0, T max ), the following equation can be obtained: Because of Young inequality, we can get According to the known conditions, the existence analysis in (25) is also obtained by combining Young's inequality.
According to (24), the above conclusions are obtained. When the boundary constraint parameter μ and separation parameter σ meet (25), it is proved that the solution of the chemotactic equation system is complete, so that the general curve of the chemotactic equation can be obtained [15]. Figures 2(a) and 2(b) show the general curve of chemotaxis equations.
According to the curve trend in the figure above, under the influence of two parameters, the result curve obtained by solving has symmetry characteristics, so it can be seen that the whole solution of chemotaxis equation system with logistic source term has complete existence [16]. Next, we will prove the global limitation of the solution.
If σ ≥ 3 + (23/3)‖z 0 ‖ 2 L ∞(Ω) is still assumed, then there is a constant Q ′ for any t ∈ (0, T max ), which makes the following equations set up: en, let h(t) � 4 Ω x 2 + Ω |∇z| 4 + Ω x|∇z| 2 ; the quantitative analysis process exists in the same reason, and the inequality can be gotten: Using Young's inequality for (27), we can get that, for any μ > 0, σ > 0, there is a constant Q > 0. It is known that Ω |∇z| 2 ≤ U, where U is the effective constant; combining (18), then according to Gagliardo-Nirenberg interpolation inequality, we can get In the above equation, e represents the natural base number; s represents the failure sequence. λ 1 > 0, which is the first nonzero eigenvalue of the operator ∇ in the bounded region Ω under the boundary condition. Because ‖x(·, s)‖ L ∞ (Ω) has been proved, it can judge the solution of chemotaxis equations with logistic source term, which has limitations [17].

Evaluating the Well-Posed Behavior of Solutions of Chemotaxis Equations.
On the basis of the above analysis, the global existence and limitations of the solutions of the chemotaxis equations are determined to evaluate the wellposed behavior of the solutions of the chemotaxis equations. According to the limitation of the solution of the chemotactic equation system for logistic source term, the linear term of the equation system is controlled. In this paper, the symmetrizer chemotaxis equation is used to simplify and form a symmetric hyperbolic function compatible with the boundary [18]. However, the difficulty in this process lies in how to select proper symmetrizer to make some linear terms disappear in the process of symmetrization evaluation. On the other hand, after choosing the proper symmetrizer and transforming the equation with logistic source term into the chemotaxis equation with symmetric hyperbolic rate, due to the particularity of the boundary, it cannot directly use the energy estimation to get the estimation of the spatial derivative and the spatialtemporal mixed derivative but needs to use the time derivative to control and then gets their estimates. e relevant Mathematical Problems in Engineering chemotaxis equation is given again, and the well-posed solution of the chemotaxis equation group with logistic source term is evaluated. e equation is the same as (19), and the quasi-transformed steady-state equation group is as follows: In the equation, ε 1 and ε 2 represent the intrinsic parameter and the well-posed parameter, respectively; ∇ε 1 and ∇ε 2 represent the gradient operator of the two parameters, respectively; Θ represents the behavior condition; ∇Θ represents the gradient operator of behavior condition; ΔΘ represents the missing function operator of the behavior condition; F(x) represents the steady-state control function. e specific description equation of behavior conditions is as follows: In the equation, κ represents the set of behavior condition constraints. In order to meet the practical application, it is necessary to consider that the solution of chemotaxis equation system meets the stability requirements, so as to ensure that the solution result will not be completely unbalanced due to the change of a certain parameter. e satisfied condition is inf x ε 1 > 0, inf x ε 2 > 0.
To solve the ordinary differential equation, combined with (29), there is a new equation: According to the satisfied conditions, the above equations are solved and the following system of equation is obtained: In the above equations, Q 1 is an unknown normal number. Taking (32) into the third equation and (3)-(2) of (29), combining with boundary conditions, the following elliptic equation can be obtained: erefore, if we want to obtain the well-posed solutions of the steady-state equations (29) with insulation boundary conditions, that is, (30), it is equivalent to getting the wellposed solutions of the elliptic equations (33) [19]. According to the elliptic equations, the corresponding variational functional can be taken as in the following equation: According to the above equation, the general curve of the elliptic system of equations is drawn. According to the maximum and minimum values of the curve changes, the extremum element of the functional is found out; that is, the function that makes this functional takes the extremum; namely, the solution of the elliptic system of equations [20]. Figure 3 is the curve diagram of the elliptic equation.
According to the intersection points marked in the above figure, it proves the minimizer of functional I in space D � Θ|( Ω |∇Θ| 2 dx) 1/2 + ( Ω |Θ| 2 dx) 1/2 < ∞, x ∈ Ω . e norms defining the space D are I(Θ n ) of the equation is a convex functional. In the next step, we need to know whether ‖Θ n ‖ D ≤ Q 1 holds, where Q 1 represents a positive constant. By using the counterproof method, supposing that there is ‖Θ n ‖ D ⟶ ∞, whether Θ n is positive or negative, there is e Θ n > Θ 2 n , so there will be I(Θ n ) ⟶ ∞. But this contradicts the fact that Θ n is a minimization sequence of I, so ‖Θ n ‖ D ≤ Q 1 holds. According to this equation, there is a weak convergence subsequence Θ ni of the minimization sequence Θ n , so that when ni ⟶ ∞, there is a limit Θ in D(Ω) that satisfies Θ ni ⟶ Θ, so that the evaluation on the well-posed behavior of the solution of the chemokine equation system is completed, and the analysis of the solution of the chemokine equation system with logistic source term is realized [21].

Test Results and Analysis
Compared with the traditional analysis method, the two analysis methods consider the calculation error when dealing with the solution of the chemotaxis equations. e test results of the proposed analysis method are set as experiment A and the test results of the traditional analysis method as experiment B. e sensitivity index and specificity index of the parameters of the experimental test object are set, and two analysis methods are used to analyze the solution of the chemotaxis equations with sensitivity and specificity parameters. Table 1 shows the sensitivity and specificity values.
Two analysis methods are used to analyze the solution of chemotaxis equations with logistic source term under the influence of sensitivity and specificity parameters, and the In the equation, R(x) is the reliability measurement result of the analysis result x; λ (a,b) is the allowable variation of the parameter within the allowable error range (a, b); S is the support parameter; C is the reliability parameter. At the same time, the error prediction model is as follows: In the equation, ω is the prediction weight coefficient; β is the prediction limit parameter; n is the number of direct predictions; k is the linear parameter; μ i is the floating value of the i − th group of test data. Figures 4 and 5 are the analysis error diagrams of the two methods when the sensitivity and specificity parameters are different.
It is known that, in the above two groups of test results, the dotted box represents the allowable error range, and the solid dot represents the error node. According to the distribution of error nodes, we can see that, under the proposed analysis method, the error nodes in group A are completely included in the allowable range of error, these errors are relatively concentrated, and the difference floating is very small. However, under the traditional analysis method, some of the error nodes in group B exceed the allowable error range, and the error nodes in the range are relatively scattered with different values. According to the test results of the two groups, the proposed analysis method is superior to the traditional analysis method in controlling the calculation error.

Conclusions
According to the solution, problem of chemotaxis equations with logistic source term, a more targeted analysis method is proposed. e sensitivity index and specificity index of the parameters of the experimental test object were set, and two analysis methods were used to analyze the solution of the chemotaxis equations with sensitivity and specificity parameters. rough quantitative analysis, the global existence and limitations of the solution of chemotaxis equations are obtained. e reliable control of the calculation error is realized. Accuracy of the solution result is improved, and a more cautious analysis method is provided for the solution problem of equations in the mathematical field. It was found that the proposed analysis method is better when compared to the traditional analysis method in controlling the calculation error.

Data Availability
e data used to support the findings of this study are available upon request.  Mathematical Problems in Engineering 9