A Dynamical Innovation Diffusion Model with Fuzzy Coefficient and Its Application to Local Telephone Diffusion in China

This paper studies the innovation diffusion problemwith the affection of urbanization, proposing a dynamical innovation diffusion model with fuzzy coefficient, and uses the shifting rate of people from rural areas stepping into urban areas to show the process of urbanization. The numerical simulation shows the diffusion process for telephones in China with Genetic Algorithms and this model is effective to show the process of innovation diffusion with the condition of urbanization process.


Introduction
Technological innovation promotes the evolution of industries by altering the competitive market structure and value chains of industries [1].The market success of an innovation is determined not only by technological performance, but also by the interaction of numerous factors.Consequently, it is understandable that many studies choose to investigate the key factors that influence the acceptance or the rejection of an innovation during the process of diffusion.Based on these key factors including price and advertising, researchers have contributed to the development of the diffusion theory by suggesting analytical models for describing and forecasting the diffusion process of an innovation in a social system [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16] such researchers include Fourt, Mansfield, and Mahajan.Fourt and Mansfield discussed the diffusion pattern in external influence such as mass advertisement or internal influence such as oral communication separately [8,12].The sales model for forecasting TV made by Bass in 1969 has settled the base for the following research on innovation diffusion with both external factors and internal factors taken into account [3].However, the potential adopters in these models are static or fixed at the time an innovation is introduced and remains constant over the diffusion process.Obviously, such an assumption is not tenable with regard to either theory or practice.
In response to this shortcoming, Mahajan and Peterson proposed a dynamical diffusion model where the potential is permitted to vary over time [13].There are many relevant factors affecting the potential, including socioeconomic conditions in the social system, increases or decreases in the population of the social system, government actions, and efforts to influence the diffusion process, such as advertising and pricing.Such a proposition was followed and developed by Lackman, Dodson, and several other researchers [17][18][19][20][21][22][23][24].Kalish, for example, conjectured that the qualitative structure of an optimal pricing policy will remain the same upon the inclusion of additional marketing variables such as advertising, in the diffusion model.Jain and Rao show that price affects the diffusion rate via the coefficients of external and internal influence [20].
Because different adopters have different attitudes of risk or decisive patterns, there will be great difference in adoption possibility and rate.The consumption habit, the living standard, and the environment around also make a great difference in the adopting rate.For this reason, a diffusion model is presented in a different section in the paper [25].The author believes that the characteristics of adopters in different patches make differences in the diffusion rate.Furthermore, there will be migration between patches, which would also have an effect on innovation diffusion.With great economic development and urbanization, we

Conversion name/way of conversion
Explanation for conversion (A) External information, oral communication non-users become users in colony 1 (B) External information, oral communication non-users become users in colony 2 (C) Give up using Users become Non-users in colony 1 (D) Give up using Users become Non-users in colony 2 (E) The members shifting between the colonies Users in colony 1 continue using the innovation after shifting into colony 2 (F) The members shifting between the colonies Users in colony 2 continue using the innovation after shifting into colony 1 (G) The members shifting between the colonies Users in colony 2 give up using the innovation after shifting into colony 1 (H) The members shifting between the colonies Users in colony 1 give up using the innovation after shifting into colony 2 (M) The members shifting between the colonies Non-users in colony 2 still remain as non-users after shifting into colony 1 (N) The members shifting between the colonies Non-users in colony 1 still remain as non-users after shifting into colony 2 should not ignore the potential adopters who have changed their economic conditions or status.They will have great consumption demands for goods which they would not have even considered before.For example, the data of telephone subscribers of China in 2002-2006 has got tremendous growth which cannot be explained by any diffusion model before.Furthermore, the rate of urbanization development is different in different stages.So it is more acceptable if we take the uncertain rate instead of the deterministic rate to simulate the process.We use the rate of people from rural areas stepping into urban to show the process of the urbanization development of this area.Consequently, this paper constructs a fuzzy innovation diffusion model and tries to find the principles of innovation diffusion with the effects of population increasing and urbanization.In this context, model uncertainty is portrayed through the fuzzy transition process from nonuser to the adopter of an innovation because of the uncertain rate of urbanization.We organize this paper as follows.In Section 2, an innovation diffusion model with the population increasing and the conversion between two different colonies is constructed.In Section 3, a stability analysis is done on the population increasing and the diffusion rate of different groups.An empirical analysis is also done with the data of telephone subscribers of urban and rural areas in China compared with the Bass model in Section 4. Some concluding remarks are finally given in Section 5.

