Добавил:

Fragga Опубликованный материал нарушает ваши авторские права? Сообщите нам.

Вуз:

МИРЭА - Российский технологический университет

Предмет:

Теория принятия решений

Файл:

Курсовики / Теория прин. реш(теория игр).doc

Скачиваний:

Добавлен:

26.05.2014

Размер:

122.88 Кб

Скачать

☆

<<< < Предыдущая 12 / 22

3. Связь с линейным программированием.

В соответствие с введенным понятием цены игры и записи его для случая применения произвольных стратегий со стороны одного из игроков, решение данной матричной игры не изменяется от прибавления ко всем a_ij одной и той же положительной константы a, так что не будет ограничением рассматривать игры с положительной ценой V>0. В этом случае оптимальные стратегии p₁,p₂,…,p_n; q₁,q₂,…,q_m и цена игры v конечной игры двух партнеров с нулевой суммой с заданной матрицей выигрышей есть

где z_min, w_max, X_i, Y_j определяются как решения двойственных задач линейного программирования

z=X₁+X₂+…+X_n=min

с ограничениями a₁_jX₁+a₂_jX₂+…+a_njX_n1 (j=1,2,…,m),

X_i0 (i=1,2,…,n),

w=Y₁+Y₂+…+Y_m=max

с ограничениями a_i₁Y₁+a_i₂Y₂+…+a_imY_m1 (i=1,2,…,n),

Y_j0 (j=1,2,…,m).

Текст программы:

#include<bios.h>

#include<iostream.h>

#include<conio.h>

#include<string.h>

#include<stdlib.h>

#include<stdio.h>

int** array;

int str,stb,f[100],f1[100];

//funkcii minmax i maxmin

int func(int n)

{array[str][n]=array[0][n];

for(int i=0;i<str;i++)

{if(array[str][n]<array[i][n])

array[str][n]=array[i][n];}

return array[str][n];}

int func1(int m)

{array[m][stb]=array[m][0];

for(int i=0;i<stb;i++)

{if(array[m][stb]>array[m][i])

array[m][stb]=array[m][i];}

return array[m][stb];}

void main()

