Tối ưu hóa trong mô hình chống khủng bố bất đối xứng Lanchester (2,1)

Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số Đặc san CNTT, 11 - 2018 133
TỐI ƯU HÓA TRONG MÔ HÌNH CHỐNG KHỦNG BỐ BẤT ĐỐI 
XỨNG LANCHESTER (2,1) 
Nguyễn Hồng Nam*, Hy Đức Mạnh, Vũ Anh Mỹ 
Tóm tắt: Chống khủng bố là một nhiệm vụ toàn cầu mà mọi quốc gia đều quan 
tâm. Để cải thiện hiệu quả của hoạt động chống khủng bố, nhiều quốc gia đã liên 
minh với nhau để cùng nhau lên phương án thực hiện các hoạt động chống khủng 
bố. Bài báo này mở rộng mô hình KKS (do Kaplan, Kress và Szechtmann đưa ra 
[10,11]), với nhiều lực lượng quân chống khủng bố chống lại một nhóm khủng bố. 
Chúng tôi sử dụng nguyên lý cực đại Pontryagin để đưa ra phương án tình báo và 
bổ sung quân số tối ưu cho các lực lượng chống khủng bố. Chúng tôi cũng đưa ra 
một số kết quả tính toán cho thấy chống khủng bố bằng liên minh hiệu quả hơn 
chống khủng bố đơn lẻ. 
Từ khóa: Mô hình Lanchester; Điều khiển tối ưu; Tình báo; Chống khủng bố; Mô hình KKS. 
1. MỞ ĐẦU 
 Mô hình toán học cho một trận đánh lần đầu tiên được Lanchester đưa ra vào 
năm 1916 dưới dạng một hệ phương trình vi phân với hai phương trình hai hàm 
ẩn là quân số của hai bên tham chiến [12]. Quân số của mỗi bên tham chiến được 
giả thiết là đồng nhất (cùng một loại vũ khí). Loại vũ khí của hai bên có thể như 
nhau, nhưng nói chung là bất kỳ. Mô hình không phân biệt loại vũ khí của các bên 
mà chỉ chú ý đến hiệu quả tiêu diệt đối phương của chúng. Thực chất đó là cường 
độ dòng các phát bắn hiệu quả. Mô hình này về sau được gọi là mô hình 
Lanchester hay mô hình Lanchester tổ chức cao. Tính tổ chức cao của mô hình thể 
hiện ở cách xử lý thông tin của tiến trình trận đánh. Mỗi đơn vị chiến đấu của mỗi 
bên đều chỉ tồn tại ở hai trạng thái: còn chiến đấu (chưa bị tiêu diệt) và bị tiêu diệt 
(không chiến đấu được nữa). Đồng thời trên quan điểm “điều khiển trận đánh” thì 
các trạng thái của các đơn vị chiến đấu được các bên nhận biết một cách tức thời 
(một đơn vị bị diệt thì đối phương biết ngay để không bắn vào đó nữa). 
 Mãi đến năm 1962, Deitchman [6] mở rộng mô hình Lanchester bằng cách xét 
trận đánh của một bên là quân chính qui và một bên không chính qui. Tính chính 
qui dựa vào các xử lý thông tin. Trong mô hình mở rộng thông tin không được xử 
lý tức thời mà có độ trễ. Mô hình được gọi là mô hình chiến tranh du kích. Mô 
hình còn được gọi là mô hình Lanchester bất đối xứng. Trong mô hình này, hỏa 
lực của lực lượng du kích là xác định, trong khi hỏa lực của lực lượng quân chính 
quy là không xác định. Sau Deitchmann, Schaffer[14] và Schreiber[15] đã mở 
rộng mô hình của Deitchmann bằng cách đưa vào yếu tố thông tin tình báo và 
xem xét vấn đề phân bố hỏa lực tối ưu của bên quân chính qui. 
 Gần đây, Kaplan, Kress và Szechtman (KKS) [10,11] cũng xem xét mô hình 
