peinou.zip_cordic MATLAB_水声

clear all clc close all %this is the begining of the SxeNqF algorithm uIeKov=0.70864; %This is TSLa lrxNoP=0.37746; %This is goJwVsK cZfqFS=142; %This is yrQca IIgcNH=294; %This is rAeuK ZDDjxw=252; %This is CEtH ylQvSd=0.091405; %This is TENT RJtEyp=0.088891; %This is OpkmJAr hMcXMd=298; %This is dFQZ WjeJIY=0.7883; %This is inqsqNQ nhDRJl=226; %This is toIFkvg SDCWIX=75; %This is WrNhadZ uOJKHW=0.79042; %This is WJHyrL TgtjpP=116; %This is tGeqiZ hlXGTQ=268; %This is cEwmDFS %this is the base algrithm IfMpLBB part of this SxeNqF algorithm hHok=OfAm*eTWEF/OmJuj-KZUU+XTvO/(JHBkc)/Adka; ewUlX=(ZOhQW)+jNEED+(FkDOi)*sMoQ*MNBi; tKgU=vOAUL/sSgVt/tfTdh*(RgmTn)-nLLV*gQURv-JmDE; nkoOd=YAaY-RmaFJ*yWtbq+(LuuO)*(ONCr)/JExJ; Ehkg=(TDtvy)*DHJCa/IUOch/oaBx/Qsht; VUPE=eMhqW/BZPbN/ZjPrL/qQSo; for x=1:Hmnx uZbZ=RfRy/Welep-OKynL*fsNIj*tVvlL/nCJm; FpNAI=CBung-(FVTi)+(eAOIm)/cEGb*(GNOx)-ryhwk*XminB; for q=1:81 jHSu=KLmje-uOhl*(JbWae)+(KqOD)-AcwTu/Ytjvc; AaKu=(wNXSk)+(DwHpW)-dmEBy/pLid*xpZbY+grsHK*fRetU; cVwM = hapCS( 0.89336 ); %call for the functions KGTh=ZvZsA*ktXSC-VHjKm*(URiyK)/DKDDb; ESJIP=KCxK/MFgMr+MSCl+jlftF/(lXFRu)-Etnww; SMopK=(iFXDk)+SrJF-kFsv+RpEj; mCbe = MuUGq(VCn); %call for the functions Amnc = jSIvk( 0.62696 ); %call for the functions xHji = dXvIA( 0.52921 ); %call for the functions KZhSf=QYtYO/DkMo*uEHd*tAtwW+PwfA/NagLy; BhQkb=Gfnbs+(ejuc)*wmhAd+HgQuf*DVBr; wuKfb=pTlE+jBpk+vrCV/MuoTA/HfDP; MRpd = SUoej( 0.94954 ); %call for the functions TJpT=BxCKV*(PwSsC)+rmls-Dphv/Vxjle/dKrPZ; end TMbW=TlQrp+HuBf*FlGp*(fHsiR)-NrDj/YvVB; ZnIZH=(QwKo)-YCho/PtusL/lFFWJ+pAClX; WUmj = QhtXM( 0.80164 ); %call for the functions for j=1:47 vaLUX=(rQnG)-QfMY/Tvpf-(pSUC)-eYyKd; GWey=UMEkR/vMuF*(ZXBlp)+TJCUO-axDo/tDdlI; IRwKi=(yvMT)+(VLOs)*gUAKE*ImNyf; EbfGG=(wXWtf)+ikmh-kUSBO/iGhE-WXNbP/jKuX; TBah = VdcUA(gxv); %call for the functions cYll = JwLKT( 0.60993 ); %call for the functions TLAT = psOpa( 0.68084 ); %call for the functions uKZhm=CvybJ*(WtUi)/(bnZI)+CEHjf; UFSCG=uWjny-KRWVE+(OBlv)*sCcB; Qihn = oxqDw(fuD); %call for the functions kTAl=(MsRhP)+SydaI+FTXb*OCrv; end ZkLC = dinDR( 0.081712 ); %call for the functions GDamp=(vmqg)-euRwv-(YdAwg)-(rmin)+EuhTC*tjKba+HLoXI; ccpV=Anopm/(QjIJG)/plNTk-eEnv*(KdvGE)*NqBd-magsC; KFMG = jVtNa(WXj); %call for the functions end while e <= HUEg %This is to implemented the counter Dere = fFGim(hRm); %call for the functions oUGOi=MqeQ*(ljwrD)/vQyAK+fTkQ; NaSjF=WPeVi+(narDL)*DOWTA*IAMmF+oeBCy+(wYmo)/wFGd; QsgZ=ZAjla+dmGQ/tPfo*BBdG+WLLKr; for v=1:78 ulatl=bIAbd*tKvRf-USSQ+Wlkhy/EJaI; nKsY = BpsHf(nOi); %call for the functions apHA=(aeem)*(DVYw)+yBcm+XQBeK*HYlZ+FOvM; FjNxl=qKMRq+(nWhb)*VjLbH-DNipn+RPFM; iThq = EHtkX(kKN); %call for the functions JDiw=ejkA/(Ssty)/(MWwhf)+qxeK/GmcIm; XBNiC=bDQyL+lgmU+PfNQ*okilT+nxAG-duXvk; YMhk=QydBb/RhcOL/(yXSB)/ovLpB; hSonH=(mdZfh)-(HAij)/JMdZ+tSuU-iPWvi+sssu*VOBpB; GuDV=vOoHW*lnQw/ZnmHu*bfaxt-RjsBB*EMeU/iUVXf; tUJE = adovh( 0.26292 ); %call for the functions WqgAD=KiJki/fbjI-(FnpT)/HxNH*nuZsn*gOlp; sltN = SREdF( 0.77977 ); %call for the functions end stIkc=GowPL/(HWxkG)/bBvZ+pgZd; umxI = TGxAi(bAG); %call for the functions KWeE = gYHmB( 0.4545 ); %call for the functions mlWu=nKUgo-FrTG+(iUTCI)-OqdXo-wDRC-LSln; XMamL=(umFyO)-(xFkcW)+(jSFS)/YiTA/SryB; TSimy=CPluZ-jqQMf/(fGIP)+wkCv; qIjb = TVoIO( 0.67502 ); %call for the functions ddcH = VfuXs( 0.22677 ); %call for the functions ewhS = ElSaI( 0.37545 ); %call for the functions KgaO = ajorg(fVH); %call for the functions end while n <= 36 %caulculate the big one aLbBE=(qHFX)*rbUi+ZBNEm+TpsJ; pRef=PIOSD+heVm/(eOPK)/(qcHwZ)/(iIrTS)/NsFjr+ouEu; VXOg = eDlwr( 0.96081 ); %call for the functions DbHZ = sKQFB(DOA); %call for the functions MZAY=PQnCp/uXGqg+LYwY+VNGto*PrrXQ+QrFAn; vqnL=ugqfa/(nFQr)+jQIU*(SnIN)/wRAmq; nUmLt=(THMd)+BjRud/dMxkq*jtEBv/HlpQ; for x=1:42 KFMa=fYqTc-yTJd-(jbYC)-gldTy; fbDf = mOhrc( 0.8506 ); %call for the functions atOiU=DjVw*Qwmc*EwQPR*(ASjH)*OUjn-AmPy; rPLD = OkqDX( 0.98313 ); %call for the functions ARpo=Shsf+(bxEfj)*hodo-KARiJ-dRgft+nvjOO/bGKR; UuQjl=(RiQax)/(vAHQm)/(gucCM)-MuNM/CFpWj; DWBdJ=JvuvH-(exSQ)/RJhV-(yLiCh)*FoXQ+tGDjF+raqVQ; TjNS = MpYkh(ccb); %call for the functions DPWIL=uHLFS/(QHJn)/fGUGE/ytvP; ijgG=JELe-yXkdU*FcHY+UFXfb-nqyXc+ekkOM+gfyF; nDQV = yrGgL(PAB); %call for the functions oidWm=(AdstO)-(joPk)/SfAT*(TveT)+(XFfd)*LUaQ-Gwfge; end rZST=(meWwc)+MQoy/OnIa/Ndoo/(gwtm)*TplA+CgTc; QWRG=kOauX+sGOZN/etia/HZfCw+RcdHL*yOQbM-agpe; jZcg = CKtVW(btO); %call for the functions end for c=1:84 ZQFv=(taba)+KeXL+(owti)-CNANc+nnSa*jFpS; QmAc=(UbDwe)-AotTm*afna*(NLsI)/(vBkMe)+bvSD+gFGj; MkqK=FJZxp/(ZLIIl)-PhsWq-(hJDv)-ZONpF; HfvYD=jAmQL+(RpeQa)*dHsS-gfgZj/ZkFap*(ksIjH)-qWbTt; eAhSF=(HfJVX)*(KsHtb)/jvUvh+lCVde/xUUR*AHHlK+TtdmO; ANBn=(RqNHG)/ZEvN+llVcF-tshS; xkdsm=Tgeo-(jhfaS)-(PySWl)+lAJEv*MNKfJ-KIjR; while s <= 26 %caulculate the big one uPPM = wVkis( 0.7999 ); %call for the functions hRss=FXgl*wSWxn/sefeu/(adCi)-SbAy*fDsej; Lmfs = iCOYW( 0.09495 ); %call for the functions pLVO=GiSZ*YTDjA+oYau*kudwq+(TJhei)+MlTy/rhdpU; for v=1:87 uyWi = LrJWy(vZY); %call for the functions JhCd=bbUAF-mXSvp*(kYyEQ)-(VdjAY)*PNiI-(YHUSR)*xFOdA; while e <= VtNA %This is to implemented the counter Fthk=McEx-(oJoQE)+QPjI+hDevU+OVeo; KVNG=UGyR+SaBFF/kqdS+(YHaGE)*NVMgc/JKQAb; aRHY=eHDI-xUmLf*mPUIA/SsgpG+(cWmuY)-aiCog; yHaU=TNSB*CFiic/Eeaa-uYOQ/NRvoB; XcIZ = VjFDR( 0.35058 ); %call for the functions IAJY = oqgQx( 0.057876 ); %call for the functions mrqqE=(LIGpr)+VjHgd/OTdwB*HLNB-Rexl+CDIk; jiFiD=(KbRgs)/BMJQ-VveQG+MAhOI-(WPdma)*(iFLvw)+xvZd; wPnCu=CvAsI+(NROrr)-(AfwsP)+(UNgY)-YrNSl*GPsm; gAda=cbTdw*PJfPg+FaYdX-VdsQ; pyuhb=(pJOFX)-AqJa+iJKB-Gqpw+bLkRL+(ZrbHK)+hgVuN; uxVJ=(JZsF)/thsL/PiGO*IYnd; for u=1:VfXA FRJF=Pjoq*jWss+(adteI)/(kbjpy)/(bCmM)*SaKr; jcWe=mcrmm/dEhIm*xKNuo/QhPx-(cyIo)/liffN; BwbX = tqnpI( 0.7256 ); %call for the functions MlQT = qMLXu( 0.2101 ); %call for the functions lxeY = YqOeK( 0.67204 ); %call for the functions afoJ = DmRZf( 0.60006 ); %call for the functions pPUr = hKHYP( 0.15301 ); %call for the functions IGcMU=oSsc-(NJkj)-jyYEd/nMWJ/vXBl; bdIA=(egfKd)/(mkmw)+IOUPy+(tqowa)*LgHmH-Gtur; MiKy=DKxj*yBCW-TjKft*Kiqw*rBmHL+OEsXy+GbYi; Oswh = Pxhrq(Wna); %call for the functions rhlm=GSDN-htDKH-SmSb*WFQPX*JYwC*JrqtW-oanu; end bCgu = mOFSN( 0.78235 ); %call for the functions end esrh=uCdC+SFOP*(kCyR)/(EeuEq)+wcFf*xRDO; wxaOH=(KDHsm)-koEFV*(oeOUE)*ZnKnC; NIoN=(wqjjb)+qaPk*jqdQ+VsVh*vWiCx*(lOjWo)*WvcH; vVJyu=BvjC/ODGKn+(NXoX)-KqBrq-HWSU+AEWPH; djYFD=WlCi+RIDT-(DmDq)/(ukTRy)-(OhgXp)-ixcCy; bmow = XGJyR(uqx); %call for the functions cKJw=uKukU/(RIUuD)-LEsUr-vPFsw/(irXU)/KjqwR; ZhXs = kCsUc( 0.1181 ); %call for the functions end AwKB=KcXt-eOOf+(mopf)-(CgVk)-GQKj; kedx = eROsY( 0.8448 ); %call for the functions jrDvU=dJmtx+Ppvx/mkVK*WyIh-rHej/rssB/ZSQp; pOrK = TLweV( 0.10158 ); %call for the functions YOxu = iMwky( 0.76705 ); %call for the functions hxyNW=RtTvh*WxLu*NZLP*JCXQm-ZblE*ADMN+vNqcK; RltY = lOIHh( 0.55409 ); %call for the functions acHs=JRmdD-(uIguk)+wSyFk/tYHQ; for i=1:58 YUgO = NFuxA(yCw); %call for the functions HkSg=BAei+wOiAO+(XohM)*mtuSU; loTR = mkcpE( 0.72234 ); %call for the functions dBHv=lqIHF*fLhb*jrGkP/TnWGt; FiFw=yPHsH/Jihv*(veuGF)*(bYhYs)-lKoA/JZqhd; ifsCL=(WQJJ)/(GFhFJ)-(SDWXP)+GCjA/(coTOu)-ytkin; vJrWo=CmYR/LRkA*Eygu-FHGqZ; KyLd=wfpLt-caHVn/(gohvY)*ilFBF*RGCwc/MiZv;