Accès libre

Using Linear Programming with Interval Coefficients in Optimization of a Military Action

  
24 juin 2025
À propos de cet article

Citez
Télécharger la couverture

Introduction

In military field is more realistic to use models of linear programming with interval coefficients method (LPIC) when we want to find the optimum solution.

In paper 1 we present the LPIC problem:

“Suppose we are given the following LPIC problem: minZ=j=1n[cj_,cj¯]xjj=1n[aij_,aij¯]xj[bi_,bi¯],i=1,,mxj0j \matrix{{\min \,Z = \sum\limits_{j = 1}^n {[\underline {{c_j}} ,\overline {{c_j}} ]{x_j}}} \hfill \cr {\sum\limits_{j = 1}^n {[\underline {{a_{ij}}} ,\overline {{a_{ij}}} ]{x_j} \ge [\underline {{b_i}} ,\overline {{b_i}} ],\,\,\,i = 1, \ldots ,m}} \hfill \cr {{x_j} \ge 0\,\forall j} \hfill \cr}

The best optimum and worst optimum are computed as follows:

For the best optimum: minz_=j=1ncjxjj=1naijxjbi_,ixj0j \matrix{{\min \,\underline z = \sum\limits_{j = 1}^n {c_j^{'}{x_j}}} \hfill \cr {\sum\limits_{j = 1}^n {a_{ij}^{'}{x_j} \ge \underline {{b_i}}} ,\,\,\,\forall i} \hfill \cr {{x_j} \ge 0\,\forall j} \hfill \cr} where cj=cj_ c_j^{'} = \underline {{c_j}} and aij=aij¯ a_{ij}^{'} = \overline {{a_{ij}}}

For the worst optimum: minz¯=j=1ncjxjj=1naijxjbi¯,ixj0j \matrix{{\min \,\bar z = \sum\limits_{j = 1}^n {c_j^{''}{x_j}}} \hfill \cr {\sum\limits_{j = 1}^n {a_{ij}^{''}{x_j} \ge \overline {{b_i}} ,\,\,\forall i}} \hfill \cr {{x_j} \ge 0\,\forall j} \hfill \cr} where cj=c¯j c_j^{''} = {\bar c_j} and aij=aij_ a{_{ij}^{''}} = \underline {{a_{ij}}} .

So the best optimum solution will be , at (x1, x2, … , xn) and with coefficient settings cj, aij and i, ∀i and j and also the worst optimum solution with (x1, x2, … , xn) and with coefficient settings cj, aiji și j.

The general solution of LPIC will be Z = [, ] i.e. optimal value of the problem will be between (the best case scenario) and (the worst case scenario).“

Optimal Distribution of Firepower on Objectives

In the military field, there are a multitude of situations that can be “optimized” using various mathematical models.

Problems of distribution of forces and means on objectives refers to: optimal implementation of defense against tanks, optimal implementation of defense against air attacks, optimal choice of combinations of means of combat against different objectives, distribution of resources on different types of weapons and others.

Inspired by Căruțașu (2014) who presented and found an optimal solution to a problem of distribution of firepower on objectives, we create a similar problem, but we used numerical intervals instead of exact real numbers:

The own unit has 5 combat means that use ammunition of the M1, M2, M3, M4 and M5 types, and the enemy has 7 types of objectives O1, O2, O3, O4, O5, O6 and O7 (meaning both equipment and personnel).

Based on the probabilities of hitting each of the 7 enemy targets, the average ammunition consumption of each type necessary to destroy them was established, because of the uncertainty the values are written as intervals (the data is given in the table below).

Each destroyed objective causes the enemy certain percentage losses given in the table below, also these values are intervals because we can't be 100% sure.

Also in the table are given the quantities of ammunition of each type that the own unit has at its disposal as well as the number of objectives of each type that the enemy has.

The problem is to determine the number of enemy objectives of each type that must be destroyed so that the losses caused to the enemy are maximum, taking into account the quantities of ammunition of each type that the own unit has at its disposal.

Problem data

O1 O2 O3 O4 O5 O6 O7
M1 [10,12] [8,10] [6,8] [7,9] [9,11] [7,9] [8,10] [690,720]
M1 [3,5] [2,4] [4,6] [3,5] [4,6] [5,7] [5,7] [360,380]
M3 [6,8] [3,6] [5,7] [8,10] [7,9] [4,6] [7,9] [520,540]
M4 [1,3] [1,2] [3,5] [2,4] [1,2] [2,4] [4,6] [200,220]
M5 [0,2] [0,1] [1,2] [1,3] [0,1] [1,2] [3,5] [80,100]
Loss [0.5%, 1%] [0.25%, 0.5%] [1.25%, 1.5%] [1.75%, 2.5%] [0.5%, 0.75%] [0.75%, 1.25%] [2.5%, 3%] -
No.objectives 12 15 8 7 13 6 7 -

Let xi be number of Oi type objectives selected for destruction, i=1,7¯ i = \overline {1,7} .

The function to be maximized has the form:

Max Z=[0.5,1]x1 + [0.25,0.5]x2 + [1.25,1.5]x3 + [1.75,2.5]x4 + [0.5,0.75]x5 + [0.75,1.25]x6 + [2.5,3]x7

For the constrains:

no more ammunition can be used than the one the unit has at its disposal

[10,12]x1+[8,10]x2+[6,8]x3+[7,9]x4+ [9,11]x5+[7,9]x6+[8,10]x7 ≤ [690,720]

[3,5]x1+[2,4]x2+[4,6]x3+[3,5]x4+[4,6]x5+[5,7]x6+[5,7]x7 ≤ [360,380]

