数学建模社区-数学中国

标题: 博弈论 / 对策论 [打印本页]

作者: 浅夏110    时间: 2020-6-11 15:53
标题: 博弈论 / 对策论
1 引言
/ }- h: U2 x6 i社会及经济的发展带来了人与人之间或团体之间的竞争及矛盾,应用科学的方法来 解决这样的问题开始于 17 世纪的科学家,如 C.,Huygens 和 W.,Leibnitz 等。现代对策论起源于 1944 年 J.,Von Neumann 和 O.,Morgenstern 的著作《Theory of Games and Economic Behavior》。/ Y2 @" A, e: h  w8 Z! U

7 `- @8 ^# g( U& ?7 I) ?8 g0 f对策论亦称竞赛论或博弈论。是研究具有斗争或竞争性质现象的数学理论和方法。 一般认为,它既是现代数学的一个新分支,也是运筹学中的一个重要学科。对策论发展 的历史并不长,但由于它所研究的现象与人们的政治、经济、军事活动乃至一般的日常 生活等有着密切的联系,并且处理问题的方法又有明显特色。所以日益引起广泛的注意。
, }2 ~! a8 F* L4 T
) e9 m! q8 a" J7 \, o" t5 N( u深度学习的生成对抗网络的目标函数就是这个原理:二人零和博弈思想,用极大极小原理来判断某个对策是否有鞍点。
' W8 x; Y# ^, W3 H
, h* N+ B* v; J7 N在日常生活中,经常看到一些具有相互之间斗争或竞争性质的行为。具有竞争或对抗性质的行为称为对策行为。在这类行为中。参加斗争或竞争的各方各自具有不同的目标和利益。为了达到各自的目标和利益,各方必须考虑对手的各种可能的行动方案,并力图选取对自己最为有利或最为合理的方案。对策论就是研究对策行为中斗争各方是否 存在着最合理的行动方案,以及如何找到这个合理的行动方案的数学理论和方法。1 _4 x4 ~3 P& ~$ _
) E& U- h& N3 t/ p
2 对策问题
$ h9 D1 Q0 `! A8 D0 h: [5 _对策问题的特征是参与者为利益相互冲突的各方,其结局不取决于其中任意一方的努力而是各方所采取的策略的综合结果。 先考察一个实际例子。 . p1 R3 {9 ^& V

1 q2 E  l2 m- k6 R' e" m+ ^例题1 囚徒的困境
+ ]9 B$ l: c' m& Y 警察同时逮捕了两人并分开关押,逮捕的原因是他们持有大量伪币,警方怀疑他们伪造钱币,但没有找到充分证据,希望他们能自己供认,这两个 人都知道:如果他们双方都不供认,将被以持有大量伪币罪被各判刑 18 个月;如果双 方都供认伪造了钱币,将各被判刑 3 年;如果一方供认另一方不供认,则供认方将被从 宽处理而免刑,但另一方面将被判刑 7 年。将嫌疑犯 A、B 被判刑的几种可能情况列于表 1。
. g$ H: S) w& ^, @( W0 T: O: L
! Q3 B* C) I. O8 E- n5 I
) a* C1 [7 }2 h$ v' t$ d
. Q# ?5 h1 }: j0 q表 1 中每对数字表示嫌疑犯 A、B 被判刑的年数。如果两名疑犯均担心对方供认并希 望受到最轻的惩罚,最保险的办法自然是承认制造了伪币。     从这一简单实例中可以看出对策现象中包含有的几个基本要素。# \9 y$ f% n1 p  t5 j& K# a

