Research Article  Open Access
ShengYuan Wang, WanMing Chen, Rong Wang, XiaoLan Wu, "Multiobjective Evaluation of Coevolution among Innovation Populations Based on Lotka–Volterra Equilibrium", Discrete Dynamics in Nature and Society, vol. 2021, Article ID 5569108, 14 pages, 2021. https://doi.org/10.1155/2021/5569108
Multiobjective Evaluation of Coevolution among Innovation Populations Based on Lotka–Volterra Equilibrium
Abstract
The collaborative evaluation of enterprise innovation populations is a hot issue. The Lotka–Volterra model is a mature method used to evaluate the interaction mechanism of populations and is widely used in innovation ecology research studies. The Lotka–Volterra model mainly focuses on the quantitative characteristics of the interactive populations. The growth mechanisms cannot explain all the synergy mechanisms of the innovative populations. The collaborative evaluation between enterprise innovation populations is a typical multiobjective evaluation problem. The multichoice goal programming model is a mature method to solve multiobjective optimization problems. This paper combines the Lotka–Volterra model and multichoice goal programming method to construct a threestage multiobjective collaboration evaluation method based on Lotka–Volterra equilibrium. An evaluation example is used to illustrate the application process of this method. The method proposed in this paper has excellent performance in computing, parameter sensitivity analysis, and model stability analysis.
1. Introduction
Innovation activity is a systems engineering issue. It is difficult for a single organization to have all the resources necessary for innovation. It has dynamic evolution, symbiosis, and concurrence formed by the connection and transmission of material flow and talent flow among various innovation populations (IPs) under a specific innovation environment in a certain region. The innovation resources in a certain region are limited. The basic problem addressed by theoretical research and social practice in the fields of economics and management is how to allocate resources optimally [1–6]. An appropriate scale, along with the development of an innovation ecosystem, can support the effective development of the IP. In this case, the resources in the innovation ecosystem can be fully utilized, and the innovation performances can be maximized [7, 8].
The population scale in a region determines the population density. Existing studies generally use the population density index to analyze the relationship between population scale and innovation performance. Some studies found that industrial clusters have a positive impact on enterprise innovation [9, 10]. However, other studies found that industrial clusters have a negative effect on enterprise innovation [11, 12]. Given these two different conclusions, it is necessary to establish a new framework to analyze the relationships between an enterprise’s population size and its innovation ability. Therefore, the topic of population ecology should be addressed. Hannan and Freeman [13] analyzed the enterprise problem based on the logic of population ecology and founded the organizational ecology field. This present study further develops and expands the population ecology framework.
The Lotka–Volterra model is often used in population ecology to analyze the cooperative or competitive relationship of populations. Studies have shown that the introduction of Lotka–Volterra, a population competition model in biology, into market competition and diffusion has produced better analysis results [14, 15]. Lee et al. used the Lotka–Volterra model to analyze the dynamic competitive relationship in the Korean stock market [16] and confirmed that the role of competitors will change over time. Kim, Lee, and Ahn studied the Korean mobile communications market through the Lotka–Volterra model and found the symbiotic relationship between competitors [17]. Kreng and Wang used the Lotka–Volterra model and the Lyapunov function to study the dynamic competition and balance point of PDP and LCD TV [18]. On the basis of previous studies, this paper selects interdependent innovation populations as the research object and uses the Lotka–Volterra model to study the level of synergy and balance between innovation populations in the innovation ecosystem. Different from previous studies, when using the Lotka–Volterra model, this study is not simply based on the scale of the innovation population to explore the synergy relationship, but on the number of innovation populations, input and output optimization, and other multiple objectives as the research basis to explore the relationship between innovation populations. The research organization of this article is shown in Figure 1.
As shown in Figure 1, the research organization of this article is as follows: (1) this paper uses population dynamics (Lotka–Volterra model) to analyze the influence mechanism between interactive innovation populations. (2) This paper builds multichoice goal programming (MCGP) and Lotka–VolterraMCGP models for innovation populations scale optimization. (3) This paper uses TOPSIS method to evaluate the collaboration of innovation populations.
In this study, the interaction model of IPs is constructed based on the perspectives of ecological theory and innovation theory, which form the theoretical basis of this study. A multichoice model is used to determine the appropriate scale of the IP under resource constraints, population synergy, and maximum output targets. To resolve the multichoice optimization problem mentioned above while taking into account the collaborative development of the IP, a comprehensive method is required. Based on the perspective of resource constraints, this paper constructs a dynamic model of the growth for the IP. This paper estimates the suitability of population size by using the multichoice goal programming method. Two proposed models are constructed to obtain the appropriate population scale. This research constructs the theoretical model from two aspects, namely, input constraint and output maximization, which is theoretically innovative. Our research also has practical significance because it provides an appropriate analysis method for various innovation subjects to analyze and plan the development of IPs.
2. Literature Review
Over the past 10 years, the concept of innovation ecosystem has become popular among the rapidly growing literature [19] that typically focuses on business and strategy. Oh et al. [20] criticized the concept innovation ecosystem with regard to its usefulness and distinctiveness in relation to existing conceptualizations of innovation systems and with regard to the biologically inspired “eco” qualifier, which is arguably a flawed analogy to natural ecosystems.
Related literature uses the theory of organizational ecology to study the coevolution process of corporate populations. Baum, Korn, and Kotha empirically studied the competitive advantages and survival rate of incumbents and new entrants in the telecommunication service industry after the technical standards have been established [21]. Low and Abrahamson studied the changes in corporate population organization in three different stages of emerging industries, growing industries, and mature industries [22]. Geroski and Mazzucato used the market size model, negative feedback model, infectious disease model, and density dependence model to empirically study the evolution process of the population of the US auto industry [23]. Hannan, Carroll, Dundon, and Torres [24] and Hannan [25] used the densitydependent model to empirically study the evolution process of the corporate population of the European auto industry.
Adomavicius, Bockstedt, Gupta, and Kauffman qualitatively studied the interdependence and mechanism of the three types of technologies: components, basic common technologies, and products and applications [26]. Adner and Kapoor empirically studied the influence of the difficulty of developing upstream components and downstream complementary products on the competitive advantage of integrated innovation leaders [27].
Almost all of the abovementioned studies are based on a certain industry, starting from the enterprise level, using the number of enterprises entering and exiting the market as the research data source, using ecological related models to verify the evolution of a single population or the coevolution relationship between two subpopulations. In reality, especially in hightech industries, it is more common to form an interinfluenced and interdependent technological innovation ecosystem around the industrial chain. There are few relevant literatures focusing on technological populations in the innovation ecosystem, and qualitative research is the main focus.
3. The Model of Innovation Population Growth
In this section, the growth model of the innovation populations (IPs) is to be constructed. Based on the Lotka–Volterra model, an innovation population relationship model is proposed, and the equilibrium point is analyzed.
3.1. Model Construction Ideas
When there are abundant resources, populations can grow at geometric or exponential rates. As resources are depleted, population growth rate slows and eventually stops. This is known as logistic population growth. The environment limits population growth by changing birth and death rates. On average, small organisms experience increases per capita at higher rates and more variable populations, while large organisms have lower increase rates per capita and less variable populations. In view of resource constraints and the need for specialization, it is difficult for any single firm to develop and commercialize a technologybased offering from start to finish [28, 29]. Increasingly complex constellations of organizations have been emerging in the form of innovation ecosystems, where actors interact with each other to create, deliver, and appropriate value.
The study of enterprise IP dynamics should consider the influence of regional constraints. The ecological sense of a population is a collection of certain organisms within a given time and space. The region where the population grows is a relatively homogeneous nonlinear region that is different from the surrounding environment. There are universal temporal and spatial constraints in natural hierarchy systems.
The spatial distribution characteristics of different ecological regions such as size, shape, boundary, nature, and distance make up different ecological zones, forming the differences of ecosystem and regulating population growth. The model of enterprise IP dynamics focuses on the quantity change in the IP. Its changing rule is based on the nonlinear growth principle of biological population quantity. The growth model of most species is nonlinear in nature. The number of innovative enterprises may change rapidly with the influence of incentivebased policies and innovation resources in a given period and in a given area.
Competition and synergy within populations are also important factors, based on the intraspecific competition principle of biological populations. Competition exists within the biological population. The larger the population scale is, the more intense the competition will be. Competition among populations has the function of population size adjustment. There is also a certain competition mechanism in the IP, and this competition mechanism will suppress the excessive expansion of the IP, to some extent.
Therefore, intraspecific competition is also one of the processes for the survival of the fittest. Thus, this mechanism should be an important component of the growth model for entrepreneurial population. There is also a competition or synergy relationship between the different IPs. Based on the points raised above, this study uses the growth dynamics model of biological population theory to investigate the development characteristics of IPs.
3.2. Deduction of the Innovation Population Relationship Model
3.2.1. Dynamics Analysis of Logistic and Lotka–Volterra Models
According to the logistic model, we construct an internal relationship model of innovation population 1 (IP_{1}) as follows:
indicates the population growth rate of phase t. indicates the number of individuals in the population in phase t. Within a certain period of time (phase t), is the maximum population scale in a constant environment. Each unit occupies resources and is defined as . reflects the promotion of the population’s growth. reflects the retardation of growth due to the consumption of limited resources by the population.
If , then . The synergistic effects are dominant effects in the population. Resources within an innovation ecosystem can support an increase in the number of individuals in an IP. Thus, the growth can be sustainable.
If , then . The competition effect is dominant in the population. Innovation resources are less able to support the increase in the number of individuals in the IP. Thus, the growth is unsustainable.
According to the logistic model, this paper constructs an internal relationship model of innovation population 2 (IP_{2}) as follows:
represents the number of individuals in the innovation population 2 in period t. Researchers should consider the impact of the IP_{2} on IP_{1}. Then, the logistic model can be modified as follows:
is the influence coefficient of population 2 on population 1. If , population 2 has a synergistic effect on population 1. If , population 2 has a competitive effect on population 1. After the formation of the dependent symbiosis system, due to the promotion of population 1, the size of population 2 will also increase. The scale change of population 2 can be described as
is the influence coefficient of population 1 on population 2. If , population 1 has a synergistic effect on population 2. If , population 1 has a competitive effect on population 2. In the system of IP_{1} and IP_{2}, the symbiosis mathematical model is
Among them, , . is the contribution of population 2 to population 1, which means that the resources that population 2 supplies to population 1 are times the resources that population 2 supplies itself. According to the dependence and independence conditions, then . Similarly, we can get .
Equation (5) is called the Lotka–Volterra model. The Lotka–Volterra model is based on the logistic model of a single species and considers the dynamic growth of competition and symbiosis of two or more entities simultaneously in the ecosystem [30], which can accurately describe the competition and symbiosis between corporate populations, and It can determine the influence of the core population in the evolution of the entire ecosystem [31]; the Lotka–Volterra model has better data fitting and prediction performance [32].
The Lotka–Volterra model of dualpopulation or multipopulation growth is a differential dynamic system to simulate the dynamic relationship between populations in the innovation ecosystem. Based on the numerical value of , the type of interaction between species can be judged [32]: ① When , , it means that the innovation populations are independent, develop independently, and do not affect each other. At this time, the Lotka–Volterra model expresses no symbiotic relationship. ② When , , it means that the two parties compete with each other. One party grows while the other party declines. There is no symbiotic relationship between the two. ③ When , or , , it means that one party is attached to the other party during the symbiotic evolution of the innovative population, showing a parasitic mode of constantly requesting resources from the other party to maintain its own growth. ④ When , or , , it means that both sides of the innovative population have obtained extra highquality resources in the evolution process, but the symbiosis coefficient of one of them is zero, indicating that it has not obtained extra resources, but no losses have been suffered, and the innovation ecosystem is now in a symbiotic mode of partial benefit. ⑤ When , , it means that the innovative population is in a mutually beneficial symbiosis mode. Among them, if , it means that the symbiotic relationship between the two parties is asymmetric and mutually beneficial symbiosis; when , it means that the innovative population has obtained equal benefits in the process of symbiotic evolution, and the resources are exchanged in equal amounts, forming a symmetric and mutually beneficial symbiosis.
In innovation activities, competition can occur between populations that use common resources. Symbiosis in the innovation ecosystem does not exclude competition. Innovative populations in completely or part of the same living space need to conduct technology, talent, and market interaction in the factor market. The competition for capital and then separating, expanding, and alliance niche occupy a more favorable living position, enhance its core competitiveness, and form new offspring to adapt to the environment of the innovation ecosystem through mass reproduction. However, when one party in the innovation ecosystem relies on another core or dominant population to obtain resources and living space, a parasitic relationship is formed. Under the parasitic relationship, the symbiotic subject has a oneway exchange of interests. Because of the oneway asymmetric exchange, this state is not extensive. Therefore, the system will gradually develop in the direction of symbiosis that is conducive to mutual dependence and mutual benefit.
3.2.2. Analysis of Equilibrium Point Stability
The equilibrium point of the evolution of the innovation populations means that the output of both parties has reached the maximum and remained stable. The following uses the stability analysis of the equilibrium point to discuss the symbiosis stability of the innovative populations 1 and 2. When the two populations reach a symbiotic stable state, the differential equations can be expressed as
Solving the equations can get the equilibrium point of the symbiotic relationship between the two IPs:
The Jacobian matrix is used to solve the equilibrium point of the symbiotic evolution model of the innovation population in the innovation ecosystem. Obtain the determinant Det (J) and trace Tr (J) of the Jacobian matrix. When Det (J) > 0 and Tr (J) < 0, the local equilibrium point is in a stable state; otherwise it is not a stable equilibrium point. The following is a comparative analysis of the determinants and traces of the Jacobian matrix (Table 1). On this basis, the stable equilibrium point and stability conditions of the symbiotic evolution model of the innovative population are obtained.