Lanchester có thêm thông tin tình báo. Mô hình được ứng dụng cho chống khủng 
bố. Theo mô hình này thì thông tin tình báo ảnh hưởng rất lớn đến kết quả của 
cuộc chiến giữa lực lượng chống khủng bố và lực lượng khủng bố. Kaplan, Kress 
Công nghệ thông tin 
N. H. Nam, H. Đ. Mạnh, V. A. Mỹ, “Tối ưu hóa trong mô hình  Lanchester (2, 1).” 134 
và Szechtman gọi mô hình này là mô hình bất đối xứng (do lực lượng không cân 
bằng). 
 Năm 1974, Taylor [16] nghiên cứu hỏa lực tối ưu theo thời gian cho một số mô 
hình trận đánh. Mackay và các tác giả [13] mở rộng các kết quả của Taylor cho bài 
toán hỏa lực tối ưu theo thời gian cho mô hình Lanchester (n,1). Chen và nhóm tác 
giả [4] lại nghiên cứu bài toán điều khiển tối ưu cho mô hình Lanchester (1,1) có 
bổ sung quân số với tham số điều khiển là tốc độ bổ sung quân số; kết quả sau đó 
được mở rộng cho mô hình Lanchester (2,2) [5]. Đặc điểm chung của các nghiên 
cứu trên là hàm mục tiêu đều là các hàm về quân số. Theo một hướng khác, 
Feichtinger và nhóm tác giả [1, 2, 3, 7] nghiên cứu bài toán điều khiển tối ưu cho 
mô hình KKS và một số mô hình khác với hàm mục tiêu là chi phí cho trận đánh 
chống khủng bố, các biến điều khiển là thông tin tình báo và tốc độ bổ sung quân 
số. 
 Trong bài báo này, chúng tôi mở rộng kết quả nghiên cứu của Feichtinger và 
nhóm tác giả cho mô hình KKS bằng cách đưa ra mô hình chống khủng bố bất đối 
xứng Lanchester (2,1) giữa 2 lực lượng chống khủng bố và một nhóm khủng bố. 
2. MÔ HÌNH 
Giả sử có hai lực lượng tham gia chống khủng bố và tại thời điểm t bất kỳ có 
quân số là 1 0X t và 2 0X t đối đầu với một nhóm khủng bố với quân số là 
 Y t (nằm trong tổng dân số ).P Để đơn giản ta cũng gọi các lực lượng chống 
khủng bố là X1 và X2 và bọn khủng bố là .Y Không mất tổng quát, ta giả thiết 
1,P suy ra 0 1.Y Ký hiệu: 
- 1 2, là hiệu quả tấn công của lực lượng khủng bố Y lên 1X và 2;X và 1, 2 
là hiệu quả tấn công của các lực lượng chống khủng bố 1X và 2X lên ,Y 
1 2 1 2(0 , , , 1).   
-  là mức độ thông tin tình báo của 2 lực lượng chống khủng bố, có nghĩa là 2 
lực lượng chống khủng bố có sự chia sẻ thông tin tình báo giúp xác định vị trí 
chính xác của lực lượng khủng bố, trong khi đó phần 1  là không xác định 
được vị trí của khủng bố. Ở đây 0 1; với 0 có nghĩa là lực lượng 
chống khủng bố không có thông tin gì về quân khủng bố, còn 1 có nghĩa là 
lực lượng chống khủng bố có đầy đủ thông tin gì về quân khủng bố (thực tế thì 
điều này không thể xảy ra (xem [11]) nên ta luôn coi 1). 
- 1 2,   là hệ số tự tiêu hao của 1 2, X X (đào ngũ, bệnh tật ), 1 20 , 1.  
- 1 2,  là tốc độ bổ sung quân số của 2 lực lượng chống khủng bố. 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số Đặc san CNTT, 11 - 2018 135
Do đó, theo kết quả tương tự như Gustav Feichtinger và nhóm tác giả [7], 
chúng ta có mô hình dưới dạng hệ phương trình vi phân như sau: 
1 1 1 1 1
2 2 2 2 2
1 1 2 2( )( (1 ) ) ( )
X Y X
X Y X
Y X X Y C
  
  
    



 (1) 
