Discrete Dynamics in Nature and Society

Discrete Dynamics in Nature and Society / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 5569108 | https://doi.org/10.1155/2021/5569108

Sheng-Yuan Wang, Wan-Ming Chen, Rong Wang, Xiao-Lan 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

Academic Editor: Maria Alessandra Ragusa
Received06 Jan 2021
Accepted27 May 2021
Published08 Jun 2021

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 three-stage 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 [16]. 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–Volterra-MCGP 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 Abra-hamson 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 density-dependent 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 above-mentioned 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 high-tech 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 technology-based 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 incentive-based 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 (IP1) 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 (IP2) as follows:

represents the number of individuals in the innovation population 2 in period t. Researchers should consider the impact of the IP2 on IP1. 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 IP1 and IP2, 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 dual-population 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 high-quality 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 one-way exchange of interests. Because of the one-way 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.


Equilibrium pointDet (J)Tr (J)Stability conditions

Unstable

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 P4. The parasitic equilibrium point may be P2, P3, and P4 (Table 2).


Influencing factor valueSymbiotic relationshipStable equilibrium point

Parasitic relationshipP2, P3, P4
Favor symbiosisP4
Asymmetric symbiosisP4
Symmetry symbiosisP4

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. Output-Oriented 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 decision-making 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 three-layer MCGP method to help forest managers obtain appropriate solutions for forest resource allocation. The multiple-choice goal programming (MCGP) method proposed by Chang [36, 37] is described as follows:

Here, indicate the value of the i-th goal exceeding and not reaching the expected value of the goal. is the function of the i-th object. X is the decision variable, representing m alternatives . is the expected level for the i-th 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–Volterra-MCGP 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 [3840]. 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 (IP1 and IP2) have a synergistic effect between them. Variable interpretation and data selection are as follows (shown in Table 3):(1)P1: population scale of IP1 (number of enterprises in IP1, unit: number)(2)P2: population scale of IP2 (number of enterprises in IP2, unit: number)(3)E1: R&D expenditure for P1 (unit: 100 million yuan)(4)E2: R&D expenditure for P2 (unit: 100 million yuan)(5)O1: innovation output of P1 (expressed by the number of granted patents)(6)O2: innovation output of P2 (expressed by the number of granted patents)(7)ET: total R&D investment (ET = E1 + E2, unit: 100 million yuan)(8)OT: total innovation output (OT = O1 + O2, expressed by the number of granted patents)(9)K1: maximum population scale for IP1(10)K2: maximum population scale for IP2


AreaP1E1O1K1P2E2O2K2OTET

A11918616101578871644811904172221412621801012027
A21887214561664451487411133452158610441880311801
A314150132813196613567998325148439831468091653
A412283123917278712658944248139387501867251487
A511133108018622011033909208123626291985821288
A67712964118919984879510894063271283251072
A722575257178153637143336726100878507858
A82159417469764260713300361890850594717
A92508684239256988693113242934226354797
A102236571146165833632115182334816439686
A1118314017835409765829016008789435691
A1216953245436331070026910788146514593

We can get the IP1 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 MCGP-achievement, 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.


AreaSample observationsMCGP modelL-MCGP model
IP1IP2MIP1MIP2LMIP1LMIP2β12β21

A11918611901506920411827514020.100.10
A2188721113175465061681511620.100.10
A3141509981227916671507410920.100.10
A412283944179870170378500.300.10
A511133909191290183107330.560.10
A67712795123610119523660.190.10
A722577146716758630411260.100.10
A8215971333841334473310080.100.10
A925086930227377643800.100.10
A1022366320141864813860.100.10
A111831658081345529750.100.10
A121695700056136779040.100.10

As seen in Table 4, there is a significant difference between the observed value and the optimized value. The population scale of IP1 and IP2 does not belong to the appropriate number of population. Innovation ecosystem can carry out self-organized evolution. Cooperation among populations and among enterprises within populations can promote a win-win situation among multiple populations. However, the negative cases in reality challenge the self-organization evolution of the innovation ecosystem. The reason lies in the lack of detailed and in-depth 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 isStep 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 positive-ideal solution and the farthest distance from the negative-ideal 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 j-th criterion.Step 3: determine the positive-ideal solution (PIS) and negative-ideal 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 positive-ideal solution (PIS) and negative-ideal 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.


AreaVi1Vi2Rank

A10.0110.0350.1650.0830.3344
A20.0240.0080.0140.0880.8661
A30.0110.0160.0870.0930.5172
A40.0560.0160.0980.0630.3913
A50.0850.0290.1630.0420.2058
A60.0500.0710.2550.0350.12012
A70.0480.0690.2490.0370.13010
A80.0310.0490.2030.0590.2256
A90.0620.0520.2160.0300.12211
A100.0500.0410.1850.0460.2009
A110.0320.0530.2120.0560.2097
A120.0240.0340.1610.0720.3095

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–Volterra-MCGP and then solved by LINGO on a PC with CPU time of 3.2 GHz. The average relative performance of MCGP and Lotka–Volterra-MCGP is measured by CPU time. The CPU times of MCGP and Lotka–Volterra-MCGP 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–Volterra-MCGP model, this paper changes the relevant constraints and model parameters to obtain the following new model:

Model optimization results are shown in Table 6.


AreaP1P2O′TETHTK′1K′2LMP′1LMP′2

