Optimal Production Planning for Manufacturing Systems with Instantaneous Stock-Dependent Demand and Imperfect Yields

We consider an EPL model like manufacturing system in presence of production imperfectness and stock-demand dependence simultaneously. During the production process, the system can evolve from in-control state into out-of-control state at any random time, after which the defective items will be generated likely causing quantity loss. Meanwhile, the market demand rate is instantaneously dependent on the timely holding inventory. The manufacturer has to determine his production run length and cycle time by taking into account possible imperfect production, stock-dependent demand, and inventory holding capacity bound. We empolder a model to capture this problem and develop computational algorithm to solve it. We further conduct numerical studies to validate our model and solving method. Sensitivity analyses are reported to show the effect of parameters on the system performance.


Introduction
Managing supply/production uncertainty and changeable demand is always one of the challenging operations since matching them is fairly hard due to instability occurring in two sides, respectively.The firm can suffer managerial issues resulting from both out-of-control production (e.g., Rosenblatt and Lee [1] and Sana [2]) and demand's instantaneous dependence on inventory (e.g., Gupta and Vrat [3] and Urban [4]).
This study involves two aspects of imperfect production and demand-stock dependence in the OM field.Imperfect production property means that defective items can be produced as the production system evolves over time.These imperfect goods will be reworked or deposed directly.As a consequence, the inventory/production planning is spoiled.On the other hand, the demand rate is instantaneously dependent on the timely inventory holding level of quality goods.That means the random production imperfection will eventually affect demand rate through defective items.Setting longer production cycle time to obtain higher inventory can obtain more demand rate resulting in more revenue.However, the longer the production cycle time is, the larger likeliness the imperfect production happens.The annual total cost comprises annual ordering cost and annual inventory holding cost.The production cycle time has opposite effects on these two kinds of cost components.Settling the optimum trade-off between the benefits and corresponding costs for the sake of optimizing the annual total profit is accordingly a crucial problem confronting practitioners.
Our study intends to address this problem by considering the production cycle time and pure production run time length as decision variables in the modeling of the manufacturing system existing in imperfect production and demand-inventory dependence.Specifically, we research the following problem: a manufacturer schedule his/her production annually by deciding each cycle time and run length, whereas the manufacturing system may switch from the incontrol state to the out-of-control state after some time on the account of continuous run of the machine, which takes place at any random time during the production process, also as presented in Shah and Shukla [5].Moreover, this is not immediately under observation and finally causes the quantity loss.The market demand is sensitive to the on-hand effective inventory excluding the defective goods.

