• 検索結果がありません。

S. Sarkar and T. Chakrabarti

N/A
N/A
Protected

Academic year: 2022

シェア "S. Sarkar and T. Chakrabarti"

Copied!
23
0
0

読み込み中.... (全文を見る)

全文

(1)

Volume 2012, Article ID 264182,22pages doi:10.1155/2012/264182

Research Article

An EPQ Model with Two-Component Demand

under Fuzzy Environment and Weibull Distribution Deterioration with Shortages

S. Sarkar and T. Chakrabarti

Department of Applied Mathematics, University of Calcutta, 92 APC Road, Kolkata 700009, India

Correspondence should be addressed to S. Sarkar,sanchita771@rediffmail.com Received 26 April 2011; Revised 28 June 2011; Accepted 15 July 2011

Academic Editor: Hsien-Chung Wu

Copyrightq2012 S. Sarkar and T. Chakrabarti. 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.

A single-item economic production model is developed in which inventory is depleted not only due to demand but also by deterioration. The rate of deterioration is taken to be time dependent, and the time to deterioration is assumed to follow a two-parameter Weibull distribution. The Weibull distribution, which is capable of representing constant, increasing, and decreasing rates of deterioration, is used to represent the distribution of the time to deterioration. In many real- life situations it is not possible to have a single rate of production throughout the production period. Items are produced at different rates during subperiods so as to meet various constraints that arise due to change in demand pattern, market fluctuations, and so forth. This paper models such a situation. Here it is assumed that demand rate is uncertain in fuzzy sense, that is, it is imprecise in nature and so demand rate is taken as triangular fuzzy number. Then by usingα-cut for defuzzification the total variable cost per unit time is derived. Therefore the problem is reduced to crisp average costs. The multiobjective model is solved by Global Criteria method with the help of GRGGeneralized Reduced GradientTechnique. In this model shortages are permitted and fully backordered. Numerical examples are given to illustrate the solution procedure of the two models.

1. Introduction

The classical EOQ Economic Order Quantity inventory models were developed under the assumption of constant demand. Later many researchers developed EOQ models taking linearly increasing or decreasing demand. Donaldson 1 discussed for the first time the classical no-shortage inventory policy for the case of a linear, positive trend in demand. Wagner and Whitin 2 developed a discrete version of the problem. Silver and Meal 3 formulated an approximate solution procedure as “Silver Meal heuristic” for

(2)

a deterministic time-dependent demand pattern. Mitra et al. 4 extended the model to accommodate a demand pattern having increasing and decreasing linear trends. Deb and Chaudhuri 5 extended for the first time the inventory replenishment policy with linear trend to accommodate shortages. After some correction in the above model 5, Dave6 applied Silver’s 7 heuristic to it incorporating shortages. Researchers have also worked on inventory models with time-dependent demand and deterioration. Models by Dave and Patel 8, Sachan9, Bahari-Kashani10, Goswami and Chaudhuri11, and Hariga 12 all belong to this category. In addition to these demand patterns, some researchers use ramp type demand. Ramp type demand is one in which demand increases up to a certain time after which it stabilizes and becomes constant. Ramp type demand precisely depicts the demand of the items, such as newly launched fashion goods and cosmetics, garments, and automobiles for which demand increases as they are launched into the market and after some time it becomes constant. Today most of the real-world decision- making problems in economic, technical, and environmental ones are such that the inventory related demands are not deterministic but imprecise in nature. Furthermore in real-life problems the parameters of the stochastic inventory models are also fuzzy, that is, not deterministic in nature and this is the case of application of fuzzy probability in the inventory models.

In 1965, the first publication in fuzzy set theory by Zadeh13showed the intention to accommodate uncertainty in the nonstochastic sense rather than the presence of random variables. After that fuzzy set theory has been applied in many fields including production- related areas. In the 1990s, several scholars began to develop models for inventory problems under fuzzy environment. Park14and Ishii and Konno15discussed the case of fuzzy cost coefficients. Roy and Maiti16developed a fuzzy economic order quantity EOQmodel with a constraint of fuzzy storage capacity. Chang and Yao17solved the economic reorder point with fuzzy backorders.

