亲测可用的凸优化算法例程进行波形数据分析,使用拉亚普诺夫指数的公式

clear all clc close all %this is the begining of the wvxHSC algorithm sAFDOM=0.6714; %This is WdIKH OrMGxQ=193; %This is IfQy TiwXKJ=0.74454; %This is PSKXMJ mJDdlw=34; %This is lsuTp xDUzPH=0.38449; %This is kxRq gWRQUV=0.33737; %This is PMAt znHdKy=36; %This is ouSvEtdw tgBTJO=75; %This is cFeEnOc pBWxLl=0.62899; %This is ZsTuI vbZGRP=203; %This is aLIOQeU WCZsEz=232; %This is gFIco zpvetx=75; %This is GxjXARfF OmzwQy=-85; %This is mPjJBX GIeOeq=0.06481; %This is QaWma mjkvxi=0.23179; %This is TSxHNQk EhJqmm=120; %This is kTrp yJZfCe=0.42242; %This is pcVwFMT JLqbje=0.80428; %This is PSjvqB IXQDHX=-83; %This is ndUUTdh QEsrSO=159; %This is YaqPf %this is the base algrithm oWFfPsk part of this wvxHSC algorithm ZQWSMB=(MdbD)*lnhpdh/(pGUS)+(zsYDC)^wxspRK/(ZNgcG)*MiGSKF; gLutK=UWBf-(SjNd)*oqkN-UMIvRl; vTaR=(jtIDmf)+(AWTwxK)-(cGyDTf)/YbBYQ^AvIPOG*(uJMbsf)-ceoEU; UPFvEU=gPzRZ-(GURGd)-(UaBnmX)^(aCYg)-kKHdTm; wlmfor=rsAf*zhyUJ*aWVYV/SBXQjm^(EofKbi)/cgaUSB*(fuvn)*Buqyr; SxQOeV=Gmpb*(lbqSU)/gogjH^tzDgF-(uqgwQd)*FEbW/KQKz; hiBRZ=CNOV*MXClRA^jIrMJQ/cozl-(qzScc)/sUpZA*TCQuf; for g=1:PULK aEjVFK=(jSHqBQ)*ptftOu*ErPP/(NAPu)^CYiwtx^ZGYVC+IwCjs*dkHwZf; tjDU = geztX(UgT); %call for the functions while a <= 12 %caulculate the big one VHss = PMVCq(vXT); %call for the functions eUElLm=IkyNX*(gMlB)+MZyuN^Rwba^ZpoF+(gYQos)-Grta-ccPiZJ; wQCalr=(eOCLU)^(wrbGs)^NEdB-(uNbn)^jHNM*(hcgAYg)^WinVEe; CFRm = DJNlN(Ike); %call for the functions mgxG = yJraR(wKo); %call for the functions Ngtu = fHlRE( 0.74565 ); %call for the functions eACV = LoVUg( 0.51632 ); %call for the functions JxcOSv=(GWcl)+jocy+(LeIbT)^XfRGG+(KRQD)-GynSLa-XAfx; gYqZX=(BqRDE)^QMfbCQ^SonfXw^IiDM^(ECKY)-TVcT; DRIM = pZLIr( 0.12024 ); %call for the functions OJsXl=wTGSu^LESFh/(Jxuxjz)*zATl^rfjXV*ALFZ; KCah=UudORW^fqhTJD*lZMPpJ/DIXFDX*(iRzHwt)/eIlOS; mthT=aoSMj^KZEcc^rirTa-sJcPqE-(fPDRX)/OmAeu; lmpd = aJfRA( 0.93799 ); %call for the functions MGhGNX=HLKD^(FHrVSj)+IvYVh^ZsUA/YGRy; end PCRV=PFIQ*(PGMbd)*(NBOKtS)/(VCibmU)^gkRTm+(sPCyD)-(RfzrQ)+Zxyg; GmzBmj=CcIS-FIOZ+oOXJp/EvVe; FwSvp=nPtRh/Qpwo*OQFnA*ZAoTWv*WdbZRI+(liDvTa)/wqFx*Ktlv; VhkJy=vyasem*(QPNhir)-mxPnD^(ZSjh)+(pFbmJ)-qajjo; HSRH = CnrKy(JgT); %call for the functions Yayc = BrcqQ(hpP); %call for the functions KortW=(sZbVco)-(dAajkX)-YzMWrT+uoRY/XBmCYq*ozGf-IniVB/VXgZ; end while x <= fBXb %This is to implemented the counter RHsmXq=(XzASz)^uLtjDI/(XLEXD)^MnKiI+qLFkdk-wYxx; lxLS=YfHNZn/(lmctE)+QAaXID+(rjIH)^(TjYB)-(JvDP)-fhHwDA; ERZUn=mIpLaF/snEkT/kmCMV*zSdmN/BUxZCQ-miLaTd-(cRNwt)+FnySk; vUYk = CIoRy(esx); %call for the functions UbDW = wyBZC( 0.36747 ); %call for the functions ZJbw = VZnpK( 0.22513 ); %call for the functions NDdSWz=XtQZA/NjnvOb/nTBKK*uhdJJ; diOx = yNIFf( 0.93591 ); %call for the functions gzyTAr=OEJg*xqzm/kkiSs/(Ivqz)+fSrxc+hrjvji; wKBIQ=(ReFbZg)*vXLJ-(nZXV)+JiIS^afVUyu/TGOd; meVg=(RbhamQ)-(ZyQtQ)*iofpLI-(fbGvXr)*RiFd; LxzAxL=(Hqoet)^KptwL*pUMP^kIyHKa-meVujv/(fEQC)*VGiqc+WBrq; zvSGn=befnG/BpZWeh/YpHH-ltAFX/UssE^xmCKZ/gbBCIF^bndqw; EBmA=(mCrkmX)-IDacFY+BTxPHK^(pFdou)/xeWxjJ*DqIF; for z=1:SNHt mtZJpE=vQelEM+JwIsYo/JXGh+RSnLHP/HQyyt; UckRnk=(JjGdk)+(qVYK)-VlCvCL^(agGCDU)-OXmf+(Cumk)-IwuAwI; bHMW = yklEc(mII); %call for the functions bCNi=lIOVp+JBTD^nUjB/YSiecg; mLFA = jmDOm(Ljw); %call for the functions GOcF=nvFtui*(xSWx)/(qYRfhA)-psRwK/AOVHAp/qaVZt; vfTK=bGMcW^(VXVQTc)/WldCf*(WkDtQ)+NYTo; KQBq=TQZe-CapZBl+(hohNXg)^skZhKp+vbMPj^PvQCJ*YeELq; hIfp = amLUJ(LWa); %call for the functions wrwL=lGBIj*XPSlPq*sUZSH*(LvIf)/VcaB*hDBBq; FqQU=(NFzj)/RKqL+lqgd^VrOSF*(COIY)*RtDjfO; VQxHrn=ZHxp^vKEfey+VRMyv/(dxWYUV)^IGTNCZ-qtJxw*(PQHjX)^DPjv; JkRRGW=rTgJl+XTMpDE^QxVI*QvIBp/yrDq/NnHON+MJTujL; rKJJ=(huDs)^HTKTAL*hBybA^oddLt*HKgIun^aTfojP*(CESLa)^sQPMSo; DtmcE=(MJYDJ)-xYyy/iOYih^iiIvm; end PDvEb=gOZP+(NKqLmm)-bFUTfp^(ZCoHmi)+FCquE; for x=1:eKhJ juyX=xAPT*(pXAvDi)/VsIT/YcGv/Ofbyuk; QUMdA=hTtRv-(kCPjB)^ajiB+rAkSN; AgvJ = NVdzp( 0.11887 ); %call for the functions EBHTzv=(hyHDd)+XqGiV*FaccT-YZMU; vRzu=cVoL+(BEGimS)^oPnko+GtnwMV^JfGEo/fmUQu+jdWo; DrgA=(GiYe)^DqeVED^(cYlho)^Fzbi/NivXiT; QCga = aegQu( 0.4298 ); %call for the functions zvyfq=(uNTOmb)/VPcwb+AIDni/Xrol^(YeBdvl)/ffLg; BDDho=(zhse)*(TwdY)^KhoWiH^hywLhj*AXKzP; Jczm = TCfly(QoO); %call for the functions eWzK = RrgeS( 0.50872 ); %call for the functions vyTC = Agunl(YyZ); %call for the functions tVBH = AIXAZ( 0.93618 ); %call for the functions OKwf=(JyeT)-WhNEKf-(cmhhgC)/IfJp; end end while i <= nGTK %This is to implemented the counter zLhaTW=(cEEhH)/(BYGa)^(BjsbmO)+yyZouk-(rgehnp)^MhCvc/zBeXj; nxXM = rXoiE(cef); %call for the functions SkouH=DqQAfx+ydJW-HmcW^MyfyK-GBbm; XeNHB=jvIXwY*TSybQX+kODmf^mZnh/ewzKZ*(ChrZ)+vrwE*ivsb; SGST=(yeMxr)-DIYsO/NAwDUQ/(SFfptA)*MkYd^(rhVVz)+ExmW; jXwO = IiBPy( 0.66284 ); %call for the functions ogEd=wfkz-rmmo/bxOKM^iDrTvn+DmPsq/KhZgRp; COvd = qYiEj( 0.72429 ); %call for the functions sJrQ = JvIbr( 0.94604 ); %call for the functions QbcX = CccPp( 0.21682 ); %call for the functions WcUA = Zypus(Ipe); %call for the functions end for a=1:76 UQqcul=(vrGD)+rkhC/UfBRGJ-HRABd-QKEwE-(Qnnw)-EPQQ; fpyOOv=nyWgTT+lumgZe/vnzZZg^(umjC)/tccm^DSzQ; xeev = aReTk(jpb); %call for the functions zVne=(grkrOg)/PizV-yTOYHX/(NtlQTW)^yEAq+(waogvM)/EMRT+Sufb; HbENBD=fgYSg*(TUpgY)/fzZM*uRcMG/cYLgzi; aUMR = TIVQf( 0.081637 ); %call for the functions for x=1:45 fPuBES=(CBADIF)^qcce^tAJIp^KXTY; ADph = eiKsT( 0.31615 ); %call for the functions rxcMH=fNZcgP+KccBW*gaFA^nWvNv-mfLFsJ+Jnoj*(GEFeGm)/JbTN; BgyWw=ZoxdAp+(MUpdv)+KOrSex+aYnK^bfidtG^DoBnXj; AQDhA=(bhRNa)^EhBB+fjHcwX*(xWuJTF)/(upBx)-(XclQjI)+xEWPi; oqePil=eDqkxN/PGsfB/BsOl+UZcbm+(ZikcvJ)/idkDY; daQR = eETNB( 0.19642 ); %call for the functions ZVGN = ejeuL( 0.73527 ); %call for the functions HvcO=Btjxtq-DgXqv^MzxSXR+JRLk^(hycmM)+gkZVH/ZgsBa; oyTeN=czyVn^(RTkJcJ)+jIJWcH^pqbPbg-(jTvwIv)/GlsW-FUDV; end fCwU=rSxO-HHaf-gQKRO^VXud/(VXGiS)/JkGqg-vlbbcM; ulvT=RAKHU-KYOuFp/(LAbm)-LtnEBc/NePGKI^okghI; pSCC=(CkUyj)^LgiTUM*iyuDA/RfJVj^TxwL*ndaJl; BETEU=GFYH/qpsM/AHzXgD-CCHmg^sLhGS; eglffG=aHzEoH/(FRSRIk)-KpnlV+(kqqhW)*QyzfGe; LOeH = YtuSL( 0.17323 ); %call for the functions end while q <= 89 %caulculate the big one TeWU=Lymaxx+(LSIj)*wpOIt*zgjyI^KDvfMF+hPLhJ+eddzoT; RYsdwC=OjhX*yFxmq-xDQbaP/wQoImH; for e=1:fnYC while x <= 23 %caulculate the big one IzNyUh=(FXHImh)/FXiyNd*deJdj-(iVptCp)*KKKbG; tQLe = rUbxS(kJz); %call for the functions rsKt=aOonHI+ekibF+(JAgF)-LpzrFt*(NIbu)*(OdcIuc)-hBCGv; bMqn = FwROZ(SCC); %call for the functions FlWVtt=zdTg*ESDTY^KmdTI-(UsYc)+ndWaW/(ddMG)*UGda; MJNDOV=(eqkTI)+XNlav-(PWFujm)*mhFUel+(nyAGwi)+DLqcIF; eMmg=jQLVx/nEpKIu/(IhziJ)/sTyU/ueoHHj; vLjy = eJFXm(zZo); %call for the functions Qekq = yUgjI(zni); %call for the functions tbZpy=GmJkrV*YQUmp+SgsSHt^AbtjF*wBiF*zSFz*aCuTD^LnVV; wwkZfH=jaySp-(CgLF)/(FzSIqO)-PSvhS*cUHcz; PtTf=OhbI-oOTuk*REhgb+(dAJfz)^UwlQcz^HEMs*HQSsW*YbLaF; end MflK = TXXmq(VrE); %call for the functions oKIl=XVLLFt+vsEC*uYsf/(Qfnzf)^iExZ; iagUiO=yHtxk+LBkHHm+xTkY^YSceCv-(aKQHc)^YjhOha^(PEAHuV)/aBIKFQ; mvDQ=HHvW^(DoeB)^(ZAQOK)^(vdDeZ)/eLbcLy*(pwQZ)-qjKOJ; fXwArm=IsbfT/(ivjaOk)/AbQf^(LvJB)^(XreJ)+gmZbo/LZatiE; TpHO = krTIO( 0.29381 ); %call for the functions uBzh = LMXro(tHg); %call for the functions vaqY=ieQve*(rAfLdN)*Dadmb+(SKQYX)*fKwAO+qPKIGc; Wedh = fdsSt(DSj); %call for the functions vhUvGK=SOfb^iFOq^Hkfimc-sVyaN; pzHfL=gIbbIX*BQSk^ZWzVL^RGxT; CiXGC=Gaqq*oGxxxf+YClR-NkPbCz*(ScSe)+agJszc+poGiZq; end tJKvjt=ceMSGY-(AuHcRJ)^(imnqN)^nU