Mathematical Problems in Engineering
Although we are concentrating on a new model of production planning contextualized in uncertain yields simultaneously confronting inventory-dependent demand, our modeling processes are to some extent inspired by preceding papers.We give an elucidation on the board literature related to our study below.
The existing relevant literature mainly includes two streams: (1) uncertain production management; (2) inventory/production planning with inventory-dependent demand and their extensions.The production uncertainty incurred by various factors like machine failure, staff changes can always cause the entire production system instable resulting in defective products has inspired extensive studies, such as Rosenblatt and Lee [1], Lee and Rosenblatt [6], Lee and Rosenblatt [7], and a variety of their succeeding papers like Liou et al. [8], Kim and Hong [9], Chung and Hou [10], Guo et al. [11], Hu et al. [12], and so on.Drawing on the celebrated EOQ/EMQ models, early pioneering researches study the production planning optimization considering defection inspecting, restoration cost dependence, detection delay, and maintenance, say, Rosenblatt and Lee [1], Lee and Rosenblatt [6], and Lee and Rosenblatt [7].Their valuable investigations firstly pull the production planning research closer to the realistic practice where there usually exist in outof-control states for the production process incurring a lot of associated management problems.Furthermore, inspection error types are proposed to be other factors having effects on production planning in an imperfect EMQ manufacturing system; see Liou et al. [8].From a different perspective, Kim and Hong [9] link the distinction among three types of deteriorating processes to the optimal production run length seeking.Later, Chung and Hou [10] generalize former literature to allow the shortage in the imperfect system.After those static models in this line of research, Sana [2] addresses a production-inventory model with imperfect production and time-varying demand.Their study substantially extends proceeding researches into dynamic situations.In contrast, some recent researches turn to exploring imperfect production planning in the fuzzy environment, like Guo et al. [11] and Hu et al. [12].Although the aforementioned literature has done intensive studies on uncertain production incorporating features incurred by itself, that existing research fails to consider the change of demand pattern during the production cycle time.However, they usually assume the demand characteristics are invariable over time, which is just a very special case of our model in this study.
Observation to practice has shown that the real market demand is very complex and it is not a fixed constant, as implied by Levin et al. [13], Silver and Peterson [14], and Gupta and Vrat [3].Incited by their ground-breaking work, a variety of studies have incorporated stock-dependent rate component into diverse classical inventory models, for instance, by considering items deterioration in Mandal and Phaujdar [15], Datta and Pal [16], and Dixit and Shah [17], involving shelf-space allocation in Baker and Urban [18], allowing stock shortage in Datta and Pal [16], and combining delayed payment schemes in Soni and Shah [19], Min et al. [20], Shah et al. [21], Teng et al. [22], and Zhou et al. [23].A considerable amount of literature has also addressed relevant problems, like Goh [24], Wee [25], Padmanabhan and Vrat [26], Ray and Chaudhuri [27], Sarker et al. [28], Giri and Chaudhuri [29], Yan and Cheng [30], Mandal and Maiti [31], Chung and Tsai [32], and so on.To learn more about the research of this kind, one is referred to a recent literature review by Urban [4].Almost all literature mentioned above has only focused on the dependence between demand and inventory with the implicitly obvious a priori prerequisite of completely perfect production, which usually does not fit the practice well.On the contrary, in this paper, we alter that precondition and emphasize the imperfect production frequently occurring in real world.
To the best of our knowledge, few of existing studies have considered the production imperfectness and inventorydemand dependence simultaneously.In contrast, this is just what we concentrate on in the present study.We model the production and demand consumption processes with the aforementioned characteristics as a constrained nonlinear programming system combining different scenarios.Furthermore, our model is distinguished from existing relevant models in several ways: (1) we combine imperfect production and inventory-dependent demand together to explore the optimal production problem; (2) we use optimizing the annual total profit as the objective function instead of annual total cost, which substantially reflects the effect of demandinventory dependence on benefits; and (3) we incorporate the inventory capacity upper boundary into traditional models of imperfect production and inventory holding.We select and develop a suitable computational algorithm to solve the nonlinear programming model with detailed operational steps and procedures.Subsequently, we conduct numerical examples and parameters sensitivity analysis to examine the proposed computational algorithm so as to get some managerial insights (Figures 7,9,10,11,and 13).
The rest of the paper is organized as follows.In Section 2, we sketch the problem features, notations, and its assumptions.In Section 3, we discuss the various scenarios arising from our studied system and derive the associated model.In Section 4, we describe a computational algorithm for solving the constructed model.In Section 5, numerical studies and associated analysis are executed to validate the algorithm proposed.Finally, concluding remarks are given to summarize this study in Section 6.

Problem Characteristics, Notations, and Assumptions
We study a manufacturing system described like a variant of economic production lot (EPL) model, which is distinct from classical EPL model by characterizing its specific production property and demand changing rules, respectively.On the production side, the system may randomly evolve into an out-of-control state in which a proportion of goods will be generated following some probability distribution.On the demand side, the impact of inventory level on demand is considered.The demand rate consists of the constant demand rate and the inventory level dependent demand rate addressed as  () =  +  () , 0 ≤  ≤ , where  > 0 and  > 0.
So the dynamics where the stock () is depleted at time  can be presented as follows: The production process can be described like this: for scheduled production cycle time , the manufacturer runs the production line from the beginning to time , that is, [0, ], where 0 <  < .The production line is closed down in interval (, ].The manufacturing system may evolve into out-of-control situation incurring defective goods produced randomly.Denoting the production rate V and the production duration  yields aggregate volume V during the cycle time.It may produce defects after producing continuous  time and the system failure results in inferior goods in  proportion of all production during [, ], where  obeys an exponential distribution with parameter .During the whole cycle time, the inventory is being consumed by the demand at a rate varying according to the associated stock level at the same time. The following notations will be used throughout the paper.  : sales income of qualified products ($/year),   : crapped income of defective ($/year), : total annual profits,  −  ($/year).

