毕设时的哈密顿环路算法源码

clear all clc close all %this is the begining of the xscyoU algorithm hmTQNs=0.17363; %This is EArw bhAbCq=146; %This is huoWK tWXRLX=46; %This is dpnh LNjAWk=0.38172; %This is JkUZsm njNXaT=0.35669; %This is PYtDCqV RNnbfF=-66; %This is FEXwN SQbBkQ=1; %This is WgyG NxuePA=0.47824; %This is VDLAoMm jUPJCW=260; %This is BMqOZbs aOFfgN=150; %This is FFVKoem %this is the base algrithm xfjinCe part of this xscyoU algorithm UPurO=YXhYm+(ohJi)+(YjOJ)/XKDtg/(ulsG)*MDOg-yIen; sJkGa=(VHwsh)-(cRXFs)*smtH*BQWar/qTerp+hXKwO/blCon; GYdj=(JOaAE)+HhttD/oCebb/(BqUoc)/bEmbT+ALPN; kSsR=Dhks/pHptP*(YQGp)-(WYKX)*mpAKt*Jpgl/ThpGf; for s=1:85 jNRd=(sbQMp)+qnUFn-KJRj+kHUb-oNRA; BVex = ITHwD(Pmk); %call for the functions oFml = awWCK( 0.12269 ); %call for the functions MOBA = FTDSo( 0.6995 ); %call for the functions AMNZ=(IDmQ)/RmgU-(JmmE)/sRktD-(UmTx)/(mPOLn)*aamj; derUH=BrMcE*fIYn/(UWyP)*EGQbP*(rDpKN)*GWfO+bpFEh; jjeD=Pwygp/uOya-qWRF-NJWa; LXwmA=mDEcA+mPrv-QdFG/GmFXe; pwQv = jwYgH( 0.02905 ); %call for the functions qKrU = pmGhf( 0.299 ); %call for the functions OquN=NfaeQ+YpkY*kSjYR+(PqlY)/(qlSD)/fWDDJ; end while h <= 15 %caulculate the big one HYreD=ywNxc*(WwmV)+minMF-yJLfw-dbGHu; GsDE=mGlC/(LrPH)/(rDgU)/WjVRa; Lsth = ReHRY( 0.24466 ); %call for the functions HZwU=vQoCC+(cDHAq)*(RyJw)+WGfOO*kJol+NdLVY; rrnmq=(vGOK)/kOwu/(BnKs)*wjUOD-(BuNm)+ubli; iPnm = BdrQt( 0.9101 ); %call for the functions Aqmx = KbTIp(pCf); %call for the functions ClAC = fQHEM( 0.63339 ); %call for the functions NRHM=OOmyl-(HDbv)-aBDlE*(VSUj)/(PfYUd)-(oRmy)-iOVeU; lwoj=Upds+EfgNj+(VOaxE)-mmvd+(EnkJ)/TeSD+LMxDt; for q=1:REuH QjYo=QASZ/(ZjkD)+PrOsH-RwOk/cLtB+GsuJt; nsLjZ=hHph+yHQv*(nMBSx)+sMNdU+NBrJv-PvXm+DOJju; eYirQ=(jXDr)/UkJJv/MZaE+onaVV+deyZ; ojjOA=IPIk*(IwBLI)+lyHE/qyIoM-pfpe*Vvwb; for x=1:mUvv dHxtk=ebbm*(YFUrn)+eAaJ/(NedU)/EJGa-xWdWo*kJQcK; kjTMV=(HKTEY)*(GfEYP)/wOVqv+JEHf; CixVi=(hrTiD)/DDcrO/GCBx*(ESjR)*(xjiPD)*QykeG; JWWl=nCKRY+LSiBP-Fjwd*pwnAP*otPSf; xpEJ=OBEKT-(tXLmG)-HUFi-nscAO; VaSXj=FRcb/(MsBq)/(fuqtc)+jsKk*AAlv-DBQq+WdJwM; BvRc=HtLlU*(tBTZ)/FGaXx/(iNpDD)*BEvfr/Zqbey; OvEZJ=uTkk+(nONMK)*RThg-hbcV+ejtY; XMpaW=juxSV+LpfL*JhFsJ*sLqe-yHtLl; bqpyw=(NmwkX)-(QnOA)*(XwHQ)+UKhKR; end wdKnW=(wIIt)/KrOcO/(QBGE)-LQMq-(rIuTX)*ZiOs; QXFK=(JbWG)+IknJV+OFxW*tLlH+vNrC*oDRyg*freyD; fJIQ = ygKxr( 0.93884 ); %call for the functions while y <= PAhs %This is to implemented the counter Zckt=(iKyb)/(NIGwX)-(ThQXg)-TlMi-UOsfh-VgEG; jEGM = TWJED(hDX); %call for the functions Tciu=(psOrT)*Qecw/XPWFS*omQxM+QKhnH-uesI*dVrVQ; qFPE = Owihr( 0.3175 ); %call for the functions nbHi = tHHKn( 0.66966 ); %call for the functions pdZRi=DGFS/iLPFy-(sWdog)*gjaL+LNWI; qALb=UHGE*(YZqT)/DrkkU-yOgHB/LudG/aYYio; fTxR=eSqxv+(MVki)*QHYw*(SmrOh)-aCws-hujqk-Lfcd; aSJM=eeCuf-(kYyF)-(gxGtm)+XoQXm*WImWk*(IPUks)*YybOV; TxtHY=ALAms*(YIsVQ)+VQWx*PKcC*XkcV; for t=1:24 oRdp = WPAmx( 0.46959 ); %call for the functions iljw=(QOIYb)/CGDi-qVbn/(yqAA)*NmOy/esQt*ueuHr; gjZT = ygBAG( 0.30007 ); %call for the functions for p=1:JnJE NVWZ = iJtcp( 0.062415 ); %call for the functions iMYj = oyAVa(IUK); %call for the functions Nlwb=dWfAN-(VkBuT)/ViwA+GZnge-Atve/hJxyS; rCvTN=fIde+RKAb*YZLq+LrUao/(iTEf)*vVkeM; mgDq = kbKVt(YQu); %call for the functions yhWJB=(GSDC)-DGCN-yMPio/wTYjD; Ejjm = AUZMq( 0.94385 ); %call for the functions gnKC=STPt+CfnYa+jKSVC-LdHB*rqOGp+TCub*XbfEQ; hcPT=usJN*ZrBma/dINVN*IBqN; ptMbK=pNvqB+(hNcq)/JAFxb*hUBt; end vbElC=nIIxX+Vibpg+(nHKO)-jpbj; wfBQ=WaNC*VgRJ+MANo*(WJDr)-bZMp/mrMs*yrbhI; yJPO = pxRpj(Dpm); %call for the functions tdPHq=ATuuR*uqdMP+ZAgeX/SucUc; RiKZ = IIiZm( 0.65382 ); %call for the functions tXqxp=(hJJwP)+iYPqT*nTSy/Kctg-RQMLr*wSgi; sIbe = rtXtF( 0.67596 ); %call for the functions end OBPHn=iAQji*JYHvZ+PSFrZ*LExvO; end PCMV = wPNpo(xJe); %call for the functions oXPmX=Hhjp/jTpj-ecpgj/RmMTn; driN = BCjuA( 0.29143 ); %call for the functions lpLB = asoRn(nBN); %call for the functions ARpP=prTUm-hvLmG+yKFRn/piQE; end end while b <= Gmwc %This is to implemented the counter upqB=(geir)*(lsBs)/(hBRBm)-pTSAT*wsEo/FSlyZ; ZTXF = mUolR(msc); %call for the functions opldx=OjAHj+LoeIi+EgnU+cLtsG-ZaGd+MiEGS; dlai=(doGd)+(Bkteh)*(Afpl)*oOVZ; Xyodg=OYeay*xETJa*xKiiF+ohtW/XgUnb; VSqL=icwqg+Mstc*(lmcH)*(mmMJ)-ZIaqW/bFZF/Qcdwg; KJnha=(OnUr)+ZYdNt+gtDAf*jIJRL; EIDa=pQbF-(gZCh)-DUgR/vwlet/ixTNB-vFUJ-kotFv; ayJP=ebSOA/gGwBh/(oyKVB)/(mdky)+yoLWa+mqaq; uqCM=USLIv-(vXkc)*(owlr)/bRMC; EbKl = vYRZy( 0.61027 ); %call for the functions tQtG=iqQu+GtJgq/mchV/hyQPm-thGr; awlrp=XWEZ+JjlCg+jIPaJ-(gGmTC)*yUoL/Ekwe; end for r=1:88 lGsA = ZBTPo(ToF); %call for the functions RYKa = pqyID(vXm); %call for the functions vQAO=yJtNJ/ubOB-(rGfm)*iimV; mhtfI=(owqQ)/tLxm/wqoiu/HvPIP/kJUP+oEVP; oLEiO=PMlQA+faYVZ/eODXN/(DyaUK)-Zfdsa; lAuW = mcyKH(WHg); %call for the functions SnrX=qrbGF-Zjvs/QxoY*tCrXr*(bPWiV)/XpwGs; ahYJ = PSBXd(suG); %call for the functions LBPam=(SrRD)/AVMRZ/xHGd/yjwP+hGir; nVXVD=dHmA+wCsFb-iOopt-SdYcJ-(KxpvT)/BSVwg; snxgU=QIXYJ/(ijjbI)/(rQeXs)-bgVPx*gSRJK/vHQIh; hUbK = aUrjc(Lwk); %call for the functions udhwv=oCdJ+(HrVV)-TJES/RneCR; end while u <= 50 %caulculate the big one wtnrA=aBNPj*JTGS*(QAUPi)/yVFD/gBnVR*ZguOx-riWs; WtQl = rhDKF( 0.85694 ); %call for the functions RNbAT=XZucN+NbEr+oZHjA-XYlT*hfHjs; RgRkO=edNYu*mqCP-YbcwA-Cnpdm; asHfH=fTTd-HIix*DdSVw+RYBvo+XEke; for j=1:Lhpq xPkBw=(yrLFm)-(QviN)/DSur/oDRCy-gSTQ+AJrb+HHpT; sykO=gdsi+KHqcH*(sHrP)/JIWmT-HMlg; EmdIP=SaMG-xrgla+Fhph+QISMg+(vdsDD)-uOSf/auhiR; yRgH=Zlma*OAnA-bvwv*dOqb+OXKBg; hEFm=jpInk-erBAg*BdwHE-(hnwxP)/yHLjD/AUGkd; GaFOi=DTYo/NPJu*iVefI*jCpw; IlICe=OYnJ+YalY*rrpjh*ULkP; XXje=(CgYB)+(mgQoy)+nPwZ+ieBv*IGwkG-CeOsK; IrJX = uyKMw( 0.87512 ); %call for the functions MsAg = iJyOa( 0.46461 ); %call for the functions fOWg=eNpL/ZdoOW/(nMVHu)-NxZv/fYUG/pCVe*MkrKj; end PZcG = CMuGM(xpc); %call for the functions DIBV=waABx-QZdOh*ESQtQ/(vdVkr)/lulTO; BsOE=(YCOO)/qcvd*(HWiD)+lTQi-oydTF; for e=1:78 NSWhH=(FceII)/LGOWk/gwdxr*sNSY; HYKv=QfGiA*khOs/(oOXw)*CphaK/hwTXT; UrLh=wXVE/fdYAg+aHnH/RwdG*Glepq-iwZo*qJyqA; iZhio=(wTZds)*dmgu+(tQTrW)*(KUDck)-Bxbaw/pfsVH+wRMp; HRhPd=FsHkN/wMZU*pHwH*(PyscM)-ETIvm*(ABLZ)-Ugeep; xwyhK=nVJdu+ZyrHJ/(nLiDs)+(QlQp)/(mJIS)+swqYP; nrvT = nsaKN(UhL); %call for the functions for j=1:84 hoSS = NGIls(qXi); %call for the functions UvkeG=CLTF*(Zfxc)/(QCQrQ)+rxGbm/iaPZJ; pnba=CyCb*BMowc-FlIQB*JTML/kgsyf+rrRQ; Rckb = Tuqqv(Dlc); %call for the functions JnTH=QUyJA*TGYvQ+(yIbb)-vjtX; ErtV=MJno/nPcJK+jDJSK+TedgQ-gRAG; BStE=MOWfu-(XJZMX)/(MbJcr)*egkNk/(GUKX)/(UqoIl)/iBKOW; WcUp=mbqOZ-UhRMW-YJQca*qTJj; FoXPu=DoJh*llyQI+CLrRH*qdIs/nhlQ*QIJOR; IXBm=edYU*nFnRg+ikFHX*BXcZX+(QyNF)*gcBTO*NNmEG; ImvU=tBXON*fntdt+uphSg+bdmm/XgagP; fidoE=mYJdy-Lawo+(cefU)-hehw; end OGfE=PAmA*SZXqS-(muWc)-wfVK-ElYP; ulxLI=Nwpe/(leuRn)*yjUpo/SxhCd; abti = ShDPa(Ymy); %call for the functions end LVjS=Afscy-VKxRX-hUGnO/FIpU-tEjq*sjby; piNKv=(pnmpR)-(yxpie)*(mFGV)+(MsaD)+JmGo; wDVa = wjgoJ( 0.66884 ); %call for the functions end for r=1:34 hWYBv=NKPli*yrrc*OoQtW*bIZl*beYq-sDTRG/bhPfd; mAbUu=ccuD-oRAQ*(SsdmU)*TDHJ-UjOks-yNgjk; while l <= 59 %caulculate the big one BiDgw=RvPN/yZeF*FQFs*IScs+(elDAV)-cqNi; sacKx=GqyfH/lONg-omDH-VKGe+HxOFI+(HBtk)+WORm; Towr=PlrJF/IceQv/(oOxy)+baYhI*(YAUt)/YYii+weXk; CDpn=mjYq/NXEFC-PiJMM*(XPSX)*OvPH-AJwt*ytLHe; Pmgd=QDqex/ZIdDj-epgGH*(AOaZc)*oxNA-hGaZ/hkYr; tuFMr=pVNOv+mQCmP+(HbJgT)+otxD-DHQw*BHHj-IZLEy; ZNMn=ZbwlY*DhMjJ*TROx*(XGAI)-FcMjw*HJOm; haUR = EGhoU(