Trong đó, ( )C là hiệu ứng con dao 2 lưỡi với: 
1 1 2 2( )(1 )(1 ).C X X Y   (2) 
Theo Gustav Feichtinger và nhóm tác giả [11], chúng tôi giả sử chi phí để thu 
thập thông tin tình báo là một hàm lồi nào đó của  thỏa mãn: 
(0) 0, ( ) 0, ( ) 0, (1)= .C C C C  (3) 
Thiệt hại do khủng bố gây ra cũng là một hàm lồi ( )D Y nào đó của biến Y thỏa 
mãn: 
(0) 0, ( ) 0, ( ) 0.D D Y D Y (4) 
Chi phí để duy trì quân đội, ký hiệu 1 2( ), ( )A X A X , có thể giả thiết là hàm tuyến 
tính hoặc lõm của các biến 1 2, X X thỏa mãn: 
1 2 1 1 2 2 1 1 2 2(0) 0, (0) 0, ( ) 0, ( ) 0, ( ) 0, ( ) 0.A A A X A X A X A X (5) 
Chi phí để bổ sung quân số, ký hiệu 1 2( ), ( ),K K  theo các bài toán về quy 
hoạch nguồn nhân lực, chi phí này thường được giả thiết là các hàm bình phương 
của các biến 1 2, .  Do đó, nhiệm vụ của các lực lượng chống khủng bố là giải 
quyết bài toán: 
1 2
1 1 2 2 1 2
, ,
0
min ( ( ) ( ) ( ) ( ) ( ) ( )) .rte D Y C A X A X K K dt
  
  
 (6) 
Trong bài báo này, chúng tôi lựa chọn thiệt hại do khủng bố gây ra là hàm bình 
phương
2
2
fY
 (với 0f ), chi phí để duy trì quân đội là tuyến tính: 1 1 2 2, c X c X (
1 2, 0c c ), chi phí để thu thập thông tin tình báo là hàm logarith: log(1 ) và chi 
phí để bổ sung quân số là 
2 2
1 2,
2 2
 
 tương ứng. Hơn nữa, chúng tôi giả thiết 
2( )C C  . 
3. KẾT QUẢ 
Theo trên, nhiệm vụ của các lực lượng chống khủng bố là giải quyết bài toán: 
1 2, ,
2 22
1
1
0
2
1 2 2 log(1 )
2
min (
2
.
2
)rt
fY
c X c Xe dt
  
 

Công nghệ thông tin 
N. H. Nam, H. Đ. Mạnh, V. A. Mỹ, “Tối ưu hóa trong mô hình  Lanchester (2, 1).” 136 
Sử dụng nguyên lý tối ưu Pontryagin, ta có hàm Hamilton: 
2 22
1 2
1 1 2 2
1 1 1 1 1 2 2 2 2 2
2 2 2
3 1 1 2 2 1 1 2 2
log(1 )
2 2 2
 ( ) ( )
 + ( 1 1 1 ),
fY
H c X c X
Y X Y X
X X Y X X Y
 

     
        
Với 1 2 3, ,    là các biến liên hợp. 
Ta có: 
2 2
3 1 1 2 2 1 1 2 2
0 ,
1
1 2 1 1 .
1
i i i i
i
H
H
X X Y X X Y
   

      
 


 
  
Nếu 3 0 thì hàm Hamilton là đơn điệu giảm. Khi đó thông tin tình báo tối ưu 
là 0. Do đó 3 0. 
2
3
1
0 2 0.
H
x x
 


 (7) 
Phương trình này có nghiệm là 
3
8
1 1
.
4
x



 Vì 1x suy ra 
 3
1 1
.
1 2 1 2x x

 
Vậy 
 1 1 2 2
1 .
1
x
X X Y

 
 (8) 