Modeling
There is a comparative independence and stability between colony 1 and colony 2. Each of the two colonies has its own main communication styles, broadcasting channels and ways, and consumption habits.However, there also exists migration between them for specific reasons such as improved living standards.The members moving into new community will be accepted into the new consumption system since they are affected by the new environment.In the interior of a community, for some related external factors and oral communication, nonusers adopt innovation and serve as users, or users give up the innovation for some reasons after experiencing it.

Conceptional Model.
Based on the foundation above, the concrete convert relations and relations between the inside and outside of the two communities are shown in Figure 1 and Table 1, respectively.The arrowheads denote the flow directions of the members.

Mathematical Model
Hypothesis 1.The members of each colony are divided into users and nonusers.The nonusers include those who adopted the innovation and have given up using it.
Hypothesis 2. The internal increasing rate of the two colonies obeys the logistic distribution functions.
Hypothesis 3. The nonusers of each colony are still nonusers when entering another colony, which can be shown in Figure 1 in the conversion way of  and , and the users of each colony can either keep on or give up using the innovation for various reasons when entering another colony.
Hypothesis 4. The number of the members in one colony shifting to another colony is proportional to the number of the total members of the colony, and the coefficient is fuzzy for the uncertainty of urbanization.
Based on the conceptional model, the hypothesis, and the conversion relations above, we build the increase equation for variables  1 ,  2 , where   () represents the number of members in colony  in time ,  = 1, 2. According to the empirical statistics, the birthrate and ratio to death are linear functions of the population of the colony.Suppose that  0 + (− 0   ) is the birthrate of colony ,  0 is the birthrate of colony  without a resource limit, and − 0   is the birthrate of colony  with a resource limit, where  0 and  0 are both positive constants.The ratio of death of colony  is  1 + 1   > 0, where  1 is the natural death ratio without a resource limit and  0   is the ratio of death increasing with a resource limit; that is to say, the birthrate will decrease with the increase of population and the death ratio will increase with the increase of population.Thus, the number of members increased in time  in colony 1 is ( 10 −  11 ) 1 − ( 10 +  11 ) 2  1 and  10 −  11 > 0 according to Hypothesis 2. From Hypothesis 4 and the conversion way of  and , the number of increased members shifting from colony 2 in time  is  2  2 +  2  2 =  2  2 , and the number of members shifting into colony 2 from conversion way of  and  is  1  1 +  1  1 =  1  1 , as seen in Figure 1.Based on the analysis above, the increase equation of  1 is Similarly, the increase equation of  2 is .
Consequently the increasing model of colonies' members is governed by Suppose that   () is the number of users in colony  at time , and   () is the number of nonusers in colony  at time .Now we will build the equation of variable  1 () in colony 1.From conversion way  in Figure 1, because of the internal oral communication and medium, the number of nonusers who have changed to become users is  1  1 +  1  1  1 , in which  1 is the probability for nonusers to become users caused by medium, and  1 is the probability for nonusers to become users caused by oral communication between users and nonusers.Owing to natural deaths, the decreasing number of users is ( 11 +  11  1 ) 1 .From the conversion way  in Figure 1, there are  1  1 users who have to give up using the innovation for some reason, because of movement of population, there are  1  1 users shifting out of colony 1 and  2  2  2 users shifting into colony 1 from colony 2.Here   is the probability for the members of colony  converting into another one and   is the probability for users to continue using the innovation in colony  after converting into another colony.Based on the analysis above, the increase equation of variable  1 is .
Similarly, the increase equation of variable  2 is .
Consequently, we can say that the dynamical equations of innovation diffusion based on members shifting between colony 1 and colony 2 are governed by This model focuses on the analysis of the affection on the customers by the process of urbanization.The urbanization process is expressed by the mutual conversion between various consumer groups.That is to say, with economic development and improved living standards, the low-consumption groups constantly improve their levels of consumption and step into consumer groups of higher levels.From the right expression of the following equation: .
( 10 −  11 ) 1 − ( 10 +  11 ) 2 1 is fatally affected by natural population increase, and the policy of family planning, and little affected by economic development, while  2  2 −  1  1 is mainly affected by economic development.Although certain regularities of economic development can be followed, the effect of this on the customers is uncertain and will vary during the diffusion process of an innovation.However, the models mentioned in the introduction, including system (6), are crisp and obviously cannot meet the actual changes in the law of development.Since these models stressed some main influencing factors, the impact of which on the consumers might be exaggerated, and the combined effects of other factors were ignored, some serious disturbance might occur.Moreover, during the process of diffusion, the impact on the diffusion rate by these factors is of great difference in different stages.For this reason, we try to define   ,  = 1, 2 as fuzzy numbers and use fuzzy theory to express the uncertainty of urbanization of China.