Some inventory problems with fuzzy shortage cost are analyzed by Katagiri and Ishii 18. A unified approach to fuzzy random variables is considered by Kr¨atschmer19; Kao and Hsu20discussed a single-period inventory model with fuzzy demand. Fergany and EI- Wakeel21considered the probabilistic single-item inventory problem with varying order cost under two linear constraints. Hala and EI-Saadani 22 analysed constrained single- period stochastic uniform inventory model with continuous distribution of demand and varying holding cost. Fuzzy models for single-period inventory problem were discussed by Li et al. 23. Banerjee and Roy24 considered application of the Intuitionistic Fuzzy Optimization in the constrained Multiobjective Stochastic Inventory Model. Banerjee and Roy25also discussed the single- and multiobjective stochastic inventory model in fuzzy environment. Lee and Yao26and Chang and Yao17investigated the economic production quantity model, and Lee and Yao 27 studied the EOQ model with fuzzy demands. A common characteristic of these studies is that shortages are backordered without extra costs.

Buckley28 introduced a new approach and applications of fuzzy probability and after that Buckley and Eslami29–31contributed three remarkable articles about uncertain probabilities. Hwang and Yao32discussed the independent fuzzy random variables and their applications. Formalization of fuzzy random variables is considered by Colubi et al.

33. Kr¨atschmer19analyzed a unified approach to fuzzy random variables. Luhandjula 34discussed a mathematical programming in the presence of fuzzy quantities and random variables.

(3)

Many researchers developed inventory models in which both the demand and deteriorating rate are constant. Although the constant demand assumption helps to simplify the problem, it is far from the actual situation where demand is always in change. In order to make research more practical, many researchers have studied other forms of demand.

Among them time-dependant demand has attracted considerable attention. In this paper exponentially increasing demand has been considered which is a general form of linear and nonlinear time-dependent demand. Various types of uncertainties and imprecision are inherent in real inventory problems. They are classically modeled using the approaches from the probability theory. However there are uncertainties that cannot be appropriately treated by usual probabilistic models. The questions how to define inventory optimization task in such environment and how to interpret optimal solution arise. This paper considers the modification of EPQ formula in the presence of imprecisely estimated parameters with fuzzy demands where backorders are permitted, yet a shortage cost is incurred. The demand rate is taken as triangular fuzzy number. Since the demand is fuzzy the average cost associated with inventory is fuzzy in nature. So the average cost in fuzzy sense is derived. The fuzzy model is defuzzified by usingα-cut of fuzzy number. This multiobjective problem is solved by Global Criteria method with the help of GRGGeneralized Reduced Gradienttechnique, and it is illustrated with the help of numerical example.

2. Preliminaries

For developing the mathematical model we are to introduce certain preliminary definitions and results which will be used later on.

Definition 2.1fuzzy number. A fuzzy subsetAof real numberRwith membership function μA:R → 0,1is called a fuzzy number if

aA is normal, that is, there exists an elementx0such thatμAx0 1 is normal, that is, there exists an elementx0such thatμAx0 1;

bAis convex, that is, μAλx1 1−λx2μAx1μAx2for allx1, x2Rand λ∈0,1;

cμAis upper semicontinuous;

dsuppA is bounded, here suppA supp{x∈R:μAx>0}.

Definition 2.2triangular fuzzy number. A Triangular Fuzzy numberTFNAis specified by the triplet a1,a2, a3 and is defined by its continuous membership function μAx, a continuous mappingμA :R → 0,1as follows:

μAx

⎧⎪

⎪⎪

⎪⎪

⎪⎨

⎪⎪

⎪⎪

⎪⎪

xa1 a2a1

ifa1x < a2 a3x

a3a2

ifa2xa3

0 otherwise.

2.1

(4)