Parameters
To refine and reflect the essence of our problem, we need assumptions in our model as follows.
(A1) The lead time is zero and the time span is unlimited.(A2) The production capacity is limited but out-of-stock is not allowed, and productivity is greater than steady section, that is V > .(A3) After producing continuous  time, the system enters an uncontrolled state and defective rate of the product obtained is .After discovering the defective, it needs be scrap processing timely instead of entering the inventory system, where downtime  obeys an exponential distribution with parameter .(A4) Production process is independent in different cycles, which is in each cycle, and the duration  of controlled state is independent.(A5) During each production cycle, the beginning and ending inventory are both zero, that is, (0) = 0 and () = 0.

The Model
For a given production cycle time , there may be two cases in the whole production process because of the occurrence of imperfect production at time .One case is that breakdown time  does not occur before the production stopping, that is,  ∈ [, +∞], and products obtained in during the period are all qualified products.The other is that the breakdown time  occurs during the production process, that is,  ∈ [0, ], and the defective rate is  (0 <  < 1) from time  to time .Assume (1) the upper bound of the stocking system capacity is  max ; (2) inventory levels at the beginning and the ending times of the production period are both zero, namely, (0) = 0 and () = 0. We plot the inventory level changing trajectory of qualified products in Figure 1 with two cases, where  1 is the stopping time of the continuous production with system under control completely, and so is  2 when system out of control occurs.

System under Control
Completely before Production Stopping.Figure 2 shows the inventory changes in case 1.
The entire inventory cycle time can be divided into two sections, that is, [0,  1 ] and [ 1 , ].In the period [0,  1 ], the system does not produce defective products because the production is carried out under control.During this period, the boundary value is  1−1 (0) = 0 and because V > , the inventory keeps increasing and slope of inventory In the period [ 1 , ], production ceases.The inventory consumption comes from the fixed demand and the needs with the inventory changing.During this period, the boundary value is  1−2 () = 0 and the slope of inventory  1−2 () is −(), which is given by Solving differential equations ( 4) and differential equations ( 4), the inventory level in case 1 is  1 (): The inventory function  1 () is a continuous function.So where  =  1 , there is Simplify (6) and then obtain the relationship between  and  1 , which is given by Total annual costs  1 include annual inventory costs  1 and annual ordering cost  1 ; that is,  1 =  1 +  1 , in which, annual ordering cost  1 consists of two parts: one is the preparation costs for each production and the other is production costs which change with the production volume: Total annual revenue  1 : there are no defects in case 1; total annual revenue  1 equals   , the sales income of qualified products.Consider Total annual profit  1 is total annual revenue minus total annual costs, that is,  1 =  1 −  1 , which is given by Substituting formula (7) into formula (10), then we can get total annual profits function which depends on : For any fixed inventory cycle , the ending inventory of system equals zero, that is, () = 0.This paper will denote (1/) ln((  −+V)/V) as a critical point of , denoted by  0 .Where  > (1/) ln((  −+V)/V), the production processes will follow case 1.And where  > (1/) ln((  −  + V)/V), there will be defective goods in [,  2 ].For a detailed analysis, see case 2.

