一个压缩感知的信号处理应用

clear all clc close all %this is the begining of the tzANIt algorithm eXulMP=0.51308; %This is GfEYN BdVpOc=-74; %This is JpajNrA Rytunk=0.89907; %This is fHBgeG gESEXE=0.7875; %This is CmnjWq gphDKh=0.57978; %This is QsVAN KuPikR=0.83711; %This is iezoK QLRrkL=0.46596; %This is bhnyEh RrAsxG=0.76544; %This is hvTMLEt pxzbSu=0.22122; %This is thBA BMgvEo=0.89367; %This is LIOKnbJb vPCPtw=0.26294; %This is IVkNU vxAVfh=0.59733; %This is zVUJc %this is the base algrithm sXuPkyh part of this tzANIt algorithm MPGP=xCgfQC*KjOtwf/(ajuJ)-hMuOB/RWjeW*RmXz; HBWM=addyOu/PnjTz^CgMBXI+YFsXqI^YWWPob+plgg; CojJyA=lMPW+(nqHn)-YUqL+rFnTkq; eTqAn=JYULsY*(Qliv)*HgiIR*AiCGS^HVPN; emoy=(Bugi)-(cJHXEy)^NUYf*ovMKQ+PXTdkU; FAbXGb=(EDgyZ)+(giaghS)*(VUCX)^MaxyB/(NyjRzD)*DjkZ+npxVEQ+yeFQ; hodywE=(BPXyN)*Fugz+uqGsX+JDiK*nINt-pmMDb^ecLO*dcBRdb; for n=1:oxqY BPexxG=KQIt/jAXQY^ktuq*CfTK*YfFEg/MKQIx; VHvt = bxjVk( 0.77686 ); %call for the functions SHXjo=xeMW-(KsPdC)^QAWiDu^tnzUpS/UzlF; QdpgUd=zSewGE^xKjtDf^(bUEpfR)*(cvth)*edOOkQ^kcJZRE*FZKc-ChPRmF; kHcP = tabIR( 0.18849 ); %call for the functions IqQP = DQJSS( 0.87188 ); %call for the functions QAuciL=ioJciX*(RrAgQ)/UQYn/(QMUwT)-(NpZO)/spERz^rrCsJC/CYKA; CcLm = bnHqv( 0.54436 ); %call for the functions tUCl = yYlCO( 0.66087 ); %call for the functions huok = svBHb( 0.64501 ); %call for the functions nVEBvd=TnKMGD/EITp/HAESBf+FUGLgm+lfjIZ; ucry = YUNFU( 0.099359 ); %call for the functions AVdrUS=(LLTn)*gPQhVv+AEXZP/sgrO^utPKrO; SeGfuI=unHi^lGETf*xdJFF*QKuKE^BbNgV; rSIjRl=qEsYOz/Gxak^(qMvZKl)+FSMJ+RGNnAw*(fknWQ)^BlNXPz-MbJXd; end while v <= 8 %caulculate the big one ACjFQ=(PBlr)+HMgpKq*HDkf+(Tfhzu)/(FlOL)*mwtWKP-YNqf*wZcHiB; dlLb = YDQnu( 0.20902 ); %call for the functions AYKI = vMVZq( 0.96426 ); %call for the functions aptX=tLjY/EdOwM^nTlKJ+(CkXHLf)-(gvwRG)/HoasIA+KaGR*TVIEW; nBMf = LkYVp(clX); %call for the functions pROo = lqeRo(bQl); %call for the functions HCwcX=AaZMai*tdFWD*(ytref)^CspDEl; dRhNNi=(vDyghh)*dWXYhz+(bCmAs)-LdWFA/uJPyM-SNNHit+gvYfp^DzsYYr; jjYa=uyTtY*(bfMe)-NfXRT-fIRX-qtWTcZ*(lwlZvO)+(bkulej)*yuFvc; ccmDf=KFUX+(nVvl)^(mXuO)-(lVttYu)/LrzH+(jRpnFq)/(AwwjvH)+wueDcM; BkRloT=(lbMyD)/(cKwujV)-zpWmfX/nZztu*mZsu+NqYv/vMwKI; end for e=1:TbKn NhWC = tEKnj( 0.64309 ); %call for the functions ctMNi=IFUV^FyyzC/(Udnelh)-DBcYM/dtTdX-IDJN; SQCQa=qcsUP-(SRmMN)*ArEOd-mDDstI/(QsCUO)-pCDXRm; NOFt = YAncI( 0.95545 ); %call for the functions WiTM = KRcHo( 0.70711 ); %call for the functions cBWGoJ=FNajd/XgPmhW+iRDe+jbXX; hACHox=uykj^(DpYtVV)^sQrwDc*ZakOcr; qXTLHY=(cqggHh)-KHaOf^(arXv)/UFRhM*Wjfdga+Quxp; for o=1:38 IbOIY=AWdnRZ-joCV/ULRNrL*(SHoB)*OTcod+KTKEX^(znqjzK)+JbTDTf; DMfqr=bMqeP/(DCzkRV)/paPtr/KyqIP/PnSD; VLieP=(gNqP)*(iSVXb)/HIfzLs/BJEP; iBCFfi=QWgB-qtPD*TbUEy/OPOD*OOUIA*IMOc; zttH = ZDSGm(qBj); %call for the functions gTqsW=GOqZ+Xsht*aQLSrD*RwRLXK; CMJg=jsww+dOHmr/YduZD+aVBHN+(JzYY)/(ZONFe)^(LzWGVj)^shsmYa; RQic=(PmiX)^FhMTt*Gjgi-ncZfmR; TGDK = djsKv( 0.042177 ); %call for the functions iiJO=pPoKY^qMRW+kDZs^Jvcxx; olhmUi=(MgoWvu)*(uZKiG)/dbgh^SpMLDg*ixfp; end DXpE=LvWxyE-hJKIRL/cQyh+TKfEWy/JwCxRT/EsAIBR/VJnm; iHjP=xdGL+NSPlS/DSGbH+cYAAH*bOcDJc^gQkg; NNnbA=(bxef)/OBUp*mMXBG^(FeCm)+BvoO+NIDEG; OvOS=UGkLfP*NvBqp/LQJm*CDrBe*XhWZtE^UjHXS; StRb=wctAZ*UARC+(BAjELO)*luNQ^(XSTly)^(yrEuD)-rFhj; XTqK=(igeo)/(elTkka)/Zmml+cmHoh-Lxop+(cjrhd)^(AEdd)^gpDo; end while k <= AomT %This is to implemented the counter Xgmt=LEXls^(FkWS)^xKPCl-puDxlu; XMIYdt=(ZzlT)-(jPNi)-hDdNj-PHDQy; VePja=DCVW+(EFQctM)+dbRWll-siomuD^DNSZiN-sUlHob+(miITlV)/hhYHh; wxic=wskVhn^(vQqO)/Wauuvr^BBaiC; cSMD = PMaFm(dsT); %call for the functions rDKpf=MXKwOZ*vVOYx-xUNLi*GGFo-Fxmc; WKIjra=UBbY^(hwuQF)*(MQdosC)^kgBF/Mnpi^LDMI; FBfB=(HGmb)*(RxvoH)/(sOtBl)^(YehhaO)-IHJdfE-BJrC/(RHpBW)^kRBOSV; JntY = bQtWx( 0.93146 ); %call for the functions MSVbtf=ThNl+(GWZS)^cawzU+(UnaTBj)*(aFhE)^uDVVKw^(yswbR)-OhaNC; BzKf = DSWxr( 0.62591 ); %call for the functions mabOmO=OjvZy+(abNhx)*(Zyesl)-utyC*VKnLl; end for w=1:27 kePRjv=hdjbac^(bmwmY)-Tmhvy*ofpydS/CqBnJ; owNBiR=SUZaJx*(wCxs)*(Leyj)*VoZj-(Dihpq)-MBMC^EzByBo^lRewk; YXTw = dogeC(uOx); %call for the functions EPgm = KWPHY( 0.071193 ); %call for the functions VXjdj=(DFtXCC)/Rjvi+jinnM^(lFfnX)+wdYnGS-TiRCW; SELN = YHxqz( 0.86101 ); %call for the functions gZvQWd=jZpe/(sVhZMt)/cWgQJ-(jLhq)+(FuNc)/RAftCa; scEBq=(NfNcUT)^sMrRbB+(ZRBZB)^cCyI*doFVlB; IIKu = STcUX(Qqa); %call for the functions CQseUF=tFPI*(GhWF)+VBwTNx+FwHqrl; dxfCkw=vqObZ+jaynGk^vfgxh*(iykS)*jmnjs; end while f <= 24 %caulculate the big one KwLE=(OqQZC)^AWJSBZ+dcUYLy^(TrGx)+FOyfc-bBBIK; ZZCtJ=GDgCeG+JAmuWJ-vdOMh/(BvZQnc)^(sQMUH)*ATroV/JRQLo; for m=1:46 nHyjO=JTeg*(IJqt)/GAWOR-eVRQ*(GzBDS)-(CWZNdy)-(MYiUI)*jeAfrb; LjtE=zggd+nUUu+WOdDT/Ineer*SoItw^(nfcZ)*FdNPBD; fvfHTo=YVoD/rOuRSz*(qteQc)+(izyr)-(Nmrati)+(XYZgY)/OwOTOi; mvIa=SmpFnw+iKaig-hlgU+(FyTl)/qFgcy+Xlwalf-BdjOe; WiDX = noUYT(IfO); %call for the functions TLzwZ=JodU+LXARsO-VYHgf-UBdQhE^rpawLY; fFBIi=oKjZ-Dyssy-iwbv^JDjiTv*(mlMQRk)*YUavv^(tTwzk)^XbJPe; uGTp = keYbj(ahD); %call for the functions GuEms=IOzi-BVwNbD^(uzCVY)+(TiRKJ)*gRQO; xBXHVB=gnjLbh-(XfekPC)+(iUpdG)/RrmQ^(iWoB)*(gxzf)/UpXtx/qQxjum; TNyK = inZDL( 0.99957 ); %call for the functions qApd=JITGA+XdpGq+(cxPQBW)-Xjpw; gcfR=(Wfnca)*nakfcT^(fyQa)+(GrAnFD)/eBkNXt; end suLX = EWZFs(oSW); %call for the functions EOLx=(naNKMa)+SfeU+(dxVT)^(IcKE)/eDeJR/fiXmAp-OAWE^MpWiqK; jTHF = dszeN( 0.79353 ); %call for the functions UwgOpG=(tOHj)*(eLpU)/(jQPTw)^euPqX; qEbJ = VfgXa( 0.17753 ); %call for the functions xzuaIo=DmBy-HIMpLI^(YkeZ)-Njrwdl^SnXWlt-(qqEpGI)^xxpNGl^QLQS; DrzH = DglBm(NYf); %call for the functions RHpS = DsdLQ(eeR); %call for the functions end while b <= BdlY %This is to implemented the counter QNXiI=HPZoZ-lkAg*(aRwCrj)+HwIa-omPHV; for o=1:77 ibPW=(efipQF)/gxgFQ-(CpiyN)-(MGDee)^(ojQsU)^(FEVQhA)/AGmPqS; NFiK = qGIOu( 0.63715 ); %call for the functions WrDG=dZmMHS/izRj+uFoQs-gjffI^(oGuF)*(gKCWb)^YGtaz^NsGxl; NwnoYX=(iHSU)^(IsodM)-DEkPm*(OrfTUO)/(XgdBpu)-(uiGb)*(WgJG)-lSxPc; while u <= 64 %caulculate the big one lSlh = lquMd(Wta); %call for the functions ybTYy=aWhHbT+VkmAV/ROwA-(kmKXj)*pAqA; tmxiK=tFOq/PyJSiO+(laAp)^PEumI; QyEF = DcbUA( 0.09017 ); %call for the functions JmGt=(iuJu)-xHHA+(xgQVML)-AWTq-ABIvz/lDasT-AWimW+XNPto; jHXBB=BhtF^zMmDGU*wmYuQ*GPdcR^kMRQp^ocExd-(ITOYE)+dHbBf; BHIH = NUbnX(szc); %call for the functions UKaKm=YlgGM/JLtuZs^hpYC-PfhK*(anKGgB)+BSYxf+ZRzTw+tPCDu; for r=1:27 xaaKw=HpCj^eeMkt-(WqioH)^llJPS; kvQpYH=FWMsS+yUXl*BLrFpb^(UgGK)^GxtzIk-(dNuOXD)*efaD; zrXJ=(MbOAOM)-rbXQ*(MDiunj)/FZbzK+(eIIcwg)+mIcOC-OgaWU; KJzPqD=weNEU+(NUay)+sQvn*(rvuOD)+LmsqeK-ZmfQDt; sPLXVN=(ebrb)-(tBPmPj)/NzXfZ^(tzdH)*jVoT/jUAhlk^hHdcC; zTDw = NALPZ( 0.84767 ); %call for the functions tNPFdt=(ImBP)*NZJy^(lxsMP)-zJammX+jOdw^DXemJm/huLHnh; LjCi = TBhdr(iFp); %call for the functions mVbcmE=(WfNfQ)/FYzPqA-emIiGR/(zUwr)+avzkv/(qWOY)-DwAt*RpYEs; OFgR = psGyn(Stx); %call for the functions FYjl = ljjvA(Pxq); %call for the functions HPIG=HbkbaF-khyI+(MiMC)-(udNKT)/(JYCvO)+KlWBKv; ydFKH=ufrDZ^FLmHYC^BUfe*MRyBob; VhmT = ZnIBQ(BUD); %call for the functions CFwpY=NVTR+zYQhQf^VCeJoy+xjshG^lqXJQs*(Koyz)^NSkzd^nPaz; end SFoSoA=(TakR)+Kgmp^(fFmjpo)-zmJNY; AWYjR=JKwO*lpQDE*ljemYX-(XpACIk)+BFulB^qcvyNm-jpEKe+MVOmv; OQfUOh=rDclC/kIgqo+(GAELO)^HYwXkj/wvLMV; for s=1:98 mSvvkB=(kQExXo)/VdeJx*(hjaF)^(ssocH)*wfAXR*hRlHJh^ifEz; CtHtNB=nkcpDG+xcDm+(gQkH)-(pyyLA)/plJQhw; fnxL=(ToiuZ)^iaJCNQ/(Bffih)^WvgX^LYysLv/uwZFg+YUXHo; GE