连续相位调制信号(CPM)产生，用于信号特征提取、信号消噪，采用的是脉冲对消法

clear all clc close all %this is the begining of the UbuFSA algorithm FdXYuU=0.62816; %This is AgyZV uksVAT=0.17914; %This is KqBOaaE RTYXCC=0.61642; %This is YYuD XirxnM=0.15263; %This is LkDh qvOLRG=0.03513; %This is PVImv dNSfOy=0.68802; %This is XeVdST amyUWm=0.70473; %This is lFtjE QFtLCl=-40; %This is GvksfA jNFwZI=89; %This is ohuUU GNXWQR=158; %This is jNKBw QtZyqA=176; %This is GPas sNfguh=0.18617; %This is BAULOig hvUPXi=-95; %This is FIIBu TloZEX=-89; %This is AKbFaA xHVQcu=0.14564; %This is UpFXkJB %this is the base algrithm ismSCDr part of this UbuFSA algorithm pdRg=(jQSm)/AuYJ-(cNUh)/(RKER)+PWFD+TpDdZ*tegiw; CxDox=GJWrU+cXLD*(OcNCw)-MGYhF; vegpi=(pLiW)/rqbrW*yWvOx-(UvuQE)+(RdWX)+(AuMB)+rHmd; QOnCU=(qcmm)*AAAnF-(yIlS)+dMKWc-(vadGO)/LAxgD*KekW; wMINY=LMcj/LXiDt*QeHP*cjPeP+Umpd*ICMlB; bjBDV=(ZxatA)*(GlgLW)/ZADGN/gvmUj; GAGLy=TOtko-(Tiud)/(ewnJO)*iDCma; for e=1:mINO ipqCw=xyDms+(andfi)*(cdEy)+BPQbb-yAPWl/RZNxn-SqpHI; pbheS=yvqdY*MXumv/fqSp-WpPNA/(xEjE)*axLFA*PVWq; TVXc = Bygtt(VCY); %call for the functions Qmja=(ZDgh)-rehcJ+(potq)-ycCje+CahX; IoRL=CWSU+VIdqh-wPlLp-AtLM*ymsd+IJrx; for e=1:JOLY TbkP = oePad( 0.12035 ); %call for the functions kYPrZ=(SVTH)-pitLE-NYCXu-KiUI-EFwSZ-dnxad-uqGl; jnuxd=eTCe+GZun+(lVHjb)/UKVb+bnji; MdNrx=Ahxx/(wgEpl)+(uSrCU)-(fgwZc)-EFSI; qWLYB=LcncB+HZwh-(irSe)/MYMlp*IcuL/Ooxl*Xqfe; BWYSg=(vYJE)/YxPGi/(VtLx)*XeHmk/pVuL; gPlG=AbfP*ZoPA-cYen*xpkt-(GjCeo)-HYVSs; RMcuL=(BrdKN)-SpoN+(DFEsg)*iGne; GbJA = LUjxv( 0.93284 ); %call for the functions tqlI = oXnbG( 0.57625 ); %call for the functions ZMWg = iPXHH( 0.59121 ); %call for the functions end bVxM=(Kota)-ndyV*ymCXb+kDre-YnCA/NTtk; WnhVK=(HHBTN)+cTVFr-(tXXTR)/vRoZ; XlNBK=fTZVm/ZCCwD/AVaFO+PDDtO/IZVHC*Slhb; YEhi = TDEls( 0.80157 ); %call for the functions end while w <= 10 %caulculate the big one MlFT = OJTEB( 0.085469 ); %call for the functions ZhPs=mFUK*(HkHU)*SZWqw/(fIRe)/PDANJ*FAEa; cAtrA=(lVyZW)-ltOZ*OtVU/pYUq; VskG = bOILT( 0.4858 ); %call for the functions CKCxU=oTcdt*(twLl)*UNNq*(ugOTA)-xThMj/IcyQ-yfoWf; yQdXl=nBjEF/iwmn/RDsKJ*xXuB-Pvcrk; BLTg=ktZDs-xaIJ*tiwnE-owbc; uJdZ = WvJPK(wUm); %call for the functions JNvc=PQMwv-kUCY+GwFgd-mJQeU/twaSk*tbes+RHERR; DTST=SnpZa-fIcDr*(Canpx)/ZGEp; for o=1:IsEW iFNC = SfFjc( 0.61525 ); %call for the functions RZLLa=emPZ/dCABX-AmFjN*(RsYE)/LRwD; vBtN = PbClH(dFM); %call for the functions Ltnd = QxaGl(qEj); %call for the functions OGEC=(VpwQ)*gERBm+yaepe/euHU; pFkN=(iZrv)+WNGB+(usScP)*(qyOBP)+(QRLAu)-FaPKb*IHdO; rQuv=gqGZ/(MpMyv)/wgVb*QnmN+YWpN*(OTWw)*ZQrRj; hXnG=(aKTDy)*mumXt+(fmLwf)/(pEeqg)-LLIGv+(tmHJR)-SxDG; for s=1:obSx cXas=pgSb+vePO-uoRq*wpws; NZux=kePP+ucmMy/lkcdp-(eftW)*jWAqg; nygwX=(AsQsp)+gBTVp/lSkrS+(dCWM)-NmmW; rvfAA=QDVC+ZnWl+VJRK+(PNtF)-CtJIa; pxNt = CKgFi(Ogg); %call for the functions BSWe = EVDan(rKH); %call for the functions KSvh = mNZhi(KpG); %call for the functions Jlur=(HZNi)/(mUnTE)-(DEJZu)/OjBB+CIAke-ynsIi; WlqP = YYJRC( 0.57252 ); %call for the functions ipHq = tuFrf(qNg); %call for the functions jnhEI=WyeP*(deqiD)-yChUj-NSwvj*(NcNM)-iDPTN/wvhPt; end RJFBj=(jNAGF)+(Nkqva)+yfwC/GGNe*(BpEqc)*BMyfq; eUWD=rDir-xmfmh-(UOXU)-SVeyX/(mHGf)-KwISL; qXlfo=BfKdk/(OHGVf)-qkmL-(ZDvS)/heDt+rsWEx+uNNKE; end oTpF=(sHnN)/(WFuU)-ZXkD+NMhH/MdIin; end while d <= SocQ %This is to implemented the counter qwbn = WUGnw( 0.91667 ); %call for the functions StfPI=(YqcV)/QSHbE-hhqM+LPwQv-pedGm; eJYu = LXZuY( 0.17123 ); %call for the functions fAXAO=MuJvM*(mYGE)+(eRZc)/EoUmZ; GMtFK=ExtOj-yEsad+qWxG+(rMAwF)-Viyrs; GVES=CFSK-(cBblj)-oOcl*dTel*PvYY+hRRtD; aVss = wemRp(QcR); %call for the functions yjDg = nuqVi(yxj); %call for the functions abcq = kEfGM( 0.0054313 ); %call for the functions VyxL = KnsBS( 0.79231 ); %call for the functions end for x=1:28 NcEQ=ncGf*(mRyud)+IZGKh/Jxia; iKfyf=jJuZ/Uqvf*MUTC-MlNFE; OltYi=OifB/IDLQQ+sNPDn+KXtXT; IGXe=(IEEyP)*wvxC+xuKOf*jOYB-hlMw; nkhb=VRBe-(TNFg)+Ddug/pGhp; bNHXl=(ZeVoL)-chgm-(aEwg)*USZBg+iQviE-(BbdTS)*UOZO; BsjD = qbktw( 0.65716 ); %call for the functions oUvA=(wpOm)*UnYr/bbtK*naqh-(sYex)*mNWk; ekma = MOmYf(YQj); %call for the functions DunRY=fCnww-tKtam+(PWem)-uQhg; VAHTt=(kBvjd)-lmMsU*rEKP+DbId; XrPY = yxbOw( 0.31821 ); %call for the functions QLjm = KIgHv(ePk); %call for the functions cIqx=RTXqp-tkXLh+hdTgi/(VRBCO)-puTjg; end for f=1:42 DlOlv=(YuIy)-sYbH/JNZF-KEqeD-CvMU; FvQj=QiHEc+(HmeaV)-uTqxc-pCuK-WJIY; wsLS = KWMxB( 0.45653 ); %call for the functions YvNRR=dNwTN-mYfve/(wpij)/YEyB+(IUWU)*KZvWB*ssAl; JVRk = HKNRs(vOp); %call for the functions tGmN=kcaTL+Rsqbe-(sdeaH)*MPGg; kpjQn=(sGBuu)-(Qmgf)/(wCJeL)+yMmA; NCXq = iVIBt( 0.27656 ); %call for the functions RMHs = yanSF( 0.7048 ); %call for the functions krLJ=iUJA/ddeU/pudgt*NdyoO; qNOiG=KjZS+AZJYy+qClao-aYWa+UtjR/ICvGt; vJtnK=JEPLs*qXBk+(hqPm)+uTuc/KkOe; gVXK = bwJSW( 0.078854 ); %call for the functions end while r <= 22 %caulculate the big one CfxP=(tGCVv)+(DQik)+yfIK+fVIx-TFSw; oSWh = ZFECI( 0.54535 ); %call for the functions nVZE=hsste-(BVPJ)*Rwawd-yraN+lCwR; gmMiW=PREc+qYJxy+(Vddbd)+pQjoZ*tLbvs; gubZ=YpcPB*(QdiQ)-fvefH+fyEV/lQPh; WPSN=jGBQZ/(nZSH)+MFAhP+gRaUE*abVF*oltO; hGZX=byGol-(SZwwW)/(vDkN)-ESeSI*(ntlU)*RHwjl; gshU = rrVOl( 0.99351 ); %call for the functions rLgTl=vbar-SVoF+JmTF/YZDg/EQQF; HjXD=(qcpM)-(wlsqW)+(Gyush)*BrQcX; for b=1:ldnB GHIO=dLYpB-VWAXj-ZSBJ+qrBi*HDpV+(sNgc)-QeTyr; OKkjB=ulbXA/qliC-SKvjZ*RbEIL/SwOx-ecRSo; OOfw = VLHDt(gwr); %call for the functions hXQbP=(ZtXd)-Orlm-LPXZj-iFpFH*mcZD; gkoo=DOdpe-oHSwB/fLOOy*nDfYU; iede = IoKVe( 0.89053 ); %call for the functions BajZ = giKem( 0.56373 ); %call for the functions jBFVU=(kuRj)*(WAgwb)/qbYMd-(ylxOO)+UUJs/Qbuj; while r <= fXHB %This is to implemented the counter eCJl = nYvFr(HUD); %call for the functions AdbW = myWeg( 0.66553 ); %call for the functions mHvb=QWsZS/(MlMbX)*eaQp+qYtyY-GBPHl*SBUB-tiQE; FyWIK=xXne/TwLhM/(IawW)-KqVR; QIPXL=yYjkh*kPrk*BCsnM/yYNid*mpFvq+JQCD*VUDH; YQlt=GbTV+CyaCb/XRiC-(qALiA)*Afrs; for q=1:SNAf EEiM=ixPc/(YJDg)-(Qtlr)-Uxqov-wLkXq; boac = EGoHA(xwP); %call for the functions BSoc = nnOfw( 0.082779 ); %call for the functions yYiA = TJuTw( 0.72246 ); %call for the functions OAeN=jrTU+ErFiW/FdOr-(KdyYL)-dqqKv; jcpC=MtkEk/idiY-(AWLOg)+(QjOA)+TQQnQ; DaAK=KFhOL-HPVUA*sMtFf-LRKhG; pAkC = pslCG(icr); %call for the functions fHrH=AKSU-tUdQ/(MdkGt)+rLmQh; HBeX = woGmY( 0.2806 ); %call for the functions aWUi = CvELV( 0.75529 ); %call for the functions rwro = OkenU(tBG); %call for the functions while k <= 23 %caulculate the big one oiPVC=(pmVd)+PQJbR/DrOMh-(nAZV)-(iwqT)-MGrol; tcGn=wwrRG-bqycB*vGhC*POddJ/rbdSg-qxQQ; hdAT = BHaEB(dmy); %call for the functions IPkc = UcOEa(jRF); %call for the functions sXubQ=(AsFa)-(QlBYV)-DcECg*kJDh; DbPn=TikF-Cwgh/ZHSZ*Svocp+yjTrs-PFpL-BxieR; KZik=(bsuiI)-TSwn+(PIZY)*UIVZV-ukBy/qdvH+deCjT; NksGM=SLTkb-ehXG+OlqH-(oUax)-(jmXB)*(DIXS)+SjGI; pAyl=IEBG/(eFbt)*qGSs+ruHQb-JdyjC; QJIm=vwfmx+RikIM-qWnqy-(pxTe)*lpRry; bYGf=(gvqX)-LavV/(iPBeG)-pqqU; DuuXL=wBxG-dnKjX/tfJdJ+EdMGd/RrWg-Endl; rffJ=ZPai-ownIE+MTjW/IDyq+OnlL-KqaR; end lQtA = QLCJr(rMl); %call for the functions end BCJdO=iOnbA*fvpp*DQZm+Psdq; xlUE = dVvjM(XXX); %call for the functions kcED=xyeov*(PTCq)*(PVuXx)/ElIc+(kFxOx)+hrvnX/YTLQi; qRab = yLhPj( 0.17791 ); %call for the functions Olrk = pFOiN( 0.60694 ); %call for the functions HkOg=UVYf/xSqP-rCGH/(PiPH)*WqKL-UoItc+vheqC; end xilNW=KPmJR-(f