[6,8]x1+[3,6]x2+[5,7]x3+[8,10]x4+ [7,9]x5+[4,6]x6+[7,9]x7 ≤ [520,540]

[1,3]x1+[1,2]x2+[3,5]x3+[2,4]x4+[1,2]x5+[2,4]x6+[4,6]x7 ≤ [200,220]

[0,2]x1+[0,1]x2+[1,2]x3+[1,3]x4+[0,1]x5+[1,2]x6+[3,5]x7 ≤ [80,100]

no more objectives of each type can be destroyed than those the enemy has:

x1 ≤ 12

x2 ≤ 15

x3 ≤ 8

x4 ≤ 7

x5 ≤ 13

x6 ≤ 6

x7 ≤ 7

Sign restrictions

-are the sign conditions imposed on decision variables:

x1 ≥ 0

x2 ≥ 0

x3 ≥ 0

x4 ≥ 0

x5 ≥ 0

x6 ≥ 0

x7 ≥ 0

We will convert the problem in to a minimization problem and change all “≤” inequalities into “≥” inequalities.

The new model is:

Min(−Z)=[−1,−0.5]x1+[−0.5,−0.25]x2+ [−1.5,−1.25]x3+[−2.5,−1.75]x4+ [−0.75,−0.5]x5+[−1.25,−0.75]x6+[−3,−2.5]x7

[−12,−10]x1+[−10,−8]x2+[−8,−6]x3+ [−9,−7]x4+[−11,−9]x5+[−9,−7]x6+ [−10,−8]x7 ≥ [−720, −690]

[−5,−3]x1+[−4,−2]x2+[−6,−4]x3+[−5,−3]x4+ [−6,−4]x5+[−7,−5]x6+[−7,−5]x7 ≥ [−380, 360]

[−8,−6]x1+[−6,−3]x2+[−7,−5]x3+ [−10,−8]x4+[−9,−7]x5+[−6,−4]x6+ [−9,−7]x7 ≥ [−540, −520]

[−3,−1]x1+[−2,−1]x2+[−5,−3]x3+ [−4,−2]x4+[−2,−1]x5+[−4,−2]x6+[−6,−4]x7 ≥ [−220, −200]

[−2,0]x1+[−1,0]x2+[−2,−1]x3+[−3,−1]x4+ [−1,0]x5+[−2,−1]x6+[−5,−3]x7 ≥ [−100, −80]

x1 ≤ 12

x2 ≤ 15

x3 ≤ 8

x4 ≤ 7

x5 ≤ 13

x6 ≤ 6

x7 ≤ 7

x1 ≥ 0

x2 ≥ 0

x3 ≥ 0

x4 ≥ 0

x5 ≥ 0

x6 ≥ 0

x7 ≥ 0

We obtain:

The best optimum:

min − =−x1−0.5x2−1.5x3−2.5x4−0.75x5− 1.25x6−3x7

−10x1−8x2−6x3−7x4−9x5−7x6−8x7 ≥ −720

−3x1−2x2−4x3−3x4−4x5−5x6−5x7 ≥ −380

−6x1−3x2−5x3−8x4−7x5−4x6−7x7 ≥ −540

x1x2−3x3−2x4x5−2x6−4x7 ≥ −220

x3x4x6−3x7 ≥ −100

x1 ≤ 12

x2 ≤ 15

x3 ≤ 8

x4 ≤ 7

x5 ≤ 13

x6 ≤ 6

x7 ≤ 7

x1 ≥ 0

x2 ≥ 0

x3 ≥ 0

x4 ≥ 0

x5 ≥ 0

x6 ≥ 0

x7 ≥ 0

The worst optimum:

min − =−0.5x1−0.25x2−1.25x3−1.75x4−0.5x5−0.75x6−2.5x7

−12x1−10x2−8x3−9x4−11x5−9x6−10x7 ≥ −690

−5x1−4x2−6x3−5x4−6x5−7x6−7x7 ≥ −360

−8x1−6x2−7x3−10x4−9x5−6x6−9x7 ≥ −520

−3x1−2x2−5x3−4x4−2x5−4x6−6x7 ≥ −200

−2x1x2−2x3−3x4−1x5−2x6−5x7 ≥ −80

x1 ≤ 12

x2 ≤ 15

x3 ≤ 8

x4 ≤ 7

x5 ≤ 13

x6 ≤ 6

x7 ≤ 7

x1 ≥ 0

x2 ≥ 0

x3 ≥ 0

x4 ≥ 0

x5 ≥ 0

x6 ≥ 0

x7 ≥ 0

We use LINDO (a linear programming solver) to find the both optimums:

Figure no. 1:

LINDO best optimum

(Source: Author's compilation)

Figure no. 2:

LINDO solution for best optimum

(Source: Author's compilation)

Figure no. 3:

LINDO worst optimum

(Source: Author's compilation)

Figure no. 4:

LINDO solution for worst optimum

(Source: Author's compilation)

So the best optimum is −=−87.25 (Figure no. 2) for the minimal problem, so for our maximum problem =87.25, and also in Figure no. 2 we have the variables. Similar the worst optimum is =43.75 for the maximum problem with the fixed variables used to get this solution.

So we can conclude that optimum solution Z of this problem lies between 43.75 and 87.25, so the enemy's losses will be between 43.75% and 87.25%.

Conclusions

In recent decades, mathematical modeling of combat operations has developed as an indispensable tool in decision-making at all levels. These models provide the possibility of analyzing different variants of the use of forces and combat means, as well as their dynamics during the conduct of military operations.

Many operational research optimization models, such as linear programming, assume that the data used to build a model is accurate, but in the real world most of these assumptions are only approximately true. So it is more realistic to use models of linear programming with interval coefficients (LPIC).