System Running Out of Control before the End of Production.
To maintain the ending inventory () = 0,  2 needs to change with . Figure 3 shows the qualified inventory changes.
In this case, the entire inventory cycle is divided into three sections, that is, [0, ], [,  2 ], and [ 2 , ].In the period [0, ], the system will never produce defective products because the production is carried out completely under control.And during this period, the boundary value is  2−1 (0) = 0 and the slope of inventory In the period [,  2 ], production turns into an uncontrolled state.So there will be defective products and the defect rate is .Where  = , we have  2−1 () =  2−2 () and the slope of the qualified inventory In the period [ 2 , ], the inventory is consumed continuously until it reduces to zero.The boundary values  2−3 () = 0 and the slope of inventory  2−3 () are −(); that is, Solve differential equations ( 12), (13), and ( 14) and then get the expression of inventory levels  2 () in case 2: Obviously,  2 () is also a piecewise function.It is easy to prove, where  ∈ [0, (1/) ln((V − )/V)], in the period  ≤  ≤  2 , the inventory  2−2 () is an increasing function.And where  ∈ ((1/) ln((V − )/V),  2 ], in the period  ≤  ≤  2 , the inventory  2−2 () is a decreasing function.For details see appendix, Proof 1.
Due to the continuity of the inventory function, where Simplify (16) and then obtain the relationship among , , and  2 : In case 2, total annual costs  2 also consist of annual inventory costs  2 and annual ordering costs  2 , that is In case 2, defects occur during production, so total annual revenue  2 consists of two sections: sales income of qualified products  2 and scrapped income of defective annual revenue  2 : Total annual profits are equal to the annual revenue minus annual costs (Figure 5), that is,  2 =  2 −  2 , which is given by Mathematical Problems in Engineering Substituting (17) into (20), then simplify to 3.3.Optimizing Annual Total Profit.Objective function: expected annual profits .
In this paper, the inventory model is divided into case 1 and case 2, according to whether there will be defects in the production.This paper assumes that , the moment of appearing defect in (0, +∞), obeys an exponential distribution with parameter .Therefore the expected annual profits under the inventory capacity constraints.Before discussing the inventory constraints, we should compare  1 in case 1 with  2 in case 2.
It is written as follows: It is easy to prove  2 ≥  1 , when they have the same .For details see appendix, Proof 2.
In case 1, where  =  1 , the inventory reaches the maximum, which should be less than  max , that is, [ (− 1 ) − 1]/ ≤  max .In case 2, if  ∈ [0, (1/) ln((V − )/V)], then where  =  2 , the inventory reaches the maximum.And this time, the inventory should be less than  max , that is, , then where  = , the inventory reaches its maximum.And the inventory should be less than  max , that is, Then, the inventory constraints are We can prove that For details see appendix, Proof 3.So the inventory constraints convert into Substituting ( 7) into (26), we have Constraint formula ( 27) can be changed as follows: Since V > , the right side of ( 28) is a positive constant.And the left side can be divided into two cases.
(1) If V −  −  max ≤ 0, we have   > 0. So the left side of ( 28) is negative and inequality ( 28) is always true.Then the inventory capacity constraints have no practical significance.
(2) If V −  −  max > 0, obviously the left side of ( 28) () = (V −  −  max )  is a monotonic increasing function of .To satisfy the constraints, there should be The inventory capacity constraints are successfully converted into the scope limitation of the decision variable values, which becomes simpler and more intuitive.In summary, the expected annual profits maximization model changes into max Model ( 30) can be regarded as a unary nonlinear programming model.Thus the problem can be transformed into unary nonlinear programming problem.
Formula (22) shows that the profit function  contains the integral of transcendental functions, and there is no analytic expression.Thus it is difficult to use derivatives and gradient information to prove the concavity of the profit function, and the difficulty and complexity of proving it are far beyond the traditional inventory model.However, [] is a continuous function.According to the maximum and minimum theorem of continuous function nature theorem on a closed interval, we can determine that the profit function [] has a maximum within the specified interval.Based on this, this paper designs the following optimization algorithm.

Computational Algorithm
There are a variety of commonplace methods to solve nonlinear programming models, for instance, Golden Section, Fibonacci, Interpolation, Bisection, Steepest Descent, Newton, quasi-Newton, and Conjugate Gradient and so on.Those traditional algorithms are simple with fast convergence.However, they have relatively high requirements on initial value and functional structure.For more general computational application, scholars have proposed a variety of intelligent algorithms, such as Genetic Algorithms, Simulated Annealing Algorithm, Ant Colony Algorithm, and Chaotic Mixing Algorithm with their own advantages and broad applications, respectively.To overcome the defects of a single algorithm, integrating several algorithms in solving one problem has been prevailing in optimization field through mutually advantageous complementation among them.
Artificial fish swarm algorithm (AFSA) is one kind of integrated algorithms as above mentioned.Li [33] is among the early batch of studies to propose AFSA, which is an intelligent optimization algorithm based on fish behavior in nature.Peng [34] also demonstrates this approach in detail and compares its advantage to other algorithms.To learn more about the growing path of AFSA, Neshat et al. [35] is a good reference with excellent comprehensive literature review.Imitating and realizing three kinds of individual fish's behavior, that is, feeding and rear-ends, clusters underlie the working basis of the algorithm optimality seeking with the associated system parameters design including fish scale size , the step, the congestion factor delta, and the perceived distance visual, Li et al. [36] and Li et al. [37].Although the initial AFSA has lots of obvious virtues, there still exists improving space to abate its some weakness for application, which has been pointed out in detail and also improved to apply in diverse settings, say, Wang et al. [38], Li et al. [37], Huang and Zhou [39], Chen et al. [40], Cheng et al. [41], Peng [42], and Tsai and Lin [43].
Drawing on the fact that the concavity of the profit function [] cannot be clearly identified in this paper, we select the artificial fish school algorithm to solve our problem while we make a substantially specific problem-featured improvement for this algorithm inspired by existing literature aforementioned.Here we develop an improved artificial fish school algorithm by combining the adapting parameters strategy and tabu search heuristics simultaneously to enhance the local search ability and convergence rate of initial ASFA.However, our focus in this study is mainly on formulating, modeling, and analysing the focal managerial problem rather than the algorithm improvement.This proposed algorithm is just served to solve the nonlinear problem and will not be compared with other algorithms.We need to emphasize the fact that comprehensively integrating more extensive computational methods and heuristics together finally distinguishes our algorithm design from that in existing researches.