The Stability Analysis on the Model
Since the market is attracted to the variety of innovation users' size and the market potential, we only discuss the variables of members  1 ,  2 and the variables of users  1 ,  2 .Before discussing the stability of the members increasing model, we will introduce some related concepts and then make the stability analysis on the innovation diffusion.

Related Concepts.
We now recall some definitions needed through the paper.
Let   denote the set of upper semicontinuous normal fuzzy sets on   with compact support.That is, if  ∈   , then  :   → [0, 1] is upper semicontinuous, supp() =  ∈   : () > 0 is compact, and there exists at least one  ∈ supp() for which () = 1.The -set of , 0 <  ≤ 1, is Definition 1.If for ,  ∈   , there exists  ∈   such that  =  + , then we say that the Hukuhara difference of  and  exists, call  the -difference of  and , and denote  −  = .
The approach of Hukuhara differentiation suffers a grave disadvantage in so far as the solution has the property that diam (()) in nondecreasing in .Consequently, this formulation cannot really reflect any of the reach behaviors of ordinary differential equations such as periodicity, stability, bifurcation, and the like and is ill suited for modeling.However, Hüllermeier suggested a different formulation of the FIVP based on a family of differential inclusions at each  level, 0 ≤  ≤ 1, where now [(⋅, ⋅)]  :  ×   →    , the space of nonempty convex compact subsets of   .The idea is that the set of all such solution Σ  ( 0 , ) would be the  level of a fuzzy set Σ( 0 , ) in the sense that all attainable sets   ( 0 , ), 0 <  ≤ , are levels of a fuzzy set ( 0 , ) on   .Considering Σ( 0 , ) to be the solution of the fuzzy DE   = (, ), (0) =  0 , thus captures both uncertainty and the rich properties of differential inclusions in one and the same technique.It has been shown that the solution set and attainability set are fuzzy sets under fairly relaxed conditions on  [26].
Theorem 3 (see [26]).Let  0 ∈ Then the attainable sets   (  , ),  ∈ [0, 1] of the family of inclusions are the level sets of a fuzzy set ( 0 , ) ⊂   .The solution sets Σ  (  , ) of ( 14) are the level sets of a fuzzy set Σ( 0 , ) defined on   (  ), where Property 2 (Lyapunov Stability of a Family of Differential Inclusion, [26]).If  ∈   is a fuzzy set and U, W ⊂   are closed subsets of   , define the distance from W and Hausdorff separation, respectively, by The significance of these definitions is that the metric space (  ,  ∞ ) of fuzzy sets,  * (, W), is the distance of  ∈   from W ⊂   and is the analogue of (, ) in   .Correspondingly,   (U, W) is the Hausdorff separation between U, W ⊂   with respect to the metric  ∞ and is the analogue of the Hausdorff separation (, ) in   .
The existence of the unstable manifold Ȗ of ( 23) is established in exactly the same way by considering the system (23) with  → −; that is, The stable manifold for this system will then be the unstable manifold Ȗ for (23).Note that it is also necessary to replace the vector  by the vector ( +1 , . . .,   ,  1 , . . .,   ) in order to determine the  − -dimensional manifold Ȗ by the above process.This completes the proof.
If the original point is replaced with the equilibrium point , fussify it and the  level set of  is called equilibrium core, denoted by   .Now we prove that ∀ * ∈   such as Theorem 4; that is, ( * ) = 0 and ( * ) has  eigenvalues with negative real part and  −  eigenvalues with positive real part, then there exists a -dimensional differentiable manifold  tangent to the stable subspace   of the linear differential inclusions (16) at  * such that for all  ≥ 0,   () ⊂  and for all then the equilibrium core   is stable, and there exists an  − -dimensional differentiable manifold  tangent to unstable subspace   of ( 16) at  * such that for all  ≤ 0,   () ⊂  and for all  0 ∈ , then the equilibrium core   is stable.For all  * ∈   corresponding to   , lim That is, ∃ 0 , if  >  0 ,  * ( * ,   ( 0 , )) < ,   (  ,   ( 0 , )) ⩽  * ( * ,   ( 0 , )) < . ( According to Property 2, set   to be stable.Similarly, the equilibrium core   corresponding to   is unstable and can be proved in the same way.

Stability Analysis on Colonies'
Members.Here we will discuss the stability of the colonies' members based on the the following different cases: (a) there is no shifting between members of different colonies; (b) one colony's members shift into another while another's members do not; (c) the two colonies' members are shifting into another, respectively.In this section, we define  as triangular fuzzy numbers, just as Figure 2. From Figure 2, And   ,  = 1, 2 in system (6) follow this definition.It is easy to achieve that the equilibrium cores (0, 0), ( 10 , 0), (0,  20 ) are unstable and equilibrium core ( 10 ,  20 ) is stable.So ( 10 ,  20 ) is stable in  2 + .This completes the proof.This indicates that if there is no shifting between the two colonies, the members of the two colonies will trend their maximum eventually and maintain stability.
]  > 0, system (3) has one and only one positive equilibrium core, which is stable in  2 + . Proof.
) . ( then the coefficient matrix of the linear system of system ( 48) is The feature values of it are both less than 0, so the equilibrium core is stable.

Now let us check Tr
and since ( 10 ,  20 ) is the equilibrium core of system (3),  10 ,  20 will satisfy the following equations, respectively: Consequently, ), and then Tr  < 0.
From (55), we can learn that From ( 56) we can conclude that   −    0 < 0 and det  > 0. As far as (Tr ) 2 − 4, det  > 0 clearly comes into existence.Therefore, the equilibrium core ( 10 ,  20 ) is stable.Consequently, the equilibrium core ( 10 ,  20 ) is stable in  2 + .This completes the proof.Accordingly, the members of the two colonies will trend their maximum eventually and maintain stability if satisfying certain conditions whether there are members shifting from one colony into another or not.

