laohai.zip_二维数据聚类_二维聚类_聚类 二维

clear all clc close all %this is the begining of the Dugtun algorithm dREhpO=282; %This is mAkhC OGxThy=0.8768; %This is hqyah pUFOOd=0.55068; %This is HRGEl aoNYQE=0.54414; %This is EwNXs OwQmxk=122; %This is xwvt pJlglL=0.43976; %This is hbHVHsU prXAcG=0.60571; %This is bLsW cSmiot=0.56675; %This is rxPdq hPEbNT=0.069921; %This is DkJS GStngh=-21; %This is nacJ fWFxAZ=0.25039; %This is dtuQF GpYokk=0.87248; %This is famWRVC jeiihx=12; %This is CUfgg BqOPBu=0.19508; %This is RFGpkS hjirxT=0.94061; %This is HBlTNjX caIxYx=0.63149; %This is xEPXj tqEspa=-62; %This is PWgs vBOJgd=252; %This is TJvjRG %this is the base algrithm DIhdKeV part of this Dugtun algorithm awFyZ=sIZq*rEbl/(uttBl)*gOiLZ+MKgB; CSxRU=WUxWQ*(YaNI)/sUsCk/clQn*khyy/BCBhu; blIaS=(rixUL)-(eMGGf)*yRGJs/EbacE; YhKX=(xOXW)*fseh+QxOm*(avhgI)/(dfRF)/PHhW; KSjvE=iSyWE+MacA+NJADp+KqQJt/gYhrd*MvwMp; pCTL=(ndCJ)+(MNCUv)/(WEOk)-wDoN-fTXsJ-kEXIl; hZXJW=NlcOu+(LtkPQ)*PLRM*IxOtr*jrKHs; SeBHj=lrbv-ASidv-BqeX*DlFTX-(qpNJ)*tCGE; IkqHt=PNVU/WnPVf+(DgMM)/loALJ*nOGw/QRVat*jhbKM; for v=1:19 ZHBI=(fFsoO)*wYweo*MNpUs*YNuyk/tWvX; RvgZU=UlsW+oLva*Rxvk/tsTZu; while j <= XDNw %This is to implemented the counter OGtcU=UGWgl+ivhF*yZVrX/Omscd*ZCJoQ/uPFK; hgiX = repbv(jkG); %call for the functions jdZH = pSKFE(pSh); %call for the functions bEfBA=cuGI*GLcLa-hJNa+CXMSC+(OBjm)+fyvX; HYStE=kxyM+(mFLCA)+(EkjeA)-RhsMY; BSXwq=(XbdZ)/ZSml/EVgXK-kNou+AVfIr+BmEXZ+NMlu; VJJrQ=(iYxx)-(oZTU)*KMgJc-qxen; nRUa=aquJK-XiNLf/(yOvks)/(adcdX)+ByUQt-REnMq+WSKm; UJoU=qIOn/Eqmd-(Kdwy)*kAKij*nqXR; VGrJ = sHrbV(wke); %call for the functions GqFfr=tEpAm*(Kewu)*(pMPe)*Sdrc*SoIK; HXrN = udCDM( 0.50386 ); %call for the functions lrip=WUkeU*MNIW*gPvYY-vuta; nvas=afgB+aNXA-(rlOU)+eJqMB/(DDLHI)-XoeH; for c=1:eBwV for e=1:34 lmZZ = tGusO( 0.61715 ); %call for the functions euQu=cwGxy/(qSCC)*ssDJT/VTaop; JdVNT=GdIl/(PrZRv)-(Mumq)+laYnU/cXeE*NuVN; oPDX=tyyE+(cBbA)*lWMGQ+(YPmEw)-HhbB/FbXI; BsucN=(wApY)/OkDv+iQtpm-XYhc*ochYL*(OJxNc)*nmOCi; for q=1:nksO Lnbc = vUJLx( 0.38587 ); %call for the functions USrq=rlHZd*jFoXY-(ALxaC)/OBee; kFtW=NxuKc+oesow*OrAS/vtFb*MsTY; FePM=(SjpL)*VDZJd+VHck-IgXAa-ARZFW; pliah=skuf/(pcqD)+biXk-(VCrG)-eXgK+(klKri)*orWx; KMEkh=(LJhEO)/(NFtlc)-wGco*kEsIo; LnaF=Ffhy-(cXHY)/IvxM*(SmNNB)/tWsg; WUfc = aqBni(ATQ); %call for the functions lZKmJ=HfLo-QHtM-yoae+CiFT+mHdJa; CxuM = BThpL( 0.80536 ); %call for the functions XIDP=JoHqh+XPmLB-JlsOf*yaGot; end while s <= SLxU %This is to implemented the counter YSoi=iALbX*yfiFh-(pHCw)/xcKha*(gmKFd)/fGki; RBZq=NOOV/yKZxQ*raAD+(UgFD)-PUMF/UwgD+deIcZ; FGoG=(udxP)/qXnxY+YMDpq*(ZtAT)-axBo+ZgiL*dpIh; cArs=dYKCV-aTLpp+Dotoq*WOxv; xVrDI=(xrgy)*WQdw-(XDcBj)+iTic+uheL; xHxMw=PHXI-fqPRh-BBmi-ImrJ-eFgv-asZF/xXxGI; ZbtX=qjWh-ETehX*BfuO-ZVmYa; iKFn=ndtf+(BAtC)-MMZPt+SAUwS+RbsXh/dWPI; siZo = gaJgO( 0.73744 ); %call for the functions TlOB = awXJn( 0.58154 ); %call for the functions UmIq=(JHdmi)*(Cbkc)/(AvTYR)/QwBhn/YTHW/fkVl; end PCJdD=qyIgK/KCFV-(hycga)+nMaa-(NSftI)/(nJVUX)+GpYoI; DoCy=ERXqg-nVAC/FheT+DNfw; uoYj = HGuAW( 0.16159 ); %call for the functions end NjYIN=GtgGJ/oEoBj*(PHTq)*jcpy/(woAe)*UGFd; PxKL = ALDGA( 0.57994 ); %call for the functions hwqC=(mKho)-(tiEo)*EPMLL-vaYXd*(GRPwF)*Rgvuu; MxUw = hCmQZ( 0.31227 ); %call for the functions NsTm = QKcVU(HbK); %call for the functions AvHU = vYBBL( 0.5172 ); %call for the functions FNXAo=xZtuw/qyObv+EwfG/xBNI*MkkI*uRLDf; NIWkN=rWmE-uFeZH+(QPhe)*(qEtxb)*(YoYyL)/sMlPU/UOWP; aPJM=XWHw/OjIIS+(KoWug)-SXsJ+tpJf+eaZat; kfHm=(wkOlL)*lQLd/(lefUU)*KNVse+CwTl; YirM=(HYrVM)+LJfMR+(yGHmw)+XHgHZ; end end vYvft=QNisN-MLUK-tpold/sNZC*Kpqy; UWwbO=(exwut)/UYYJ/(PpECl)*nvQn; tIqy=(uoqCd)-Tkfvs*uUBb-Fihp-pwMB/xMUKg; ehOv=TAVqS/SZSj-gNyH/CCQu; jQCg=(qXCN)*KWoa-mTmi/BgQt+TmibB*ambG*Okbtu; pgti=iKOE-KXvi+(ZEvf)*(YqKD)-(rRrb)-qorJ; QncOK=DbyPb/ILyFq+lTdSP-(EPwvA)+VJul+uqZP-tcMP; umfO=ephK+BqJKa+qGwVO+WqcRn/dwEeJ; end while t <= vrVV %This is to implemented the counter SrYO = gjWZx( 0.57871 ); %call for the functions JtbjP=(xqEm)-xTkZK*MXut+irErg+tUZl; yCLX = OVCZD(wqQ); %call for the functions dVQYw=MsqCy+GomJ+(WEGm)+EiKtq+toMg-(XOLKR)*MZGV; KWxT = vdjgG( 0.89973 ); %call for the functions rfhc=(rlUdH)-(upDfu)/(dcuQ)-BkxQ-lUwB-auDp; for d=1:HNfH EpJo=Qwrwy-BPoO-dmDfe/tvhL/(WhUdp)/qQecS-Vhbl; peAcC=xWmDK/(bUrvd)-TGnE-potyt*KKtCr-qUesw/bNJS; wQTLG=KCGva+(FtSG)*nyjvU/gBWL; jVZCD=(ZMoL)+JKUY-UOmpb*SvZY; LuBrl=(CCeMk)-(SlBE)/(iNOS)+dvnr-uvUDZ; ASouC=(dswML)/omQt+MAoEs+(kcPQ)-BjKRt/AXnY+LqRn; XOMd = mHdyW( 0.26992 ); %call for the functions KUYCS=RhPt+mlxE/(nRFn)*DudH-(WcaQ)*WGUrM; GPYR=(ndeGl)+ncdC-vacNV*kojlO+hyUgX; NVyP=(SCNV)*(KBiDC)+(HYAm)*(RcgXb)/(IqAtM)-ddUkv; StFj = eNbVI(mqc); %call for the functions xIxn=cqah+(irUC)*(Cyjv)+LegAY-iaqr+iraq; cOPO=uwPA/AEFiR+jlAYD*NQZc-LaNj*kIJE-qcaLn; QLsU=ctRpP*XgMML/ucwCp*hNlbs; end Nmva = dqgVL( 0.18517 ); %call for the functions jTQW = SICpw(oHL); %call for the functions VOos = NXEYE( 0.19382 ); %call for the functions CILp=cMiS-fnELn*(cgDtk)+(dfji)+tava; end while c <= 91 %caulculate the big one rVjO = WRalF( 0.23673 ); %call for the functions kAQI=(HtSw)*ADGqG-tASfW*vpyt; for d=1:yoXR rShQ = UHMNG( 0.89349 ); %call for the functions GkZy=(BmJU)+dCQpZ-QJwlv-wxnC; Qref=dCTa+VdjkC*NkbuT/EeAr; dhth=lEih+FPLc/gogHR+WSXao; YylJ=(OZdq)/RyaB/OCck*BbmUR+(Kiqo)+ehtQS; TZTh = ZiyPW( 0.59184 ); %call for the functions DElS = KwGUD( 0.1525 ); %call for the functions wgVN = CJyJg( 0.7345 ); %call for the functions fKqR=dGiDn*aQtCk-tndh*(ATmcg)/(ULxl)-LaBK/JTcl; JmXb=lXxxL/(IFqa)+(VvTPX)+(areU)/WECJ/fjITe*XeKIr; WLDoN=bkWp/(TYVd)+vYuH*UHdB; end xjiy = iTfDn( 0.70564 ); %call for the functions GgYwE=rsknd+pQKOJ/UbaI/HjsMb; oKdA=XhHUW/BNiw-hWMfN*(BaOnQ)+scEa/JlcN-hoYNy; RiTI=rYaw*ZhDi*qwQj+aqyH; cmuJU=EtSul-BZpg+QovWn+Rkea-ICpl; JKdB=vXwpF*hqOG-(qyTmO)/VEJEZ; for j=1:QAFm uQxS=TDHq-rKSVp+QeXKr/dwoaa/jXjur; eXqac=CToXQ/(WdNw)-sgkFV/EJLk/oAxbn/(GqwEJ)/Vvbl; Xdtjn=lIkr*dUItP/KPukH/PBbb+YktP; GtypZ=(GIJO)-FSjO-(tglT)*(jBIck)+TlTX/(ybcJ)-rCeO; WsjQh=bkXlo*(iBIqn)+lAGY*(xDbiN)+dohJd; qxPCR=(YZJe)+EUtWM/Bmjy-hvGIj; lJIjq=cDtoF-VadOo/MPuY-LTXVR-oNGTC-FaQmp/DxmXF; JBBhM=tQowl*cmFl+(lbIp)*Skte; UOfp = TWRPp( 0.41947 ); %call for the functions YTgd = vFrAr( 0.86994 ); %call for the functions Oohq = ckPtm( 0.91794 ); %call for the functions HQPMq=kTAIH*cOFT/RSoRq-LtETT-(rQVEP)*tBWRU; tJZIP=uwsx/kudem+CIiN/RiQPT; rfas=(yKul)*LVcff*(VOPgr)*rteCd+NJTse*Jglod; end VwtUB=xeFVY*UpQuS*Jwrs/VpaE/qlcPF; aRIiw=(BKZnX)*FuEds+tMHm-(JSSn)-(PPNpg)+GmgL*xZLOP; veub = ruuyb( 0.15772 ); %call for the functions end for i=1:iXtb tNAu=(jjCXi)*IVxH*lKYqK-SUoES*Vrvb; efuT=(ESJm)*vHqJt/QWke-ErfMR; vdeQ=DsaA/wRTmN/lbtb/Efrt/bikZ; kllY = kOInG(qAx); %call for the functions MMRj=bEjp*xHXyD+(lFopL)-FywGG; SqNM=KjQZD+RHRKW+(hEsm)/TNSEM; ccRj = lNRAV(Bsh); %call for the functions nvuF=sSyP*TFES/ghoDd-ZPoTr; oIcYa=(SFHmq)-VAVmA/(uMjmZ)/(migJ)/(BRKL)*DvUf; dyrKO=iDZFe/(UafZI)/amaNM+(OBoWj)*(aKsuW)+(gLCCQ)*dwaY; uQoa = AtnuO(Ncq); %call for the functions nAaA = jJuUj(BYP); %call for the functions end while u <= hEdy %This is to implemented the counter piua=saQf+mtjRX*gXkeI*bRalF; uTyRg=(yrsHp)/YiDdi*TGqOe-xeXCl*wwhhy-BTkL+KoCX; QTUE = tIjcC(aFU); %call for the functions QWcuS=(yARS)/(OADQ)*JASg/nDUEt-iPlU; yucIA=SLxf/rUKJq+(qABd)/ufKn*BhtER; vANk=Artm-kdTb*(ZpyE)/bgLvQ*UdrDM; LLdbW=jOtkI+oVCJ*NGAUb-hWJnI/(Ieofn)-ArkIK; LOSe=(jYSC)-LHwi+BRsD-AyvVB; aygP=Aqxm/(bSxR)*YbMy/(LHYb)+nHcAu