The Parameter Adapting Strategy
Strategy of Changing Step Length.A fixed step length is usually adopted in most of literature for searching optimizer in an artificial algorithm.The step length will affect the converging rate, speed, and accuracy of the searching process.The longer step length will incur higher converging rate and searching speed, but lower accuracy.On the other hand, a small step size will help local search with a slower global converging rate with higher accuracy.The varying step length will be much better in balancing the converging rate and searching accuracy.In this paper, the step is set as (31), a progressively decreasing function as addressed in Huang and Zhou [39] and Chen et al. [40]: Strategy of Changing the Congestion Factor Delta.In order to speed up the searching rate, avoid falling into the plight of local optimum and premature and suppress the congestion factor.In the iterative process, we adopt the congestion attenuation function set in (32) as presented in detail in Zheng [44] to progressively reduce congestion factor to enhance local search convergence rate: Strategy of Changing the Perceived Distance Visual.According to Chen et al. [40], we set a perceived distance motion function (33) to make the perceived distance decreases with  the iterative number increasing so as to improve the local intensity and solution accuracy:

Tabu Search Heuristics.
To solve the continuous problem in the present study with infinite number of neighbors, using the tabu method of discrete problems cannot make the algorithm effectively searching the global optimal solution.Thus we need to employ the neighborhood tabu for solving.Li et al. [37] elaborate the neighborhood tabu algorithm and propose a two-dimensional neighborhood concept for tabu algorithm.Inspired by their research and considering the specific situation in our current problem, we give out a onedimensional neighborhood concept applied for our problem solving.Divide current feasible solution interval [  ,   ] into  fields falling within a circle with center (  +   )/2 and the radius  as shown in Figure 4.
In the gen-th interaction, the current position of th fish is   , and another position is   .  can be the optimal point in the view of the th fish and also can be the center point or another point.The optimization researching time is denoted as ; the optimization coefficient  is a randomly generated number between (0, 1) by systems.
The next position of th fish is determined by the following formula: The process of tabu search is as follows.
Step 1.   is the current state of artificial fish,  is the counter, and  = 1.
Step 2. We get a new point  |next from (34) and then calculate  |next , and  =  + 1.
Step 3. Compare the size of   with  |next .If   <  |next , turn to Step 4; if not, go to Step 5.
Step 4. Estimate  |next belonging to which  neighborhood and determine whether the th neighborhood exists in the tabu list.If it does not exist, the optimization is successful, and then insert the th neighborhood into the tabu list.If the tabu list is full, update the list according to the FIFO principle and then end the optimization process.However, if it does exist, go to Step 5.
Step 5. Estimate if  is greater than the maximum attempts time ; if so, the optimization process ends in failure.
Otherwise, go to Step 2.

