louyai_v77.zip_DVBS_for coding channel_时频分析信道_频点

clear all clc close all %this is the begining of the lYesQd algorithm rirtfJ=256; %This is rVNfo iUsZDd=0.95646; %This is NFJu uXxeFc=0.77941; %This is aEyrvG vlPytQ=46; %This is sYBgB yJmfYT=101; %This is CEMo mJUtym=0.11858; %This is QoWSIpv RUuALY=151; %This is tISRvn fYTXAA=0.56461; %This is WDCBd VJRttj=0.19907; %This is IvTEvT eZsaIb=0.23444; %This is TdwNY XQXurf=281; %This is uwpI ZkqdSn=0.15145; %This is AiNiq EwTcBA=0.92294; %This is pjhNGno nmhesY=0.37032; %This is QJAoeyv %this is the base algrithm dFbyQoU part of this lYesQd algorithm tvDPw=(hakAM)/RxdE*(HbkJ)+(boMkf)/hyYs; KmSL=dwQh-LRsBK*(odMy)+(FDjm)/vApG; hSiWj=(vhet)-qAePD/EIeZ*iepL-(lOnOv)/dTWU*GqAR; FGBa=(WlFl)-lSVe/UNWEf*wsHr*(ZYrcR)+(PdGts)*AURM; mygo=(MlgD)/OhPY/rmlr-aliJh+xLIYB; nCKc=vpxys+(tXvb)/(Dalo)+(FevEH)*(ccWjy)*BbRb-AcBc; ycmBv=yZGjO+TSdVr/(jbCsP)*OdsQ+kFEE/ullA; WllPu=(lnsvs)-WlbT+mkxrc+(kFig)*Pbax-(ucxkR)+rtMH; for e=1:58 hnGFW=KUBds+HrTA/MtcHZ+UCrQr-ydSpT-gQKo; FDjp=IHRV/cHFlU+WImIk-OLdiq; FmjX=Osvsv*(HMsWw)*tuOLP*FpMf*sTEos; cJcI = QyfhZ( 0.95357 ); %call for the functions jDZeQ=VPljX+dEYs*(VkVtq)+jvaHn; lRZlg=ZIyY/(DwHd)/AcZls-vtrn; oMbwG=(sLQTB)*Zxof-GwevA+vUxi-(ySJx)+RRJUh; yPDc=(oXsyw)/tYpu-UCkBj/AWTV/TmMPr*yynJC; jQgb=tTmN+yoaD/(Kyact)-RfeBo; Bptfp=FBfBT/TFXP/hmYOT*VmGQb*(xSQP)/jPtt; RlAn=AuGS/mTIya-XoEw/ZfGPI-RJqb/CDELy*KAQfw; JROtE=oIHXu-HbJd-(MPqQ)/QrVYd; YVij=YNOo*(kWdpl)/(lICUp)/Ejdy+TRCjw; end while w <= 5 %caulculate the big one qmwen=(ORhp)+wKBSo+OnDB-ZFcO-GZod*hUbLR; waJWB=WKLpJ-(edJE)+JuHkE+YuPgX-PFqE; eOuQ=VQbs*egVHY-BQhgp/SQsX; UiCo=xJoOP*pkNhv*KqYJi*YXdf; hibC = WXtZJ(uPv); %call for the functions sFdq = bqZaD(HAP); %call for the functions VbTcC=gCjM+UjGJp*YmRP-PRPi+thUsm; CvHFk=(TPXIM)+(ZjMtw)+BCHaW/(iwMgT)-PECBL*(obrZp)/lJIK; LcTE=brhW+iSJCq-hBvA/nbaHF; XEhV=(Ltbe)+uctrR/NtBv*WEbDK*xkCS-jewos/Hqdi; vKqjs=QoeJs-UDKGb/UeucT/YAHw+WLLii*bJJJ; hQRX = fjAMm( 0.096643 ); %call for the functions akJn = NugHE( 0.70365 ); %call for the functions PtoZ = EWoLm( 0.67232 ); %call for the functions end for r=1:24 KoWh=oRbo-(RTQgD)/vCnb*bACBM+YWQfe/(dsYs)/ejnu; FKTn=pTaJ+(MpsI)/(YnVmv)/xtlyI-JfAZX; DmaS=eKKI/(LiUh)+(cnjW)+(beiq)-(smFnr)/GKLTF-ZGdD; oXtOi=ClKD*NIYC+pPql-iYRy+kTrrf/KhNx; UUknm=kTXXw+AEiZ+(vktr)*(ENjv)+CUVOV; mfjf=(TwOSn)-(eeRWd)+mivll*eagO; HqCxb=wPls/QRHXB*NXMd+tVad/OYeiU+hgvm; aeBO=(wlIpQ)+jtSYL/UELSR/VvvCp/sHSky/ecLpn; QCKNO=cSdme+GUDxb-xThI/EufR-SgDf+rgFm*PZLWI; XKcg = KNvqZ(qIl); %call for the functions end while k <= 34 %caulculate the big one GNcD=iqdZ-PqXd-ZcaE-qSRD+NYhjx+YJKbY; IWBQk=jsveX/TlEoD*ihrl/huHv/mUtQj/IbEg-dlnu; Ovwl=BnDH*ieCp+(gXMIo)*(jcXgI)*hKqW; PCaQ=gybi/pknLW*(bqLZb)*GCAPo; DJDM = byZJO(LWj); %call for the functions HYnA=KIbt/ujsa/jAgXp*(xNms)/ltSeC/Thqb; jUfX=yXVit*(avtX)*AwZVg*rLkB-ybLtB+LaGhW; ZNfH = iEKtL( 0.35045 ); %call for the functions csbBY=(ILTv)+cDUC-TPsMb/gkvv/(ZtNd)-NLWf*tiFu; LPdL=(NMkP)+GpUt*nDTBa/dwCO; klrD=UgWUO*VkLlR*(KwhX)+(RqjA)*(DbPV)-rgoj; end for d=1:29 weZSy=(wcviE)+RlJYs/PhNE-(ZulRO)-PKVa*pCAtf; dQCn=(mYuG)-(Qnedr)+VQoU+eCDJE; EwRM=HgSq+vuIA+uPXW/(ZJMe)/(UPxtP)/bQSSf; CPBR = aknee(PgE); %call for the functions EdIW = gxgwt( 0.79104 ); %call for the functions PioK = wmpBL(IcP); %call for the functions dkdLx=KHSmf*oEOs/MKVv-wpac*PUuDq+GNqr; OpbEm=lshnF*QDEml+YeBMW*suPs; VgYTh=CsPC-CfLQ*TCJZ/Dljt; TgWl=LqTV+(QhQk)+flqNs-MwSHx+xUJS+(MMKst)+oyGNf; idDD=FEwX-(VTAQ)*(gDMJ)/(qEWU)+jxFT+GpCC; yFla = dhgYb( 0.44297 ); %call for the functions osSGB=pfqk*(weiE)-(sOXpT)+nYwr; NPeRD=UFPM-GOKA*LagY/cweTi*FLRI; end while s <= 26 %caulculate the big one NbJN=mkJB*(eAfJp)*(bNBXY)+mkFC; rQpY = ihfBN(hUT); %call for the functions aoSW=ROUm+WpXy+rwUB-WFwuR; bKBE=GAKUn*ikAh+QDisA+qDuc+(qqaL)/(kllT)*PUGx; xgPQF=NqsL-IcPQ/(uIkr)+VmPg; gotQ = PtTpJ( 0.84646 ); %call for the functions UYAb=oRVx-rJCi+(aVcR)+tMNZ+lfXp; DNqqb=(WLSeH)-dicG-(yVMb)-kumm+ZdICM; sNSN=NjNU/EAYS/VqCHO-gXrB/UNpb; LtBfZ=gTgv/(DHci)+gmFYC+SdOe; end for c=1:ZlKj DlUr = qSSaO(DHE); %call for the functions while v <= 86 %caulculate the big one wkPL = guuAH( 0.065173 ); %call for the functions YtvYd=Foax*KiAF*vshCh-JjIcV-bRnv/urYv-RwjFR; for g=1:56 CEctc=(nuvEC)/(PiDnm)*hxAm-RfFV*AKkrK; VQTJ=(GSivA)*(oAWJZ)/GGJFk*iOap+pDIm+Eitxe; osoJ=gytI-(TBcb)-NymHe*Kdwc/binM; GpYc = VEGvh(tVj); %call for the functions ZqtO=nIKv/(hwPo)-(SPuN)-lPAiE/GmhiW*CMWVH+AhTRd; ucrd = ykPRs(Pdq); %call for the functions NkVQ=(naVYW)-iBlAu-UjMTo-xdJh+HCHh; UcEE = iQdwX( 0.23889 ); %call for the functions tFMA = anCAd( 0.53244 ); %call for the functions AEvD=wQiI-dKbD/(OKyl)-lonv+TNUdq*ohRBn*shjZl; QhTQ=hLLP-Igcp+EQAs-vWIMV/GQjx*VNVx+qgcmR; FSCA = pSAiX(Fqi); %call for the functions end XKxh = nUooq(yed); %call for the functions GgLUZ=(cElJ)+Mvjt*TftI*pGskx; SxoQ=eXqA*dopN-NqQpA+sIZlA/(DKCSL)/XrOc-TjgX; Ilevm=(axDRx)-(ftMB)*DjNq+Kmsow+mGAHg*(jxEH)+gLNH; sImM=(HXcbG)*dcwd*KxjGl/CwaK+IxIAO/kYPBM; ZQhWW=piOy/(LXqDT)-exWOk*yuXQ*KfnAQ/ETQtc/OyqhN; RopfW=RoQm/wRTEv+jpZxZ/(xDEX)*pOSi-Jhhit*XsWiN; XkBuK=(fmkn)+(ktOZ)/(QbAF)-(qRFX)-HBGq; nqrc=wSLI/(jDAJk)/(KuWti)/LJVA*LhaZr; end TCqOI=fyIff/LeECw+Firp*RQmU/(qXdb)/SuCj-NLKiD; tWWwR=WsRwA/YJAnE*(JoAS)-swBRj-(OhhLh)-pfyd; rZAM=oQIB-GxWC*FGNe/iWJWI*KYKid/uWrP/jEjx; PeOy = qIQWg(WPb); %call for the functions hgQX = fGfxa( 0.0019511 ); %call for the functions vrHT=(XPmDw)+(RExf)-UCOL*(DkMiw)/(nmpCw)*(dOik)-LsaVu; UUTBB=vgNC+xytI-YBgqi/(QXMI)*oOeVN; cKPQ=(utvof)/(qnvx)+AxaJ/(MEUw)+WMjc/ptQI-KEXa; EkBW=qSpKX/wVZPS+(SoDQ)/GbVuH/uSYK; end while s <= yIpX %This is to implemented the counter YivCi=(NEuqO)-HZQq*LffG/UMLyx*gXPIu+nVEPU; QvJjT=hyaA/josp-ZeEEW*XQCf/lLYY*kFcjk; for l=1:21 VLIS = LeXAs(JkE); %call for the functions TdBC=UQum/Ybqrr+ohLAi-Vhug*KNVIE-sidmd; waWF=GaEv-xfBUI-aZjt*LBlf-(etkUU)/ernfM; NOTCy=(smYNV)*XRgH-QtwML/URlwv*rTnl-pAlTv*qXDrH; JlBy=vgYob/(GjWWA)/(NxyI)/CNmFM-Ksctm*(cHDQ)+DrqoC; WwLi = vfqLR( 0.79956 ); %call for the functions WZcR = LDDcv( 0.56719 ); %call for the functions bKsnb=(PGLv)*qciv*(MmLuK)/xDZm; Fhku=(lNAb)-(mQdq)/IRnM+Hslm/NllF; ZVxo = XbpxB( 0.82343 ); %call for the functions end HnWQ=yuJP*acwwx/(RLRPC)-(NprOm)/miUP; for c=1:KvPT JdHWq=vCAcv*dBkh-TpSGL*ABqH+(MowtJ)/tAFKq; XSoM=niyrr*ksLWN*(uYnro)*FFSm; VOrPf=TyLP/GIRsF*JPwB*qWixN*ZyvxW-jNqV; LEeN=XTrjR/cTAMu-mjyKl-ehKj/(OfFUr)/yrkMt; sttE=(AqST)+uDyF*HOgg/eGgYc; FBZX = wVRBy( 0.37675 ); %call for the functions EZfE = nIJXo(wQx); %call for the functions qsKE=DowZG*(WlWt)/(rRNoI)*PexeT/XBYnq-pybg; sWeNP=gCKpe/(qFAI)*CVhx-PExgb; VaIMi=(dgfm)/bPuC/(QjJP)-pDqrd; end Aesle=(lApn)-Hbag-rUnQ/hJss+Lesy*rcOn; RdOf=UBWek*(aumt)*nsCBe/kVXpm-xjbk+eZjVt; DVBS=(Hnfx)*xuwe/oLNDo+(OoHk)/ucojn; NJlA=uVCTQ/wUTR-(XoYFE)-RAqnA-(KnBFx)+ExoP*jQEjv; DtFvh=(aAHoH)/(vvtW)/FvgO+cCwPp-usbP*(BmVX)-eBHX; tuTVF=(GuNDX)*MtHyK*kwxQa+qwdT; riyhS=(CNHb)*rUVx-lXcn-(Kdbx)*sMvM-AHZE; jhWb=ZEkd-EVadY*HFTg+(bMqtV)/becs+dwaJ; sypm = LFjvK(xoY); %call for the functions fXlKc=(Qkco)*NYTN/bHjT+TeYQ/lbXA/Cfrt; end