% \* \( j" g6 A' I) z8 m2.1  对策的基本要素 9 E; D& z. R6 p8 F" S
(i)局中人 5 q) n, B1 n1 d; b+ `2 n
/ R- o6 @+ a. o# s

, L# T& Y7 P2 Q5 m, A" ]
  c& E" z$ T1 o- V& j' Y$ j' F(ii)策略集
9 o* B" U$ S5 h, J2 ]7 A& \; R' Y; _$ r5 E( q. m

9 D0 ~$ N2 P8 c" \, ?' x2 Z4 u2 }" [. i2 T6 y6 e5 r! f' k
(iii)赢得函数(支付函数)
2 \8 m. I6 V, W8 ?/ a& c. Y% Z4 }. Q3 Z9 m7 w9 @5 e3 e+ o
: V* Q/ @! K& e& R  a: m

) {. t; ?* u+ S1 ~+ p6 Y' |本节我们只讨论有两名局中人的对策问题,其结果可以推广到一般的对策模型中 去。 6 U4 U- h& X& h5 O* c7 X
1 f5 X: R/ l- X2 W9 Y
2.2 零和对策(矩阵对策)
4 d  Q9 |1 K) q0 B: ?零和对策是一类特殊的对策问题。在这类对策中,只有两名局中人,每个局中人都 只有有限个策略可供选择。在任一纯局势下,两个局中人的赢得之和总是等于零,即双 方的利益是激烈对抗的。 4 g5 [3 o2 D5 c7 t; s4 n$ j0 E& F
4 l) \% F# C4 b$ d
& `$ I/ m  f! i

; r$ H2 y" ?- w& S0 \例题2  ; K9 ]! R4 Z; G) V9 ?; ?

8 ~( m4 b& o5 F: t, L2 Z
3 ~9 ^+ C) O0 B6 L" ?2 b7 _# ^8 s