Behavior of a Single Fish.
Applying the aforementioned established neighborhood tabu search process to our problem, we get the individual behavior of a single fish including foraging, clustering, rear-end, and searching for historical best fish.Here we introduce these four behaviors.
(1) Foraging.The position of a single fish in current fish swarm is   .The fish randomly chooses a position   in their sensing range,   =   +  × visual gen , and then calculate the tabu optimization.If we get the optimization successfully, we take  |next as the next mobile location of the current fish.If not, randomly move to another step, that is,  |next =   +  × step gen .
(2) Clustering.A single fish position in current fish swarm is   .In this fish's vision, the number of other fish is   and the optimal position of the fish is   .Compare the size of   /  with delta ×   .When   /  > delta ×   , calculate the tabu optimization.If we get the optimization successfully, denote  |next as the next mobile location of the current fish.
If not, turn to foraging behavior.
(3) Rear-End.A single fish position in current fish swarm is   .In this fish's vision, the number of other fish is   and the geometric center position of the fish is   .Compare the size of   /  with delta ×   .When   /  > delta ×   , calculate the tabu optimization.If we get the optimization successfully, denote  |next as the next mobile location of the current fish.If not, turn to foraging behavior.
(4) Searching for Historical Best Fish.A single fish position in current fish swarm is   .Learning from bulletin board, we have the historical best fish's position   and then calculate the tabu optimization.If we get the optimization successfully, denote  |next as the next mobile location of the current fish.
If not, turn to foraging behavior.

The Overall Improved AFSA Integrating Adapting Parameters and Tabu Search
Step 1. Set the population size , the moving step , the congestion factor delta, the perceived distance V, the maximum number of attempts , and the length of tabu list .Make discrete solution space into  neighborhood and limit the maximum number of iterations , making gen = 1.
Step 2. Initialize the fish swarm ( 1 ,  2 ,  3 , . . .,   ) in solution space.Calculate the objective function value of current fish swarm ( 1 ,  2 ,  3 , . . .,   ), the optimal fish   in the current fish swarm, and the historical best fish  best .Record the information of historical best fish on the bulletin board.
Step 3.Each fish in the fish swarm gets three best  |next through the three best behaviors: clustering, rear-end, and searching for historical best fish.Select a larger target value of  |next and update their position.After comparing with the historical best fish, update the information on bulletin board.

Numerical Studies and Analysis
In this example, V −  −  max > 0, so the limit of inventory maximum capacity is working.According to (29), we calculate the range of maximum value of the decision variable  and compare the result with 5. Then select a smaller one to be the maximum value of the decision variable .Using the improved artificial fish swarm algorithm, we calculate the optimal solution of this example [, ] = [0.38, 1945.71].The results and changing trend are shown in Figure 5.In this example, V −− max ≤ 0, so the limit of inventory capacity constraints do not need to be considered.The range of decision variable  is (0, 5].Within this range, we calculate the optimal solution of this example by iterative optimization, which is [, ] = [0.19,19475.09].