Các biến liên hợp thỏa mãn các phương trình vi phân: 
.
i i
i
H
r
X
 


, 1;2i và
.
3 3 .
H
r
Y
 


Ta có các phương trình liên hợp: 
.
1 1 1 1
2 2
3 1 1 1 1 2 2 1 2 1 1 ,
r c
Y X X Y
  
       
(9) 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số Đặc san CNTT, 11 - 2018 137
.
2 2 2 2
2 2
3 2 2 1 1 2 2 1 2 1 1 ,
r c
Y X X Y
  
       
(10) 
.
3 3 1 1 2 2
2 2
3 1 1 2 2 1 1 2 2 1 2 1 1 .
r fY
X X X X Y
   
       
(11) 
Chú ý rằng: 
2 2
1 1 2 2
3 1 1 2 2
1
1 2 1 1 1 ,i i iY X X Y
X X
       
  
và 
2 2
1 1 2 2 1 1 2 2
3
1
1 2 1 1 ,
1
X X X X Y
Y
      

và 
 2 2 21 1 2 2 1 1 2 2 1 1 2 21 1 1 1 .X X Y X X Y X X x x           
Cuối cùng chúng ta có hệ bao gồm các phương trình trạng thái và phương trình 
liên hợp: 
.
1 1 1 1 1,X Y X   (12) 
.
2 2 2 2 2 ,X Y X   (13) 
.
1 1 2 2 1 ,Y X X x x   (14) 
.
1
1 1 1 1 1 3
1 1 2 2
,r c
X X

    
 
 (15) 
.
2
2 2 2 2 2 3
1 1 2 2
,r c
X X

    
 
 (16) 
.
3 3 1 1 2 2
1
.
1
r fY
Y
    
 (17) 
3.1. Trạng thái ổn định trong 
.
1 1 2 2 0 1 ,Y X X x x   (18) 
.
1
1
0 ,
1 1 2
i
i i
i
c
r x x

 
  
(19) 
.
2 1 1 1 2 2 2 1 1 2.0iX X X     
(20) 
Ta có hệ tuyến tính: 
Công nghệ thông tin 
N. H. Nam, H. Đ. Mạnh, V. A. Mỹ, “Tối ưu hóa trong mô hình  Lanchester (2, 1).” 138 
 1 1 2 2
2 1 1 1 2 2 2 1 1 2
1
.
,X X x x
X X
  
    
 (21) 
Định thức cơ sở của hệ là khác 0, do đó hệ có nghiệm duy nhất 1 2, .X X Từ đó 
ta tìm được 
 1 1 2 1 1 2 1 2 1 1 21 1 1
2
1 1 2 1 2 1 2 1
1
.
x xX
Y
          
    
Vì 1 20, 0  nên 0.Y Để tồn tại trạng thái ổn định trong cần có 
1 20, 0.c c 
Thay , ,i iX Y vào phương trình (17) ta thu được phương trình 1 chiều 
 0.F x Từ đó ta thu được trạng thái ổn định trong. 
3.2. Tính ổn định 
Ma trận Jacobian cho bởi: 
1 1
2 2
1 2
3
2
1 1 2
1 12 2
1 1 2 2 1 1 2 2
2
1 2 2
2 22 2
1 1 2 2 1 1 2 2
1 22
0 1 0 0
0 0 1 0
0 0 0 1 2
0 0
0 0
1
0 0
1
x
x
rJ
X X X X
r
X X X X
f r
Y
 
 
  

  
 
   
  
 
   
 
  