0 ?, B3 [0 P7 U( W" ?. q8 \% g2 @
5 p( i. {  i( @
' P8 v' O6 A7 |& E& v什么是鞍点
4 g6 l: s( j0 W8 a, x, a1 l1 V. V' E
) s* @5 o+ x5 z5 ^3 _
4 a& E+ A9 ^0 t1 [
最优纯策略
/ P; Q) d* F" a5 M8 j
+ c! V) x# M! D4 X- K, z
2 ?: ~! H- D+ s  F# w" r, i( }- D4 [5 Y; S- K  e; j
极大极小原理-----判断某个对策是否有鞍点
( ?& C4 ?- J0 y' {" b给定一个对策G ,如何判断它是否具有鞍点呢?为了回答这一问题,先引入下面 的极大极小原理。 1 p( i) V# Z: e/ G4 t( G# ]2 p

! @! q% A9 Z; t( N6 Z8 J( R
3 ~' S0 `2 r1 w" B/ Q
/ C8 ~" Q. F% z+ S9 O, Z
6 f& _1 ?3 R6 g+ O) g3 q5 @
* d: s+ m1 Q( o9 x  v' l
+ o) x$ b+ S4 Y7 Z# s    上述定理给出了对策问题有稳定解(简称为解)的充要条件。
) J: i; H. c# I+ O) Q+ M. M* T1 t& p" i& j' c
对策问题的多个解之间的关系(两条性质). \" m* t  Q& M) r, c9 M2 |! i
  当对策问题有解时, 其解可以不唯一,当解不唯一时,解之间的关系具有下面两条性质:   o( ]$ s1 Q% i6 @! z( ~6 F

% ?+ s5 v! y; d4 h" d  h  n
; @1 [3 @4 W+ s$ l
: M6 j  B5 b0 r+ F7 @6 Y3 零和对策的混合策略 ( q8 h7 N  R4 n. `& Q; }
具有稳定解的零和问题是一类特别简单的对策问题,它所对应的赢得矩阵存在鞍 点,任一局中人都不可能通过自己单方面的努力来改进结果。然而,在实际遇到的零和对策中更典型的是 μ + ν ≠ 0 的情况。由于赢得矩阵中不存在鞍点,此时在只使用纯策略的范围内,对策问题无解。) e1 L4 {! j7 Y

4 U7 t- W  k3 {7 ?4 [3.1 零和对策的混合策略6 k& r5 W. d, J/ m& Q' h  R' k

  M& ~9 e) M0 _5 |, L# v/ @  y+ z8 r) R! D5 d

3 ^0 p; ]" v/ `; ^0 P) z: Q3.2 混合策略对策问题的鞍点
  N% v3 T. t+ U, H$ n0 e" J4 A& D* I! V6 M6 q- I3 |/ G
; K" c) t/ ~2 D: Y
) i4 q9 B2 l. K' W
( \" a' r% K0 H, ?
+ j/ x/ j& Y0 l0 r7 n6 F: B
9 [# y8 {3 ?4 \/ t6 C4 q, T* B
0 R! U3 U, t8 Q
" a. X! N% U; @" g7 a
使用纯策略的对策问题(具有稳定解的对策问题)可以看成使用混合策略的对策问 题的特殊情况,相当于以概率 1 选取其中某一策略,以概率 0 选取其余策略。
5 u; d+ o) U$ c# N! w9 Y* X+ N: u! n3 Q( I6 @' b
例题3
/ `# d5 X( b9 x& ^$ O0 t4 A  ?( R* X3 A& r) G
5 c$ m9 e5 S9 {& }4 i9 V% _

1 [0 M8 K2 H0 v  b; k  m5 `6 @  a- J; A$ `7 K" i
, B9 z* L1 d+ h$ u+ Z4 l# M, h- g! J

5 ]. Z* \' k; O0 U
! @4 H! `/ M3 `+ u+ F' m2 `, J* Y  v) E
: d+ ?# }2 u& R7 j! c) B3 V

& n" J  D4 s) K1 G2 W5 y6 s. w3.3 关于对策解集性质的主要结果的三个定理% b- W7 v, K# R3 s1 `
# b* M: j# I: _8 G: Q: Q) V

8 h; a5 h) M  i1 z9 }6 f/ |
, e6 h/ \' }& w; U4 零和对策的线性规划解法8 Y) X+ u% O3 T% \) R, j/ ^

0 B3 b$ m0 }4 G0 F+ }' r+ h( }
# Z3 |: X% F! x& G& x4 C, y
$ w1 Q9 Z) `. H7 O' K5 C/ ]: r0 ]: n, L

" R' W' a/ t4 F& E  \' Q# y) O# u& v; H1 T

1 Q2 R* y0 ]5 K. H
! C6 G  z& `/ b& u- X3 p" S例题4
# n- W' V6 P( H4 G3 d( l" k' r8 L5 ~/ j+ n7 ~- {' [
/ p+ d- }8 o2 l
6 Q# z  v* h7 d: h) p
编写程序如下:
# W. o  ~, L1 ~- o4 o% Z
! ^$ A& D9 B/ M( ?0 `/ @clear
/ ^7 a4 t9 i7 u0 Y' `a=[1/3,1/2,-1/3;-2/5,1/5,-1/2;1/2,-3/5,1/3];b=10;
# V2 T& o( l% x% ]4 _! K3 Pa=a+b*ones(3);   %把赢得矩阵的每个元素变成大于0的数 $ O) p. O4 q" {* z8 K
[x0,u]=linprog(ones(3,1),-a',-ones(3,1),[],[],zeros(3,1)); 7 H0 ~& x1 F' v: ]$ n2 S
x=x0/u,u=1/u-b ) s& o8 y# J; B8 J; s; N  \1 C
[y0,v]=linprog(-ones(3,1),a,ones(3,1),[],[],zeros(3,1)); ) S) \, P; {7 ~* n* F0 C
y=y0/(-v),v=1/(-v)-b
1 T! Q1 W7 E+ B6 C+ s+ ?$ B# p
( A$ T; }) S5 t- p2 s4 d; E  t7 m* ~
/ K! B$ a- t; q. j8 {: X* K. R* Q- g- |6 j

0 l) g7 T$ c" ~! v下面我们使用式(2)和(3),利用 LINGO 编程求例 4 的解。LINGO 程序如下:
& F* G6 E( m' b2 N& Mmodel: 0 u4 T6 z& V0 e& _
sets: ; \9 c; c1 v) X; v& w
player1/1..3/:x;
+ x: w! W( e$ Y7 U8 D1 u# dplayer2/1..3/:y;
  q1 J0 B2 ^+ ^0 P0 U' xgame(player1,player2):c; + d/ z. |. x( P: r4 Q9 h) H( D
endsets
. P% p4 \# z, S5 G9 f) zdata:
* a* I4 P* G2 L( Actrl=?; !ctrl取1求局中人1的策略,ctrl取0求局中人2的策略;
) k6 v$ o9 i8 B8 b9 a' yc=0.3333333 0.5 -0.3333333
9 y6 [* O. Y5 ~  D* Y# _3 E. v-0.4 0.2 -0.5 ( Y# t, Q# [/ C. g4 p, _# C4 t
0.5 -0.6 0.3333333;
7 ?8 o* Y7 F3 M, Genddata
) D' L2 u- s( T* J; V/ X% h6 Nmax=u*ctrl-v*(1-ctrl);
: f& N. `& w& V9 P@free(u);@free(v);
" K' u: u5 R' f@for(player2(j)sum(player1(i):c(i,j)*x(i))>u); ; W/ m8 R8 I& r7 y
@for(player1(i)sum(player2(j):c(i,j)*y(j))<v);
& J6 p  G' c4 V1 L+ s4 P% o@sum(player1:x)=1;
1 S* D( o  f" o( L@sum(player2:y)=1;
, I6 E4 }) @# Eend 4 l: o9 s; D! \- B: l4 n7 a6 r/ Z
1 b1 n2 j- H9 q! X2 ~! N5 q: z' I
由定理4知,混合对策问题的求解问题可以转化为求不等式约束的可行点,而 LINGO软件很容易做到这一点。我们编写如下Lingo程序求解上述问题。 , w, }8 n9 O, _3 ^+ P$ m, P

3 e1 o2 J5 R6 umodel: 8 G# v( s9 n! T9 Y9 a5 |
sets:
1 u# F' N7 C: E4 tplayer1/1..3/:x;
  |; r! p0 v5 Splayer2/1..3/:y;
9 y. [% B' V7 a6 v/ Ugame(player1,player2):c; ) @2 a* l! R  |5 L7 M8 C% I
endsets 9 p) a8 g; j/ P. W
data: 0 g& j) {7 t2 j- W
c=0.3333333 0.5 -0.3333333
, |/ Y* @. |9 h; N: T  -0.4 0.2 -0.5
7 x9 G- c) J1 S; \+ l% z  0.5 -0.6 0.3333333; ; G, s6 \* \& F+ n5 K. x  U
enddata
2 @( _  `  J9 q9 d! x+ W/ L3 N: a5 v& w@free(u); & J: o2 ~" X) f& D  i4 i, l% A+ Z
u=@sum(game(i,j):c(i,j)*x(i)*y(j));
) S: D& I" a+ I+ I  @8 N6 g@for(player1(i)sum(player2(j):c(i,j)*y(j))<u); @for(player2(j)sum(player1(i):c(i,j)*x(i))>u); 7 f0 t2 k" q- N; L. q) [" S
@sum(player1:x)=1; 3 U0 ?' r9 d+ L; `" }
@sum(player2:y)=1;
4 `! K* h/ W; a( Fend
( e2 n$ _! K& v* V0 x7 x/ P
6 f5 t! ~8 q9 U' S& \- p 5 二人非常数和对策 + M  ~  I# r7 `8 V1 k$ R7 v
5.1 常数和对策
1 R# |  D6 `4 }) X5 Z所谓常数和对策是指局中人I和局中人II所赢得的值之和为一个常数。显然,二人零和对策是二人常数和对策的特例,即常数为零。 对于二人常数和对策,有纯策略对策和混合策略对策,其求解方法与二人零和对策是相同的。 二人非常数和对策也称为双矩阵对策。也有纯策略对策和混合策略对策两种策略。
1 G0 b# }. X1 x( C
' c, n3 y4 c3 ]" E1 k5.2 纯策略问题 ---分析囚徒困境5 I+ [/ q3 M1 S. M# A3 X) N% N+ X
例1给出了典型的二人非常数和对策,每人的赢得矩阵是不相同的,因此称为双矩阵对策。) y: [+ G* F5 P1 a" s6 f6 |
9 c6 I3 j# s1 t5 G% S+ k# t9 j
问题分析: 这是一个二人非常数和对策问题。从表面上看,两犯罪嫌疑人拒不供认,只能被判18个月徒刑,结果是最好的。但仔细分析,却无法做到这一点。因为犯罪嫌疑人A如 果采用不供认策略,他可能被判刑的刑期为18个月或7年,而犯罪嫌疑人B 可能判的刑 期为0或18个月。而 A选择供认,他被判的刑期为0或3年,此时,犯罪嫌疑人B 可能判 的刑期为3年或7年。因此,犯罪嫌疑人 A一定选择供认。基于同样的道理,犯罪嫌疑 人B 也只能选择供认。 选择供认是他们最好的选择,各自被判3年。3 O' |1 ?2 f. Y, s5 [: m
  x3 }- g# Q: u' @; j9 ?+ [' X  R