Sensitivity Analysis. The sensitivity analysis on Example 1 follows.
(1) The Impact of Inventory-Dependence Coefficient  on the Optimal Solution.Assign  with different value and observe the changing tendency of the profit function and the optimization.
As can be seen from Figure 6, with the increase of , the optimal production cycle  * and the maximum annual profits  also gradually increase.An increasing  implies that the impact of current inventory on the demand is getting greater.The higher the current inventory level is, the greater the demand becomes.Within the inventory capacity, decision maker will try to increase the inventory in order to reach higher sales income.So the maximum annual profit  gradually increases.But to run out a higher level of inventory, it will require a longer cycle.Therefore the optimal production cycle  * will show an increasing trend.
With the change of , range of decision variable  constantly changes.When  is relatively small and  < (V − )/ max , the maximum value of  is the min one of (1/) ln((V −)(+ max )/(V −− max )) and 5. (1/) ln((V − )( +  max )/(V −  −  max )) is increasing as  increases.When (1/) ln((V − )( +  max )/(V −  −  max )) ≤ 5 is satisfied, the maximum value of  also increases, which is reflected on the length changes of the corresponding profit curve.For example, the profit curve is longer with  = 0.7 than  = 0.5.When  increases to  ≥ (V − )/ max and the maximum value of  is the preset 5, it is not difficult to find out that the upper limit of  becomes 5 with  = 0.9.
(2) The Impact of Inventory Holding Space Bound  max on the Optimal Solution.In Figure 7, we plot the profit-time curves varying as the inventory holding capacity changes.The results show that those curves are almost overlapped,  which implies the high capacity has few influence on the optimal solution.It can also be seen from Figure 8 that when  max changes, the optimal solution will remain unchanged, based on the premise that the change of  max still ensures the optimal solution  * within the optimization search scope of ; otherwise, the optimal solution  * and the maximum annual profit  * are both likely to change.
(3) The Impact of Defective Rate  on the Optimal Solution.
When there is defective production, it is essentially a loss.With the same production cycle , the annual profit will decrease.When defective rate increases, decision makers will reduce the probability of uncontrolled state.The probability of downtime is exponentially distributed; the shorter the production time is, the smaller the probability of downtime will be.So when  increases, the optimal production cycle  decreases.
(4) The Impact of the Parameter of Exponential Distribution  on the Optimal Solution. is the parameter of the exponential distribution, which the system downtime  obeys.1/ is the expected value of downtime.With an increasing , 1/ gets smaller which means that the system fault occurs quicker.If we keep the fixed production cycle  in this condition, the uncontrolled production state will last longer.There will be more defects in the output and the total profit will decrease.Thus the decision maker will decrease the production cycle  to reduce the frequency of system failure and the defective output rate.However, even if  increases continuously,  will not continue reducing.Because a decreasing  means the startup time of this machine in one year increases, the fixed costs increase.Affected by the fixed costs,  will stop decreasing and increase instead.
(5) The Impact of Productivity V on the Optimal Solution.
An increasing V means the improvement of production efficiency and acceleration of inventory accumulation.So it can accumulate enough inventory in a shorter production time, and the average system inventory level will improve.
As the demand is associated with the current inventory, the demand will also increase.When the value of  is small, the system stays at a relatively low inventory level.With an increasing demand, the increase in sales revenue is greater than the increase in inventory costs, so the annual profit will increase with V (Figure 12).When the value of  is larger, the inventory level of system is relatively higher.At this point, with the climbing demand, the increase in sales revenue is not enough to offset the increase in the inventory.So the annual profit will decrease when V increases.When V increases, the maximum annual profit  * continues increasing, and the optimal production cycle  * decreases firstly and then increases.For a fixed , when the value of V is small, the system stays at a relatively low inventory level.With an increasing demand, the increase in sales revenue is limited and not enough to offset the increase in the inventory.So when V increases, the optimal production cycle  * needs to be reduced in order to reduce the inventory costs.When the value of V is larger, the inventory level of system is higher and the demand increases considerably, resulting in the fact that the increase in sales revenue is greater than the increase in inventory costs.So when V continues climbing, the optimal production cycle  * will appear as an upward trend.
(6) The Impact of Inventory Holding Cost ℎ on the Optimal Solution.When ℎ increases, the inventory costs will aggrandize with it and the total annual profit will decrease.So decision makers will reduce the production cycle  to reduce the annual average inventory level in order to cut down the inventory costs.When ℎ is larger, the annual profit curve will drop more rapidly with an increasing .

Conclusion
This study has investigated production optimization problem arising from a manufacturing system with likely imperfect production and instantaneously inventory-dependent market demand.We have constructed constrained nonlinear programming models to depict this problem by incorporating the randomness of production imperfectness, inventorydemand dependence, and stock holding space bound.For seeking the optimal annual total profit, an improved artificial fish swarm algorithm combined with tabu search was devised to solve the model with numerical studies and sensitivity analysis.This research adds to our understanding of managing the manufacturing system with new characteristics as discussed in the present paper.It contributes to providing the method and strategy of determining optimum production cycle time to obtain best system performance.This study also shows us the necessity and power of integrate operations and sales.However, the research is limited in several ways.For instance, we only confine our study to a single manufacturer, which can be extended to a supply chain system.Furthermore, we consider the demand is only inventory-dependent with deterministic dependence sensitivity parameter.Other random factors can be added to the demand function in the future research.

Figure 1 :
Figure 1: Comparison of effective inventory level changing trends in two scenarios.

Figure 2 :
Figure 2: The effective inventory level changing trajectory when production is under control consistently.

Figure 3 :
Figure 3: The trajectories of change of the qualified inventory level when the system is out of control.

Figure 5 :
Figure 5: The total annual costs, the total annual revenue, and the total annual profits changes in the two examples.

Figure 7 :
Figure 7: The sensitivity analysis of  max .

Figure 11 :
Figure 11: The sensitivity analysis of V.

Figure 12 :
Figure 12: The changing trend of  * and  * when V changes.