A1191861190308922027560464230231547181341218
A2188721113269091801547288226461447178471140
A314150998258941653415350169801297131531004
A412283944216851487386207147401227120401002
A51113390919630128835620713360118210974970
A67712795165101072269328925410347885881
A722577141032685810930827089282142734
A8215971380257179499625919271821652
A92508693724879711035230109012306690
A10223663254956869838426838222034623
A11183165835976918056421978551731673
A12169570029475939399020349101602716

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 self-organized evolution, and cooperation among populations and among enterprises within populations can promote a win-win situation among multiple populations. However, the fact that there are negative cases in reality challenges the self-organized evolution of the innovation ecosystem. The reason lies in the lack of detailed and in-depth 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 Ching-Ter Chang of Chang Gung University for his suggestions on this paper.

References

  1. 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
  2. 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
  3. A. Razani and Z. Goodarzi, “A solution of Volterra-Hamerstain integral equation in partially ordered sets,” International Journal of Industrial Mathematics, vol. 3, no. 4, pp. 277–281, 2011. View at: Google Scholar
  4. M. R. Mokhtarzadeh, M. R. Pournaki, and A. Razani, “An existence-uniqueness theorem for a class of boundary value problems,” Fixed Point Theory, vol. 13, no. 2, pp. 583–592, 2012. View at: Google Scholar
  5. M. R. Mokhtarzadeh, M. R. Pournaki, and A. Razani, “A note on periodic solutions of Riccati equations,” Nonlinear Dynamics, vol. 62, no. 1-2, pp. 119–125, 2010. View at: Publisher Site | Google Scholar
  6. 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
  7. 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
  8. 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
  9. 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. 3-4, pp. 223–247, 1995. View at: Publisher Site | Google Scholar
  10. 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
  11. 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
  12. J. Wang and J. Wang, “An analysis of new-tech 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
  13. 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
  14. S. A. Morris and D. Pratt, “Analysis of the Lotka-Volterra 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
  15. 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
  16. S.-J. Lee, D.-J. Lee, and H.-S. Oh, “Technological forecasting at the Korean stock market: a dynamic competition analysis using Lotka-Volterra model,” Technological Forecasting and Social Change, vol. 72, no. 8, pp. 1044–1057, 2005. View at: Publisher Site | Google Scholar
  17. 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
  18. 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
  19. 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
  20. 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
  21. 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
  22. 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
  23. 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
  24. 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
  25. 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
  26. 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
  27. 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
  28. 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
  29. 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
  30. V. Volterra, “Fluctuations in the abundance of a species considered mathematically1,” Nature, vol. 118, no. 2972, pp. 558–560, 1926. View at: Publisher Site | Google Scholar
  31. G. L. Zhang, A. Daniel, and V. Mc Adams, “Technology evolution prediction using Lotka-Volterra Equations,” Journal of Mechanical Design, vol. 140, no. 6, pp. 61–101, 2018. View at: Publisher Site | Google Scholar
  32. 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
  33. C.-T. Chang, C.-Y. Ku, and H.-P. Ho, “Fuzzy multi-choice 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
  34. A. H. I. Lee, H.-Y. Kang, C.-Y. Yang, and C.-Y. Lin, “An evaluation framework for product planning using FANP, QFD and multi-choice goal programming,” International Journal of Production Research, vol. 48, no. 13, pp. 3977–3997, 2010. View at: Publisher Site | Google Scholar
  35. Y.-T. Chen, C. Zheng, and C.-T. Chang, “3-level MCGP: an efficient algorithm for MCGP in solving multi-forest management problems,” Scandinavian Journal of Forest Research, vol. 26, no. 5, pp. 457–465, 2011. View at: Publisher Site | Google Scholar
  36. C.-T. Chang, “Multi-choice goal programming,” Omega, vol. 35, no. 4, pp. 389–396, 2007. View at: Publisher Site | Google Scholar
  37. C.-T. Chang, “Revised multi-choice goal programming,” Applied Mathematical Modelling, vol. 32, no. 12, pp. 2587–2595, 2008. View at: Publisher Site | Google Scholar
  38. 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
  39. 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
  40. 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
  41. 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
  42. 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
  43. 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
  44. L. Schrage, LINGO Release 8.0, LINDO System, Inc, Chicago, IL, USA, 2002.
  45. 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
  46. 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
  47. J. Ye, “Fuzzy decision-making 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
  48. O. Mohsen and N. Fereshteh, “An extended VIKOR method based on entropy measure for the failure modes risk assessment-a case study of the geothermal power plant (GPP),” Safety Science, vol. 92, pp. 160–172, 2017. View at: Publisher Site | Google Scholar
  49. Ü. Ş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
  50. 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
  51. 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
  52. 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
  53. 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
  54. N. Chitsaz and M. E. Banihabib, “Comparison of different multi criteria decision-making models in prioritizing flood management alternatives,” Water Resources Management, vol. 29, no. 8, pp. 2503–2525, 2015. View at: Publisher Site | Google Scholar
  55. 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
  56. A. S. Chakrabarti, “Stochastic Lotka-Volterra 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
  57. S. Baskaran and K. Mehta, “What is innovation anyway? youth perspectives from resource-constrained environments,” Technovation, vol. 52-53, pp. 4–17, 2016. View at: Publisher Site | Google Scholar
  58. 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
  59. T. Wei, Z. Zhu, Y. Li, and N. Yao, “The evolution of competition in innovation resource: a theoretical study based on Lotka-Volterra model,” Technology Analysis & Strategic Management, vol. 30, no. 3, pp. 295–310, 2018. View at: Publisher Site | Google Scholar

Copyright © 2021 Sheng-Yuan 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.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views167
Downloads263
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.