Stability Analysis on Innovation Diffusion.
Since the equilibrium core in system ( 3) is stable and we are interested in the asymptotic behavior of system ( 6), we only discuss the stability of system ( 6) on the base of system (3) in it's equilibrium core.That is to say, we only discuss the stability of innovation diffusion when the members of the two colonies are in their stability.Suppose that the global equilibrium core of system ( 3

𝑁 2𝛼
) have one and only one positive equilibrium core, which is stable in  2 + .
If [ 1 ]  = 0, [ 2 ]  > 0, system (58) can be reduced as ) .(59) Take the second equation of system (6) into consideration. 2 is the unique globe equilibrium core of this equation in  2 + , where Assume that  1 is the positive solution of the following equation: then follow Define  0 =   −   , ( = 1, 2), then system (59) can be rewritten as The coefficient matrix of the linear system of system ( 63) is Since the feature values of  are both less than 0, the equilibrium core ( 1 ,  2 ) is stable.So this equilibrium core is stable in  2 + .If [ 1 ]  > 0, [ 2 ]  > 0, the demonstration process is similar to the process of Theorem 4. We can conclude that there exits a unique positive equilibrium core ( 1 ,  2 ) and the equilibrium core is stable.Consequently, the innovation will reach its maximal potential and keep stable.This completes the proof.

The Empirical Analysis
Owing to the household management system of planned economy in history, China has divided the population into urban population and rural population.Owing to the strict control, rare shifting, huge gaps between industry and agriculture, and the different environment and living level, there exists great difference between urban population and rural population in communication channels, information spreading, and consumption habit.However, with the reform and exposure to the outside world, the economic development converts China into urban society from village society, and the inhabitants begin to move.According to the related news, there is 1% rural population moving into town per year since the 1980s.The ratio between urban population and rural population may be 7 : 3 after 60 years.Thus, it is not enough to discuss their consumption systems, respectively, without taking the population moving between the cities and countries into consideration.Based on this, we will try to build the dynamical model for telephone diffusion under the condition of the existing large population moving between town and country.
To make it easier, we suppose that this model satisfies following hypotheses: (1) since the standard of living in town is superior to that in the country of China, we suppose  2 = 0, which indicates that there are no people in China shifting from town into country; (2) since the communication by telephone is popular and necessary, we suppose  1 = 1,   = 0 which indicates that the telephone users in country will continue using it after shifting into town and users in town will not give up using it; (3) the population in China will reach 1.6 billion and maintain stable, which is obtained according to some related research on population in China [28].

