
A decentralized supply chain network under uncertain cost is studied to obtain the optimal decisions of the enterprises in a situation in which the cost is uncertain. The supply chain network members adopt a chance-constrained approach to make decisions. The second-order cone-constrained variational inequality problem is used to construct the chance-constrained supply chain network equilibrium model. Then, the existence and uniqueness properties of the proposed equilibrium model are discussed under some mild assumptions. For the discontinuous functions in the feasible region of the model, the proposed model is converted to a second-order cone complementarity problem. The numerical results show that the uncertainty and risk attitude of retailers and manufacturers have different effects on supply chain network members. When the risk attitude is high, a small change in the risk attitude will significantly change all decisions of supply chain members. If the supply chain member is affected by the uncertainty positively, its profit will increase as its risk attitude increases. Moreover, it is appropriate to adopt a chance-constrained approach when the supply chain members can estimate the distributions of the competitor's strategies.
Citation: Shuai Huang, Youwu Lin, Jing Zhang, Pei Wang. Chance-constrained approach for decentralized supply chain network under uncertain cost[J]. AIMS Mathematics, 2023, 8(5): 12217-12238. doi: 10.3934/math.2023616
[1] | Yayun Fu, Mengyue Shi . A conservative exponential integrators method for fractional conservative differential equations. AIMS Mathematics, 2023, 8(8): 19067-19082. doi: 10.3934/math.2023973 |
[2] | Yong-Chao Zhang . Least energy solutions to a class of nonlocal Schrödinger equations. AIMS Mathematics, 2024, 9(8): 20763-20772. doi: 10.3934/math.20241009 |
[3] | Tingting Ma, Yuehua He . An efficient linearly-implicit energy-preserving scheme with fast solver for the fractional nonlinear wave equation. AIMS Mathematics, 2023, 8(11): 26574-26589. doi: 10.3934/math.20231358 |
[4] | Karmina K. Ali, Resat Yilmazer . Discrete fractional solutions to the effective mass Schrödinger equation by mean of nabla operator. AIMS Mathematics, 2020, 5(2): 894-903. doi: 10.3934/math.2020061 |
[5] | Erdal Bas, Ramazan Ozarslan . Theory of discrete fractional Sturm–Liouville equations and visual results. AIMS Mathematics, 2019, 4(3): 593-612. doi: 10.3934/math.2019.3.593 |
[6] | Dengfeng Lu, Shuwei Dai . On a class of three coupled fractional Schrödinger systems with general nonlinearities. AIMS Mathematics, 2023, 8(7): 17142-17153. doi: 10.3934/math.2023875 |
[7] | Mubashir Qayyum, Efaza Ahmad, Hijaz Ahmad, Bandar Almohsen . New solutions of time-space fractional coupled Schrödinger systems. AIMS Mathematics, 2023, 8(11): 27033-27051. doi: 10.3934/math.20231383 |
[8] | Xiaojun Zhou, Yue Dai . A spectral collocation method for the coupled system of nonlinear fractional differential equations. AIMS Mathematics, 2022, 7(4): 5670-5689. doi: 10.3934/math.2022314 |
[9] | Zunyuan Hu, Can Li, Shimin Guo . Fast finite difference/Legendre spectral collocation approximations for a tempered time-fractional diffusion equation. AIMS Mathematics, 2024, 9(12): 34647-34673. doi: 10.3934/math.20241650 |
[10] | Xiao-Yu Li, Yu-Lan Wang, Zhi-Yuan Li . Numerical simulation for the fractional-in-space Ginzburg-Landau equation using Fourier spectral method. AIMS Mathematics, 2023, 8(1): 2407-2418. doi: 10.3934/math.2023124 |
A decentralized supply chain network under uncertain cost is studied to obtain the optimal decisions of the enterprises in a situation in which the cost is uncertain. The supply chain network members adopt a chance-constrained approach to make decisions. The second-order cone-constrained variational inequality problem is used to construct the chance-constrained supply chain network equilibrium model. Then, the existence and uniqueness properties of the proposed equilibrium model are discussed under some mild assumptions. For the discontinuous functions in the feasible region of the model, the proposed model is converted to a second-order cone complementarity problem. The numerical results show that the uncertainty and risk attitude of retailers and manufacturers have different effects on supply chain network members. When the risk attitude is high, a small change in the risk attitude will significantly change all decisions of supply chain members. If the supply chain member is affected by the uncertainty positively, its profit will increase as its risk attitude increases. Moreover, it is appropriate to adopt a chance-constrained approach when the supply chain members can estimate the distributions of the competitor's strategies.
Fractional calculus is a popular subject because of having a lot of application areas of theoretical and applied sciences, like engineering, physics, biology, etc. Discrete fractional calculus is more recent area than fractional calculus and it was first defined by Diaz–Osler [1], Miller–Ross [2] and Gray–Zhang [3]. More recently, the theory of discrete fractional calculus have begun to develop rapidly with Goodrich–Peterson [4], Baleanu et al. [5,6], Ahrendt et al. [7], Atici–Eloe [8,9], Anastassiou [10], Abdeljawad et al. [11,12,13,14,15,16], Hein et al. [17] and Cheng et al. [18], Mozyrska [19] and so forth [20,21,22,23,24,25].
Fractional Sturm–Liouville differential operators have been studied by Bas et al. [26,27], Klimek et al.[28], Dehghan et al. [29]. Besides that, Sturm–Liouville differential and difference operators were studied by [30,31,32,33]. In this study, we define DFHA operators and prove the self–adjointness of DFHA operator, some spectral properties of the operator.
More recently, Almeida et al. [34] have studied discrete and continuous fractional Sturm–Liouville operators, Bas–Ozarslan [35] have shown the self–adjointness of discrete fractional Sturm–Liouville operators and proved some spectral properties of the problem.
Sturm–Liouville equation having hydrogen atom potential is defined as follows
d2Rdr2+ardRdr−ℓ(ℓ+1)r2R+(E+ar)R=0(0<r<∞). |
In quantum mechanics, the study of the energy levels of the hydrogen atom leads to this equation. Where R is the distance from the mass center to the origin, ℓ is a positive integer, a is real number E is energy constant and r is the distance between the nucleus and the electron.
The hydrogen atom is a two–particle system and it composes of an electron and a proton. Interior motion of two particles around the center of mass corresponds to the movement of a single particle by a reduced mass. The distance between the proton and the electron is identified r and r is given by the orientation of the vector pointing from the proton to the electron. Hydrogen atom equation is defined as Schrödinger equation in spherical coordinates and in consequence of some transformations, this equation is defined as
y′′+(λ−l(l+1)x2+2x−q(x))y=0. |
Spectral theory of hydrogen atom equation is studied by [39,40,41]. Besides that, we can observe that hydrogen atom differential equation has series solution as follows ([39], p.268)
y(x)=a0xl+1{1−k−l−11!(2l+2).2xk+(k−l−1)(k−l−2)2!(2l+2)(2l+3)(2xk)2+…+(−1)n(k−l−1)(k−l−2)…3.2.1(k−1)!(2l+2)(2l+3)…(2l+n)(2xk)n},k=1,2,… | (1.1) |
Recently, Bohner and Cuchta [36,37] studied some special integer order discrete functions, like Laguerre, Hermite, Bessel and especially Cuchta mentioned the difficulty in obtaining series solution of discrete special functions in his dissertation ([38], p.100). In this regard, finding series solution of DFHA equations is an open problem and has some difficulties in the current situation. For this reason, we study to obtain solutions of DFHA eq.s in a different way with representation of solutions.
In this study, we investigate DFHA equation in Riemann–Liouville and Grü nwald–Letnikov sense. The aim of this study is to contribute to the spectral theory of DFHA operator and behaviors of eigenfunctions and also to obtain the solution of DFHA equation.
We investigate DFHA equation in three different ways;
i) (nabla left and right) Riemann–Liouville (R–L)sense,
L1x(t)=∇μa(b∇μx(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t), 0<μ<1, |
ii) (delta left and right) Grünwald–Letnikov (G–L) sense,
L2x(t)=Δμ−(Δμ+x(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t), 0<μ<1, |
iii) (nabla left) Riemann–Liouville (R–L)sense,
L3x(t)=∇μa(∇μax(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t), 0<μ<1. |
Definition 2.1. [42] Falling and rising factorial functions are defined as follows respectively
tα_=Γ(t+1)Γ(t−α+1), | (2.1) |
t¯α=Γ(t+α)Γ(t), | (2.2) |
where Γ is the gamma function, α∈R.
Remark 2.1. Delta and nabla operators hold the following properties
Δtα_=αtα−1_,∇t¯α=αt¯α−1. | (2.3) |
Definition 2.2. [2,8,11] Nabla fractional sum operators are given as below,
(i) The left fractional sum of order μ>0 is defined by
∇−μax(t)=1Γ(μ)t∑s=a+1(t−ρ(s))¯μ−1x(s), t∈Na+1, | (2.4) |
(ii) The right fractional sum of order μ>0 is defined by
b∇−μx(t)=1Γ(μ)b−1∑s=t(s−ρ(t))¯μ−1x(s), t∈ b−1N, | (2.5) |
where ρ(t)=t−1 is called backward jump operators, Na={a,a+1,...}, bN={b,b−1,...}.
Definition 2.3. [12,14] Nabla fractional difference operators are as follows,
(i) The left fractional difference of order μ>0 is defined by
∇μax(t)=∇n∇−(n−μ)ax(t)=∇nΓ(n−μ)t∑s=a+1(t−ρ(s))¯n−μ−1x(s), t∈Na+1, | (2.6) |
(ii) The right fractional difference of order μ>0 is defined by
b∇μx(t)=(−1)n∇n∇−(n−μ)ax(t)=(−1)nΔnΓ(n−μ)b−1∑s=t(s−ρ(t))¯n−μ−1x(s), t∈ b−1N. | (2.7) |
Fractional differences in (2.6−2.7) are called the Riemann–Liouville (R–L) definition of the μ-th order nabla fractional difference.
Definition 2.4. [1,18] Fractional difference operators are given as follows
(i) The delta left fractional difference of order μ, 0<μ≤1, is defined by
Δμ−x(t)=1hμt∑s=0(−1)sμ(μ−1)...(μ−s+1)s!x(t−s), t=1,...,N. | (2.8) |
(ii) The delta right fractional difference of order μ, 0<μ≤1, is defined by
Δμ+x(t)=1hμN−t∑s=0(−1)sμ(μ−1)...(μ−s+1)s!x(t+s), t=0,..,N−1, | (2.9) |
fractional differences in (2.8−2.9) are called the Grünwald–Letnikov (G–L) definition of the μ-th order delta fractional difference.
Definition 2.5 [14] Integration by parts formula for R–L nabla fractional difference operator is defined by, u is defined on bN and v is defined on Na,
b−1∑s=a+1u(s)∇μav(s)=b−1∑s=a+1v(s)b∇μu(s). | (2.10) |
Definition 2.6. [34] Integration by parts formula for G–L delta fractional difference operator is defined by, u, v is defined on {0,1,...,n}, then
n∑s=0u(s)Δμ−v(s)=n∑s=0v(s)Δμ+u(s). | (2.11) |
Definition 2.7. [17] f:Na→R, s∈ℜ, Laplace transform is defined as follows,
La{f}(s)=∞∑k=1(1−s)k−1f(a+k), |
where ℜ=C∖{1} and ℜ is called the set of regressive (complex) functions.
Definition 2.8. [17] Let f,g:Na→R, all t∈Na+1, convolution of f and g is defined as follows
(f∗g)(t)=t∑s=a+1f(t−ρ(s)+a)g(s), |
where ρ(s) is the backward jump function defined in [42] as
ρ(s)=s−1. |
Theorem 2.1. [17] f,g:Na→R, convolution theorem is expressed as follows,
La{f∗g}(s)=La{f}La{g}(s). |
Lemma 2.1. [17] f:Na→R, the following property is valid,
La+1{f}(s)=11−sLa{f}(s)−11−sf(a+1). |
Theorem 2.2. [17] f:Na→R, 0<μ<1, Laplace transform of nabla fractional difference
La+1{∇μaf}(s)=sμLa+1{f}(s)−1−sμ1−sf(a+1),t∈Na+1. |
Definition 2.9. [17] For |p|<1, α>0, β∈R and t∈Na, Mittag–Leffler function is defined by
Ep,α,β(t,a)=∞∑k=0pk(t−a)¯αk+βΓ(αk+β+1). |
Theorem 2.3. [17] For |p|<1, α>0, β∈R, |1−s|<1 and |s|α>p, Laplace transform of Mittag–Leffler function is as follows,
La+1{Ep,α,β(.,a)}(s)=sα−β−1sα−p. |
Let us consider equations in three different forms;
i) L1 DFHA operator L1 is defined in (nabla left and right) R–L sense,
L1x(t)=∇μa(p(t)b∇μx(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t), 0<μ<1, | (3.1) |
where l is a positive integer or zero, q(t)+2t−l(l+1)t2 are named potential function., λ is the spectral parameter, t∈[a+1,b−1], x(t)∈l2[a+1,b−1], a>0.
ii) L2 DFHA operator L2 is defined in (delta left and right) G–L sense,
L2x(t)=Δμ−(p(t)Δμ+x(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t), 0<μ<1, | (3.2) |
where p,q,l,λ is as defined above, t∈[1,n], x(t)∈l2[0,n].
iii) L3 DFHA operator L3 is defined in (nabla left) R–L sense,
L3x(t)=∇μa(∇μax(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t), 0<μ<1, | (3.3) |
p,q,l,λ is as defined above, t∈[a+1,b−1], a>0.
Theorem 3.1. DFHA operator L1 is self–adjoint.
Proof.
u(t)L1v(t)=u(t)∇μa(p(t)b∇μv(t))+u(t)(l(l+1)t2−2t+q(t))v(t), | (3.4) |
v(t)L1u(t)=v(t)∇μa(p(t)b∇μu(t))+v(t)(l(l+1)t2−2t+q(t))u(t). | (3.5) |
Subtracting (16−17) from each other
u(t)L1v(t)−v(t)L1u(t)=u(t)∇μa(p(t)b∇μv(t))−v(t)∇μa(p(t)b∇μu(t)) |
and applying definite sum operator to both side of the last equality, we have
b−1∑s=a+1(u(s)L1v(s)−v(s)L1u(s))=b−1∑s=a+1u(s)∇μa(p(s)b∇μv(s))−b−1∑s=a+1v(s)∇μa(p(s)b∇μu(s)). | (3.6) |
Applying the integration by parts formula (2.10) to right hand side of (18), we have
b−1∑s=a+1(u(s)L1v(s)−v(s)L1u(s))=b−1∑s=a+1p(s)b∇μv(s)b∇μu(s)−b−1∑s=a+1p(s)b∇μu(s)b∇μv(s)=0, |
⟨L1u,v⟩=⟨u,L1v⟩. |
The proof completes.
Theorem 3.2. Eigenfunctions, corresponding to distinct eigenvalues, of the equation (3.2) are orthogonal.
Proof. Assume that λα and λβ are two different eigenvalues corresponds to eigenfunctions u(n) and v(n) respectively for the equation (3.1),
∇μa(p(t)b∇μu(t))+(l(l+1)t2−2t+q(t))u(t)−λαu(t)=0,∇μa(p(t)b∇μv(t))+(l(l+1)t2−2t+q(t))v(t)−λβv(t)=0, |
Multiplying last two equations to v(n) and u(n) respectively, subtracting from each other and applying sum operator, since the self–adjointness of the operator L1, we get
(λα−λβ)b−1∑s=a+1r(s)u(s)v(s)=0, |
since λα≠λβ,
b−1∑s=a+1r(s)u(s)v(s)=0,⟨u(t),v(t)⟩=0, |
and the proof completes.
Theorem 3.3. All eigenvalues of the equation (3.1) are real.
Proof. Assume λ=α+iβ, since the self–adjointness of the operator L1, we have
⟨L1u,u⟩=⟨u,L1u⟩,⟨λu,u⟩=⟨u,λu⟩, |
(λ−¯λ)⟨u,u⟩=0 |
Since ⟨u,u⟩r≠0,
λ=¯λ |
and hence β=0. So, the proof is completed.
Self–adjointness of L2 DFHA operator G–L sense, reality of eigenvalues and orthogonality of eigenfunctions of the equation 3.2 can be proven in a similar way to the Theorem 3.1–3.2–3.3 by means of Definition 2.5.
Theorem 3.4.
L3x(t)=∇μa(∇μax(t))+(l(l+1)t2−2t+q(t))x(t)=λx(t),0<μ<1, | (3.7) |
x(a+1)=c1,∇μax(a+1)=c2, | (3.8) |
where p(t)>0, r(t)>0, q(t) is defined and real valued, λ is the spectral parameter. The sum representation of solution of the problem (3.7)−(3.8) is given as follows,
x(t)=c1((1+l(l+1)(a+1)2−2a+1+q(a+1))Eλ,2μ,μ−1(t,a)−λEλ,2μ,2μ−1(t,a))+c2(Eλ,2μ,2μ−1(t,a)−Eλ,2μ,μ−1(t,a))−t∑s=a+1Eλ,2μ,2μ−1(t−ρ(s)+a)(l(l+1)s2−2s+q(s))x(s). | (3.9) |
Proof. Taking Laplace transform of the equation (3.7) by Theorem 2.2 and take (l(l+1)t2−2t+q(t))x(t)=g(t),
La+1{∇μa(∇μax)}(s)+La+1{g}(s)=λLa+1{x}(s),=sμLa+1{∇μax}(s)−1−sμ1−sc2=λLa+1{x}(s)−La+1{g}(s),=sμ(sμLa+1{x}(s)−1−sμ1−sc1)−1−sμ1−sc2=λLa+1{x}(s)−La+1{g}(s), |
=La+1{x}(s)=1−sμ1−s1s2μ−λ(sμc1+c2)−1s2μ−λLa+1{g}(s). |
Using Lemma 2.1, we have
La{x}(s)=c1(sμ−λs2μ−λ)−1−ss2μ−λ(11−sLa{g}(s)−11−sg(a+1))+c2(1−sμs2μ−λ). | (3.10) |
Now, taking inverse Laplace transform of the equation (3.10) and applying convolution theorem, then we have the representation of solution of the problem (3.7)−(3.8), |λ|<1, |1−s|<1 and |s|α>λ from Theorem 2.3., i.e.
L−1a{sμs2μ−λ}=Eλ,2μ,μ−1(t,a),L−1a{1s2μ−λ}=Eλ,2μ,2μ−1(t,a), |
L−1a{1s2μ−λLa{q(s)x(s)}}=t∑s=a+1Eλ,2μ,2μ−1(t−ρ(s)+a)q(s)x(s). |
Consequently, we have sum representation of solution for DFHA problem 3.7–3.8
x(t)=c1((1+l(l+1)(a+1)2−2a+1+q(a+1))Eλ,2μ,μ−1(t,a)−λEλ,2μ,2μ−1(t,a))+c2(Eλ,2μ,2μ−1(t,a)−Eλ,2μ,μ−1(t,a))−t∑s=a+1Eλ,2μ,2μ−1(t−ρ(s)+a)(l(l+1)s2−2s+q(s))x(s). |
Presume that c1=1,c2=0,a=0 in the representation of solution (3.9) and hence we may observe the behaviors of solutions in following figures (Figures 1–7) and tables (Tables 1–3);
x(t) | μ=0.3 | μ=0.35 | μ=0.4 | μ=0.45 | μ=0.5 |
x(1) | 1 | 1 | 1 | 1 | 1 |
x(2) | 0.612 | 0.714 | 1.123 | 0.918 | 1.020 |
x(3) | 0.700 | 0.900 | 1.515 | 1.370 | 1.641 |
x(5) | 0.881 | 1.336 | 2.402 | 2.747 | 3.773 |
x(7) | 1.009 | 1.740 | 3.352 | 4.566 | 7.031 |
x(9) | 1.099 | 2.100 | 4.332 | 6.749 | 11.461 |
x(12) | 1.190 | 2.570 | 5.745 | 10.623 | 20.450 |
x(15) | 1.249 | 2.975 | 6.739 | 15.149 | 32.472 |
x(16) | 1.264 | 3.098 | 7.235 | 16.793 | 37.198 |
x(18) | 1.289 | 3.330 | 8.233 | 20.279 | 47.789 |
x(20) | 1.309 | 3.544 | 9.229 | 24.021 | 59.967 |
x(t) | q(t)=1 | q(t)=t | q(t)=√t |
x(1) | 1 | 1 | 1 |
x(2) | 7.37∗10−17 | 4.41∗10−17 | 5.77∗10−17 |
x(3) | −0.131 | −0.057 | −0.088 |
x(5) | −0.123 | −0.018 | −0.049 |
x(7) | −0.080 | −0.006 | −0.021 |
x(9) | −0.050 | −0.003 | −0.011 |
x(12) | −0.028 | −0.001 | −0.005 |
x(15) | −0.017 | −0.0008 | −0.003 |
x(16) | −0.015 | −0.0006 | −0.0006 |
x(18) | −0.012 | −0.0005 | −0.002 |
x(20) | −0.010 | −0.0003 | −0.001 |
x(t) | λ=0.1 | λ=0.11 | λ=0.12 |
x(1) | 1 | 1 | 1 |
x(2) | 1 | 1.025 | 1.052 |
x(3) | 1.668 | 1.751 | 1.841 |
x(5) | 3.876 | 4.216 | 4.595 |
x(7) | 7.243 | 8.107 | 9.095 |
x(9) | 11.941 | 13.707 | 12.130 |
x(12) | 22.045 | 26.197 | 25.237 |
x(15) | 36.831 | 45.198 | 46.330 |
x(16) | 43.042 | 53.369 | 55.687 |
x(18) | 57.766 | 73.092 | 78.795 |
x(20) | 76.055 | 98.154 | 127.306 |
We have analyzed DFHA equation in Riemann–Liouville and Grü nwald–Letnikov sense. Self–adjointness of the DFHA operator is presented and also, we have proved some significant spectral properties for instance, orthogonality of distinct eigenfunctions, reality of eigenvalues. Moreover, we give sum representation of the solutions for DFHA problem and find the solutions of the problem. We have carried out simulation analysis with graphics and tables. The aim of this paper is to contribute to the theory of hydrogen atom fractional difference operator.
We observe the behaviors of solutions by changing the order of the derivative μ in Figure 1 and Figure 5, by changing the potential function q(t) in Figure 2, we compare solutions under different λ eigenvalues in Figure 3, and Figure 7, also we observe the solutions by changing μ with a specific eigenvalue in Figure 4 and by changing l values in Figure 6.
We have shown the solutions by changing the order of the derivative μ in Table 1, by changing the potential function q(t) and λ eigenvalues in Table 2, Table 3.
The authors would like to thank the editor and anonymous reviewers for their careful reading of our manuscript and their many insightful comments and suggestions.
The authors declare no conflict of interest.
[1] |
A. Nagurney, J. Dong, D. Zhang, A supply chain network equilibrium model, Transport. Res. E-Log., 38 (2002), 281–303. https://doi.org/10.1016/S1366-5545(01)00020-5 doi: 10.1016/S1366-5545(01)00020-5
![]() |
[2] |
A. Nagurney, Z. G. Liu, M. G. Cojocaru, P. Daniele, Dynamic electric power supply chains and transportation networks: an evolutionary variational inequality formulation, Transport. Res. E-Log., 43 (2007), 624–646. https://doi.org/10.1016/j.tre.2006.03.002 doi: 10.1016/j.tre.2006.03.002
![]() |
[3] |
M. Yu, A. Nagurney, Competitive food supply chain networks with application to fresh produce, Eur. J. Oper. Res., 224 (2013), 273–282. https://doi.org/10.1016/j.ejor.2012.07.033 doi: 10.1016/j.ejor.2012.07.033
![]() |
[4] |
S. Saberi, J. M. Cruz, J. Sarkis, A. Nagurney, A competitive multiperiod supply chain network model with freight carriers and green technology investment option, Eur. J. Oper. Res., 266 (2018), 934–949. https://doi.org/10.1016/j.ejor.2017.10.043 doi: 10.1016/j.ejor.2017.10.043
![]() |
[5] |
Y. X. Yang, X. Xu, A differential game model for closed-loop supply chain participants under carbon emission permits, Comput. Ind. Eng., 135 (2019), 1077–1090. https://doi.org/10.1016/j.cie.2019.03.049 doi: 10.1016/j.cie.2019.03.049
![]() |
[6] |
A. Nagurney, E. A. Flores, C. Soylu, A Generalized Nash Equilibrium network model for post-disaster humanitarian relief, Transport. Res. E-log., 95 (2016), 1–18. https://doi.org/10.1016/j.tre.2016.08.005 doi: 10.1016/j.tre.2016.08.005
![]() |
[7] |
A. Nagurney, Attracting international migrant labor: investment optimization to alleviate supply chain labor shortages, Operations Research Perspectives, 9 (2022), 100233. https://doi.org/10.1016/j.orp.2022.100233 doi: 10.1016/j.orp.2022.100233
![]() |
[8] |
D. Bertsimas, M. Sim, The price of robustness, Oper. Res., 52 (2004), 35–53. https://doi.org/10.1287/opre.1030.0065 doi: 10.1287/opre.1030.0065
![]() |
[9] |
Y. T. Guo, G. J. Ye, W. Liu, D. F. Zhao, S. Treanţă, On symmetric gH-derivative: applications to dual interval-valued optimization problems, Chaos Soliton. Fract., 158 (2022), 112068. https://doi.org/10.1016/j.chaos.2022.112068 doi: 10.1016/j.chaos.2022.112068
![]() |
[10] |
S. Treanţă, On a class of constrained interval-valued optimization problems governed by mechanical work cost functionals, J. Optim. Theory Appl., 188 (2021), 913–924. https://doi.org/10.1007/s10957-021-01815-0 doi: 10.1007/s10957-021-01815-0
![]() |
[11] |
S. Treanţă, Efficiency in uncertain variational control problems, Neural Comput. Appl., 33 (2021), 5719–5732. https://doi.org/10.1007/s00521-020-05353-0 doi: 10.1007/s00521-020-05353-0
![]() |
[12] |
S. Treanţă, Saddle-point optimality criteria involving (p, b, d)-invexity and (p, b, d)-pseudoinvexity in interval-valued optimisation problems, Int. J. Control, 95 (2022), 1042–1050. https://doi.org/10.1080/00207179.2020.1837960 doi: 10.1080/00207179.2020.1837960
![]() |
[13] |
H. Deng, Y. X. Li, Z. M. Wan, Z. Wan, Partially smoothing and gradient-based algorithm for optimizing the VMI system with competitive retailers under random demands, Math. Probl. Eng., 2020 (2020), 1–18. https://doi.org/10.1155/2020/3687471 doi: 10.1155/2020/3687471
![]() |
[14] |
S. H. Deng, Z. Wan, Y. J. Zhou, Optimization model and solution method for dynamically correlated two-product newsvendor problems based on Copula, Discrete and Continuous Dynamical Systems-S, 13 (2020), 1637–1652. https://doi.org/10.3934/dcdss.2020096 doi: 10.3934/dcdss.2020096
![]() |
[15] |
C. M. Hu, Z. M. Wan, S. H. Zhu, Z. Wan, An integrated stochastic model and algorithm for constrained multi-item newsvendor problems by two-stage decision-making approach, Math. Comput. Simulat., 193 (2022), 280–300. https://doi.org/10.1016/j.matcom.2021.10.018 doi: 10.1016/j.matcom.2021.10.018
![]() |
[16] |
Z. M. Wan, S. H. Zhu, Z. Wan, An integrated stochastic model and algorithm for multi-product newsvendor problems, Int. J. Model. Simul. Sc., 11 (2020), 2050027. https://doi.org/10.1142/S1793962320500270 doi: 10.1142/S1793962320500270
![]() |
[17] |
A. H. Golsefidi, M. R. A. Jokar, A robust optimization approach for the production-inventory-routing problem with simultaneous pickup and delivery, Comput. Ind. Eng., 143 (2020), 106388. https://doi.org/10.1016/j.cie.2020.106388 doi: 10.1016/j.cie.2020.106388
![]() |
[18] |
X. T. Shang, G. Q. Zhang, B. Jia, M. Almanaseer, The healthcare supply location-inventory-routing problem: a robust approach, Transport. Res. E-log., 158 (2022), 102588. https://doi.org/10.1016/j.tre.2021.102588 doi: 10.1016/j.tre.2021.102588
![]() |
[19] |
Y. Shi, Y. J. Zhou, W. H. Ye, Q. Q. Zhao, A relative robust optimization for a vehicle routing problem with time-window and synchronized visits considering greenhouse gas emissions, J. Clean. Prod., 275 (2020), 124112. https://doi.org/10.1016/j.jclepro.2020.124112 doi: 10.1016/j.jclepro.2020.124112
![]() |
[20] |
Y. T. Guo, G. J. Ye, W. Liu, D. F. Zhao, S. Treanţă, Optimality conditions and duality for a class of generalized convex interval-valued optimization problems, Mathematics, 9 (2021), 2979. https://doi.org/10.3390/math9222979 doi: 10.3390/math9222979
![]() |
[21] |
S. Treanţă, LU-optimality conditions in optimization problems with mechanical work objective functionals, IEEE T. Neur. Ne. Lear., 33 (2022), 4971–4978. https://doi.org/10.1109/TNNLS.2021.3066196 doi: 10.1109/TNNLS.2021.3066196
![]() |
[22] |
S. Treanţă, Characterization results of solutions in interval-valued optimization problems with mixed constraints, J. Global. Optim., 82 (2022), 951–964. https://doi.org/10.1007/s10898-021-01049-4 doi: 10.1007/s10898-021-01049-4
![]() |
[23] |
J. Dong, D. Zhang, A. Nagurney, A supply chain network equilibrium model with random demands, Eur. J. Oper. Res., 156 (2004), 194–212. https://doi.org/10.1016/S0377-2217(03)00023-7 doi: 10.1016/S0377-2217(03)00023-7
![]() |
[24] |
Z. G. Liu, A. Nagurney, Supply chain networks with global outsourcing and quick-response production under demand and cost uncertainty, Ann. Oper. Res., 208 (2013), 251–289. https://doi.org/10.1007/s10479-011-1006-0 doi: 10.1007/s10479-011-1006-0
![]() |
[25] |
Y. X. Xiao, R. Q. Zhang, Supply chain network equilibrium considering coordination between after-sale service and product quality, Comput. Ind. Eng., 175 (2023), 108848. https://doi.org/10.1016/j.cie.2022.108848 doi: 10.1016/j.cie.2022.108848
![]() |
[26] |
Y. Zhou, C. K. Chan, K. H. Wong, A multi-period supply chain network equilibrium model considering retailers' uncertain demands and dynamic loss-averse behaviors, Transport. Res. E-log., 118 (2018), 51–76. https://doi.org/10.1016/j.tre.2018.06.006 doi: 10.1016/j.tre.2018.06.006
![]() |
[27] |
Z. Wan, H. Wu, L. Dai, A polymorphic uncertain equilibrium model and its deterministic equivalent formulation for decentralized supply chain management, Appl. Math. Model., 58 (2018), 281–299. https://doi.org/10.1016/j.apm.2017.06.028 doi: 10.1016/j.apm.2017.06.028
![]() |
[28] |
Y. Kou, Z. Wan, A new data-driven robust optimization approach to multi-item newsboy problems, J. Ind. Manag. Optim., 19 (2022), 197–223. https://doi.org/10.3934/jimo.2021180 doi: 10.3934/jimo.2021180
![]() |
[29] |
S. Zokaee, A. Jabbarzadeh, B. Fahimnia, S. J. Sadjadi, Robust supply chain network design: an optimization model with real world application, Ann. Oper. Res., 257 (2017), 15–44. https://doi.org/10.1007/s10479-014-1756-6 doi: 10.1007/s10479-014-1756-6
![]() |
[30] |
Z. M. Liu, Data-driven two-stage sparse distributionally robust risk optimization model for location allocation problems under uncertain environment, AIMS Math., 8 (2022), 2910–2939. https://doi.org/10.3934/math.2023152 doi: 10.3934/math.2023152
![]() |
[31] |
T. Hirano, Y. Narushima, Robust supply chain network equilibrium model, Transport. Sci., 53 (2019), 1196–1212. https://doi.org/10.1287/trsc.2018.0843 doi: 10.1287/trsc.2018.0843
![]() |
[32] | P. Rosenblum, B. Kilcourse, Supply chain management 2018: In service of the customer, Retail Systems Research, 2018. Available from: https://www.rsrresearch.com/research/supply-chain-management-2018. |
[33] |
S. Ahmadi, S. H. Amin, An integrated chance-constrained stochastic model for a mobile phone closed-loop supply chain network with supplier selection, J. Clean. Prod., 226 (2019), 988–1003. https://doi.org/10.1016/j.jclepro.2019.04.132 doi: 10.1016/j.jclepro.2019.04.132
![]() |
[34] |
M. A. Quddus, S. Chowdhury, M. Marufuzzaman, F. Yu, L. K. Bian, A two-stage chance-constrained stochastic programming model for a bio-fuel supply chain network, Int. J. Prod. Econ., 195 (2018), 27–44. https://doi.org/10.1016/j.ijpe.2017.09.019 doi: 10.1016/j.ijpe.2017.09.019
![]() |
[35] |
T. Dinh, R. Fukasawa, J. Luedtke, Exact algorithms for the chance-constrained vehicle routing problem, Math. Program., 172 (2018), 105–138. https://doi.org/10.1007/s10107-017-1151-6 doi: 10.1007/s10107-017-1151-6
![]() |
[36] |
P. B. Chen, G. H. Lin, X. Zhu, F. S. Bai, Smoothing Newton method for nonsmooth second-order cone complementarity problems with application to electric power markets, J. Glob. Optim., 80 (2021), 635–659. https://doi.org/10.1007/s10898-021-00993-5 doi: 10.1007/s10898-021-00993-5
![]() |
[37] |
Z. J. Hao, Z. P. Wan, X. N. Chi, Z. F. Jin, Generalized lower-order penalty algorithm for solving second-order cone mixed complementarity problems, J. Comput. Appl. Math., 385 (2021), 113168. https://doi.org/10.1016/j.cam.2020.113168 doi: 10.1016/j.cam.2020.113168
![]() |
[38] |
S. Hayashi, N. Yamashita, M. Fukushima, A combined smoothing and regularization method for monotone second-order cone complementarity problems, SIAM J. Optimiz., 15 (2005), 593–615. https://doi.org/10.1137/S1052623403421516 doi: 10.1137/S1052623403421516
![]() |
[39] | S. Hayashi, Manual of ReSNA, 2013. Available from: http://optima.ws.hosei.ac.jp/hayashi/ReSNA/manual/manual_ReSNA.pdf. |
1. | Erdal Bas, Funda Metin Turk, Ramazan Ozarslan, Ahu Ercan, Spectral data of conformable Sturm–Liouville direct problems, 2021, 11, 1664-2368, 10.1007/s13324-020-00428-6 | |
2. | Tom Cuchta, Dallas Freeman, Discrete Polylogarithm Functions, 2023, 84, 1338-9750, 19, 10.2478/tmmp-2023-0012 | |
3. | B. Shiri, Y. Guang, D. Baleanu, Inverse problems for discrete Hermite nabla difference equation, 2025, 33, 2769-0911, 10.1080/27690911.2024.2431000 | |
4. | Muhammad Sulthan Zacky, Heru Sukamto, Lila Yuwana, Agus Purwanto, Eny Latifah, The performance of space-fractional quantum carnot engine, 2025, 100, 0031-8949, 025306, 10.1088/1402-4896/ada9de |
x(t) | μ=0.3 | μ=0.35 | μ=0.4 | μ=0.45 | μ=0.5 |
x(1) | 1 | 1 | 1 | 1 | 1 |
x(2) | 0.612 | 0.714 | 1.123 | 0.918 | 1.020 |
x(3) | 0.700 | 0.900 | 1.515 | 1.370 | 1.641 |
x(5) | 0.881 | 1.336 | 2.402 | 2.747 | 3.773 |
x(7) | 1.009 | 1.740 | 3.352 | 4.566 | 7.031 |
x(9) | 1.099 | 2.100 | 4.332 | 6.749 | 11.461 |
x(12) | 1.190 | 2.570 | 5.745 | 10.623 | 20.450 |
x(15) | 1.249 | 2.975 | 6.739 | 15.149 | 32.472 |
x(16) | 1.264 | 3.098 | 7.235 | 16.793 | 37.198 |
x(18) | 1.289 | 3.330 | 8.233 | 20.279 | 47.789 |
x(20) | 1.309 | 3.544 | 9.229 | 24.021 | 59.967 |
x(t) | q(t)=1 | q(t)=t | q(t)=√t |
x(1) | 1 | 1 | 1 |
x(2) | 7.37∗10−17 | 4.41∗10−17 | 5.77∗10−17 |
x(3) | −0.131 | −0.057 | −0.088 |
x(5) | −0.123 | −0.018 | −0.049 |
x(7) | −0.080 | −0.006 | −0.021 |
x(9) | −0.050 | −0.003 | −0.011 |
x(12) | −0.028 | −0.001 | −0.005 |
x(15) | −0.017 | −0.0008 | −0.003 |
x(16) | −0.015 | −0.0006 | −0.0006 |
x(18) | −0.012 | −0.0005 | −0.002 |
x(20) | −0.010 | −0.0003 | −0.001 |
x(t) | λ=0.1 | λ=0.11 | λ=0.12 |
x(1) | 1 | 1 | 1 |
x(2) | 1 | 1.025 | 1.052 |
x(3) | 1.668 | 1.751 | 1.841 |
x(5) | 3.876 | 4.216 | 4.595 |
x(7) | 7.243 | 8.107 | 9.095 |
x(9) | 11.941 | 13.707 | 12.130 |
x(12) | 22.045 | 26.197 | 25.237 |
x(15) | 36.831 | 45.198 | 46.330 |
x(16) | 43.042 | 53.369 | 55.687 |
x(18) | 57.766 | 73.092 | 78.795 |
x(20) | 76.055 | 98.154 | 127.306 |
x(t) | μ=0.3 | μ=0.35 | μ=0.4 | μ=0.45 | μ=0.5 |
x(1) | 1 | 1 | 1 | 1 | 1 |
x(2) | 0.612 | 0.714 | 1.123 | 0.918 | 1.020 |
x(3) | 0.700 | 0.900 | 1.515 | 1.370 | 1.641 |
x(5) | 0.881 | 1.336 | 2.402 | 2.747 | 3.773 |
x(7) | 1.009 | 1.740 | 3.352 | 4.566 | 7.031 |
x(9) | 1.099 | 2.100 | 4.332 | 6.749 | 11.461 |
x(12) | 1.190 | 2.570 | 5.745 | 10.623 | 20.450 |
x(15) | 1.249 | 2.975 | 6.739 | 15.149 | 32.472 |
x(16) | 1.264 | 3.098 | 7.235 | 16.793 | 37.198 |
x(18) | 1.289 | 3.330 | 8.233 | 20.279 | 47.789 |
x(20) | 1.309 | 3.544 | 9.229 | 24.021 | 59.967 |
x(t) | q(t)=1 | q(t)=t | q(t)=√t |
x(1) | 1 | 1 | 1 |
x(2) | 7.37∗10−17 | 4.41∗10−17 | 5.77∗10−17 |
x(3) | −0.131 | −0.057 | −0.088 |
x(5) | −0.123 | −0.018 | −0.049 |
x(7) | −0.080 | −0.006 | −0.021 |
x(9) | −0.050 | −0.003 | −0.011 |
x(12) | −0.028 | −0.001 | −0.005 |
x(15) | −0.017 | −0.0008 | −0.003 |
x(16) | −0.015 | −0.0006 | −0.0006 |
x(18) | −0.012 | −0.0005 | −0.002 |
x(20) | −0.010 | −0.0003 | −0.001 |
x(t) | λ=0.1 | λ=0.11 | λ=0.12 |
x(1) | 1 | 1 | 1 |
x(2) | 1 | 1.025 | 1.052 |
x(3) | 1.668 | 1.751 | 1.841 |
x(5) | 3.876 | 4.216 | 4.595 |
x(7) | 7.243 | 8.107 | 9.095 |
x(9) | 11.941 | 13.707 | 12.130 |
x(12) | 22.045 | 26.197 | 25.237 |
x(15) | 36.831 | 45.198 | 46.330 |
x(16) | 43.042 | 53.369 | 55.687 |
x(18) | 57.766 | 73.092 | 78.795 |
x(20) | 76.055 | 98.154 | 127.306 |