Linear Programming Problems On Time Scales
Linear Programming Problems On Time Scales
1. INTRODUCTION
192
Linear Programming Problems on Time Scales 193
The theory of time scales, on the other hand, was first introduced by Stefan
Hilger in 1988 in his PhD dissertation, see [22]. The purpose of this theory is to unify
discrete and continuous analysis and to offer an extension to cases “in between”.
Many applications in mathematical modelling exist for this theory, e.g., to optimal
control [2, 15–18, 25], population biology [6], calculus of variations [5, 7, 10], and
economics [1, 8, 11, 12, 26].
In this work, it is therefore meaningful to consider linear programming prob-
lems on arbitrary time scales T as an extension of the continuous-time linear pro-
gramming problems and to extend the results in [21, 24, 27–29] to arbitrary time
scale. The paper is organized as follows. In Section 2, some examples related to
time scales calculus are given. In Section 3, the basic structures of the primal lin-
ear programming model as well as the dual model on time scales are formulated.
The weak duality theorem and the optimality condition theorem are presented and
proved in Section 4. Section 5 states and proves the strong duality theorem on
isolated time scales. Examples are presented in Section 6 in order to demonstrate
our theoretical results. In Section 7, some conclusions are given.
and Z b Z b
f (t)∆t = f (t)dt, where a, b ∈ T with a < b,
a a
is the usual Riemann integral from calculus.
Example 2.2. If T = {tk ∈ R : k ∈ N0 } with tk < tk+1 for all k ∈ N0 consists only
of isolated points, then
f (tk+1 ) − f (tk )
σ(tk ) = tk+1 , µ(tk ) = tk+1 − tk , f ∆ (tk ) = for k ∈ N0 ,
tk+1 − tk
and
Z tn n−1
X
(2.1) f (t)∆t = µ(tk )f (tk ), where m, n ∈ N0 with m < n.
tm k=m
194 Rasheed Al-Salih and Martin Bohner
f (t + h) − f (t)
σ(t) = t + h, µ(t) ≡ h, f ∆ (t) = for t ∈ T,
h
and
b
Z b h −1
X
f (t)∆t = h f (kh), where a, b ∈ T with a < b.
a a
k= h
and
Z b b−1
X
f (t)∆t = f (k), where a, b ∈ T with a < b.
a k=a
f (qt) − f (t)
σ(t) = qt, µ(t) = (q − 1)t, f ∆ (t) = for t ∈ T,
(q − 1)t
and
Z qn n−1
X
f (t)∆t = (q − 1) q k f (q k ), where m, n ∈ N0 with m < n.
qm k=m
I = [0, T ] ∩ T.
A feasible solution of (P) (or (D)) is any one that satisfies the given constraints. An
optimal solution to (P) (or (D)) is a feasible solution with the largest (or smallest)
objective function value.
In this section, we state and prove the weak duality theorem and the opti-
mality condition theorem. Based on the objective functions and the constraints in
>
(P) and (D), we define the bilinear form h·, ·i1 on Em × Em by
Z σ(T )
(4.2) hz, gi1 = z(t)g(t)∆t,
0
and
(
Minimize hz, gi1
(4.7)
such that A0 z ≥ f and z ≥ 0,
respectively.
196 Rasheed Al-Salih and Martin Bohner
Proof. Using standard properties of the time scales integral (see [13, Section 1.4]),
we calculate
Z σ(T )
(4.2)
hz, Axi1 = z(t)Ax(t)∆t
0
Z σ(T ) Z t
(4.4)
= z(t) B(t)x(t) − K(t, s)x(s)∆s ∆t
0 0
Z σ(T ) Z σ(T ) Z t
= z(t)B(t)x(t)∆t − z(t)K(t, s)x(s)∆s∆t
0 0 0
Z σ(T ) Z σ(T ) Z σ(T )
(∗)
= z(s)B(s)x(s)∆s − z(t)K(t, s)x(s)∆t∆s
0 0 σ(s)
" #
Z σ(T ) Z σ(T )
= z(s)B(s) − z(t)K(t, s)∆t x(s)∆s
0 σ(s)
Z σ(T )
(4.5)
= A0 z(s)x(s)∆s
0
(4.3)
= hA0 z, xi2 ,
where in (∗) we have used the time scales change of integration order formula
presented in [23, Lemma 1]. This shows (4.8) and completes the proof.
Theorem 4.7 (Weak Duality Theorem). If x and z are arbitrary feasible solutions
of (P) and (D), respectively, then U (x) ≤ V (z).
and
Hence,
Z σ(T ) (4.11)
Z σ(T )
(P)
U (x) = f (t)x(t)∆t ≤ A0 z(t)x(t)∆t
0 0
Z σ(T )
(4.3) 0 (4.8) (4.2)
= hA z, xi2 = hz, Axi1 = z(t)Ax(t)∆t
0
Linear Programming Problems on Time Scales 197
(4.10)
Z σ(T )
(D)
≤ z(t)g(t)∆t = V (z),
0
U (x) ≤ V (z ∗ ) = U (x∗ ),
V (z) ≥ U (x∗ ) = V (z ∗ ),
where
0 = t0 < t1 < t2 < . . . < tN = T < tN +1 = σ(tN )
in both (P) and (D).
Theorem 5.9 (Strong Duality Theorem). If (P) has an optimal solution x∗ , then
(D) has an optimal solution z ∗ such that U (x∗ ) = V (z ∗ ).
i.e.,
N
X
Maximize U (x) = µ(tk )f (tk )x(tk )
k=0
subject to x(tk ) ≥ 0 for k = 0, 1, . . . , N
x(t0 ) g(t0 )
(5.12) x(t ) g(t )
1 1
x(t ) g(t )
and A 2 ≤ 2 ,
.. ..
. .
x(tN ) g(tN )
where x : I → Rn and
B(t0 ) 0 ··· 0
−µ(t0 )K(t1 , t0 ) B(t1 ) ··· 0
A = −µ(t0 )K(t2 , t0 ) −µ(t1 )K(t2 , t1 ) ··· 0 .
.. .. .. ..
. . . .
−µ(t0 )K(tN , t0 ) −µ(t1 )K(tN , t1 ) · · · B(tN )
By the standard strong duality theorem of discrete linear programming [20, Theo-
rem 5.1 on page 58] applied to (5.12), there exists an optimal solution y ∗ of
N
X
Minimize W (y) = y(tk )g(tk )
k=0
(5.13) subject to y(tk ) ≥ 0 for k = 0, 1, . . ., N
and y(t0 ) y(t1 ) y(t2 ) · · · y(tN ) A
≥ µ(t0 )f (t0 ) µ(t1 )f (t1 ) µ(t2 )f (t2 ) · · · µ(tN )f (tN ) ,
(5.14) U (x∗ ) = W (y ∗ ).
Now we put
y ∗ (tk )
(5.15) z ∗ (tk ) = for k = 0, 1, . . . , N.
µ(tk )
N
X N
X
(5.16) U (x∗ ) = W (y ∗ ) = y ∗ (tk )g(tk ) = µ(tk )z ∗ (tk )g(tk ) = V (z ∗ ).
k=0 k=0
Linear Programming Problems on Time Scales 199
i.e.,
Z σ(T )
(5.18) z ∗ (t)B(t) ≥ f (t) + z(s)K(s, t)∆s for t ∈ I.
σ(t)
6. EXAMPLES
In this section, three examples are given in order to illustrate our duality
theorems on isolated time scales.
Example 6.10. Let T = Z and I = {0, 1, 2, 3, 4}. Then, we consider the isolated
time scales linear programming primal model
4
Z σ(4)
X
Maximize U (x) = tx(t)∆t = tx(t)
0
t=0
Z t t−1
X
subject to 6x(t) ≤ t + x(s)∆s = t + x(s), t ∈ I
0 s=0
and x(t) ≥ 0, t ∈ I,
200 Rasheed Al-Salih and Martin Bohner
where we have used σ and the integral given in Example 2.4. Using MATLAB
command linprog or LINDO solver, we have
On the other hand, the isolated time scales linear programming dual model is
Z σ(4) X4
Minimize V (z) = tz(t)∆t = tz(t)
0
t=0
Z σ(T ) X4
subject to 6z(t) ≥ t +
z(s)∆s = t + z(s), t ∈ I
σ(t) s=t+1
and z(t) ≥ 0, t ∈ I,
where we have used again Example 2.4. Using MATLAB command linprog or
LINDO solver, we have
confirming U (x∗ ) = V (z ∗ ).
Example 6.11. Let T = 5Z and I = {0, 5, 10, 15, 20}. Then, we consider the isolated
time scales linear programming primal model
Z σ(20) 4
X
Maximize U (x) = tx(t)∆t = 25 kx(5k)
0 k=0
t
Z t 5 −1
X
subject to 6x(t) ≤ t + x(s)∆s = t + 5 x(5k), t ∈ I
0 k=0
and x(t) ≥ 0, t ∈ I,
where we have used σ and the integral given in Example 2.3. Using MATLAB
command linprog or LINDO solver, we have
On the other hand, the isolated time scales linear programming dual model is
Z σ(20) 4
X
Minimize V (z) = tz(t)∆t = 25 kz(5k)
0 k=0
Z σ(T ) X 4
subject to 6z(t) ≥ t + z(s)∆s = t + 5 z(5k), t ∈ I
σ(t) t
k= 5 +1
and z(t) ≥ 0, t ∈ I.
Linear Programming Problems on Time Scales 201
where we have used again Example 2.3. Using MATLAB command linprog or
LINDO solver, we have
confirming U (x∗ ) = V (z ∗ ).
Example 6.12. Let T = 2N0 and I = {1, 2, 4}. Then, we consider the isolated time
scales linear programming primal model
Z σ(22 )
X2
Maximize U (x) = tx(t)∆t = 4k x(2k )
1 k=0
Z t log2 t−1
X
subject to 8x(t) ≤ 6t + 3 x(s)∆s = 6t + 3 2k x(2k ), t ∈ I
1 k=0
and x(t) ≥ 0, t ∈ I,
where we have used σ and the integral given in Example 2.5. Using MATLAB
command linprog or LINDO solver, we have
On the other hand, the isolated time scales linear programming dual model is
Z σ(22 ) X2
Minimize V (z) = 6 tz(t)∆t = 6 4k z(2k )
1 k=0
Z σ(4) X 2
subject to 8z(t) ≥ t + 3 z(s)∆s = t + 3 2k z(2k ), t ∈ I
σ(t)
k=1+log2 t
and z(t) ≥ 0, t ∈ I,
where we have used again Example 2.5. Using MATLAB command linprog or
LINDO solver, we have
confirming U (x∗ ) = V (z ∗ ).
7. CONCLUSIONS
the dual linear models on time scales has been established. Furthermore, the weak
duality theorem and the optimality condition theorem were given for arbitrary time
scales, while the strong duality theorem has been established for arbitrary isolated
time scales. This ensures that our formulation of primal and dual models on time
scales presented in this paper is indeed the “correct” formulation. Since all efficient
algorithms in the literature can only solve continuous-time linear programming
problems by a discrete approximation method (i.e., the optimal solution obtained
by this method is a numerical solution with an estimate of the error bound), another
contribution of this paper is to solve linear programming problems on arbitrary
isolated time scales by efficient exact methods (i.e., the optimal solution obtained
by these methods is exact).
REFERENCES
1. Ferhan M. Atici, Daniel C. Biles, and Alex Lebedinsky. An application of time scales
to economics. Math. Comput. Modelling, 43(7-8):718–726, 2006.
2. Zbigniew Bartosiewicz, Ülle Kotta, Ewa Pawluszewicz, and Malgorzata Wyrwas. Con-
trol systems on regular time scales and their differential rings. Math. Control Signals
Systems, 22(3):185–201, 2011.
3. Richard E. Bellman. Dynamic programming. Princeton Landmarks in Mathematics.
Princeton University Press, Princeton, NJ, 2010. Reprint of the 1957 edition, With a
new introduction by Stuart Dreyfus.
4. Richard E. Bellman and Stuart E. Dreyfus. Applied dynamic programming. Princeton
University Press, Princeton, N.J., 1962.
5. Martin Bohner. Calculus of variations on time scales. Dynam. Systems Appl., 13(3-
4):339–349, 2004.
6. Martin Bohner, Meng Fan, and Jimin Zhang. Periodicity of scalar dynamic equations
and applications to population models. J. Math. Anal. Appl., 330(1):1–9, 2007.
7. Martin Bohner, Rui A. C. Ferreira, and Delfim F. M. Torres. Integral inequalities and
their applications to the calculus of variations on time scales. Math. Inequal. Appl.,
13(3):511–522, 2010.
8. Martin Bohner, Gregory Gelles, and Julius Heim. Multiplier-accelerator models on
time scales. Int. J. Stat. Econ., 4(S10):1–12, 2010.
9. Martin Bohner and Svetlin G. Georgiev. Multivariable dynamic calculus on time
scales. Springer, Cham, 2016.
10. Martin Bohner and Gusein Sh. Guseinov. Double integral calculus of variations on
time scales. Comput. Math. Appl., 54(1):45–57, 2007.
11. Martin Bohner, Julius Heim, and Ailian Liu. Solow models on time scales. Cubo,
15(1):13–31, 2013.
12. Martin Bohner, Julius Heim, and Ailian Liu. Qualitative analysis of a Solow model
on time scales. J. Concr. Appl. Math., 13(3-4):183–197, 2015.
13. Martin Bohner and Allan Peterson. Dynamic equations on time scales. Birkhäuser
Boston, Inc., Boston, MA, 2001. An introduction with applications.
Linear Programming Problems on Time Scales 203
14. Martin Bohner and Allan Peterson. Advances in dynamic equations on time scales.
Birkhäuser Boston, Inc., Boston, MA, 2003.
15. Martin Bohner and Nick Wintz. The linear quadratic regulator on time scales. Int.
J. Difference Equ., 5(2):149–174, 2010.
16. Martin Bohner and Nick Wintz. The linear quadratic tracker on time scales. Int. J.
Dyn. Syst. Differ. Equ., 3(4):423–447, 2011.
17. Martin Bohner and Nick Wintz. Controllability and observability of time-invariant
linear dynamic systems. Math. Bohem., 137(2):149–163, 2012.
18. Martin Bohner and Nick Wintz. The Kalman filter for linear systems on time scales.
J. Math. Anal. Appl., 406(2):419–436, 2013.
19. R. N. Buie and Jaromı́r Abrham. Numerical solutions to continuous linear program-
ming problems. Z. Operations Res. Ser. A-B, 17(3):A107–A117, 1973.
20. Vašek Chvátal. Linear programming. A Series of Books in the Mathematical Sciences.
W. H. Freeman and Company, New York, 1983.
21. Richard C. Grinold. Continuous programming. I. Linear objectives. J. Math. Anal.
Appl., 28:32–51, 1969.
22. Stefan Hilger. Analysis on measure chains—a unified approach to continuous and
discrete calculus. Results Math., 18(1-2):18–56, 1990.
23. Başak Karpuz. Unbounded oscillation of higher-order nonlinear delay dynamic equa-
tions of neutral type with oscillating coefficients. Electron. J. Qual. Theory Differ.
Equ., pages No. 34, 14, 2009.
24. Norman Levinson. A class of continuous linear programming problems. J. Math.
Anal. Appl., 16:73–83, 1966.
25. Ewa Pawluszewicz. Observability of nonlinear control systems on time scales. Internat.
J. Systems Sci., 43(12):2268–2274, 2012.
26. Christopher C. Tisdell and Atiya Zaidi. Basic qualitative and quantitative results
for solutions to nonlinear, dynamic equations on time scales with an application to
economic modelling. Nonlinear Anal., 68(11):3504–3524, 2008.
27. William F. Tyndall. A duality theorem for a class of continuous linear programming
problems. J. Soc. Indust. Appl. Math., 13:644–666, 1965.
28. William F. Tyndall. An extended duality theorem for continuous linear programming
problems. SIAM J. Appl. Math., 15:1294–1298, 1967.
29. Ching-Feng Wen, Yung-Yih Lur, and Hang-Chin Lai. Approximate solutions and error
bounds for a class of continuous-time linear programming problems. Optimization,
61(2):163–185, 2012.
204 Rasheed Al-Salih and Martin Bohner
Martin Bohner
Missouri S&T,
Department of Mathematics and Statistics,
Rolla, MO 65409-0020,
USA.
E-mail: bohner@mst.edu