解方程.txt

用简单遗传算法SGA实现的关于控制图的程序

eq0 := pi00*(p00 - 1) + pi01*p10 + pi02*p20 + pi10*p30 + pi11*p40 + pi12*p50 +
         pi20*p60 + pi21*p70 + pi22*p80 == 0;
eq1 := pi00*p01 + pi01*(p11 - 1) + pi02*p21 + pi10*p31 + pi11*p41 + pi12*p51 +
         pi20*p61 + pi21*p71 + pi22*p81 == 0;
eq2 := pi00*p02 + pi01*p12 + pi02*(p22 - 1) + pi10*p32 + pi11*p42 + pi12*p52 +
         pi20*p62 + pi21*p72 + pi22*p82 == 0;
eq3 := pi00*p03 + pi01*p13 + pi02*p23 + pi10*(p33 - 1) + pi11*p43 + pi12*p53 +
         pi20*p63 + pi21*p73 + pi22*p83 == 0;
eq4 := pi00*p04 + pi01*p14 + pi02*p24 + pi10*p34 + pi11*(p44 - 1) + pi12*p54 +
         pi20*p64 + pi21*p74 + pi22*p84 == 0;
eq5 := pi00*p05 + pi01*p15 + pi02*p25 + pi10*p35 + pi11*p45 + pi12*(p55 - 1) +
         pi20*p65 + pi21*p75 + pi22*p85 == 0;
eq6 := pi00*p06 + pi01*p16 + pi02*p26 + pi10*p36 + pi11*p46 + pi12*p56 + 
        pi20*(p66 - 1) + pi21*p76 + pi22*p86 == 0;
eq7 := pi00*p07 + pi01*p17 + pi02*p27 + pi10*p37 + pi11*p47 + pi12*p57 + 
        pi20*p67 + pi21*(p77 - 1) + pi22*p87 == 0;
eq8 := pi00*p08 + pi01*p18 + pi02*p28 + pi10*p38 + pi11*p48 + pi12*p58 + 
        pi20*p68 + pi21*p78 + pi22*(p88 - 1) == 0;
eq9 := pi00 + pi01 + pi02 + pi10 + pi11 + pi12 + pi20 + pi21 + pi22 - 1 == 0;
st0 := p00 + p01 + p02 + p03 + p04 + p05 + p06 + p07 + p08 == 1;
st1 := p10 + p11 + p12 + p13 + p14 + p15 + p16 + p17 + p18 == 1;
st2 := p20 + p21 + p22 + p23 + p24 + p25 + p26 + p27 + p28 == 1;
st3 := p30 + p31 + p32 + p33 + p34 + p35 + p36 + p37 + p38 == 1;
st4 := p40 + p41 + p42 + p43 + p44 + p45 + p46 + p47 + p48 == 1;
st5 := p50 + p51 + p52 + p53 + p54 + p55 + p56 + p57 + p58 == 1;
st6 := p60 + p61 + p62 + p63 + p64 + p65 + p66 + p67 + p68 == 1;
st7 := p70 + p71 + p72 + p73 + p74 + p75 + p76 + p77 + p78 == 1;
st8 := p80 + p81 + p82 + p83 + p84 + p85 + p86 + p87 + p88 == 1;

Solve[{eq0, eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9, st0, st1, st2, st3, 
    st4, st5, st6, st7, st8}, {pi00, pi01, pi02, pi10, pi11, pi12, pi20, pi21,
     pi22}]