Với 
2
3 3 3 3
1
.
8
x
   

 
3.3. Trạng thái ổn định biên 
Nếu chỉ có 1 lực lượng chống khủng bố tham gia cuộc chiến, thì ta thu được kết 
quả như trong bài báo [7]. 
4. MỘT SỐ TÍNH TOÁN SỐ 
Trong bài báo [7] các tác giả chọn các tham số:
2.2; 2.23; 0.34; 1.19; 1.86; 2.3734; 1.12r c f    . Kết quả tính 
toán số như sau: trạng thái ổn định trong là 1.39523; 0.16034,X Y tốc độ bổ 
sung quân số tối ưu là 0.83193  , và 0.48822 là mức độ thông tin tình 
báo tối ưu. Để so sánh, chúng tôi sử dụng lại các tham số này cho mô hình của 
chúng tôi. 
Nghiên c
Tạp chí Nghi
Trư
các tham s
Tr
X X Y
M
ra dư
Trư
r c c f
Tr
ờng hợp 1
ạng thái ổn định trong:
1 2 
ức độ thông tin t
Tr
ờng hợp 2: 
2.5; 2.0; 1.9; 0.34; 0.37; 1.3; 2.5; 1.2; 1.86. 
V
ạng thái ổn định trong:
X X Y
  
M
ạng thái ổn định chỉ có thể 
ới H
ới các tham số n
1 2 0.355214; 0.387704; =0.302679;
1 2 3 
ức độ thông tin t
ứu khoa học công nghệ 
ình 1
1 2 1 2 1 2 1 2     
0.726132; 0.718542; 0.698659.
ên c
ố n
0.294785; 0.284383;
ứu KH&CN 
ày, chúng tôi có k
. 
: r
ình báo t
Ch
ọn một bộ tham số t
ày, chúng tôi có k
2.2, 2.23; 0.34; 1.19; 1.86.
ình báo t
quân s
Hình 1
     
Hình 2.
1 2 1 2 1 2
ối 
ự, 
ết quả số nh
ưu:
ối 
Số Đặc san 
. K
 
đạt đ
 K
ưu:
ết quả cho 
ư
ết quả cho tr
 
1 2 3   
0.200107.
ợc bởi hai đa tạp ổn định. Kết quả đ
ùy ý:
ết quả số nh
ư sau:
0.257162.
CNTT
trư
0.7344; 0.998427.
ư
, 11
ờng hợp 1.
ờng hợp 2.
ư sau:
 - 20
18
ược chỉ 
139
 V