Definition 2.3integration of a fuzzy function. Letfx be a fuzzy function froma, b ⊆R toRsuch thatfx is a fuzzy number, that is, a piecewise continuous normalized fuzzy set onR. Then integral of any continuousα-level curve offxovera, balways exists and the integral offx overa, bis then defined to be the fuzzy set

Ia, b b

a

fαxdx b

a

f αxdx, α

. 2.2

The determination of the integralIa, b becomes somewhat easier if the fuzzy function is assumed to beLRtype. We will therefore assume thatfx f 1x, f2x, f3xLRis a fuzzy number in LRrepresentation for allx ⊆ a, b, wheref1x,f2x, and f3x are assumed to be positive integrable functions ona, b. Dubois and Prade have shown that under these conditions

Ia, b b

a

f1xdx, b

a

f2xdx, b

a

f3xdx

LR

. 2.3

Definition 2.4α-cut of fuzzy number. Theα-cut of a fuzzy number is a crisp set which is defined asA α{x∈R:μAx≥α}.According to the definition of fuzzy number it is seen thatα-cut is a nonempty bounded closed interval; it can be denoted by

A

α ALα, ARα. 2.4

ALαandARαare the lower and upper bounds of the closed interval, where A

Lα inf

xR:μAx≥α , A

Rα sup

xR:μAx≥α .

2.5

Definition 2.5Global Criteria methodRao,nd. In the global criteria methodxis found by minimizing a preselected global criteria,FX, such as the sum of the squares of the relative deviations of the individual objective functions from the feasible ideal solutions. TheX is found by minimizing

FX k

i1

fi Xi

fiX fi

Xi p

subject to gjX≤0, j 1,2, . . . , m,

2.6

(5)

wherepis a constantas usual value ofpis 2andXiis the ideal solution for theith objective function. The solutionXi is obtained by minimizingfiXsubject to the constraintsgjX≤ 0,j1,2, . . . , m.

3. Notations and Assumptions

3.1. Notations

iT is the length of one cycle.

iiItis the inventory level at timet.

iiiRt is the induced fuzzy demand aroundRwhich is a function of time.

ivK1andK2are two production rates varies with demand.

vC1is the holding cost per unit time.

viC2is the shortage cost per unit time.

viiC3is the unit purchase cost.

viiiC4is the fixed ordering cost of inventory.

ixT is the cycle time.

xGC is the global criteria.

xiθ is the deterioration rate of finished items.

3.2. Assumptions

iThe demand functionRt is taken to be fuzzy function of time, that is,Rt ae bt. iiK1 andK2 are two demand depended production variables, that is,K1 β1Rt

andK2β2Rt, where β1andβ2are constants.

iiiReplenishment is instantaneous.

ivLead timei.e., the length between making of a decision to replenish an item and its actual addition to stockis assumed to be zero. The assumption is made so that the period of shortage is not affected.

vThe rate of deterioration at any time t > 0 follows the two-parameter Weibull distribution:θt γβtβ−1, whereγ 0 < γ <1is the scale parameter andβ>0is the shape parameter. The implication of the Weibull ratetwo parameteris that the items in inventory start deteriorating the instant they are received into inventory.

The rate of deterioration-time relationship is shown inFigure 1.

(6)

Time β= 1 β <1

Rateofdeteroratiion

1< β <2

Figure 1

Whenβ1,zt αa constantwhich is the case of exponential decay.

Whenβ <1 rate of deterioration is decreasing with time.

Whenβ >1, rate of deterioration is increasing with time.

viShortages are allowed and are fully backlogged.

4. Mathematical Modeling and Analysis

Here we assume that the production starts at timet 0 at the rateK1 and the stock attains a levelP at t z. Then at t z the production rate changes to K2 and continues up to ty, where the inventory level reaches the maximum levelQ. And the production stops at tyand the inventory gradually depletes to zero attt1mainly to meet the demands and partly for deterioration. Now shortages occur and accumulate to the levelSat timet t2. The production starts again at a rateK1attt2and the backlog is cleared at timetT when the stock is again zero. The cycle then repeats itself after timeT.