Model Conversion.
In this section, we first use the definition of subduction defined by Zadeh, which means And then we use the -difference to simulate the diffusion process of local telephone in China.
According to the hypothesis defined above, system (58) can be rewritten as According to the definition of subduction defined by article [29], system (65) can be rewritten as That is, Then we used Genetic Algorithms (GA) to model the coefficient of model (67), which were shown in Table 4, where error is the permitted maximum error ratio between the data predicted by model and history data.
If we take the coefficient value as shown in Table 2 to the minimum of the target values in Table 3, the plots of potential curve of telephone in urban area or in rural area governed by model (67) are shown in Figure 3 or Figure 4, respectively.
In Figure 3, the red line is made up of history data and the green solid line and the green dot line represent the upper limit and the lower limit of the telephone consumers in urban area, respectively.
In Figure 4, the red line is made up of history data and the blue solid line and the blue dot line represent the upper limit and the lower limit of the telephone consumers in rural area, respectively.
From Figures 3 and 4, the fuzzy simulation effect is not perfect, and the fuzzy intervals of the consumer numbers are too big to predict the diffusion process of innovation.We could only find the time of inflection point and the maximal potential market from Figures 3 and 4.So we try to use the definition of Hukuhara difference, which,  −  = (() − (), () − ()), comes into exist.Consequently, we have the following system: We give the predicted figures in Figures 5, 6, and 7, respectively, in which the error control ratios are 45%, 50%, and 100%, respectively.
Since the the effect of simulation under other control simulation error ratio is not good, we only give the potential market of local telephone under the control simulation error ratio which is 45%, 50%, and 100%, respectively.From Table 5 from 54203.57to 73615.08 in respect to different control error ratios.

Contrastive Analysis.
In this section, firstly we consider the cases of system (58) without any fuzzy coefficients, then take the Bass model to simulate the process of local telephone diffusion in China.
We also suppose that this model satisfies the following hypotheses: (1) since the standard of living in town is superior to that of country areas in China, we suppose  2 = 0,  1 = 1%, which means that the population shifting from the country into the town is 1% of rural area inhabitants at that time while there are no people in China shifting from the town into the country;  (2) since communication by telephone is popular and necessary, we suppose  1 = 1, which indicates that the telephone users in the country will continue using it after shifting into town and users in town will not give up using it; (3) the population in China will reach 1.6 billion and maintain stability, according to data obtained from some related research on population in China.
Based on the hypotheses above, we can conclude that the population of the country will be stable in 0.48 billion and that of the town will be 1.12 billion, which is concluded by the nonlinear regression on system (3) without any fuzzy coefficients.According to the 14 groups of historical data, and the software SPSS11.0,we also make nonlinear regression on system (58) and obtain the following model: In which the coefficient of external influence in town  1 = 0.0007934, the coefficient of the internal influence  1 = 0.22536, and those in the country, respectively, are 0.00001 and 0.399696.With the calculation on the system (59), the market potential of telephone users in the town is 375684000 and that in country is 133090000.The  squared of the first equation of (69) is 0.95748 and that of the second is 0.91051.We use the software Mathematic 4 to forecast the numbers of telephone users in the town and the country, which are shown in Figures 8 and 9, respectively.And the origin point in these figures denotes the year 1989.
Consequently, we arrive at the following conclusions.
(1) When the population of the town and that of the country trend to stability, the users among them also trend to stability, the number of which is 13309000 and 375684000, respectively.
(2) Since the medium system in town is highly superior to that in the rural areas, the necessity and popularity of communication in town are also better than those in the rural area and the coefficient of external affection in town (0.0007934) is obviously larger than that in rural area (0.00001).(3) The coefficient of internal affection in town (0.22536) is much less than that in the rural area (0.3399696), which means that the inhabitants of the rural area are more reliable on oral communication than those in town.The main reason is that there are few communication channels in the rural area and the communication information mainly depends on oral communication.However, the network of communication in town is much developed and the necessity of oral communication is not so necessary.(4) The population of telephone users in town will reach its maximal market potential in about 30 years, while in rural area it will be in about 15 years later.In fact, this dose not mean that the diffusion rate in rural area is faster than that in town.In fact, the diffusion rate in town is much faster than that in rural area.The main reason for this is that the market potential in town (375684000) is much higher than that in rural area (133090000).
We also use Bass model to describe and forecast the numbers of the telephone users in the town or in rural areas, respectively, without taking the people shifting between them into account.With the nonlinear regression, another forecasting model is governed by in which the coefficient of external affection in town  1 = 0.000313699, the coefficient of the internal affection  1 = 0.233073, and those in rural area, respectively, are 0 and 0.42672.The market potential of telephone users in town is 644431789 and that in rural area is 134464950.The  squared of the first equation of (70) is 0.94766 and that of the second is 0.9047.Furthermore, we also use the Bass model to describe and forecast the number of the telephone users without considering where the users are from.The model is as follows: .  = (0.000412046 + 3.70996 × 10 −6 ) × (71216.171268− ) . (71) The market potential of telephone users in China is 712161700 and the prediction curve of this model is shown in Figure 10.The origin points in these figures denote the year 1989.
Comparing model (71) with model (70) we can see that if we distinguish where people settle down, the total market potential of telephone users in China will be 778896700 or, if not considering the place of setting, it will be 712161700.That is to say, there will be another 66 million in the market potential of China if distinguishing where people settle down.In order to distinguish which prediction model is more accurate, we have drawn the plots of the prediction curves of different models and compared them with those of the historical data, which are shown in Figures 11,12, and 13, respectively. ,  ,   denote the diffusion rate of telephones in town, in the rural area, or the total in China, respectively, and  ,  ,   denote the prediction numbers of those, respectively, by Bass model.  ,    denote the prediction numbers of the diffusion rate of telephones in town or in the rural area by model (69).From Figure 11, there is no difference in accuracy between model (70) and model (71).However, the prediction numbers of model (69) are more accurate than those of model (70), as deduced by Figures 12 and 13.
We also drew the prediction curves of telephone users in town and in rural areas governed by model (69) or model (70), which are shown in Figures 15 and 14, respectively.The curves are almost identical in Figure 15, which means that the potential of telephone usage in rural areas is virtually the same whether governed by model (69) or model (70) at any given point time , while there is a wide difference between those in town, as shown in Figure 14.With the course of China's stepping into an urban society from a rural society,   more people are moving into town from rural areas and have a stronger ability and demand for communication.So the potential predicted by model (69) should be bigger than that predicted by model (70).However, on the contrary, the potential predicted by model ( 69) is much smaller than that predicted by model (70) as shown in Figure 14.
The cause may be that the Bass model uses historical data to predict the market potential, and the historical data of