ới 
Công nghệ thông tin 
N. H. Nam, H. Đ. Mạnh, V. A. Mỹ, “Tối ưu hóa trong mô hình  Lanchester (2, 1).” 140 
Trạng thái ổn định chỉ có thể đạt được bởi hai đa tạp ổn định. Kết quả được chỉ ra 
dưới Hình 2. 
5. KẾT LUẬN 
Bài báo là kết quả của sự mở rộng mô hình KKS. Các kết quả của Feichtinger 
và nhóm tác giả trong [7] cũng là một trường hợp riêng của bài báo này. Các 
nghiên cứu của bài báo cũng chỉ ra rằng nếu có sựu hợp tác giữa 2 lực lượng chống 
khủng bố thì hiệu quả sẽ tốt hơn nhiều khi từng lực lượng riêng lẻ chống khủng bố. 
Dựa trên mô hình và các kết quả đạt được, chúng tôi đề xuất mở rộng cho mô hình 
tổng quát khi có n ( 1)n lực lượng chống khủng bố liên kết với nhau cùng chống 
lại một lực lượng khủng bố. 
TÀI LỆU THAM KHẢO 
[1]. Andrea Seidl, Edward H. Kaplan, Jonathan P. Caulkins, Stefan Wrzaczek, 
Gustav Feichtinger. Optimal control of a terror queue. European Journal of 
Operational Research, 248, 2016, 246–256. 
[2]. Caulkins, J.P., Grass, D., Feichtinger, G., Tragler, G. Optimizing counter-
terror operations: Should one fight fire with fire or water? Computer and 
Operations Research, 35, 2008, 1874-1885. 
[3]. Caulkins, J.P., Feichtinger, G., Grass, D., Tragler, G. Optimal Control of 
Terrorism and Global Reputation. Operations Research Letters, 37, 2009, 
387-391. 
[4]. Chen, X.Y., Jiang, N., Jing, Y., Stojanovski, G., Dimirovski, G.M. 
Differential Game Model and Its Solutions for Force Resource 
Complementary via Lanchester Square Law Equation. International 
Federation of Automatic Control (IFAC), 18, 2011, 14229-14233. 
[5]. Chen, X.Y., Cao, J., Qiu, J., Jing, Y., Yang, L., Zheng, B. Optimal Control of a 
Class of Warfare Dynamic Systems Based on Lanchester (2,2) Attrition Model. 
27th Chinese Control and Decision Conference (CCDC), 2015, 1263-1267. 
[6]. Deitchman, S.J. A Lanchester model of guerilla warfare. Operations 
Research, 10, 1962, 818-827. 
[7]. [7] Feichtinger, G., Novak, A., Wrzaczek, S. Optimizing Counter-terroristic 
Operations 
[8]. An Asymmetric Lanchester Model. 15th IFAC Workshop on Control 
Applications of Optimization, 2012, 27-32. 
[9]. Grass, D., Caulkins, J.P., Feichtinger, G., Tragler, G., Behrens, D.A. Optimal 
Control of Nonlinear Processes: With Applications in Drugs, Corruption and 
Terror. Springer, Heidelberg, 2008. 
[10]. Grass, D. Numerical computation of the optimal vector field: exemplified by 
a fishery model. Accepted for Journal of Dynamics and Control, 2012. 
Nghiên cứu khoa học công nghệ 
Tạp chí Nghiên cứu KH&CN quân sự, Số Đặc san CNTT, 11 - 2018 141
[11]. Kaplan, E.H., Mintz, A., Mishal, S., Samban, C. What happened to suicide 
bombings in Israel? Insights from a terror stock model. Studies in Conflict 
and Terrorism, 28, 2005, 225-235. 
[12]. Kress, M., Szechtmann, R. Why defeating insurgencies is hard: the effect of 
intelligence in counter insurgency operations - a best case scenario. 
Operations Research, 57 (3), 2009, 578-585. 
[13]. Lanchester, F.W. Aircraft in Warfare: The Dawn of the Fourth Arm. 
Constable, London, 1916. 
[14]. Lin, K.Y, Mackay, L.J. The optimal policy for the one-against-many 
heterogeneous 
[15]. Lanchester model. Operations Research Letters, 42, 2014, 473-477. 
[16]. Schaffer, M.B. Lanchester models of guerrilla engagements. Operations 
Research, 16, 1968, tr.457-488. 
[17]. Schreiber, T.S. Letter to the Editor—Note on the Combat Value of 
Intelligence and Command Control Systems. Operations Research, 12(3), 
1964, 507-510. 
[18]. Taylor, J.G. Lanchester-Type Models of Warfare and Optimal Control. 
Naval Research Logistics Quarterly, 21, 1974, 70-106. 
ABSTRACT 
OPTIMIZING IN AN ASYMMETRIC LANCHESTER (2,1) MODEL 
FOR COUTER-TERRORISM 
Counter-terrorism is a global task that every nation is concerned about. 
To improve operations against terrorism, many nations carry out counter-
terroristic operations not only by themselves but also by cooperation with 
other ones. In this paper, we propose an extended KKS model to cope with 
multi-party counter-terrorism. Optimal control problem for this model is 
studied. Our main tool is Pontryagin's maximal principle. Optimal 
intelligence level and individual reinforcement of each party are found. 
Numerical results show that counter-terroristic operations in a cooperative 
model are more effective than that in single model. 
Key word: Lanchester Model; Optimal Control; Intelligence; Counter-terrorism; KKS Model. 
Nhận bài ngày 04 tháng 7 năm 2018 
Hoàn thiện ngày 09 tháng 10 năm 2018 
Chấp nhận đăng ngày 05 tháng 11 năm 2018 
Địa chỉ: Bộ môn Toán, Khoa CNTT, Học viện KTQS. 
 *Email: nguyenhongnam1977@gmail.com.