【正文】 141預(yù)報(bào)負(fù)荷為156對應(yīng)的線路朝流量限值（）165150160155132162從表中可以看出，分配預(yù)案都會發(fā)生輸電阻塞。最后得出：，可以消除輸電阻塞。對求解結(jié)果的改進(jìn)可以發(fā)現(xiàn)其中線路1已經(jīng)達(dá)到了超出限值的百分比的最小值，且線路5超出限值的百分比相當(dāng)小。在求解時(shí)考慮到有正有負(fù)，選取松弛變量進(jìn)行代換，求和是只求所有超出限值的百分比，最后得到調(diào)整后的各線路潮流絕對值超過限值的百分比，見表八：表八 調(diào)整后各線路潮流絕對值超過限值的百分比線路123456調(diào)整后%00000得到的分配方案見表九：表九 調(diào)整后的分配方案機(jī)組12345678總出力值（）228152155117從程序中得到阻塞費(fèi)用為：。問題(3)的結(jié)果是唯一的，可以利用計(jì)算機(jī)編程或者手工進(jìn)行計(jì)算，而它的結(jié)果對問題(4)起著決定性作用。在輸電阻塞無法消除時(shí)，采用兩種模型進(jìn)行優(yōu)化求解，然后選擇較好的答案作為最終結(jié)果。九、 模型評價(jià)雖然把問題(1)和問題(2)看作是建模前的準(zhǔn)備工作，但是能否正確、合理的解決它們，將直接影響到后期的建模及求解。這種方法不同于題中所給的選取方法，但效果是一樣的。每個模型中都是把安全問題放在首位，然后考慮使得阻塞費(fèi)用最小。模型的缺點(diǎn)：在模型二中當(dāng)輸電阻塞無法消除，就有輸電的安全和經(jīng)濟(jì)兩個原則，但二者又是矛盾的，不可能同時(shí)達(dá)到最優(yōu)?？梢韵纫猿龅陌俜直葹閮?yōu)化目標(biāo)，得到一組解，然后以這組解為限值條件，阻塞費(fèi)用最小為目標(biāo)函數(shù)再進(jìn)行優(yōu)化，這樣就可以得到盡可能安全且經(jīng)濟(jì)的分配方案。uplimit=current+speed.*15。nextprice=zeros(1,length(current))。sum=0。 next(i)=next(i)+volumn(i,j)。 if price(i,segment(i))currentprice currentprice=price(i,segment(i))。[minprice,j]=min(nextprice)。 else sum=sum+uplimit(j)next(j)。 else sum=sum+uplimit(j)next(j)。 max=remain(j)。if sumdestination for i=1:length(remain) remain(i)=volumn(i,segment(i))。 if sum+remain(j)=destination if next(j)+remain(j)=uplimit(j) next(j)=next(j)+remain(j)。 end else if next(j)+destinationsum=uplimit(j) next(j)=next(j)+destinationsum。 end end if segment(j)length(segment) segment(j)=segment(j)+1。 end endendnextcurrentprice:120 73 180 80 125 125 9070 0 50 0 0 30 0 0 0 4030 0 20 8 15 6 2 0 0 8110 0 40 0 30 0 20 40 0 4055 5 10 10 10 10 15 0 0 175 5 15 0 15 15 0 10 10 1095 0 10 20 0 15 10 20 0 1050 15 5 15 10 10 5 10 3 270 0 20 0 20 0 20 10 15 5505 0 124 168 210 252 312 330 363 489560 0 182 203 245 300 320 360 410 495610 0 152 189 233 258 308 356 415 500500 150 170 200 255 302 325 380 435 800590 0 116 146 188 215 250 310 396 510607 0 159 173 205 252 305 380 405 520500 120 180 251 260 306 315 335 348 548800 153 183 233 253 283 303 318 400 800 1 2 附件[2]：，求解調(diào)整后的方案的Lingo程序model:min=(70X11)*(303+505)+(50X13)*(303124)+(30X16)*(303252)+X110*(489303)+(30X21)*(303+560)+(20X23)*(303182)+(8X24)*(303203)+(15X25)*(303245)+(6X26)*(303300)+X27*(320303)+X210*(495303)+(110X31)*(303+610)+(40X33)*(303152)+(30X35)*(303233)+X37*(308303)+X38*(356303)+X310*(500303)+(55X41)*(303+500)+(5X42)*(303150)+(10X43)*(303170)+(10X44)*(303200)+(10X45)*(303255)+ABS()*(303302)+X47*(325303)+X410*(800303)+(75X51)*(590+303)+(5X52)*(3030)+(15X53)*(303116)+(15X55)*(303188)+(15X56)*(303215)+X58*(310303)+X59*(396303)+X510*(510303)+(95X61)*(303+607)+(10X63)*(303159)+(20X64)*(303173)+(15X66)*(303252)+X67*(305303)+X68*(380303)+X610*(520303)+(50X71)*(303+500)+(15X72)*(303120)+(5X73)*(303180)+(15X74)*(303251)+(10X75)*(303260)+X76*(306303)+X77*(315303)+X78*(335303)+X79*(348303)+X710*(548303)+(70X81)*(303+800)+(20X83)*(303183)+(20X85)*(303253)+X88*(318303)+X89*(400303)+X810*(800303)。x11040。x2515。x31110。x3840。x4310。x4715。x5315。x5910。x6420。x61010。x7510。x793。x8520。x8105。*(x11+x13+x16+x110)*(x21+x23+x24+x25+x26+x27+x210)+*(x31+x33+x35+x37+x38+x310)*(x41+x42+x43+x44+x45+x46+x47+x410)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+*(x81+x83+x85+x87+x88+x89+x810)。x11+x13+x16+x11087。x31+x33+x35+x37+x38+x310132。x51+x52+x53+x55+x56+x58+x59+x51098。x71+x72+x73+x74+x75+x76+x77+x78+x79+x710。end附件[3]：，用原模型求解調(diào)整后的方案的程序model:min=pp6pp60。x21+x23+x24+x25+x26+x27+x21058。x41+x42+x43+x44+x45+x46+x47+x410。x61+x63+x64+x66+x67+x68+x61095。x81+x83+x85+x87+x88+x89+x81063。x1630。x248。x2108。x3720。x425。x4610。x525。x5810。x6310。x6820。x7415。x7810。x8320。x8915。p3p30=(*(x11+x13+x16+x110)*(x21+x23+x24+x25+x26+x27+x210)+*(x31+x33+x35+x37+x38+x310)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x81+x83+x85+x87+x88+x89+x810)+)/pp2pp20=((p1p10)+(p2p20)+abs((p1p10)(p2p20)))/2。pp6pp60=((p6p60)+(pp5pp50)+abs((p6p60)(pp5pp50)))/2。p4p400。(x11+x13+x16+x110)+(x21+x23+x24+x25+x26+x27+x210)+(x31+x33+x35+x37+x38+x310)+(x41+x42+x43+x44+x45+x46+x47+x410)+(x51+x52+x53+x55+x56+x58+x59+x510)+(x61+x63+x64+x66+x67+x68+x610)+(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+(x81+x83+x85+x87+x88+x89+x810)=。x21+x23+x24+x25+x26+x27+x21088。x41+x42+x43+x44+x45+x46+x47+x410。x61+x63+x64+x66+x67+x68+x610155。x81+x83+x85+x87+x88+x89+x810117。x11040。x2515。x31110。x3840。x4310。x4715。x5315。x5910。x6420。x61010。x7510。x793。x8520。x8105。p4p40=(*(x11+x13+x16+x110)*(x21+x23+x24+x25+x26+x27+x210)+*(x31+x33+x35+x37+x38+x310)*(x41+x42+x43+x44+x45+x46+x47+x410)*(x51+x52+x53+x55+x56+x58+x59+x510)+*(x71+x72+x73+x74+x75+x76+x77+x78+x79+x710)+*(x81+x83+x85+x87+x88+x89+x810)+)/155。pp3pp30=((p3p30)+(pp2pp20)+abs((p3p30)(pp2pp20)))/2。end38 / 39