As shown in Table 1, suppose and that the determinant Det (J) and trace Tr (J) can be determined according to the different values of and , and the condition that the local equilibrium point is a stable equilibrium point is obtained. Through comparison, the symbiosis evolution mode under different values corresponds to different stable equilibrium points, but the equilibrium points of partial benefit, asymmetric reciprocity, and symmetric reciprocity symbiosis are all the same P_{4}. The parasitic equilibrium point may be P_{2}, P_{3}, and P_{4} (Table 2).

As shown in Table 2, the symbiosis evolution equilibrium point of the innovation population is related to the symbiosis model and is affected by the size of the mutual competition coefficient between the innovation populations and the maximum capacity of the respective populations.
Because of the interdependence between the two populations, the population size cannot be zero, so the points P1, P2, and P3 are discarded. P4 point corresponds to the scale of populations 1 and 2, respectively . At the same time, the conditions for P4 to be meaningful are
The only nonnegative solution can be obtained by solving the equations above. That is, the equilibrium point is . This equilibrium point represents the equilibrium state of innovation resources occupied by the IP of the enterprise and the IP of the scientific research institution.
4. OutputOriented Population Size Optimization Model
4.1. Objective Planning Method
Objective programming is an effective method for solving the multiobjective programming problem. Its basic idea is to determine a desired value (objective value or ideal value) for each objective function of the multiobjective programming problem. However, due to the limitations of various conditions, these objective values are often impossible to achieve. Therefore, positive or negative deviation variables are introduced into each objective function to represent the situation where the objective value is either exceeded or not reached. To distinguish the importance of each objective, the priority and weighting coefficient of the objective are introduced. Then, constraint equations are established for all objective functions. From this new set of constraints, the scheme to minimize the combination deviation is obtained. The foundations of the objective programming model are simple and easy to understand, and the model and its hypothesis are in line with reality. Compared with other methods, the objective programming method has more flexibility, effectiveness, and convenience in use and implementation when dealing with multiobjective problems.
4.2. Multichoice Goal Programming
In recent years, multichoice goal programming (MCGP) has been widely used to resolve many practical decisionmaking problems. Chang et al. [33] integrated MCGP and fuzzy mathematics methods, according to different strategic directions of LCD and acrylic plate manufacturers. They considered the multiobjective expectation level and fuzzy relationship, which helped decision makers select the best supplier. Lee et al. [34] solved the problem of engineering technology selection in product design by combining MCGP, AHP, and QFD. Chen et al. [35] proposed a threelayer MCGP method to help forest managers obtain appropriate solutions for forest resource allocation. The multiplechoice goal programming (MCGP) method proposed by Chang [36, 37] is described as follows:
Here, indicate the value of the ith goal exceeding and not reaching the expected value of the goal. is the function of the ith object. X is the decision variable, representing m alternatives . is the expected level for the ith goal.
and are close to positive and negative deviation variables of . and are the lower and upper limits of the target for .
MCGP is a linear form of objective programming, which can be solved by some common linear programming software. This is because, in minimizing the objective function, the objective function can be infinitely close to the value of the objective. In the same way, in minimizing the objective function, the objective value can also approach the upper bound of the objective infinitely.
4.3. Multichoice Goal Programming Embedding with Lotka–Volterra Equilibrium
Embed the Lotka–Volterra model and the symbiotic population equilibrium point as constraints into the MCGP model to obtain the Lotka–VolterraMCGP model:
Multichoice goal programming embedded with Lotka–Volterra equilibrium is a linear form of objective programming, which can be solved by some common linear programming software.
5. An Example
Based on the classic Cobb Douglas production function, this paper chooses the accumulated assets of R and D investment as the main investment indicator and adds R and D human capital to the model robustness test part to measure enterprise innovation investment. For the research on enterprise innovation output, many studies use the number of patents to measure [38–40]. This article uses the number of patents to measure the output of the innovation population. However, the technical characteristics of different industries are different, and the tendency to apply for patent protection is also different. The method of using patent numbers to measure the innovation performance of enterprises is mostly used in the research of a single industry [40]. At this time, the company’s new product sales can be used to measure innovation output [41, 42]. As an auxiliary verification, this article adds new product sales to measure the output of enterprise innovation populations in the model robustness test part.
5.1. The Relationship between Sample Data and Variables
The two related innovation populations (IP_{1} and IP_{2}) have a synergistic effect between them. Variable interpretation and data selection are as follows (shown in Table 3):(1)P_{1}: population scale of IP_{1} (number of enterprises in IP_{1}, unit: number)(2)P_{2}: population scale of IP_{2} (number of enterprises in IP_{2}, unit: number)(3)E_{1}: R&D expenditure for P_{1} (unit: 100 million yuan)(4)E_{2}: R&D expenditure for P_{2} (unit: 100 million yuan)(5)O_{1}: innovation output of P_{1} (expressed by the number of granted patents)(6)O_{2}: innovation output of P_{2} (expressed by the number of granted patents)(7)E_{T}: total R&D investment (E_{T} = E_{1} + E_{2}, unit: 100 million yuan)(8)O_{T}: total innovation output (O_{T} = O_{1} + O_{2,} expressed by the number of granted patents)(9)K_{1}: maximum population scale for IP_{1}(10)K_{2}: maximum population scale for IP_{2}