The model is represented byFigure 2.

The changes in the inventory level can be described by the following differential equations:

dIt

dt γβtβ−1It β1−1

aebt 0≤tz, dIt

dt γβtβ−1It β2−1

aebt zty, dIt

dt γβtβ−1It −aebt ytt1, dIt

dtae bt t1tt2, dIt

dtβ1−1

aebt t2tT

4.1

with initial conditionsIz P,Iy Q,It1 0,It2 S, andIT 0.

(7)

I(t)

P Q

z y t2 T

S t1

Figure 2

The differential equations4.1are fuzzy differential equations. To solve this differen- tial equation at first we take theα-cut then the differential equations reduces to

dI1t

dt γβtβ−1I1t β1−1

a αebt, 0≤tz, 4.2

dI1t

dt γβtβ−1I1t β1−1

aαebt, 0≤tz, 4.3

dI2t

dt γβtβ−1I2t β2−1

a αebt, zty, 4.4

dI2t

dt γβtβ−1I2t β2−1

aαebt, zty, 4.5

dI3t

dt γβtβ−1I3t −aαebt, ytt1, 4.6

dI3t

dt γβtβ−1I3t −a αebt, ytt1, 4.7 dI4t

dt −aαebt, t1tt2, 4.8

dI4t

dt −a αebt, t1tt2, 4.9 dI5t

dtβ1−1

aαebt, t2tT, 4.10

dI5t dt

β1−1

a αebt, t2tT, 4.11

(8)

where

Ii sup

xR:μI

ix≥α

, i1,2,3,4,5, Ii inf

xR:μIix≥α

, i1,2,3,4,5.

4.12

Similarlya α,aαhave usual meaning.

The solutions of the differential equations4.2to4.11are, respectively, represented by

I1t

β1−1 b

ebt−1

{a3αa3a2} −γβtβ 1 β1−1

a3

1 β , 0≤tz,

I1t

β1−1 b

ebt−1

{a1 αa2a1} −γβtβ 1 β1−1

a1

1 β , 0≤tz,

I2t P

1 γzβγtβ

β2−1{a3αa3a2}

ebtebz

bγtβ

β2−1 a3 b

ebtebz

γ{a3αa3a2}

β2−1 b tβ

ebt−1

γ{a3αa3a2}

β2−1 zβ

ebz−1 b

{a3αa3a2}γβ β2−1 1 β

zβ 1tβ 1

, zty,

I2t P

1 γzβγtβ

β2−1{a1 αa2a1}

ebtebz

bγtβ

β2−1 a1

b

ebtebz

{a1 αa2a1}γ β2−1

b tβ ebt−1

− {a1 αa2a1}γ β2−1

zβ ebz−1 b

{a1 αa2a1}βγ β2−1 1 β

zβ 1tβ 1

, zty,

I3t

{a1 αa2a1} b

ebt1ebt

a1γtβ b

ebt1ebt

{a1 αa2a1} b

γ

tβ1

ebt1−1

tβ

ebt−1

− {a1 αa2a1}γβ

tβ 11tβ 1

1 β , ytt1,

(9)

I3t

{a3αa3a2} b

ebt1ebt

a3γtβ b

ebt1ebt

{a3αa3a2} b

γ

tβ1

ebt1−1

tβ

ebt−1

− {a3αa3a2}γβ

tβ 11tβ 1

1 β , ytt1,

I4t −{a1 αa2a1}

ebtebt1

b , t1tt2, I4t −{a3αa3a2}

ebtebt1

b , t1tt2, I5t

β1−1

{a1 αa2a1}

ebTebt

b , t2tT,

I5t

β1−1

{a3αa3a2}

ebTebt

b , t2tT.

4.13

Therefore the upperα-cut of fuzzy stock holding cost is given by

HC

C1 T

z 0

I1tdt α1t1

z

I2tdt t1

α1t1

I3tdt

C1

T

A1 B1t1 C1tβ 11 D1tβ 21 E1tβ 31 F1ebt1t1 G1ebt1tβ 11 H1e1t1 I1ebt1 . 4.14 Calculations are shown inAppendix A.

