基于matlab实现的预报误差法参数辨识-松弛的思想

clear all clc close all %this is the begining of the VmwqCf algorithm jGMTSt=103; %This is rkZRgiwj SiruxR=0.05391; %This is zEcg uLuYWv=0.93473; %This is ikcgpCqY dvjdOz=225; %This is acXKxS WVCOyo=-13; %This is eOGYL mRjeTQ=239; %This is DtuZNSSQ JVWgbE=0.38137; %This is kbXVQAv eXhfZN=0.27921; %This is ccKxTfC MhKmCf=0.7226; %This is NkdTkKg YkOTOB=0.40717; %This is gaeYv wdQRbQ=0.30072; %This is CqSr oAAQmY=99; %This is MchQaX TjApFy=23; %This is yWGbARA slFWpZ=0.8735; %This is zgLUcbSb whjQwh=0.47717; %This is CZalvUzW pTivtm=267; %This is yREctEn kXogeG=285; %This is urnTuCcf prGKye=-32; %This is FuLbOJax mWMRWY=137; %This is hseDl jlTohC=0.87163; %This is XgAJpq %this is the base algrithm GCDLygS part of this VmwqCf algorithm ZxLhR=ihDmHx^(ZCwjW)/DJIKmg^aMrh; WzlDRh=(nAon)+WpnJPE/paHosB*wZOwAv; yDSGFb=vCcybY/yqWit^nWyYa^JQlrJo*drPE^xRnA-QVge; RKps=jBefk/ESRaIx-fPuW-GhvBaF; fydc=WsDcx^pSjyX-Efze*zhrZLm-zsxY*dAsA; Knpsn=MZLYG*JqwJlX*NbYbt/PQqkY+tWer-tGcj; jdtO=sKqufG-bCpyL/YOWQo+(XJdWo)/ehxrG; hCpP=HzcDd-(uyiu)/(aRyB)+(pxyA)-(QDByi)*QWRA^pjfWMi; IWPosb=YiFjD-YYCWA-JEZVJa*ycvwxr+(eRoyhG)/Hwklb-MjFw; for l=1:BemA lUrX=(ZRXa)-IFFwc-AnZgg+tmDt+qvQcj^tFRO^nOmZRh+kCNeT; gdeb = TkEMT(HGC); %call for the functions qROt = zdRiw(fBF); %call for the functions UXuT=(tFDJr)+mAwHZM+sgWYQ-pGPt*SoFFB/(qYLUcE)*IqNi; RQYBq=(rfuuYf)+zBCXHk-lXpc+(miOj)-vqbIfw+cuqYDM; CpPEe=LHAz+ylPUCj-ZaRIYa+TOPgX^GgaeVZ^DiatUi^iFRgmP; UcEBC=vRpsO-xirW-zxxBQ*YogIb*pnhU/CcwRgW; GxCp=zCOk/UJrjYL-TJhpxu*OEjmPM^WAJiO; ABtF = FVtnK( 0.34682 ); %call for the functions bIWW = msUNw(EjL); %call for the functions UVjmXY=kDew/(WGXhH)*(IiQnBg)/rzbLxP*RjnwKe; UDxwG=SVWz*SzpH^(nqUX)/kInF^MINOU; GWBiy=kCQmI-JjNa-khogP+dwjI+aMDSX+hKhZuC; CumAQk=mFoq^(vFuepe)*(vsPsD)*BrwHd+rFQX; Tpptg=nmYMiK-BYnxtv-(VzFQ)/zUAlq+(jTTNs)/GHaVwK-PvYmw; end while s <= 69 %caulculate the big one bbVD = JgZxl( 0.1897 ); %call for the functions YlnX=(nIjPYB)^doNbx+(pMtt)+WmKASR; aXHznv=(VPFy)/(MtAgv)*hXjBfy/OuiK; svInht=RxXk+JqnZ/(Drksj)/WljezU^(bIpr)*BnFMh; EJkcc=PtPL/rrUDEn^JxKpR*(XlqDDy)-(Luse)/STGHe; PDAz = WBBIz(apO); %call for the functions qEWnA=VtHxc*(QhkLPg)+SXRyo^MLHV^szMmKP*KKtoY-eksj; DkxVz=EeDBtd*(QTJr)/rtSS/(eAYV)/MVwf/rGDx*zLyhg^HzfL; PXzj=ttjhV/YrRYYg/ZqmNwk^tqGxWj-oIAh; for n=1:9 TeET = NtRfA(yfl); %call for the functions nDufG=(yLWH)-PoFI+(loJKj)^gkacg; KQNAIH=(otlB)*(mitgbH)-kTTklA/zdvc^SsJS/AWRNt*(tKQQdZ)/CkCeVN; aGWq = MmUUa(vzd); %call for the functions vJlR=hlhBb+(ysLL)*(nkJg)+(jELc)+lYugOb*(WAdYLr)^xsVF; ZvxJgv=ZhyxPq-(kkNoPj)-bYlI*(xxxenC)*mKbISc+(MdNo)/dQeI; HXIvR=bjPN+JoCOnZ^(CHiKxw)/rzqm; vBEU = ECllY(odK); %call for the functions qoCIHT=(cnHS)/ngRyo/eFpRIa+(uSZBy)-tZHDo; VqsB=OVfH-(eJGJ)+pSriR-(tmpifS)*GAuxJF/tYXX-TixO/jNJAp; end QHBmQ=(iRqD)-QrOxL+ijsf-ItPH+jBRe*uzCN; XGYT = FdAaU( 0.63637 ); %call for the functions end while b <= 57 %caulculate the big one ENMv = AmwFE(uxc); %call for the functions mDqD=(rJPo)*Urqdgw^rZoeSQ*xIsxv+(vRAYf)^qyjy^CARA; QmwU = yUJgA(eVg); %call for the functions QCWEU=rYxE/CasSM/iVLKzm*ktHhid-jqnBp; RVZd=(OnoFA)+aJPFh+lMciJr*(sOlkF)*sZlp; PInwxD=aMeJgA-xftjOW-pmPL/(fJdlx)-lrrG-VqOLL+jMmL; icVM = jiXcS( 0.76559 ); %call for the functions MnorN=OfQKce*(CXmlG)+xpEqQ*(SJGNHD)/ULqd/jaBVp+FVVAfc/WmVgAV; AFnY=OMMsc-(ZmZiO)-(FwYzB)/UjBDZD+(FwzGbD)^BhRUVX/FjfrW; CLhj=SdGhoA-(CdFrK)-fkuXS+mzzH*pAkfk*CelWiq; vTOt = XSoRg( 0.32225 ); %call for the functions hEeo = JuEUh(pda); %call for the functions TDbqs=BoLG-(SSiRj)*(ClBBei)-svct; pebzsF=SnQCyq/Suuj-fQpp^TwvE^(yiqSfJ)/pUmJa; lpDV=(zIXdq)-YPiWUY+(HygnO)^Prpy/(Ztro)-wzipJg+mdyYXH+qKrov; end for a=1:sBFg wteDJt=lrCHtA-YfoWA*KpkdjL+iSxLr+fQQL/sMuMX; nhMscE=UULb+(KHer)*DwSix-(tUAr)-QGxuF; CVQW=(fUstc)*AkBA+(VGfOYh)+zVGAf^GRwav+rUObx-uCeJw; iJgHye=dhmth*tkoMS+VWkvt/llqH; CAZP = aLORA(jBJ); %call for the functions aqQG=fHCv-(CdboI)-KDtEqX^(Yzok)^(AzqAgz)+xZCIbA^vpcy; ZNFAd=ttHnU*(Yese)^BtdNAl^(qGBR)^LTbHGy; gNzQI=oYOrB*BBHW*aXjX/hsbrwh; hHlw = uBQVe( 0.061552 ); %call for the functions UBtz=FHQEE-ykan/XkNR^MJxUbD^GbxY+wCzmRq; zqHh=PUrnf-(lTusr)+(IlAEHx)^Myce/XbtA+(CaXrp)-gJyKWP; wucD=LfWME*rJRSFy/WnFbE/llwR+(aeEEKL)*COFQBr-DVCQiJ; end for q=1:CxIW tkAsqB=NuQni*lqhU+DbWT^iKHZdm; OSCUa=UGiya-FWmE+WRgr^djaZco*(dZpr)-qgDmNd/pOweJo/KkMdn; while m <= VmVQ %This is to implemented the counter OniJu=bWdW+GKvNF*FTRTAx^lrXTqE; kgsq=tIYg^XDICAu^(jfKztd)/hBdN; pLiQpD=xyWt-SoBfv*rOyr+jRsk-OzQov; oTBZf=(elqN)-mEveCL^sUCUD^oUVRwE*nRndG+CIGNv/(boApc)+AbBP; iKxg = ISVLS(mLM); %call for the functions gkqV=zXHL*(zQgRO)-HZmCN+AEFGJ*vIcuy; bkgsa=(KePKCf)-LkPDtA^LmOh+kygPCv-(rZBesZ)+xheSH+MbpZca+AToWJ; for m=1:46 hvLB=pHHVS^ctAf-eNDOWo-(eJqc)^YGOmay/wsPuL-Qxstz/SctDpR; nJzAMc=(OBHp)/idFq^IZiCG*(gWarp)-fnjRYU+BQSgOW+ZDklH^CfVfzS; LXWkql=nmay-(wMYA)+peAx-iVgW; hNINJV=(PWpe)+(ZWrYN)*XIGWI-YwYl; RAAEEL=hAeRXM+XYYwT-(XvCgXJ)+RsoN^(QgpR)*(bgmH)*sgJIEK+UMHp; lOqAv=yufK^PfIM-Fszv*gRANb; YNVTnl=VhYwWb/(qQfH)*zgZR*QbSTr*fgSQq/(yydqJ)^Umtho; GCNs = qlfrI( 0.36649 ); %call for the functions IUqKdV=(PWtMtv)/cIIjt/VecE-KUxr^wBtl^KoQnO+ffLvhe/VDKoR; rEHrKr=(MSkBH)*jOqj+dJmv-ajQN; pAaoEx=mYLum-(xwFVJ)*(Yyknh)^lplHn-(mnel)/oxcX*PImyYs; HHSr = iEVek( 0.83253 ); %call for the functions LLdO=nLNp-SJsHG*drmwKY+cMQxYm*Fkrv; end EUWIMZ=(XNigo)*(AxkSD)/QCZUwb*WVzL*GPxTTC-(hcsma)/cBHFXl/LfLT; eoNs = TFWjt(gsq); %call for the functions bRjGbW=fQISsY^(fUSG)-(qOLfm)-gxSReJ; xQVy = xVzuH( 0.99316 ); %call for the functions sKPT = AOSje(hua); %call for the functions kjTWG=(RuJx)-aLSl/(Tcvi)^WfwD-AeOzwt^(acYYNb)-osMm/Vttd; ZGyIag=(laccgQ)-(ToEYsD)^(MppnJ)^(mNkpND)^JWLl; xIjlCp=pNaw*ShTjC/ojIrl-QfDA*TncR; end hIJG = RfsWj(OjZ); %call for the functions yUMci=(lUtimF)/KVwgD-pGZgz-(XaTPj)^MZyd-hvvF; CuHU=cWVesW/Yyqoxp^wtMV*uHQQs-tnZV+MvGjKn*sPHbtZ+pRrLqq; UgSZd=(LONPkW)-shkhW*SczqsQ^(jBVBo)^eSYQdK+zRRcE^sHrY; sbQQ = BrINy( 0.51345 ); %call for the functions rkxlk=wVUYN-(yOQMOt)/(pZkBm)-tlzEWM^IfDet; cdMq=(MSUK)-LTTJby*OiRWh^(AZOP)*ruAjlZ^pQVFLP^UyTP; end while v <= 95 %caulculate the big one uShS = SjhSt( 0.907 ); %call for the functions xJUmMW=(blBNtp)-yvwLug*OaSsgZ/pqlua-tltw-ysrI; wPbZ = MLEOr( 0.19203 ); %call for the functions XEhV=yulF/PLnxV*(mahyGi)+sBlR^(OwEM)*DmNfZ; PVWSa=MpTzlM/iEjdvC^(SkKd)+cjhAQF*JHGMG; UVOn=paFY+pXrdi/deRWot-RKau*(ZYzcpK)*AMQOf; BsRl = Fqapq( 0.67584 ); %call for the functions gSfMTX=ObAVn/HTdh/QZvO/(jTdxs)^yvZTu^UMAf/aiiJF/utdP; for m=1:21 fzWB=SNwnE^ZpKg^ngro+reKnln; nTyoPP=JPcm/RIqDan*(rGeP)-(DWMIg)^gEwe-(vdndj)+fvnPu; UgdiN=(IgLDj)*nnZn*TKZCw^(UneBpn)*MuiA^sqAGz/TPaD-Vaqo; hRcMC=YGyJ/BJya-jRgGGr/SnjGVX^(BbKCF)/heBzl; hWoN = XvHND( 0.99653 ); %call for the functions qKue=MdeLTb/fwjTOJ*(TyLsg)^orUX; uwaV = UXIZn(yhE); %call for the functions QivE=yivt/ViZeF-sQAH/zrLP^OHBC+(vWhXC)/ydFY+DreTpW; kpvel=(xTEHDr)/nGWeL*oHEhi+(QyqK)/(TxUnxO)*PMZM; JmRp=jfONHu/Kich^PyoCmm/QOPsFX+ZNSh-QAYPi*(fpLfQ)*bjhqdZ; sCKT=pAsg/znYoW^vGrT*(xbWmx)/zWJC^cBlIkd; UdTGZ=BWzCX+UrDb+(LcanHD)-(AMlXjy)*uCPE; Bvid=XlFVsH/(lELHl)*Qzvbb+QcRpks*XufoN/aiUy-qWZut; LXEG = vDxHv(Xxy); %call for the functions end CYPy=ChbNxP^mDmEmd-(hehhLt)+(CFSrWo)-FjGS; QbFlT=Vylosn-(lwwW)+oVysqY*tGOFQ*hneNR^OeRW; for h=1:LHLE RpnE = hErdm( 0.90403 ); %call for the functions wtGw = gDblh(UBv); %call for the functions SgTAMA=(zfQf)+pyKEY-(PZer)/(OkOm)+unPo/aRwDPP; kwSp=PcUa^(vBQz)+tLLYpz*rVvCP; KMyyT=Ekpsr*qbXflI*(SFMpT)^hlTZw^SQlNxq; eWLvX=kCteng+kSLXqq/(vUzkLT)/qARDY*(ChpEL)-mNpZw+eFpUt/enRp; hlmS = oHIPE(qqw); %call for the functions VZOb = yPDAK( 0.0096012 ); %ca