We can get the IP_{1} and IP2 interpopulation relationship model (Figure 1) based on the symbiotic relationship.
Figure 2 shows the relationship model, influence path, and variable coefficients between symbiotic innovation populations. The small system composed of two populations can be regarded as a niche in the innovation ecosystem.
In the innovation ecosystem, multiple related populations affect each other. There are also interactive behaviors within the IP of enterprises. This kind of mutual influence behavior can be expressed as coordination or competition. Due to the constraints of the total resources in the innovation environment, the collaboration between populations or within populations may not necessarily promote innovation. The total resource variable is set in the model, which fully reflects the resource constraint mechanism.
5.2. Optimization Model Interpretation
In this case, the number of authorized patents is taken as a measure of innovation output. We use the objective solution of MCGP and the solution of equilibrium value to construct the suitability of population.
The related functions and parameters are listed below: (innovation output goal, the more the better) (R&D expenditure goal, the less the better)
In basis of MCGPachievement, this problem can be formulated as follows:
Considering the symbiotic relationship [43], new constraints need to be added to the above model:
The problem is solved using the LINGO [44] software and is shown in Table 4.

As seen in Table 4, there is a significant difference between the observed value and the optimized value. The population scale of IP_{1} and IP_{2} does not belong to the appropriate number of population. Innovation ecosystem can carry out selforganized evolution. Cooperation among populations and among enterprises within populations can promote a winwin situation among multiple populations. However, the negative cases in reality challenge the selforganization evolution of the innovation ecosystem. The reason lies in the lack of detailed and indepth analysis of the mechanisms at play in innovation ecosystem. The innovation ecosystem is a complex system, and its operation is inevitably affected by the interaction between the subsystems.
5.3. Evaluation
This paper evaluates the level of synergy among innovative populations in different regions based on the similarity between sample values and optimized values. The higher the similarity, the better the synergy. The evaluation matrix is A.
This study uses the entropy method to determine the weight of the two populations in the collaborative evaluation. Use the technique for order preference by similarity to an ideal solution (TOPSIS) method to evaluate the similarity, and sort the regions in the sample at the level of population coordination. The ideal scale of the two populations can be regarded as two criteria for evaluating the synergy of the innovation system.
5.3.1. The Calculation of Evaluation Weight
Entropy weight is an objective weight method. The advantage of entropy weight method reduces the subjective impact of decision makers and increases objectivity [45]. Entropy was originally a concept in thermodynamics to calculate the degree of confusion. Shannon applied it to information theory, making it one of the ways to deal with uncertainty [46]. The less the entropy value is, the more the information that can be provided will be. Therefore, the criterion can be assigned a bigger weight [47]. The concept of entropy weight has been widely used in several fields.
For example, Mohsen [48] combined entropy with fuzzy VIKOR for the risk assessment of equipment failure. Sengül et al. [49] combined entropy with fuzzy TOPSIS approach for ranking RE supply systems. Shad et al. [50] combined entropy with AHP and GIS in green building assessment. Hafezalkotob [51] integrated entropy and subjective weight for the engineering design. They compared the ranking results with other methods and found that the entropy method is suitable for application in similar problems. The calculation of Shannon’s entropy weight is presented as follows [52].
Assume m alternatives (A1, A2, Am) and n criteria (C1, C2, Cn) for a decision problem. Then initial decision matrix is Step 1: normalize the evaluation matrix: Step 2: compute entropy: Step 3: the weights of each criterion are calculated:
Calculated based on the data in the example: = 0.262, = 0.738.
5.3.2. Order Preference by Similarity to an Ideal Solution
Technique for order preference by similarity to an ideal solution (TOPSIS) is a popular method proposed by Hwang and Yoon [53] to determine the best alternative. The main rule of TOPSIS is that the best alternative should have shortest distance from the positiveideal solution and the farthest distance from the negativeideal solution [54]. This method has been widely adopted to solve problems in many different fields. The algorithm of the TOPSIS method is presented as follows. Step 1: construct the normalized decision matrix R: Step 2: construct weighted normalized decision matrix V: is the weight of the jth criterion. Step 3: determine the positiveideal solution (PIS) and negativeideal solution (NIS), denoted, respectively, as A^{+} and A^{−}, which are defined in the following way: where are sets of benefit and cost criteria, respectively. Step 4: calculate the distances of each alternative from positiveideal solution (PIS) and negativeideal solution (NIS): Step 5: calculate the closeness coefficient and rank the order of alternatives:where with . The best alternative can therefore be found according to the preference order of . The value is the more the better. If is close to 1, it indicates the alternative is closer to the PIS.
As shown in Table 5, this paper can successfully get the collaborative evaluation ranking of innovation populations in different regions.