And the lowerα-cut of fuzzy stock holding cost is given by

HC C1

T

z 0

I1tdt α1t1

z

I2tdt t1

α1t1

I3tdt

C1

T

A2 B2t1 C2tβ 11 D2tβ 21 E2tβ 31 F2ebt1t1 G2ebt1tβ 11 H2e1t1 I2ebt1 . 4.15 Calculations are shown inAppendix B.

1As demand is fuzzy in nature shortage cost is also fuzzy in nature.

(10)

Therefore the upperα-cut of shortage cost is given by

S.C C2

T

ηT t1

I4tdt T

ηT

I5tdt

C2

T

{a1 αa2a1} b

−ebηT

b ηTebt1 ebt1 bt1ebt1

C2

T

{a1 αa2a1} β1−1 b2

ebTTebT

bηTebT ebηT b

.

4.16

Also the lowerα-cut of shortage cost is given by

S.C C2

T

ηT t1

I4tdt T

ηT

I5tdt

C2

T

{a3αa3a2} b

−ebηT

b ηTebt1 ebt1

bt1ebt1

C2

T

{a3αa3a2} β1−1 b2

ebTTebT

bηTebT ebηT b

.

4.17

2Annual cost due to deteriorated unit is also fuzzy as demand is fuzzy quantity.

Therefore deterioration cost per year is given by

D.C C3 T

z 0

K1Rdt α1t1

z

K2Rdtt1

α1t1

R dt

C3 T

z 0

a

β1−1 ebtdt

α1t1

z

a

β2−1

ebtdtt1

α1t1

aebtdt

C3

T

z 0

f1tf2tf3tdt α1t1

z

φ123tdt− t1

α1t1

ϕ123tdt

,

4.18

where

f1t a1 β1−1

ebt, f2t a2 β1−1

ebt, f3t a3 β1−1

ebt, φ1t a1

β2−1

ebt, φ2t a2 β2−1

ebt, φ3t a3 β2−1

ebt, ϕ1t a3ebt, ϕ2t a2ebt, ϕ3t a1ebt.

4.19

(11)

Therefore upperα-cut of deterioration cost is given by

D.C C3 T

a3

β1−1 ebz−1

bαa3

β1−1 ebz−1 b

a2α

β1−1 ebz−1 b

a3

β2−1 b

e1t1ebz

αa3

β2−1 b

e1t1ebz αa2

β2−1 b

e1t1ebz

a1

b

ebt1e1t1

a1α b

ebt1e1t1 a2α b

ebt1e1t1 .

4.20 The lowerα-cut of deterioration cost is given by

D.C C3

T

a1

β1−1 ebz−1 b

αa2

β1−1 ebz−1 ba1α

β1−1 ebz−1 b

a1 β2−1

b

e1t1ebz αa2 β2−1

b

e1t1ebz

αa1

β2−1 b

e1t1ebz

a3 b

ebt1e1t1

a2α b

ebt1e1t1 a3α

b

ebt1e1t1 .

4.21

3The annual ordering costC4/T.

Therefore total variable cost per unit time is a fuzzy quantity and is defined by TVC

TVC TVC

, 4.22

where

TVC sup

xR:μTVCx≥α , TVCinf

xR:μTVCx≥α .

4.23

The upperα-cut of total variable cost per unit time is

TVC

C1 T

A1 B1t1 C1tβ 11 D1tβ 21 E1tβ 31 F1ebt1t1 G1ebt1tβ 11 H1e1t1 I1ebt1

C2

T

{a1 αa2a1} b

−ebηT

b ηTebt1 ebt1 bt1ebt1

(12)

C2

T

{a1 αa2a1} β1−1 b2

ebTTebT

bηTebT ebηT b

C3

T

a3

β1−1 ebz−1 bαa3

β1−1 ebz−1 b

a2α

β1−1 ebz−1 b

a3 β2−1

b

e1t1ebz