5.3 双矩阵对策
$ R& L. C. D5 W1 i+ k! F
' w2 T* |: V6 @! J2 X. {& e% @+ T
: m8 r/ F- ?- h. z
( L* }/ O* L( ^! i- U$ _( u
9 x; i* r! }: Z0 _  ^* d7 E2 ]9 e6 f' j$ q. }
* i" q, p+ N. C# Q6 V! \
5.4 Nash平衡点 / 纳什均衡
; ?& z* E7 C* `( S. g" O) ]5 P# C7 @" d, b$ E1 r

  z/ Y. A: H& U1 t& u
: x4 I& X/ t, G' [$ v5.5 混合对策问题
' |5 k8 w4 b. f5 C7 n如果不存在使式(4)成立的对策,则需要求混合对策。类似于二人零和对策情况, 需要给出混合对策的最优解。
# }! @* _2 N  F4 }& [
* O6 h1 f9 X6 f; e(1) 混合对策问题的基本概念:赢得值9 P* a. B# W2 I1 G
2 c% s( }$ P' @& ~% z

& y# |" M8 a  Y) h% M  N8 B6 S! R. i2 h' T0 e

9 ]+ a8 W0 i; b5 u# B对于混合对策问题有如下定理' A, s& }, H+ L9 `3 g
/ T4 D6 N/ E& d! I2 U$ y+ B! J( `

5 _" I2 Y5 t% Z
) H5 u, }5 z9 U! q(2)混合对策问题的求解方法 , I+ I8 Q) d$ @- C3 f
由定义6可知,求解混合对策就是求非合作对策的平衡点,进一步由定理8得到, 求解非合作对策的平衡点,就是求解满足不等式约束(5)的可行点。因此,混合对策问题的求解问题就转化为求不等式约束(5)的可行点,而LINGO软件可以很容易做到 这一点。( r  t( w/ ?* p5 k! N  X

7 s! h( g. _! Q$ @, U9 O7 J例题5  
9 c9 E! n0 N% }3 b  r$ G有甲、乙两支游泳队举行包括三个项目的对抗赛。这两支游泳队各有一名健 将级运动员(甲队为李,乙队为王),在三个项目中成绩都很突出,但规则准许他们每 人只能参加两项比赛,每队的其他两名运动员可参加全部三项比赛。已知各运动员平时 成绩(秒)见表 3。
$ P/ {- n. [- F* I6 u9 n$ B' p+ |3 _8 X4 y. U' s

- L- C9 J3 n: w, I8 y% N
$ e) H+ c2 R) |6 _
; s2 y4 C0 ^( B7 `5 ?. R( p3 gclc,clear
( z4 u- Y9 i8 Y' V) q* Wa=[59.7 63.2 57.1 58.6 61.4 64.8* E  J# y6 q' O: u
67.2 68.4 63.2 61.5 64.7 66.5
2 k5 `4 Q( i8 t* P$ J' o5 ?7 b' n74.1 75.5 70.3 72.6 73.4 76.9]; 1 T! K7 u0 o  P! i
m=3;n=3;kk=3;T=1000; 1 m, ?7 E% p  p9 ?- B$ u3 s; K
sc1=[5:-2:1,zeros(1,3)]; %1-6 名的得分 . b) f* Y, F$ U7 {! X6 e
sc2=repmat(sc1,kk,1); $ l- A' Y" `$ n+ M# ]8 ]0 f
for i=1:m     ( g# r' C. @$ D" X
    for j=1:n         5 {0 ~$ l( y- \6 v7 z, U; z
        b=a;         - m: O  n" I5 N7 |' ?. C
        b(i,3)=T;b(j,4)=T; %不参加比赛,时间成绩取为充分大         
: \; T* i+ F1 r& S: {2 |$ ?0 w' t# \        [b,ind]=sort(b,2); %对 b 的每一行进行排序         . V. F" B6 Y$ n9 x+ }" l( s
        for k=1:m             ) S/ o$ x; r8 H' B! w
            sc2(k,ind(k,)=sc1; %计算得分         0 [# R6 E( G; H; }3 \
        end         & E$ v1 I" {( L6 B
        A_sc(i,j)=sum(sum(sc2(:,1:m)));  %统计得分
0 V. B' J# W5 E( S2 I% e2 t        B_sc(i,j)=sum(sum(sc2(:,m+1:end)));     
$ x$ _% L% ]3 p1 w( S    end   
- e: ?6 U  a; ~3 O( t5 ?6 m' dend ! h. B- D$ A+ [$ A
A_sc,B_sc ) ^/ [6 L2 g4 R3 V% t- V
fid=fopen('txt2.txt','w');   w4 a2 k+ Q# s5 I- b4 G
fprintf(fid,'%f\n',A_sc'); 7 K; F) x! t! K1 l/ s' _: {* s
fwrite(fid,'~','char');       %往纯文本文件中写 LINGO 数据的分割符 fprintf(fid,'%f\n',B_sc');
7 G1 J% W2 G1 j! F) I2 f% m$ d: wfclose(fid); $ L; F9 x4 c# h1 h0 i+ m1 N( \5 `

0 U. D1 I; C+ B" l2 B% y# D. ^
$ a' g( w" P; P
0 ^6 y$ S7 J1 q0 J; C$ `5 B, Q* z7 b4 T3 v- d( _4 ^7 d- z1 \* b
& S2 P8 C" i9 W  {) |
按照定理8,求最优混合策略,就是求不等式约束(5)的可行解,写出相应的LINGO 程序如下:
- e+ P1 S& C. e; O: l7 x
* Q3 M+ @7 h6 x8 K# ~2 ?/ Gmodel:
2 v/ Q- v$ a0 x: `2 W1 wsets: : s* ?0 }! p8 q0 [  a9 C
pa/1..3/:x; % |7 t. V$ k8 k5 P" Z! ~1 R
pb/1..3/:y;
) v: }: d5 C0 x) |3 dlink(pa,pb):c1,c2;
7 w6 \5 e# m, K& ~- t* r1 ?3 Sendsets
7 ~, M. L, `" z! M; V6 qdata:
" b; e5 N3 j( z* d0 V3 b6 Mc1=@file(txt2.txt);
, ]9 J8 z3 C6 q: d( h0 xc2=@file(txt2.txt);
" \4 U; _0 r+ {2 x* Q$ senddata . L# m1 f0 p! g
v1=@sum(link(i,j):c1(i,j)*x(i)*y(j)); 0 z$ V! S6 Z6 E1 y0 O
v2=@sum(link(i,j):c2(i,j)*x(i)*y(j));
* X4 b* S+ N# \# ~4 }; @@for(pa(i)sum(pb(j):c1(i,j)*y(j))<v1); - R, K; V+ m$ j  h$ W! m& |
@for(pb(j)sum(pa(i):c2(i,j)*x(i))<v2);
) Y$ }8 k& s! l/ E@sum(pa:x)=1;@sum(pb:y)=1;
) P$ c3 q) k4 t9 Y@free(v1);@free(v2);
! W; y9 H4 h+ d) ]7 b7 n/ ^end
7 l! a- Y) w( `, G
1 W7 {5 a+ p/ ~& @& b' q0 n+ |, M/ x
0 a6 \, h1 F3 Y% B& p; {. ]$ p# H; c! l2 B2 x

' I3 A. Y: s8 a2 d: N
3 W$ J( u* Z6 A$ N# l————————————————
" h3 z3 [( L+ W6 h版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。8 G5 l7 I$ N) r& _% M2 a! J
原文链接:https://blog.csdn.net/qq_29831163/article/details/89454176& W% Q0 O0 |1 L) k* w0 c8 }2 ]

7 D0 W; y+ x4 P1 e3 i: e$ T, U! g. \& d' r+ Y" _2 |" x





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5