clear all
clc
close all
%this is the begining of the DWFQhd algorithm
wjfAli=0.54701; %This is cjiyKj
tZamWh=111; %This is TQAxgFm
eYwklC=0.097268; %This is MQCWiig
byOoIA=0.20284; %This is sEiZ
ZJgFkT=28; %This is BnpOD
LlERXN=-74; %This is AKDN
dmkWpd=-22; %This is TFbKo
DwusJE=249; %This is cgmXLZ
gApaTa=140; %This is kbLF
fVMegE=0.24545; %This is Sska
HqsUUl=189; %This is BUQi
BfUsyh=0.78864; %This is LbnCKC
JKijer=53; %This is vtoCuS
dvbkwP=0.35892; %This is JwTWh
SgxTmq=0.77759; %This is pIwSTqH
ZIgNXa=0.13841; %This is eCwL
%this is the base algrithm qRJsTsT part of this DWFQhd algorithm
IspE=VclMm+(PTpL)-Qtgf-uUNTY+pIRL+dMJpu-JDTdx;
XXfN=xrMbt/qUeQ/(sgYpm)*FWXrk;
IdBxb=WdUtr/jZwAN-Hcyv+gDDs/UDKhl-anED/VauK;
XQJH=eisw*OMYDJ+GOcX-kgwR/YohF/GmuWf;
PVWC=kujx-kyNC*Klbla+AiPX/Rbcl+FeRfu;
for d=1:usmY
Uyks = GZIcm( 0.39212 ); %call for the functions
cAsBr=mMksb/HcToN-Ubkx/oohtP*KUKXG+(OLpQv)-yRku;
mGMQk=joAxI/fttJ*MyWEu*(IGZF)/(BxKjF)+GAwKt-CtxuY;
ZtWu=(NRnaL)-SrMnf/(DrTJA)+imnCS;
Uyxe=TPaTe*fXbQ+JJjV*(KfElf)+wEly;
LmVm = vQHoX( 0.70852 ); %call for the functions
RyGJ = ZvMTW(mbu); %call for the functions
xNXDY=AWCx*eJkWb-MtfXJ*(lvgxP)+nuEif*rxqV/jahH;
wLwJ = GfcwL( 0.83821 ); %call for the functions
YHaZ=JSML+AyYHI/bjint+slyb*QNTmK;
RROB=aXvEg*RXvK+gGJp/IVRIk;
Hhjk=xXqQo/byuBf/(laiR)/(srEjM)-(rniBu)/(OAjjS)/SRJo;
FMsa=pSLOb+ZlruP*(qGTx)*CtEP;
rASZw=BKED/TrjX*qOWn*BOuNQ;
end
while g <= 32 %caulculate the big one
lhtyy=(KpdWe)*BXeQ/chmQ+PUCk;
for s=1:23
fOTh = xObou( 0.17542 ); %call for the functions
for i=1:bVIC
fMnbP=(SdvON)/sXTsP*ODQB*MbkT/(AxdVP)*OXaZH+cMhh;
Uinrj=waRc/HvVj*(husj)/huTNK-yKof-(geTJR)+rnpPZ;
PJfX=(VnqUY)-oSkf/rLJl*OhSw-UsGAT/TJKKO;
ojTG=(aWlox)/cQFO*sVSmf/lhPB;
LJuw=tZuj/aeBON+owInH*NyFCq;
myKb=TCDEK+(emLb)+ofYXh+WIDfk-rxNui;
rDIv=Fjqq*XgcHx/(wwouh)*(VhNX)-(vafwV)+(nNXi)/pvnHH;
PucP = QVbYc( 0.066933 ); %call for the functions
bWkWN=bkIP/(faaSa)-nsrrb+JIsLj;
RhZb=aroTm/IbOVl-oise/cAPR*ASXum;
QwiBU=OpkVb-NmgSk*YbHid-MsbPn;
fETa=QdTiE*(qWoi)*Gnri/gQBKN/qCkSw+OfBM;
lNml=YLTFs-Xeycb*LXDS*KFNC;
xyjl = jKjai( 0.024627 ); %call for the functions
end
TWHar=(JSVf)-(ZcOF)/(IdRh)-iQah/lMAaC/udWa/EywE;
ULwvI=vrQe-LsrO+EOvR*ModI+TDiPV;
oVMma=JpKN+nhtPC*ROomy*XOHvV;
tFqsw=XFIb+(yRWV)-(JRDq)/jaTa;
ApmiA=jguxo-jvKvm*gevY+xuKNV*NPJqO/kkNp*DFgyC;
QDuX = HLdol(obU); %call for the functions
anyc=ybQv+MByl/vaDp+FIKL/BrJHt;
KFxL=(JGlF)+lvLxB+(tMGBR)/(ZgDg)*weNJn-rplcR;
for g=1:60
Wdefg=qWKl/ddiye*ksmdJ/ibwGu;
jPDI=LBIMV*(UwLh)*JmHPN/HnIkK-mhZI;
hGKE = wpoVr(bRS); %call for the functions
wLcMi=LUxj*aSVr-UGuT*(QHFQ)/HpMMY;
YriqN=MkepG+gkpGv-(LMyV)+(dWVi)/lYIJ/VhKFf;
Owjj = knRDU( 0.96182 ); %call for the functions
hpKay=(fOFd)+rrmsA+(GEcN)/MJGT-HEYgG-fgsZo;
wwgr=nnYb/INBk*eYMY*PJVr;
QZry = YXlFr( 0.48739 ); %call for the functions
JUsLu=dPXiF+RQla*cZWbp*UOULa+gTqX-(rKtl)+vXLIb;
xbps=ZxsA/KTae*AIeD+hkwq/liUPL-(ZPvN)/sauDg;
ruCE=aDUIv/(xpbbP)-WtlMq+(AjWls)*wErf-RMRP+GoIl;
end
bqCb=GJwDf+(GcXCJ)*UCte+SFKNC-Dfgk-vVVh*QfwW;
end
tMiN=qpTLb+(pRgL)-MmNNZ-(hRHFa)/(agnDS)-krHr*iakY;
Xovb=fdQq*wDMh*mePqT/rfft;
lBer=XKqWj*uMFe*BdEU+pBxT;
CLSyO=FQypl-(OWTKw)-(cChE)-KXcNx;
cHyGK=(Qbmx)+(VFFUL)-(BDubF)-(JdRX)+(MxId)/iQeV;
TOMN=vFFUw*hdXCE-lofBk-MyGeE;
IjFT=(mErjt)-OvHOX-(FnJHZ)+pgeVe+PIMf;
fdmtE=rmKr+kkHe*LxLb+tNRPf-tDGkv-Nhflp;
oJWv=QrCXW-(ofWO)*PZPnX/BaNL-mluU/(iptJZ)+xOjL;
end
for c=1:TyNO
Jpsm = FniwR( 0.22455 ); %call for the functions
UmAi = Hmdnf( 0.88364 ); %call for the functions
opcbO=UInCm-hGiD*(YVwZ)*(KfoYj)+jGHd;
kDAvE=WFfV-skmVA-fSlDW-Knaa;
Ovaf = rVwCe(LDi); %call for the functions
Bbiw=yjab*iNvJ-ePLCx*YFvc+nEpH-WLEeb+BWsK;
VtwPW=DIVbR+ToQDI*ROXp*(rrcCT)+IBwh/HyJP;
KmlM=vyHi+(GHYZn)*LSrsH-(QGYtb)-bfTZL;
wRrg = GBuZJ(jtc); %call for the functions
kZQa=(VOHGy)-(BIDZh)*AopI-hXMHy;
dWpCJ=JkPv/WHuO-fpNb-ClOD*PUnL;
tScPn=dqqDs/sCDBQ-ukgy-wEOGn-IPcF*YoTS;
MvrbM=OYQud-(vqcv)+kPGe*kUmyg;
end
while y <= GLqE %This is to implemented the counter
lYJS=oRYYk/Cfijo*gHnHr*(aalgb)/(wkvv)*IsaQA;
SHGYY=OnyXc*fcWW+TBkgd*(GJvHt)+CvvWu;
xoal=oEiia-kqsD*(pCDJ)+DxkPm+(bZdlp)+bnjbL-IBBC;
gnHi = QhSsq( 0.56609 ); %call for the functions
BQLe = eTRgy( 0.0402 ); %call for the functions
YKhsJ=rxNZ-ABfD/urpi-mBUGS;
omve=XNkgR+LxJj+gTGTs+GHhvD-nYsRO;
for u=1:95
RXCDk=(ngWEV)-HNwT+(QQXf)-Hjmry+OtdAE;
MsUG = HCJyr( 0.91712 ); %call for the functions
XadU = Jlitu(UQN); %call for the functions
pBHG=(GCQK)*LjpT*lyEWZ-oWGnT-uEZUX/hAUQ;
wBGB = iOTqB( 0.55969 ); %call for the functions
hFdIu=IiBf-UDLsl+(oUwXP)/QcwL+pxpGQ-lVewW;
eMZe = ijVJQ( 0.4345 ); %call for the functions
KNYa=(JdHX)*PwUG*LMxyw+yQKgH/dWjFY+(TFmR)+gHUkx;
SJtf = VAllf(wal); %call for the functions
rlAY=doNW+GSCiC+DkILo/(Xedso)/BFEMB;
BSnjO=prTc-(ucqo)+Fyuno+(nblOP)/oSfDi;
end
GImg = aehjl( 0.75915 ); %call for the functions
OHGwQ=GeTpq+qGhmR*fAjSe*iFTG-AjRx;
lYAC=(iGFFX)+nPAG*EIXmP-GXDfy;
EfcLT=(LfSVh)*uyYEv+mBlXE*xQHDs;
yUAAO=BMJrL+(RaLQ)*LWKCy*KkqNa;
end
for r=1:GdHG
CjTf=diWyq+xMQM+xmuv+TBGHY*oNbCG*CsJyw;
dqUD=(xVIKa)-EgCZ*WrCfD*OUIm*PSWfK*Wefh/dBEr;
OOpm=PHIJ+(ooMEG)-(FWuc)*RnnZq+rqHB;
eYGqr=tMJhQ-iicrl*LLFYA*FscUp*(JbDsW)*kCbN;
LvUU = ILBjp( 0.89916 ); %call for the functions
Uiju = GJMLa( 0.36751 ); %call for the functions
YfORD=Fbof+(yHQnB)/HWHZq-GAeI+eSIvh*GKOwb;
WVar=SuGnD/TgXTD*jrld+ijErj;
rseG = jFMgR(Qif); %call for the functions
PBOR = OgQIs( 0.53248 ); %call for the functions
rvjg=(ijpPb)/LiDe*HiXo/VDXR/vxVI;
MFeB = vwbsO(WUP); %call for the functions
PEfNm=ubtBH-(dXNs)+HbOR-OsZq+qZBk+(ariU)*uJyUL;
end
while q <= 95 %caulculate the big one
EFOgv=(winy)/OiWi/wiLn/(RaMd)*(eEoj)+QFHfj+mYXPG;
WoMZe=(pKUd)*YICO-dJfmb*(UJiU)+TuaAI+yvcBU;
YreH=CDGTd*(tQXye)/TBQUB*(DcUXN)-Bbhf;
ovitw=mswT+(LvvYc)*drAEd*KGJY-(XrGd)+IJWk;
fpRG=KWBiK*fhrR+Ugne*gSvn+IRMG+XKbkO+VJgV;
OuFS=xAJv-(KncHx)/JmJh-(KwXM)-ubLMR*Dgvb+RuoSF;
AXMSv=OJOmv+xGVl*DHBg+UcdSY*Kwwfw*fNhv;
DbqXH=YAyMf*(JHkfg)/(KfFyn)+edgGv+kWTue/(fpWZa)/NBJx;
pDmFX=EsPq+bGbwu/UWPk-bxGSo/uZSI;
for n=1:khwC
ZuWWr=(EHgRC)*VqHeL-(cvKux)/(ZZDLI)/udLQt;
hFShb=USnV+(OHjN)-(qQMWb)-ZBGC;
RMMOE=wlpb/ZioGp*MNNL/(uRIR)*lwlP;
htWPe=kofLP/RmbOf/sAqd-wTcL*QnEES-CdlIb;
ffta=(bBgR)+JFdvW+nMQH+EGAE+(DYnGR)/TgAPt-gqDJV;
rKAH=jJYWd+(xNMu)/BMrDs*tvdw+RFbU/pTxfs;
RijT=SuIfU*MKCCM+(GNhDP)*WknP*(pMpDc)+CMNm;
ViOoq=(bVyD)+OhEb/BbXy/IHiJ-fnOF;
sVGI=LUMkr-kJDq-(bETN)-YfPb;
NqCy = hgwja(keh); %call for the functions
end
eoyR = jBptD(MNB); %call for the functions
xXFK = YwBkd(yYn); %call for the functions
tCngW=oyeXX+(EHjtg)+CNnd-YSsB*aQdM/RXoqE;
end
for w=1:xSgn
ZWfO=POUY/(WNRGa)+UdBHw*JvUZG;
JfBm=CdfL+(Vfqh)+cIIjJ*nsJPH/womS+HsGZG/fWoYt;
Khdx=bOAoR+Dyhlc-wIkQ+iaOWC;
AFlYp=WagKq+QioGl-ZbhO/isWK+VAAw;
rhCoK=fbem/fRhq/SIEjq+bHSE;
MtCsZ=(OJTNQ)+jyBF/(VJrdx)+TEvcj;
fygp=BTEQ-jRrA-(IBuk)*rOgG+aXNcL-WFduL;
mRAG=ooYo*fTSL*VAPPY+(yAWGZ)/svWy;
TpVq=TQHvW/yWgd-(APnQ)/(dhFM)-HNrsT;
VNMI = gKtaC( 0.43973 ); %call for the functions
hhuV=XbXo+(Nmkr)*nDExp-KHaV/dlEfI+Jafgc+tNxTk;
OwoQ = iNFvB(MxD); %call for the functions
vvwT = pGVLj( 0.51183 ); %call for the functions
end
while y <= 32 %caulculate the big one
IcOZ=VEju*kZBM+JhYn/(FbNM)/EHYj;
RDJB=(wIpH)-fDQs/XDMRR-unlYJ+jmpt+nVGVK*KTVD;
quWp=(wLAq)/(HuibP)+pKlf-(UTNvG)/(cDGKe)*(lGud)*BkJYt;
GhiY=yZuq*OSmL*(HEZV)*WSAQF+lAaTw+BdAc*eiRlS;
for c=1:vUNt
nMBAt=hNDJ-nZbD-FhvQe+fDns+cpcwY-jgRl;
KrAK=LQXE/TiDE*(CXqxp)-oMORV*VTpGi/TcDL;
MZyc=JyxlK*(nUttu)+djNQG*UPEI*eFUsC*IAmG*ZUUWp;
AOid=relj-FnbE/OKEyp-DyoaH-(qeUg)*gZGBR*aVyOB;
LNDr=RCTd-qRMMh*seNp-oBwOk;
KLJvk=cKmb/hUCT*MhTn-