αa3 β2−1

b

e1t1ebz αa2

β2−1 b

e1t1ebz

a1 b

ebt1e1t1

a1α b

ebt1e1t1 a2α

b

ebt1e1t1 C4 T .

4.24

The lowerα-cut of total variable cost per unit time is

TVC C1

T

A2 B2t1 C2tβ 11 D2tβ 21 E2tβ 31 F2ebt1t1 G2ebt1tβ 11 H2e1t1 I2ebt1

C2 T

{a3αa3a2} b

−ebηT

b ηTebt1 ebt1 bt1ebt1

C2 T

{a3αa3a2} β1−1 b2

ebTTebT

bηTebT ebηT b

C3 T

a1

β1−1 ebz−1 b

αa2

β1−1 ebz−1 ba1α

β1−1 ebz−1 b

a1

β2−1 b

e1t1ebz αa2

β2−1 b

e1t1ebz

αa1 β2−1

b

e1t1ebz

a3

b

ebt1e1t1

a2α b

ebt1e1t1 a3α

b

ebt1e1t1 C4

T

.

4.25 The objective in this paper is to find an optimal cycle time to minimize the total variable cost per unit time.

Therefore this model mathematically can be written as

Minimize

TVC ,TVC

Subject to 0≤α≤1. 4.26

(13)

Therefore the problem is a multiobjective optimization problem. To convert it as a single- objective optimization problem we use global criteriaGCmethod.

Then the above problem reduces to

Minimize GC

Subject to 0≤α≤1. 4.27

5. Global Criteria Method

The model presented by4.26is a multiobjective model which is solved by Global Criteria GCMethod with the help of Generalized Reduced Gradient Technique.

The Multiobjective Nonlinear Integer ProgrammingMONLIPproblems are solved by Global Criteria method converting it to a single-objective optimization problem. The solution procedure is as follows.

Step 1. Solve the multiobjective programming problem4.26as a single-objective problem using only one objective at a time ignoring others.

Step 2. From the results of Step 1, determine the ideal objective vector, say TVC min, TVCminand the corresponding values of TVC max, TVCmax. Here, the ideal objective vector is used as a reference point. The problem is then to solve the following auxiliary problem:

MinGC Minimize

⎧⎨

TVC −TVC min TVC max−TVC min

Q

TVC−TVCmin TVCmax−TVCmin

Q

1/Q

, 5.1

where 1≤Q <∞. This method is also sometimes called Compromise Programming.

6. Numerical Example

We now consider a numerical example showing the utility of the model from practical point of view. According to the developed solution procedure of the proposed inventory system, the optimal solution has been obtained with the help of well-known generalized reduced gradient methodGRG. To illustrate the developed model, an example with the following data has been considered.

Leta1 100 units/month,a2 200 units/month,a3 100 units/month,C1 $.1 per unit,C2 $.5 per unit, C3 $.2 per unit, C4 $100 per order,b .2, α .25,β 1.9, β10.9,β21.9,Υ .001T 7 hrs.

Substituting the above parameters, Global CriteriaGCis obtained as

GC0.0926. 6.1

The compromise solutions are TVC $548.62, TVC$503.74.

(14)

7. Conclusion

In this paper the multiobjective problem is solved by Global Criteria method. In reality, in different systems, there are some parameters which are imprecise in nature. The present paper proposes a solution procedure to develop an EPQ inventory model with variable production rate and fuzzy demand. In most of the real life problem demand of a newly launched product is not known in advance. This justifies the introduction of fuzzy demand.

The technique for multiobjective optimization may be applied to the areas like environmental analysis, transportation, and so forth. In this paper we have taken two rates of production but this can be extended tonnumber of production ratespiduring the time when the inventory level goes fromi−1QtoiQi1,2, . . . , n, whereQis prefixed level.

Appendices

A. The Upper α-Cut of Fuzzy Stock Holding Cost

The upperα-cut of fuzzy stock holding cost is given by:

HC

C1 T

z 0

I1tdt α1t1

z

