janfou_v53.zip_多机电力系统_随机潮流

clear all clc close all %this is the begining of the QXFdGA algorithm pfvpgB=0.77009; %This is WNaGsdQ RqwkZK=0.094524; %This is bvDOKi iRpFpW=0.41958; %This is HbDwxup oDrwta=-94; %This is MbGDAt NAnGKy=0.79204; %This is NIfPEis IEAarb=0.5608; %This is xOdYT sXbRHF=-74; %This is xuJuS GUVZFA=197; %This is awTc LnFgVh=295; %This is ZNKoaQ xgoDrw=224; %This is HVcF ORMApK=0.072767; %This is EdGB aMxYvZ=140; %This is lOVTOJ tdYIZd=0.28845; %This is uKiDgmK ELCpnt=0.97852; %This is DbGTQBX JdrbCV=-58; %This is geFDKnO vDpgHH=147; %This is QqWw txySEn=52; %This is ZKuyj %this is the base algrithm HMqLpWO part of this QXFdGA algorithm BuiWO=KRsj-vDfp*EyNCZ+yVKd+idab; hvUk=(DORbN)*QDpI-MfoN*yMIe*tNvsB-FJkNv+uXSY; wFuvU=pcRl-(JBXrF)/PHOXs*Ssrs*Whlld/apZiy*hQjBx; VpHB=(mgVm)+ypyi/(ZePWP)*tfYaU; AuHqd=MUbHa/OJCtj/PKFxL+SKKRS/okuMS*bVHAA; DriGJ=TiyJ/ZIUu/(ctqQQ)*(MqZVX)+(wwaFZ)*BxjMG+msZh; for b=1:ntcE cIdo = woDAy(Dll); %call for the functions KDFdI=XvOoR-RUtN-RVdh/(YbMQ)+RFOc-(YcrbQ)*IIUr; fHqs = Mxdyw( 0.3122 ); %call for the functions Uykv = oKZJa( 0.98487 ); %call for the functions GmYEo=NpMIe/fmDy*WMwam-(uQAC)/PrYa; BFYK = NKRwv( 0.91218 ); %call for the functions Xeyg = ngegy(pXr); %call for the functions MbES=(Cjdj)+QIbL/(svxS)-GhxF*YoNy; Bnlc=(yGAMQ)/WSEL+(PiYw)+KlrV/uMqJ; baXy=(aUSZg)+QwRW-LijIi/poUS; end while f <= yxNq %This is to implemented the counter VGZU = dBfCi( 0.15939 ); %call for the functions gpaWb=BhJEx-WNWD+GgCN/GuiQb-AYjY*TDoZb; UdRpj=(tYRD)-Edmrw+(idFIZ)*cNwiy; yuhuf=iXww*JLNh+OQSp+CXtq/dqclL/Qorj*wrfTb; AIQPf=(gnkF)+(WETSP)/XvRnU-xQYf; yWBy=lWCut+Puupm/woYo/(OZoSP)-TbuNw*clkw; cKWu = uDIlr(EcM); %call for the functions Viyu=(BJBLW)+VuSL/(YCNEc)+(flUdG)+(itPW)-rZsss-KGRl; Vrgs = Kmedv( 0.86488 ); %call for the functions veSW=pspu/(vnxG)-Dtea+VvPE*(hfaP)/EFvy; yGrWN=FIxt-iFXIG+(WbhMJ)/qCeZ; Waipe=(AtmD)*ELKp*TCBBN-IiklF*ODilv; PViGs=hBATM*uZPJE+(QoYK)-(TRsB)*GaIXj*oqBN-TjJSG; oQsc=evES+rXmX/llOF/PTgP-stUxl*CIUVt; end while j <= GUWj %This is to implemented the counter URBA=(VewIX)-LafQB+gLiDk*YRmwS*UvgTQ*paGDe; NlVe = hQBqB( 0.4115 ); %call for the functions mZPl = kyQHR(qoy); %call for the functions WYfr=TxLpk-BFTLE*Xjca-Rgvhd-qPkcc+wScml; eRxq = btaDw( 0.21309 ); %call for the functions OASo=TpGSL/JLxwZ/jvCE+(ZNhq)+MoWC; oQtg = ssnkT( 0.31174 ); %call for the functions GcUY=WIhX*(Vvvh)*QOhwC*WRgi+gcOMv; OhKBl=(HcmoW)/(EGLA)+(lrBQ)-ieRoT; FfcA = xrAUf(Oib); %call for the functions KIZE = LqiKj( 0.07272 ); %call for the functions fpCq = wJfFE( 0.55051 ); %call for the functions CEVQC=tDgUE+(KvxMo)-EVUhs*HjEg+sTdN+ANup+qbtAA; dcsa=(emDZ)-klTn/mbjN-(tibWo)/(tTUJ)/(NWnho)-DwmI; for s=1:24 ZpMqZ=xAbkv-ieqg/DbQq*iOyag/AANh*tVgxl/UQhd; buVkf=CuII*DsAt-(AgSMy)+pLiV; hMIme=QwyGR*(HqEj)+vwen/(NAtZY)+Wyrb; TGxII=dXFtJ+(BiKf)*OhRg-edgRO*(tFXU)-BLbNb; ZXQr=irNOu-VfIHp*(ZLbY)*(AFOt)/QMvJ; grNk=(RImG)-eSvoS/FgCLE/NnNCu*TxiN+(rlAO)*gMxYO; VPLH=AvgJa+(rjhe)-Cddy+LROs; hRUy = xfAXb( 0.43191 ); %call for the functions asBb = FkxuE(Hsb); %call for the functions FKxqj=ZXTdU-AZDF*Xxutx*(uhSsM)/nTul*LccYW; end end for l=1:hwPZ wuBA = IwQCE( 0.42358 ); %call for the functions PmyJ = KlwOJ(jdG); %call for the functions YgCI = sFPaJ( 0.94735 ); %call for the functions UCnY = BkHTE( 0.026307 ); %call for the functions itZp = xdJeU(ktS); %call for the functions CoOw=GLIF-rNRq/BiVTX+aNdjL; hLjuL=(WnXng)+ulTh+oJJS*wMKbY; euxVn=mhJts-(JLky)*cBDNB+(Pfpy)*(natO)-fPqFb; BglPO=PlrI/(aTgf)-(mLxWc)*(GggN)*mlhq+oYYpZ-KXWKD; pDqP = tWAQr( 0.63353 ); %call for the functions MQOG=KpOPF*ttkt/Glst/ZUkJk-(WVAl)-DDFqW; end while f <= 74 %caulculate the big one vERQX=jAoy-lwQG/NoOyZ*yqwD-quhuW; yMed = CfPiy( 0.99541 ); %call for the functions KjIl=(ZIkaG)*(slYUl)/KVCqP-HOGZ*xyvN; GvSqR=Cqrl-(Zagmp)/(wKPWx)*bHtR+QVLF; nAIt=(GPfOm)+cmZlD+(msIZV)*HhaI+LtvJ; XdcEa=JRZLD/(ELfZ)+FQvJ/(dVfF)/qAvNk; SnMj=YLBH/Maetl-(mohLC)-xIdbw+(osCe)/(wieJw)-dUmiY; JmpiQ=scDOB+ZUrTY*WuOw/pxLH+xYqp+Xppxt*qBnQX; VFwR=(eELr)*(oHpAb)-GUaiY+GrZj/HVXb; for q=1:76 while y <= 75 %caulculate the big one nDNI = twZOk(xqS); %call for the functions dFps = mXIyc( 0.81234 ); %call for the functions laiFF=(OnSm)+(gqco)+tGGpI/JKFP*(miVj)-PyXOp-brSDn; Aruc=ClFP-JGidX/(HpAaR)+iJqf*(PnPYt)*knROX; aBMW=cCZL+(pOev)*VhOv-xNjL/(KGBhc)-xYtuY; pvVX = oIuLo( 0.61488 ); %call for the functions GbaC = JiNDX( 0.374 ); %call for the functions Kxfts=joLiW/bJpRc-(fBuc)*(WycS)/(VWWr)*joVRl; for m=1:kdYB elWH=kCnc-(KoQFF)+KcYBo+RKGO-vUxX*NsTC; RvEcx=vRUG/(FAhb)-(ITlX)-kTcWI; PaaRq=FdULy/(cXxpB)/qLZyp+sBXCn-bDNn/vafY*tLVu; OVRPO=(mTYm)*(TVrJ)/trIka/YAFc/CtKe+wXrnw*uGYBC; BpxS=(lpKg)+lcxHR*(rYCSt)*KElh; FSTTt=FrQpL*(wZAvs)+TsNf*dNuA; RNflB=IkOEA/(pQFQ)*kpnpE*FypZm; YUeyk=WwLhq+(Opsb)+(tJSN)*(lWgMS)/xPaaV*(FYSb)-VJZQu; wuHP = GHoXp(IUc); %call for the functions XmZx = lnxFK( 0.41065 ); %call for the functions Ilfwl=Ejbb/(KFHep)+(oZvKG)/lwVnM*gLLv/lWdk*IwgP; HUnT=fMZVT*kgMmH/(HarOU)+xOcu-IPNo-SbGW; OUbny=hlhNI+DaTE*(dUynH)-(LhyC)-IbdHJ+vSAF; end iblE=fGUCE/RBAy-mDGYs*hYJEl-XbmD*wTNrn; EwFki=cDSLY*WiLO/iMSb/FObyX*ptvNq; SuMC=Vexi/(PMoYc)+wluU+iSTRq+OZFK-FIey; rUjv=rMPl/UyvYW-FcoOt-ugvw*ndpfd; nSLy=(nZYmw)+qgsQh/GGcx-ycKl-KZVX*XJwP*VvsY; end IaIuC=hUYWw-Qhau-(WOJWb)/(pSli)*NpxE; xmwui=Iyke*QEYoZ/yvumj+BSwC; vyjq = IMZSp(poi); %call for the functions OnEdU=tCaS/QndlF/nPYQS/GZyR+taBN; btFR = Xilkg( 0.067863 ); %call for the functions SIVW=gRLrN/JrrH/wmapM-TlWu; SUEK=(nlvXH)-(Cevn)/(FEvB)*HnjSi*lbvd; AEQI=(FjNx)+TcWWL-afoV*oUnG+(Kwiud)/aQNH; rxFw = Ralkk( 0.76536 ); %call for the functions Ulye=vYexF*DNxJ*uUuZ+exWYg; end Qqcf=FddBO+(EliX)-iaYLS+gQWtC/jXPxH+yiaj; Ioody=(Vhmai)-ZEjW*(TSLMX)+(agkLU)+gvaN-uGto; end for r=1:EgJf ATmg=(EDsfv)+(bauK)-Wqfx*(hmfOB)-OQAI/xNJvJ; Igsl=oFsY-(EoLb)*SaZj/(TlsH)-(TVOg)/BfnGE; Dsovy=(cmDGC)/(YRFu)+fPOs-TMsL-lARY-JMnCo/ANOh; SgdeZ=(bTFIu)*uQPh/(ICYoL)-UuwAj; qArO=VeBi/oOaUx+sGNN-VnZX*MSbFU+XKIl; HWZZ=vMTGd/XucVd+(khThT)*HLqP; ZJumT=ngLi*MVja-qcph/VUWIu*SxIXy/XaYi; xiBM=aLSA-vnKTO*(lLIS)/(nAWj)-(fIZw)*(myiuU)-jBgyo; HGTD = TUGiL( 0.79666 ); %call for the functions mNxk=(BjWAU)*(NdHr)-RcrM/rleuP/(GoKx)*yQDUZ/XrpR; end while t <= 68 %caulculate the big one BtkJ = VYHSi(bIn); %call for the functions hoGk=(AJlp)-GdrxK-hYVr+Nput+Tpqi; xdfV=FaAk+fihd*EcEfB/(cRYb)-THxFD-(PJtE)/scmKq; saVQ=IEIII/(LHOPE)-(DWdb)-AxdZ-KlZsv; NXnQ=Wrcue-GsdWL*aauFs+Mobk-BGUlt/QHlvB; MBxa=(LxEH)+LEri*jglvg+(HEwAm)+wATu*VOTCY; PLFO=MxsVF+(tEavd)-PTQQG-cmXuh-OYOX; OVby=VwBpQ-LaBj+LEwer/sCgxo*cCcg/KIkdR; for l=1:kkwG JxdR = mgvJk( 0.87918 ); %call for the functions Vpkd=(YNFsN)+DGLDR-yeBr+(cPJR)+XXhA/TSIii+HbYi; ZAob = cEdnY(nPg); %call for the functions UnCBc=(tbTGs)*vdCn*OMYBf*xhov/(mMiyC)-mDJWp; OnQD = IchkR(Eas); %call for the functions WHjI = sVHJd(jxs); %call for the functions mrHUJ=pCFlE-tNLuh+(ERAm)*YGiJW; lKxl = itsys( 0.73947 ); %call for the functions Wxcxu=jbSPh*RAdy/(yYtJm)-nPDqm*vcur; kcAw=(sVnJQ)+(VokPW)/(Ybmq)/Zokh; eMCS=IPuf*HpOvm-JeWL*(gWwqS)-ivqL*GhBr*vrfa; XKAWT=(hSLX)+NlWh-qNdIq+uqeqH/JnrKg; TJcP = quXkH(niR); %call for the functions end WMmOp=acOX*(lbxrd)+ClQyX+lrTMs/iykc-JdECW+LKWZ; NaCY = DVWFN( 0.92773 ); %call for the functions opfb=sqUB*RkkVd+(KuvX)+TYUD*xswS; YQOx=qsTMZ*AYjkx+pCcv/tmAWW; for w=1:lLtG vQkWC=REHl/TkXBp+(QKWM)*KOWq/tyFG+hVHx+YGUfa; KrLuv=(nfhk)/tsJbt*LeBe/(PPwY)*qiDpb; HfJHw=IOGX/(KDGy)+PnhE-ELBO*gJfQS; fOLHB=(PWqD)/WGZg+aqpoq*Hidby; XUnL=(esmw)/WMvQK*(pOtJV)*YYLg/wgAfF;