Conclusion
In this paper, we constructed an innovation diffusion model with migration between two different colonies.With the analysis of the innovation diffusion model based on the members of the two colonies shifting between each other, and on the empirical analysis on the telephone users of China, we conclude that the subscribers will trend to their maximal market potential whether or not there are members shifting between town and rural areas.During the process of urbanization and the relaxation of household management in China, such phenomena should not be ignored since its influence on the economic and social system may be vital and significant.Since obtaining data is very difficult, the population data is directly obtained from related information and empirical estimation.Moreover, the proposed model is acceptable without competing technologies taken into consideration.The simulation of data is not accurate enough since we have not taken the economic environment, competing technologies, and the effect of decreasing costs in communication into consideration, all of which will be included in future research.

Figure 1 :
Figure 1: Customers flow diagram for innovation between two colonies.

Figure 2 :
Figure 2: The definition of triangular fuzzy number.

5 Figure 3 :Figure 4 :
Figure 3: Plots of potential curve of telephone in urban area governed by model (67).

Figure 7 :
Figure 7: Plots of potential curve of telephone in urban area and rural area governed by model (68) (the control error rate is 100%).

Figure 8 :
Figure 8: Plot of the diffusion of telephone users in urban area of China governed by model (69).

Figure 9 :
Figure 9: Plot of the diffusion of telephone users in rural area of China governed by model (69).

Figure 10 :
Figure 10: Plot of the diffusion curve of telephone in China governed by model (71).

Figure 11 :
Figure 11: Plots of the diffusion rate of telephone by Bass model.

Figure 12 :
Figure 12: Plots of the diffusion rate of telephone in urban area by model (69) or model (70).

Figure 13 :
Figure 13: Plots of the diffusion rate of telephone in rural area by model (69) or model (70).

Figure 14 :
Figure 14: Plots of potential curve of telephone in urban area governed by model (69) or model (70).

Figure 15 :
Figure 15: Plots of potential curve of telephone in rural area governed by model (69) or model (70).

Table 1 :
The explanation for conversion relation.

Table 2 :
The data of telephone subscribers in the year 1989-2011.

Table 3 :
The error control target value of history data.

Table 4 :
The error control target value of history data.
, the number of local telephone consumers in China ranges

Table 5 :
The predicted consumer numbers of China local telephone under different control simulation errors.