I2tdt t1

α1t1

I3tdt

C1

T

⎡⎣

β1−1

{a3αa3a2} b

ebz

bz−1 b

γβ β1−1

a3zβ 2 1 β

2 β

P

⎧⎨

α1t1 γα1zβt1γαβ 11 tβ 11

β 1 −zγzβ 1 γzβ 1 β 1

⎫⎬

β2−1

b {a3αa3a2} e1t1

bα1t1ebzebz b zebz

γa3

β2−1

βαβ 11 tβ 11 b

1 βγa3

β2−1 βzβ 1 b

1 β

γa3

β2−1 b

α1t1β 1zβ 1 b 2

α1t1β 2zβ 2

b2 6

α1t1β 3zβ 3 γa3

β2−1 ebz b

1 β

α1t1β 1zβ 1

γ{a3αa3a2} β2−1

βαβ 11 tβ 11 b

β 1 γ{a3αa3a2} β2−1

βzβ 1 b

1 β γ{a3αa3a2}

β2−1 b

α1t1β 1zβ 1 b 2

α1t1β 2zβ 2

b2 6

α1t1β 3zβ 3

(15)

γ{a3αa3a2} β2−1 b

β 1

α1t1β 1zβ 1

γ{a3 αa3a2} β2−1

zβebzα1t1z b

γ{a3αa3a2} β2−1

zβα1t1z b

γ{a3αa3a2β2−1 1 β

2 β

β 2

zβ 1α1t1αβ 21 tβ 21zβ 2 β 2

zβ 2

{a1 αa2a1} b

ebt1t1ebt1

bα1t1ebt1 e1t1 b

a1γebt1

tβ 11αβ 11 tβ 11 b

β 1

a1γβ b

1 β

tβ 11αβ 11 tβ 11 a1γ b

tβ 11αβ 11 tβ 11 b 2

tβ 21αβ 21 tβ 21

b2 6

tβ 31αβ 31 tβ 31

−{a1 αa2a12b

tβ 11 ebt1tβ 11α1tβ 11 ebt1 α1tβ 11 βγ{a1 αa2a12}

b

1 β

tβ 11αβ 11 tβ 11

γ{a1 αa2a1} b

tβ 11αβ 11 tβ 11 b 2

tβ 21αβ 21 tβ 21 b2

6

tβ 31αβ 31 tβ 31 γ{a1 αa2a1}

tβ 11αβ 11 tβ 11 b

1 β

−{a1 αa2a1} γβ 1 β

⎧⎨

tβ 21tβ 21

β 2−α1tβ 21 αβ 21 tβ 21 β 2

⎫⎬

HC

C1

T

z 0

I1tdt α1t1

z

I2tdt t1

α1t1

I3tdt

C1

T

A1 B1t1 C1tβ 11 D1tβ 21 E1tβ 31 F1ebt1t1 G1ebt1tβ 11 H1e1t1 I1ebt1 , A.1

参照

関連したドキュメント

We define the basic model for serial systems as follows: each stage controls its inventory by an installation base-stock policy; external demand follows Poisson process; the

We obtained the condition for ergodicity of the system, steady state system size probabilities, expected length of the busy period of the system, expected inventory level,

In this paper, for the first time an economic production quantity model for deteriorating items has been considered under inflation and time discounting over a stochastic time

It is suggested by our method that most of the quadratic algebras for all St¨ ackel equivalence classes of 3D second order quantum superintegrable systems on conformally flat

Massoudi and Phuoc 44 proposed that for granular materials the slip velocity is proportional to the stress vector at the wall, that is, u s gT s n x , T s n y , where T s is the

The theory of generalized ordinary differential equations enables one to inves- tigate ordinary differential, difference and impulsive equations from the unified standpoint...

If in the infinite dimensional case we have a family of holomorphic mappings which satisfies in some sense an approximate semigroup property (see Definition 1), and converges to

In [13], some topological properties of solutions set for (FOSPD) problem in the convex case are established, and in [15], the compactness of the solutions set is obtained in