{clrscr();

int i,j;

//Vvod

cout<<"\nVvedite razmer matrici [Nstolb*Mstrok]:\n";

cin>>stb>>str;

array=new int*[str+2];

for (i=0;i<str+2;i++)

{array[i]=new int [stb+2];}

cout<<"\n\nViberite vid vvoda(1_Manual. Other_Random.):";

for (i=0;i<str;i++)for(j=0;j<stb;j++)

{array[i][j]=0;}

cin>>i;

if(i==1)

{for (i=0;i<str;i++)

for(j=0;j<stb;j++)

{printf("Vvedite element[%d][%d]: ",i+1,j+1);

scanf("%d", &array[i][j]);}

}

else {randomize();

cout <<"\n\t Sluchaino!\n";

for (i=0;i<str;i++)

for(j=0;j<stb;j++)

{array[i][j]=(rand()%8+1);}}

//vivod

cout<<"\n";

printf("\n\t\tMatrica:\n");

for (i=0;i<str;i++)

{for (j=0;j<stb;j++)

{printf(" %d",array[i][j]);}

printf(" \n");}

for(i=0;i<str;i++){f[i]=0;} //massivi 0

for(j=0;j<stb;j++){f1[j]=0;}

while(1)

{int g=0; //flag

for(i=0;i<str;i++) //cikl po vsem strokam

{if(f[i]!=1) //cikl po strokam s prediduchei

for(int k=0;k<str;k++)

if(k!=i)

{if(f[k]!=1){

f[k]=1;g++;

for(int zz=0;zz<stb;zz++)if(f1[zz]!=1)

if(array[i][zz]<array[k][zz]){f[k]=0;g--;break;}}

}}

for(i=0;i<stb;i++)

{if(f1[i]!=1)

for(int k=0;k<stb;k++)

if(k!=i)

{if(f1[k]!=1)

{f1[k]=1;g++;

for(int zz=0;zz<str;zz++)if(f[zz]!=1)

if(array[zz][i]>array[zz][k]){f1[k]=0;g--;break;}}

}}

if(g==0){break;}

}

cout<<"Posle mazhorirovanija:\t\t\n";

for (i=0;i<str;i++)

if(f[i]!=1)

{for (j=0;j<stb;j++)

if(f1[j]!=1)

{printf(" %d",array[i][j]);}

printf(" \n");}

int maxmin=0,minmax=0;

cout<<"\nai ";

for (i=0;i<str;i++)

if(f[i]!=1)

{func1(i);cout<<"\n"<<func1(i);}

int n=func1(0);

for (i=1;i<str;i++)

if(f[i]!=1)

{if(n<func1(i))

{n=func1(i);maxmin=i;}}

cout<<"\n\nbi ";

for (j=0;j<stb;j++)

if(f1[j]!=1)

{func(j);cout<<func(j)<<" ";}

n=func(0);

for (j=1;j<stb;j++)

if(f1[j]!=1)

{if(n>func(j))

{n=func(j);minmax=j;}}

cout<<"\n"<<"Verhnaja cena igri: minmax= "<<array[str][minmax];

cout<<"\nNiznaja cena igri: maxmin= "<<array[maxmin][stb];

if(array[maxmin][stb]==array[str][minmax])

{cout<<"\n"<<"Sedlo!!!";}

int min=array[0][0];

for(i=0;i<str;i++)

for(j=0;j<stb;j++)

if(array[i][j]<min){min=array[i][j];}

if(min<0)

{for(i=0;i<str;i++)

for(j=0;j<stb;j++)

array[i][j]+=((-1)*min+1);}

else {min=0;}

int cs=0; //skol'ko strok ostalos'

for(i=0;i<str;i++)

if(f[i]!=1)cs++;

int cs1=0; //skol'ko stb ostalos'

for(i=0;i<stb;i++)

if(f1[i]!=1)cs1++;

double ox=0,oy=0,oc=0;

if(cs==2)

{int y1,y2;

for(i=0;i<str;i++)

if(f[i]!=1){y1=i;break;}

for(i=i+1;i<str;i++)

if(f[i]!=1){y2=i;break;}

for(i=0;i<stb;i++)

if(f1[i]!=1)

for(int k=i+1;k<stb;k++)

if(f1[k]!=1)

{double opred=array[y1][i]*array[y2][k]-array[y2][i]*array[y1][k];

if(opred==0) continue;

double lx=(array[y2][k]-array[y2][i])/opred;

double ly=(array[y1][i]-array[y1][k])/opred;

if((lx+ly)>oc)

{ox=lx;oy=ly;oc=lx+ly;}}

if(min==0)

{cout<<"\nPa1="<<ox/oc<<" \nPa2="<<oy/oc<<"\nCena igri V="<<1/oc;}

else cout<<"\nPa1="<<ox/oc<<" \nPa2="<<oy/oc<<"\nCena igri V="<<1/oc+min-1;}

cout<<"\n";ox=0,oy=0,oc=0;

if(cs1==2)

{int y1,y2;

for(i=0;i<stb;i++)

if(f1[i]!=1){y1=i;break;}

for(i=i+1;i<stb;i++)

if(f1[i]!=1){y2=i;break;}

for(i=0;i<str;i++)

if(f[i]!=1)

for(int k=i+1;k<str;k++)

if(f[k]!=1)

{double opred=array[i][y1]*array[k][y2]-array[k][y1]*array[i][y2];

if(opred==0) continue;

double lx=(array[k][y2]-array[i][y2])/opred;

double ly=(array[i][y1]-array[k][y1])/opred;

if((lx+ly)>oc)

{ox=lx;oy=ly;oc=lx+ly;}}

if(min==0)

{cout<<"\nPb1="<<ox/oc<<" Pb2="<<oy/oc<<"\nCena igri W="<<1/oc;}

else cout<<"\nPb1="<<ox/oc<<" Pb2="<<oy/oc<<"\nCena igri V="<<1/oc+min-1;}

for(i=0;i<str+1; i++)

{delete[] array[i];}

delete[] array;

while(1)

{switch(bioskey(0))

{case 0x11b :goto swq;}}

swq:}

Запишем матрицу выигрышей, найдем минимальный элемент (a=-5), прибавим его ко всем элементам матрицы.

2	-3	4		7	2	9
-3	4	-5		2	9	0
4	-5	6		9	0	11

Матрица симметрична, m=n=3.

Выигрыш первого игрока будет составлять:

7P₁+2P₂+9P₃V

2P₁+9P₂+0P₃V

9P₁+0P₂+11P₃V

Поскольку V>0, то можно поделить на V с сохранением знаков неравенств. Обозначим P_i/V=X_i, тогда нас интересует min X₁+X₂+X₃ при X_i0 (i=1,2,3) и

7X₁+2X₂+9X₃1

2X₁+9X₂1

9X₁+11X₃1

Приводим задачу к каноническому виду.

7X₁+2X₂+9X₃-Y₁1

2X₁+9X₂-Y₂1

9X₁+11X₃-Y₃1

Вводим m дополнительных переменных для поиска допустимого начального приближения.

Z₁=1-7X₁-2X₂-9X₃+Y₁

Z₂=1-2X₁-9X₂+Y₂

Z₃=1-9X₁-11X₃+Y₃

Z₁+Z₂+Z₃=3-18X₁-11X₂-20X₃+Y₁+Y₂+Y₃ Решаем задачу линейного программирования:

Найти min Z₁+Z₂+Z₃ при Z0, Y0, X0 и

Z₁=1-7X₁-2X₂-9X₃+Y₁

Z₂=1-2X₁-9X₂+Y₂

Z₃=1-9X₁-11X₃+Y₃

Симплексная таблица								До деления на опорный элемент
	1	X₁	X₂	X₃	Y₁	Y₂	Y₃		1	X₁	X₂	X₃	Y₁	Y₂	Y₃
Z₁	1	-7	-2	-9	1	0	0	Z₁	-2	-4	22	-9	-11	0	9
Z₂	1	-2	-9	0	0	1	0	Z₂	-11	22	99	0	0	-11	0
Z₃	1	-9	0	-11	0	0	1	Z₃	-1	9	0	1	0	0	-1
K	3	-18	-11	-20	1	1	1	K	-13	18	121	-20	-11	-11	9

	1	X₁	X₂	Z₃	Y₁	Y₂	Y₃		1	X₁	X₂	Z₃	Y₁	Y₂	Y₃
Z₁	2/11	4/11	-2	9/11	1	0	-9/11	Z₁	-2/11	-4/11	1	-9/11	-1	0	9/11
Z₂	1	-2	-9	0	0	1	0	Z₂	-4/11	80/11	-9	81/11	9	-2	-81/11
X₃	1/11	-9/11	0	-1/11	0	0	1/11	X₃	-2/11	18/11	0	2/11	0	0	-2/11
K	13/11	-18/11	-11	20/11	1	1	-9/11	K	-4/11	80/11	-11	59/11	9	-2	-81/11

	1	X₁	Z₁	Z₃	Y₁	Y₂	Y₃		1	X₁	Z₁	Z₃	Y₁	Y₂	Y₃
X₂	1/11	2/11	-1/2	9/22	1/2	0	-9/22	X₂	-4/11	2/11	1	-9/11	-1	-2/11	9/11
Z₂	2/11	-40/11	9/2	-81/22	-9/2	1	81/22	Z₂	-2/11	1	-9/2	81/22	9/2	-1	-81/22
X₃	1/11	-9/11	0	-1/11	0	0	1/11	X₃	-2/11	-9/11	81/22	-59/22	-81/22	9/11	59/22
K	2/11	-40/11	11/2	-59/22	-9/2	1	81/22	K	0	-40/11	-40/11	-40/11	0	0	0

	1	Z₂	Z₁	Z₃	Y₁	Y₂	Y₃
X₂	1/10	-1/20	-11/40	9/40	11/40	1/20	-9/40
X₁	1/20	-11/40	9/40	-81/80	-9/40	11/40	81/80
X₃	1/20	9/40	-81/80	59/80	81/80	-9/40	-59/80
K	0	1	1	1	0	0	0

X₁+X₂+X₃=1/10+1/20+1/20=1/5 => V=5. P₁=5/20=1/4 P₂=5/10=1/2 P₃=5/20=1/4

В силу симметричности матрицы выигрышей Q₁=1/4, Q₂=1/2, Q₃=1/4.

<<< < Предыдущая 12 / 22

Соседние файлы в папке Курсовики

#
26.05.2014122.88 Кб22Теория прин. реш(теория игр).doc
#
26.05.201474.75 Кб17Теория принятия реш.DOC