5.4. Computational Experience
The superiority of the model can be tested by the model operation time. Sample data of each year is formulated as MCGP and Lotka–VolterraMCGP and then solved by LINGO on a PC with CPU time of 3.2 GHz. The average relative performance of MCGP and Lotka–VolterraMCGP is measured by CPU time. The CPU times of MCGP and Lotka–VolterraMCGP are both the same of 00 : 00 : 00 (hh : mm : ss). The two models have no significant difference in computing time, and both can be solved quickly.
5.5. Sensitivity Analysis of β_{12} and β_{21}
β (β_{12} and β_{21}) has significant sensitivity in the model. When β_{12} and β_{21} change in the same direction and the same amount, the model optimization result remains unchanged. The asynchronous changing of β_{12} and β_{21} significantly affects the model optimization results. The optimization result changed even more because of the changing of β_{21}. Therefore, β_{21} is more sensitive than β_{12}. The changing rate of product population 2 is significantly higher than that of product population 1 . Product population 2 is more sensitive to β.
5.6. Model Robustness Analysis
In order to test the robustness of the Lotka–VolterraMCGP model, this paper changes the relevant constraints and model parameters to obtain the following new model:
Model optimization results are shown in Table 6.

As shown in Table 6, this paper can successfully get the optimization results of innovation populations with robustness test model.
6. Discussion
Enterprise population interaction is a common topic [55]. This kind of relationship is ubiquitous in the innovation system [56]. Scholars have studied similar problems from resource constraints [57] to organizational ecology [58]. And, it is popular to use the Lotka–Volterra model to analyze the competition of innovation resource between two enterprises [59]. In this paper, Lotka–Volterra and MCGP are combined to construct a multiagent and multiobjective optimization model. Based on the optimization results, the population collaborative evaluation is carried out. This method takes into account the mathematical characteristics of the two models, integrates the functions of the two models, and expands the application fields of the two models. It is more suitable for analyzing practical problems.
7. Conclusions
This paper combines the Lotka–Volterra model and the multiobjective decision model to build an analysis path that takes both input and output perspectives into account. At the same time, the paper used optimization values to evaluate the population size suitability, which makes the evaluation more operable. The results show that the combination of Lotka–Volterra model and multichoice goal programming model can better evaluate the scale suitability of IPs within a region.
In the enterprise innovation ecosystem, multiple related populations affect each other. There are also interactive behaviors within the IP of enterprises. This kind of mutually influential behavior can be expressed as coordination or competition. Due to the constraints of the total resources in the innovation environment, the collaboration between populations or within populations may not necessarily promote innovation. The total resource variable is set in the model, which fully reflects the resource constraint mechanism. The innovation ecosystem can carry out selforganized evolution, and cooperation among populations and among enterprises within populations can promote a winwin situation among multiple populations. However, the fact that there are negative cases in reality challenges the selforganized evolution of the innovation ecosystem. The reason lies in the lack of detailed and indepth analysis of the mechanisms at play in the innovation ecosystem. The innovation ecosystem is a complex system, and its operation is inevitably affected by the interaction between the subsystems.
In this study, the different life cycles of IP development are not considered, and the demand for resources of IPs in different life cycle are different. Future studies could consider the characteristics of population life cycle development and add life cycle factors into the special analysis model.
Data Availability
The experimental data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Social Science Foundation of China (no. 20BGL203). The authors thank Professor ChingTer Chang of Chang Gung University for his suggestions on this paper.
References
 U. Lichtenthaler, “Open innovation: past research, current debates, and future directions,” Academy of Management Perspectives, vol. 25, no. 1, pp. 75–93, 2011. View at: Publisher Site  Google Scholar
 H. Garriga, G. Von Krogh, and S. Spaeth, “How constraints and knowledge impact open innovation,” Strategic Management Journal, vol. 34, no. 9, pp. 1134–1144, 2013. View at: Publisher Site  Google Scholar
 A. Razani and Z. Goodarzi, “A solution of VolterraHamerstain integral equation in partially ordered sets,” International Journal of Industrial Mathematics, vol. 3, no. 4, pp. 277–281, 2011. View at: Google Scholar
 M. R. Mokhtarzadeh, M. R. Pournaki, and A. Razani, “An existenceuniqueness theorem for a class of boundary value problems,” Fixed Point Theory, vol. 13, no. 2, pp. 583–592, 2012. View at: Google Scholar
 M. R. Mokhtarzadeh, M. R. Pournaki, and A. Razani, “A note on periodic solutions of Riccati equations,” Nonlinear Dynamics, vol. 62, no. 12, pp. 119–125, 2010. View at: Publisher Site  Google Scholar
 M. R. Pournaki and A. Razani, “On the existence of periodic solutions for a class of generalized forced Liénard equations,” Applied Mathematics Letters, vol. 20, no. 3, pp. 248–254, 2007. View at: Publisher Site  Google Scholar
 R. Pouder and C. H. St. John, “Hot spots and blind spots: geographical clusters of firms and innovation,” Academy of Management Review, vol. 21, no. 4, pp. 1192–1225, 1996. View at: Publisher Site  Google Scholar
 A. Cuspilici, P. Monforte, and M. A. Ragusa, “Study of Saharan dust influence on PM 10 measures in Sicily from 2013 to 2015,” Ecological Indicators, vol. 76, pp. 297–303, 2017. View at: Publisher Site  Google Scholar
 K. Head, J. Ries, and D. Swenson, “Agglomeration benefits and location choice: evidence from Japanese manufacturing investments in the United States,” Journal of International Economics, vol. 38, no. 34, pp. 223–247, 1995. View at: Publisher Site  Google Scholar
 T. Frost and C. Zhou, “The geography of foreign R & D within a host country,” International Studies of Management & Organization, vol. 30, no. 2, pp. 10–43, 2000. View at: Publisher Site  Google Scholar
 J. A. C. Baum and S. J. Mezias, “Localized competition and organizational failure in the Manhattan hotel industry, 1898–1990,” Administrative Science Quarterly, vol. 37, no. 4, pp. 580–604, 1992. View at: Publisher Site  Google Scholar
 J. Wang and J. Wang, “An analysis of newtech agglomeration in Beijing: a new industrial district in the making?” Environment and Planning A: Economy and Space, vol. 30, no. 4, pp. 681–701, 1998. View at: Publisher Site  Google Scholar
 M. T. Hannan and J. Freeman, “The population ecology of organizations,” American Journal of Sociology, vol. 82, no. 5, pp. 929–964, 1977. View at: Publisher Site  Google Scholar
 S. A. Morris and D. Pratt, “Analysis of the LotkaVolterra competition equations as a technological substitution model,” Technological Forecasting and Social Change, vol. 70, no. 2, pp. 103–133, 2003. View at: Publisher Site  Google Scholar
 C. Watanabe, R. Kondo, N. Ouchi, and H. Wei, “A Substitution orbit model of competitive innovations,” Technological Forecasting and Social Change, vol. 71, no. 4, pp. 365–390, 2004. View at: Publisher Site  Google Scholar
 S.J. Lee, D.J. Lee, and H.S. Oh, “Technological forecasting at the Korean stock market: a dynamic competition analysis using LotkaVolterra model,” Technological Forecasting and Social Change, vol. 72, no. 8, pp. 1044–1057, 2005. View at: Publisher Site  Google Scholar
 J. Kim, D.J. Lee, and J. Ahn, “A dynamic competition analysis on the Korean mobile phone market using competitive diffusion model,” Computers & Industrial Engineering, vol. 51, no. 1, pp. 174–182, 2006. View at: Publisher Site  Google Scholar
 V. B. Kreng and H. T. Wang, “The interaction of the market competition between LCD TV and PDP TV,” Computers & Industrial Engineering, vol. 57, no. 4, pp. 1210–1217, 2009. View at: Publisher Site  Google Scholar
 L. A. Gomes, A. L. F. Facin, M. S. Salerno, and R. K. Ikenami, “Unpacking the innovation ecosystem construct: evolution, gaps and trends,” Technological Forecasting & Social Change, vol. 136, pp. 30–48, 2016. View at: Google Scholar
 D.S. Oh, F. Phillips, S. Park, and E. Lee, “Innovation ecosystems: a critical examination,” Technovation, vol. 54, pp. 1–6, 2016. View at: Publisher Site  Google Scholar
 J. A. C. Baum, H. J. Korn, and S. Kotha, “Dominant designs and population dynamics in telecommunications services: founding and failure of facsimile transmission service organizations, 1965–1992,” Social Science Research, vol. 24, no. 2, pp. 97–135, 1995. View at: Publisher Site  Google Scholar
 M. B. Low and E. Abrahamson, “Movements, bandwagons, and clones: industry evolution and the entrepreneurial process,” Journal of Business Venturing, vol. 12, no. 6, pp. 435–457, 1997. View at: Publisher Site  Google Scholar
 P. A. Geroski and M. Mazzucato, “Modelling the dynamics of industry populations,” International Journal of Industrial Organization, vol. 19, no. 7, pp. 1003–1022, 2001. View at: Publisher Site  Google Scholar
 M. T. Hannan, G. R. Carroll, E. A. Dundon, and J. C. Torres, “Organizational evolution in a multinational context: entries of automobile manufacturers in Belgium, Britain, France, Germany, and Italy,” American Sociological Review, vol. 60, no. 4, pp. 509–528, 1995. View at: Publisher Site  Google Scholar
 M. T. Hannan, “Inertia, density and the structure of organizational populations: entries in European automobile industries, 1886–1981,” Organization Studies, vol. 18, no. 2, pp. 193–228, 1997. View at: Google Scholar
 G. Adomavicius, J. C. Bockstedt, A. Gupta, and R. J. Kauffman, “Technology roles and paths of influence in an ecosystem model of technology evolution,” Information Technology and Management, vol. 8, no. 2, pp. 185–202, 2007. View at: Publisher Site  Google Scholar
 R. Adner and R. Kapoor, “Value creation in innovation ecosystems: how the structure of technological interdependence affects firm performance in new technology generations,” Strategic Management Journal, vol. 31, no. 3, pp. 306–333, 2010. View at: Publisher Site  Google Scholar
 B. Clarysse, M. Wright, J. Bruneel, and A. Mahajan, “Creating value in ecosystems: crossing the chasm between knowledge and business ecosystems,” Research Policy, vol. 43, no. 7, pp. 1164–1176, 2014. View at: Publisher Site  Google Scholar
 R. Kapoor and N. R. Furr, “Complementarities and competition: unpacking the drivers of entrants’ technology choices in the solar photovoltaic industry,” Strategic Management Journal, vol. 36, no. 3, pp. 416–436, 2015. View at: Publisher Site  Google Scholar
 V. Volterra, “Fluctuations in the abundance of a species considered mathematically^{1},” Nature, vol. 118, no. 2972, pp. 558–560, 1926. View at: Publisher Site  Google Scholar
 G. L. Zhang, A. Daniel, and V. Mc Adams, “Technology evolution prediction using LotkaVolterra Equations,” Journal of Mechanical Design, vol. 140, no. 6, pp. 61–101, 2018. View at: Publisher Site  Google Scholar
 T. Modis, “Technological forecasting at the stock market,” Technological Forecasting and Social Change, vol. 62, no. 3, pp. 173–202, 1999. View at: Publisher Site  Google Scholar
 C.T. Chang, C.Y. Ku, and H.P. Ho, “Fuzzy multichoice goal programming for supplier selection,” International Journal of Operations Research and Information Systems, vol. 1, no. 3, pp. 28–52, 2010. View at: Publisher Site  Google Scholar
 A. H. I. Lee, H.Y. Kang, C.Y. Yang, and C.Y. Lin, “An evaluation framework for product planning using FANP, QFD and multichoice goal programming,” International Journal of Production Research, vol. 48, no. 13, pp. 3977–3997, 2010. View at: Publisher Site  Google Scholar
 Y.T. Chen, C. Zheng, and C.T. Chang, “3level MCGP: an efficient algorithm for MCGP in solving multiforest management problems,” Scandinavian Journal of Forest Research, vol. 26, no. 5, pp. 457–465, 2011. View at: Publisher Site  Google Scholar
 C.T. Chang, “Multichoice goal programming,” Omega, vol. 35, no. 4, pp. 389–396, 2007. View at: Publisher Site  Google Scholar
 C.T. Chang, “Revised multichoice goal programming,” Applied Mathematical Modelling, vol. 32, no. 12, pp. 2587–2595, 2008. View at: Publisher Site  Google Scholar
 G. Ahuja, “Collaboration networks, structural holes, and innovation: a longitudinal study, structural holes and innovation: a longitudinal study,” Administrative Science Quarterly, vol. 45, no. 3, pp. 425–455, 2000. View at: Publisher Site  Google Scholar
 R. C. Sampson, “R & D alliances and firm performance: the impact of technological diversity and alliance organization on innovation: the impact of technological diversity and alliance organization on innovation,” Academy of Management Journal, vol. 50, no. 2, pp. 364–386, 2007. View at: Publisher Site  Google Scholar
 B.J. Park, M. K. Srivastava, and D. R. Gnyawali, “Walking the tight rope of coopetition: impact of competition and cooperation intensities and balance on firm innovation performance: impact of competition and cooperation intensities and balance on firm innovation performance,” Industrial Marketing Management, vol. 43, no. 2, pp. 210–221, 2014. View at: Publisher Site  Google Scholar
 A. Fosfuri and J. Tribó, “Exploring the antecedents of potential absorptive capacity and its impact on innovation performance,” Omega, vol. 36, no. 2, pp. 173–187, 2008. View at: Publisher Site  Google Scholar
 Z.L. He and P.K. Wong, “Exploration vs. exploitation: an empirical test of the ambidexterity hypothesis,” Organization Science, vol. 15, no. 4, pp. 481–494, 2004. View at: Publisher Site  Google Scholar
 S.Y. Wang, W.M. Chen, and Y. Liu, ““Collaborative product portfolio design based on the approach of multi choice goal programming,” Mathematical Problems in Engineering, vol. 2021, Article ID 6678533, 16 pages, 2021. View at: Publisher Site  Google Scholar
 L. Schrage, LINGO Release 8.0, LINDO System, Inc, Chicago, IL, USA, 2002.
 H.C. Lee and C.T. Chang, “Comparative analysis of MCDM methods for ranking renewable energy sources in Taiwan,” Renewable and Sustainable Energy Reviews, vol. 92, pp. 883–896, 2018. View at: Publisher Site  Google Scholar
 Z.H. Zou, Y. Yun, and J.N. Sun, “Entropy method for determination of weight of evaluating indicators in fuzzy synthetic evaluation for water quality assessment,” Journal of Environmental Sciences, vol. 18, no. 5, pp. 1020–1023, 2006. View at: Publisher Site  Google Scholar
 J. Ye, “Fuzzy decisionmaking method based on the weighted correlation coefficient under intuitionistic fuzzy environment,” European Journal of Operational Research, vol. 205, no. 1, pp. 202–204, 2010. View at: Publisher Site  Google Scholar
 O. Mohsen and N. Fereshteh, “An extended VIKOR method based on entropy measure for the failure modes risk assessmenta case study of the geothermal power plant (GPP),” Safety Science, vol. 92, pp. 160–172, 2017. View at: Publisher Site  Google Scholar
 Ü. Şengül, M. Eren, S. S. Eslamian, V. Gezder, and A. B. Şengül, “Fuzzy TOPSIS method for ranking renewable energy supply systems in Turkey”,” Renew Energy, vol. 75, pp. 617–625, 2015. View at: Google Scholar
 R. Shad, M. Khorrami, and M. Ghaemi, “Developing an Iranian green building assessment tool using decision making methods and geographical information system: case study in Mashhad city,” Renewable and Sustainable Energy Reviews, vol. 67, pp. 324–340, 2017. View at: Publisher Site  Google Scholar
 A. Hafezalkotob and A. Hafezalkotob, “Extended multimoora method based on Shannon entropy weight for materials selection,” Journal of Industrial Engineering International, vol. 12, no. 1, pp. 1–13, 2016. View at: Publisher Site  Google Scholar
 F. H. Lotfi and R. Fallahnejad, “Imprecise Shannon’s entropy and multi attribute decision making,” Entropy, vol. 12, no. 1, pp. 53–62, 2010. View at: Publisher Site  Google Scholar
 C.L. Hwang and K. Yoon, “Methods for multiple attribute decision making,” in Multiple Attribute Decision Making, pp. 58–191, Springer, New York, NY, USA, 1981. View at: Publisher Site  Google Scholar
 N. Chitsaz and M. E. Banihabib, “Comparison of different multi criteria decisionmaking models in prioritizing flood management alternatives,” Water Resources Management, vol. 29, no. 8, pp. 2503–2525, 2015. View at: Publisher Site  Google Scholar
 R. Cerqueti, F. Tramontana, and M. Ventura, “On the coexistence of innovators and imitators,” Technological Forecasting and Social Change, vol. 90, pp. 487–496, 2015. View at: Publisher Site  Google Scholar
 A. S. Chakrabarti, “Stochastic LotkaVolterra equations: a model of lagged diffusion of technology in an interconnected world,” Physica A: Statistical Mechanics and its Applications, vol. 442, pp. 214–223, 2016. View at: Publisher Site  Google Scholar
 S. Baskaran and K. Mehta, “What is innovation anyway? youth perspectives from resourceconstrained environments,” Technovation, vol. 5253, pp. 4–17, 2016. View at: Publisher Site  Google Scholar
 Y. Cui, J. Jiao, and H. Jiao, “Technological innovation in Brazil, Russia, India, China, and South Africa (BRICS): an organizational ecology perspective,” Technological Forecasting and Social Change, vol. 107, pp. 28–36, 2016. View at: Publisher Site  Google Scholar
 T. Wei, Z. Zhu, Y. Li, and N. Yao, “The evolution of competition in innovation resource: a theoretical study based on LotkaVolterra model,” Technology Analysis & Strategic Management, vol. 30, no. 3, pp. 295–310, 2018. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 ShengYuan Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.