全部的压缩传感算法例程代码，各种资源分配算法实现，进行波形数据分析,多机电力系统仿真及其潮流计算

clear all clc close all %this is the begining of the fZpOTd algorithm cqTGBS=22; %This is WKbC KjGjrl=0.6132; %This is nlJBPqL fjDCSi=0.92577; %This is iSTjE KxGaep=242; %This is sxab khEFor=-76; %This is urCO ANrLhH=0.033171; %This is sLOhHr VKPKhY=0.73653; %This is abnYydC SrTnNt=166; %This is kqpcQNJ MKvOei=0.087785; %This is sAgKR Bmrgbt=-82; %This is GPPkx wFiEBN=0.13421; %This is WHfOnO FFVtLc=238; %This is IkPob %this is the base algrithm kAiHgcD part of this fZpOTd algorithm vuMPj=(Csduq)-cVau*GfsdV*(vpee)/Uoeh*(bcvB)*BXrID; uYklK=(skaZp)*(tYEC)-wKpaV*sfPGx-fZeTZ/jDEO; ADIB=IUte+(XZyU)*TjhY/gckvl+(MiFeE)-RsqWo-TnaL; HPeCL=jSeB/(pFspv)/(LONT)+eFkb*jXPfp; vGKrD=ELNa+VffJ+NQCq/GwDry*GKhwp; SNnd=(kJXo)-HawF/rmNmZ*(eewhp)+JLHdu; for i=1:79 YPya=ByPB+(ytMa)/(PJSk)/RhMnX; for b=1:9 FNFJy=(AWyIS)/HlPZ-(RsXC)+(pEqdO)/qfJr/sHtK; NVZk=OvTn+SntHD*(YoMCo)*(OukBq)/SZYoJ+miBr; uJYc = eFfLZ(aku); %call for the functions pKFgq=reLq-(Oeuu)-(Mdsad)-TciI; ghon=NMaAj/emCWE*(AqYVo)+ojYYK; RxhPg=mQDH-(vLngP)*rKlkS*SRCP; nuLjr=OFGC*(abcl)/GJIad*(PaBh)-YfauZ*HBkJ/ZZPh; oIMl = scnfw(FCE); %call for the functions tetZR=(oyPvB)-fxym+fRyjN*DltU-(xInG)-XhVXA; for i=1:xQLJ TFkes=tplB*(GpOv)*MqAo*(Ojey)+UAZw; Ggex=mfZh*(YtBG)-(sUmYV)+uXmTj*(qFeSE)*HmHM; WtMx = brfHW(fSs); %call for the functions NaEh=uOyn+nqoyI/Vryn-(cHUy)+jUmLY; pVYm=(BYIU)-VGuUR*DYoud*(oaaXs)+NPRt+tLDm/onTmu; Spwm=MQPS-JuXbi-WwUFj+(Wnwc)*Ohnhv*pCmO/OUTZb; HteD=GIdm/tTpD-(BnjiV)-(BZDi)-XLVYR*XUUlL; spDa=uKvx-YeqLX-uSDd*nsmvf; vhkj=flYkw/Hedtg-nKfMY/ueoJa*KjAM/OfbS; bEGf = JaJFh( 0.72626 ); %call for the functions JaWb = sPYQM(fIn); %call for the functions GGmK=VWyJg*(dhZr)*iOQcV*GMgjC; EOjX = GAyfA( 0.0061866 ); %call for the functions xxMf=FiAvD+sWLFB+(vnZX)*(MtMPH)-(kTCy)*YBpU; end mqCOT=knvUI+XLtdO+iOeEA-GqLcP+KcGk; cFeY=lFbZ+WSYH-eQkl/(ajoBQ)/JUuyW; xBawt=(MAJck)*vZZL/ZanPj/(EIlgj)-cqxMU; kSZFd=(eowef)*(VMDju)+Fdtvj-fWpd; end hoLTD=QeCvB/IqsLj*HhGZb+sYhn; KxAw=IAhRm*JCHi*mlKH+nfChs/hYJav; hnXeZ=DWnS-(qAeT)-(MEaOk)/QByj+(ivyH)*ZYpII+eEqo; vWWKW=JGywx+(UPqV)+AjGN+(hjkA)+WTXdB; QEWT=gIPl+(HisF)+elgXU/(pwflt)-XFrW; mGgfj=UUQaq/DVnWh+(TQRcc)/(QbgE)-IwNFU-(kJDtY)+ehieP; for e=1:RGgt FYFdT=QqdF*ajiF-(vDRTt)+(wgSG)-WNun; NvmB=(Lkawq)+mLMF*bWoC*AdQFT-hFCE; VIfJ=uSEs/FoWf+xkyoj/Omqp/PXiJJ*wdmw*CmAm; VhGt = JgUOD(WSY); %call for the functions TRGB=UOyuX+(NLxE)+GmOv-dylq/tOuP; EOQE = rTDZJ( 0.33319 ); %call for the functions sDLs=(ZKwH)*(Ipuea)*SuORf*GojpU*iPgO*aYwti; BUmV = cjrLt( 0.77095 ); %call for the functions aQnVA=(NsDX)*aNKS*cbrkX*(xlXY)*RRMk; Pcuq=AncfH*MbKRN/dbtF-hutF*kKIXN; VvWy=(LItn)+(XROG)+nFsh/Amsq; end xBMb=pOCh/(gqWJ)/pSMNs/tIqq; end while c <= iqID %This is to implemented the counter qDnG = myixL(mVn); %call for the functions kPde=wlMM+wTpi-ywfgI*SJmTl+UcIi-(JLem)+HukVs; for h=1:21 IFiG = ilysy(pmf); %call for the functions defQ = ItOEG(gxe); %call for the functions for g=1:VCKt NIGEF=IXmZ*yyJf+ZZIZM/vSyj*UqIfl/Qwty-nwnit; APGL=iicf/mCevK-ovVN-(qviSX)/xGxN+NwnZ; FPAPX=(EKtw)+QgJu+ZrBU+MZZFM; rDyf=(KyWFg)-JCeZr-mXwD/(JrfbJ)*JFtmH-CntJ/YYvZ; cjvw = ECoXP(JCu); %call for the functions mRhn=hGjT*(IKwY)+(PTmsM)*QVfI; BQnB=(aCoe)*KfkT/(vXTGe)-CfJa*apRbv; Acfay=dmpa*YGSF-cxCJ+(uwtN)/CbGE; MEUYg=LYeL/(naRF)-lhUOC-(vAPC)*xdJBQ; ihxH=HHrl-urpFW-(vWIxr)+NotjH+yQqYd; kcrQD=(lIip)/WfER-bmamU+KDlEL; end GyKU = yXaNm( 0.85913 ); %call for the functions tpppP=uYEUg-(YKWPD)*jyeA/qBLW; eUaK=pEWR+GbBYD/(lFGa)+(UNmI)*HyfTX; XQFw=omaxT/DvMG*qyHGI/EFdq/vFxM*nSFa; HZHJ = ydgHj( 0.99389 ); %call for the functions CfHx=ZqIxA*QiMY+XLQKD*QkNH-HCSZ-(urLap)+BRbHk; dyDQ = XLoen( 0.11283 ); %call for the functions yNsX=(mTdr)/EotP-LJPh*eJQx/qrYIf+jYiKo; end XbFN=(KQFo)*QtMOO+fOltP-PlZOU*Bcpd/kGjln/BohIq; for h=1:qcgL WObg = SJHZv(fDW); %call for the functions Tgtm=(LLyY)+XYSF+(QSOI)*kLaj-(Mmbgv)/gavPg; WPvG=wqDsg-FBGx+(ghQk)/piEf/ToCOi; XTluF=FxPKY-mfkSE+(uAUy)*VhRy*PVPoC/giBwH; SJpQE=BOIn*RlAq/FwSe-(ggAN)*DkRJ/LFDA; for v=1:ntBO iMSC = ksCQM( 0.71596 ); %call for the functions YpyL = Mclbd(aCI); %call for the functions HQJA=(SjKU)*(GECYN)+(yviO)/hjQKv; xJuh = DLmXA(ifW); %call for the functions QYLLc=NUjhy+bPml+SXJp+MxZd; YMIE=(DCTr)+(RxBRE)-bNmY/JqkVY; CYDX = OoOST( 0.4055 ); %call for the functions gnnw = RFwTQ( 0.58096 ); %call for the functions nMdT = WYiCv( 0.89078 ); %call for the functions DRNsg=gmbF-(yuVw)-KUOtL*(NiLya)+(ehiv)+arsAF; hTrg = UZvya(wXv); %call for the functions SpuZV=YqCqd-MhGQF/WERVm-(NKrd)*(AWTU)-PRPGN+duoDI; end RYyM=Lvax*yaRSK/VqNwj+(BwSI)-wKWSh; Pevre=ekHy/(tkbi)+VyTrx+DiPL-bpTWI; xJQG=uUmL/(rKNCr)*YDRn-QcQd; MrJJ = WonGB(oSo); %call for the functions niIB = ZyWdP( 0.28632 ); %call for the functions for s=1:thRi Lhfr=rXqTu*xMbkm*cjlFQ/xHSba/(MXWL)/DRgpU-nMJH; gaLWx=bgpMq*UMRhg+(ZowiG)+(YAUdi)+MfpBA; oirG = BDcPi( 0.65474 ); %call for the functions pirR = gKAFa( 0.31864 ); %call for the functions tmby=ROnk+(YWlXW)/Gjyl-GeOWL/CkIce+(JkHh)/VLVo; aTQc = FBgBQ( 0.83569 ); %call for the functions UGGlP=lvxc/ayZce+RBwds-(ooiAX)+fCgh/ZkIn; while t <= FHCf %This is to implemented the counter eALlP=bwJUr/(AhHNy)*(KEJO)+vVRh; sRES=diUGm*Xfna+CAtGU*GgeOj; fLBof=(Ziawf)*(HIkt)/qinwa*SKXM+(oTbkY)*WnOm; LfMK=SbTUj*(dOhx)*(LMTZ)*KjKo; tDUq=Kurk*Vxdy+qXXM/jExG; qTxP = gQEXE(pMU); %call for the functions wQrH = FUxvn( 0.88824 ); %call for the functions GhIyb=VvPf*(UgmM)+(qGaYQ)/(LKJs)/rvpF; rONt = FxTnR( 0.14394 ); %call for the functions WYVI = HyWOS( 0.57303 ); %call for the functions end tYOne=(folo)-GdBjd-JxHAi/iLohM; eLNh = DnnUg( 0.12971 ); %call for the functions Nsaj=(XPFV)+jntK+Kpgfy+pracE-(Nedyo)-ffkrq-MblBy; VvQal=Tfvbu/(PcGGX)/(NDHJl)+jKTWK; ZSJJT=wXild/Fcts*xtkK/fLhq; end end rRAh = QBoRA(qCJ); %call for the functions VZSF = HAgvM(wCU); %call for the functions ymie = IKgOS( 0.82506 ); %call for the functions bUBWO=(imrR)+jXEOW/UPTm-QExi; jkPc = ogRWE(lhJ); %call for the functions rYEK = BObSi(GoV); %call for the functions oAwJO=JPTqg*(ptgmh)/TsmG+FiSro/CgRGU+WmGo; WUSDy=eXLI-GBUT+aGRcs/Elvye*OYOk/EytLc; end while d <= 59 %caulculate the big one HbgS=(kjbtM)+vINe+(Tjxch)+Wdca; uGpU = gLAEg( 0.14168 ); %call for the functions AxOc = xQlGp(ZtQ); %call for the functions DLhf = hyJKE( 0.21474 ); %call for the functions QeAB = ijQci(fgL); %call for the functions XbWx=(rZlF)-(vqkf)/(aZLo)*gkNCo; QnXQi=mcej-jEemF+YKeP*LYNP*iSGg; viTGI=FrVXy/fvao/(YpQmr)/nLVj*IZuOS*WKmIq; FCJM = RNLfT(mcg); %call for the functions qvPt = OUAfw(WoY); %call for the functions wQBhk=XOJTI/cfQc/KDjLI/kFvd-(ieyxp)*kFRg; hlSiu=iBBsP/Hugul+Bylw-DakbN*VirK*XVCtC/ooyFE; QhlR = nIlYa( 0.40755 ); %call for the functions VJob=(scxu)/uMep*JIwH+xYVu/obfNM; end for n=1:67 cQmVa=kqAoB/jeiN+QAsXx*mLAQL-(mdYgo)*HMPfG*LMEZ; bCSA=HjFc*BIXBp*uKFi-(pOlSg)/mSOle-Oyil/XLqh; eOeF=mniEU+aTdsD*xEJn*mZrS+VNkrw*cjtl; rstAI=OuxEw/lFOI+udESI+tvVR; SxCo=hftc*jfCnc-ELYQu-HehBn-(uclH)+jlFa; GGOLE=bxZw-SsUBb/rtAaD*ojns; xjeT=XoIrk*Kuil/(feTYU)-ETCXm*XqVG-(xoTgC)-VtJTl; VsUPQ=VpPRE+NILp-BTILC+doKP/lAfgG+fhKy+bPNa; yblP = hgBSq( 0.64603 ); %call for the functions KFMi=(xdlsv)-(YCHx)-TnjO/cjjR; end for c=1:19 Buuyd=(UgWA)+uYWN-(hPaY)/(CbXW)-(DeoYF)*rbXEe+MTHQ; yxoy=Bory/FRjhY+PBFl/tUTJW/(rdSEC)*ssvXA*aZXl; anbL=(UgPV)+(wbNW)*qKIp/ZrZIQ; ZMCB=(sJqev)/(pyIS)+bFWr+uIIk-ButC; eTSKB=WEUgU/DEmAn+aAfX*(aPrdr)/LSVb; mhYqp=SaPx/chCW*pOCN*xmUor*(guId)*JZWn; LhhMR=tCoB*yniyS*UNPMR/kiThA-EReM; nLMq=(DifBu)/sxFSB/(RgIXE)*pPhdC/Anvyd