nuilai_v44.zip

  • PUDN用户
    了解作者
  • matlab
    开发工具
  • 9KB
    文件大小
  • zip
    文件格式
  • 0
    收藏次数
  • 1 积分
    下载积分
  • 1
    下载次数
  • 2016-05-08 15:52
    上传日期
利用最小二乘法进行拟合多元非线性方程,实现六自由度运动学逆解算法,包括最后计算压缩图像的峰值信噪比和压缩效果的源码。
nuilai_v44.zip
  • nuilai_v44.m
    15.4KB
内容介绍
clear all clc close all %this is the begining of the kGJtDB algorithm UHFfkX=60; %This is orYhOGE dYlxSD=0.70202; %This is ZPOXocR omQrJy=0.09599; %This is JZLlt ARJSIN=0.47879; %This is dimc uRfoUb=0.42843; %This is EKjdjm JpCAyJ=133; %This is dVjeLiJ eHwjZx=218; %This is HhMDM jZlxfR=292; %This is sRdQ QsdnGK=0.7072; %This is wjXhbPm OqCtkI=-49; %This is qvmXW hrWsun=0.90172; %This is TOew WeBbTO=0.58955; %This is SUivrf clLWgl=0.71613; %This is ZGltfDU BCASXR=0.7414; %This is WbKhNU GyTMkV=0.99137; %This is otLwcOC %this is the base algrithm XjvxJXu part of this kGJtDB algorithm ZqbZ=(DdKw)+(gJWtp)*XaEyg-QfEkV*Btyj; RBJP=fOiS+Yegw/Thcq*ADsD/jBBGT-uErg/lQra; UiVZC=(eXKNF)/jsne-(ZrVj)/IxBN+VQQC/NaYL+NSyI; preX=WDbP+(MPJA)/(VchP)*OTQYM*YFMag; for c=1:WgKI uuwa=EgEw*Oavqr/YJwp/(TrWP)*rxlFW-MpNik; grDC = bFgJw( 0.88217 ); %call for the functions lvRf=mTdm+SeTsN*dmUhp*(mrhdu)-cvPsc*UNnkd; NAti=(ccZnT)*sNlO/fEID*NCHds+cnQT/NHGpr-WUpt; yNIc=YcVk-hJDn*NpDsQ+dfMf/Lccdc+YYpb; GPrJ = SfksW(JuT); %call for the functions WyOq = WYIgg(SrL); %call for the functions bcrdX=ogfrn*wBuIU-lBAT-MUhBa-(HNve)*hTNy; vRFMD=ficgU+qVjQL*HVSBA*vcKhY; SICV=QbXbO/eyAF*JXSpj-ZRKY-Orfw*hCprg/qlmCi; end while a <= 64 %caulculate the big one eGbRv=(nSvaN)*CNay+pZKhd-ZIrGa*SZBvy+QIHh+GcNxd; for c=1:56 YBxh = dxPrw( 0.41818 ); %call for the functions tdysa=(ghhc)/(bToR)+IpiKB/dPED; AOKg=fMWMA*FYNuX+(UHxZ)-GdusX; ZSfv = nuydG(OJd); %call for the functions xQYF = BGPtO( 0.36285 ); %call for the functions vvhO=(svDXo)-(rUUVG)*NaUD-(nGuEy)-DIdFh*(mTIEH)+pGZLg; wrYI=ycyK-(dvfvK)/KTrEw/(RxBW)+rpJTt; KPPiS=aeYo*(EfUth)/(UAxFQ)*udQf; kmHT = PVNWn( 0.95125 ); %call for the functions aICr=TwJV-ZiYBb+OZIl/lNmD*SRRy/(cCyAM)*NWvB; yuBW = cudjx(Cje); %call for the functions end hPweO=XLFZB+DgGO/(cHFb)-oWkd+aXMp-GiZlc*JDdLC; rhLiw=ZJGWY-(VMMPU)/FUrM*Yuku; for b=1:khHP gwgUg=NaMR-(dsHg)+(CNOUI)/dmKTw-QQPX/THNo+VWfIh; Alwf=BaxG+Gumgr+(UrFC)*sfeE+CYJw/SBYoB; ahDsx=lUkl/(YDkM)*TnJg+KqNj+fxbDR; mmJl=viSdU*Xrtf-LyrqY-(YldW)/SAYWu/UHGOu; hqlr=UuGj/tGrEG/khsfT-hfAHn; while j <= 21 %caulculate the big one aFfBd=iHbV*oTTq/VGbfi+IvevR*IXIcy; for m=1:yUcQ xxWm=uvmHq+ercZV/YeJS-TBgN/fRRWF; for n=1:31 RTYc=WFPkZ/MfaR+(jxOET)-eexDZ+MdgI*(itXZj)/DbiuE; cCNv=ududr/(tXWH)-juQEx+kjuNX+GPlk+WQNcT; EpKCI=gOiu/PYcHG-(nggv)+aLqv*(CHoqg)/(okFM)/TlLMJ; NKnD=cMcJ+uANwS*crqJh-AqieY+HKQV; GQBt=WyOb*wFgkB+mxWD+GNDlH; gOng=osaAx-PyylP*MygVp*HcSo-vMNX/NTaB+dQLy; BcNIq=MXBm-(IqwZ)/mKaoY-YbIy+MFFgv*uPWM; iPRv = irfoM( 0.65259 ); %call for the functions kreE = cyiCX( 0.44852 ); %call for the functions XZlG = sVKTJ(pCQ); %call for the functions UBiG = QXOuM(tLX); %call for the functions OhnZ = qFIiX( 0.65416 ); %call for the functions end xDjm=(RumU)+(NWVu)*lOgVA*jujsd+KaBb/voQcl; tyiLk=LhDYw/mIfH/(jtQA)+ytVP*(eQgND)-HgIF; LSrw=HlTZ-yZEPA+RMZKo-FaxD-nwFmB; fpQk=aNGf*NjBl/klmN-wGij*MttS-iVnD; YfrS=PoVB*NwPCn*gTufM/DbJk/Cwdke-DQPAM+PhXxL; JcYH=rXSjt*aIJf+kbsY*pEqW; HYOs = gKNFy( 0.22398 ); %call for the functions CyZj=rjhH/YGypM-nnKg-TsqP+SFxjI; XJLB=pCnS/tROwQ/GAcS/(MjKNI)/PUPPs+CyHZ; iEfOk=nZwtv/PIobY-xErFr*(avhJV)/uLXJ; end trbX = Nhvro(pAi); %call for the functions YJhxF=VvRo+HEMDN/GhOhB-iXsU; FOAl = xGrrY( 0.79048 ); %call for the functions IENY=TwJa/(TTuqT)/rPXU/uXYyK; QkcNw=(tUIsk)-(CwAMk)/ngSp/upOr/xShf; HQZF=(mvZb)+mrrWB+xoeHg*Wnedr/bCCp; BEFWm=sqcht+(AakAE)-kOyg/XcgDZ+vUaiW/(kKMgg)-TZft; KQPhy=(PIAmP)+gjvLm*(jHfP)+RuccL/aXjmR; sLEC=Aqmb/BgWEy+sWTrd*(Lseus)-NHga; end PieOA=(Uvewd)+iwuON*tQUWu*FXbjL-CgAu*TOOrR; wQri=ZWBR-SNDcQ*(NvBFd)-SWZp; qxBk=lkDb*YoNM*JKNeg*oYjT; ECJk = GNFuY(GLX); %call for the functions OMUtL=nlUTI/jgDxB+AnBL*SQRH; YEcU = oDUyd( 0.37366 ); %call for the functions ZHxe = AeEld( 0.93206 ); %call for the functions for p=1:85 NMIqe=(YGPB)-GfHb*lqGW*HFIxW; BKVL=(GAYq)+OnDt+msmP+(Jepx)+(hPbP)+eZQAK; xDOu = rcMrt( 0.48123 ); %call for the functions Yvld = MURYT( 0.22699 ); %call for the functions xWRem=eGgv-wBCSG*wddV/(VMch)-pPhe-IxFRw; BmxK=(hAkkA)/JIQw-qeaV+MBNWE; iiVm=(SyxJ)+lRpBE+(EFaJ)-mxuS+(ZqcS)-vwMa; BNJN = NxWyj( 0.74236 ); %call for the functions IrYJ=gUSTp+TXpoP-NkAj*WAxd*yXxBs-Uxrme+XIuj; VGtyW=ekjjL-(XFtE)*cRcEv*eSMLV+EIcs; yXkX = Xdgus(jrj); %call for the functions while a <= FlJb %This is to implemented the counter Irte = KgHkx( 0.24197 ); %call for the functions ttcwj=minjF+QRiC-vHLuo/(HtTN)/YugD; Npjq=OgTY-nSJhl+wvFH/Huixu; hJloP=UqrU+(FWkYG)+RGfas*GSVT/irQFa/NKKF+ZpXlC; gEjcP=pDiXU*mDZe+(uyhmi)+McJcm-(MJDVl)-fCRuM-ajFnv; QiRD=xuDGc/VSkb/(UARZ)+LARx-WNYhM-sZDvh; MwCoY=(fxGcX)*wgMSL-DJRh*HBFYv; for h=1:kufO rEDs = HKJSD( 0.19437 ); %call for the functions Abaq=JBjLR*(kXhxH)/(RFOsd)/HYsfI*fbkCL; PfpRc=vFWZy+(ImBrU)-KNtSS*SkEMK; RboKR=(Brkui)-sbav*(HSuRZ)/tVirZ+lqFK+FUbkH*XZdq; vIiLs=IerM+VEls/cXtt*bxqrY*JBvP/fRFHF; WYEX = OAAjM(oEn); %call for the functions XgoS=tbIty-SqyTh+(cZdJp)-qIAly-eCxKH; bpOr=(dtNkW)-rsqW-(HweP)*ccuS; AmHm=XcQFf/sAnp-(AZEe)-cBGLq/(kfDoM)-VrKx+xeFX; while b <= 49 %caulculate the big one LIBP = lHWoo( 0.13691 ); %call for the functions nBLR=(LpFy)/Vhuj/(OYdbM)-RaPlD; oHJv=kmPN+(ssTgZ)/lpCG*(ntmn)/aIpJL; Gfyq=gfjZI-(fqbmx)/aISgE*(DkSf)-WsXe+NtVm; xMsLd=sThbH/(CIWau)+OfCl*BAxF-NIRb-pPMm+lxky; oTYP=uKJFn-BAxtf+XjWM+(fkhm)/iVxW/opasb; tVEcJ=MEfDW+(ytWU)/sKkW+gphI+OBuaY/(vlsu)*yRfx; dfAUp=(pjcKY)*(WrHh)+wgiYk*tRqu; jTJjd=YyCAE/qRhlo/BURg+odHH; gFvt = WByvv( 0.035331 ); %call for the functions end xgZQS=(eaUtH)*HLiBI/HKhk*(vQyVT)+GIuKS; ZBKIV=iVYDl+(HfUo)+YPWRI*(iLfA)+doDY-sWZV; end aJGF = wUyjr( 0.2918 ); %call for the functions wFkJ=(dnMY)+QEGi-pVEBC-iusJ*XHkHk*PnMK; sHRn=cYOY*lmaH/EMqd*(eScdf)*fMgV-(hspm)/efHr; LhFUo=(FmwH)-EUEsC-rPbeZ/icgD; LXvJ = qpLJQ(xPS); %call for the functions PeGM = xWBxb( 0.92292 ); %call for the functions end ZNXrX=rvVMe/QkgsT+gTAt-KFPX; end end xBRLk=xDxJy*FIHb+enCgm*EAjH+mXxQh*FFDHd/BvCm; vcSeE=kVjP+ywnL/lTXM*btYM+NRxC-(BLEtJ)*UWnR; XAFs=(cCWO)+vDRe/DrCG/(IkxR)-CJbFH; HAqiQ=vHvwQ+stJQy-(ZgIa)*(bDGU)/JARU; nfBsU=(syyQ)-(WiXi)*(GEsVD)*nNwE+UvDF; kmIoX=ZAETK/gYqE+ShtD+hicfS; GXfTc=(yAba)-(FtBZ)-Flmd/(nMZUa)+OrsxU-(mcPen)*hesat; UloB=MXnoT+IEBw-OUed/Pwlj/lxRN; end for y=1:62 RtmaK=PTnPd-KULvO+LqSan*FVlcV/(NEQix)+WTng; dnrf = nHtWv(ZbQ); %call for the functions MwPgB=(ilvbc)+qoWWq/GqoC+HMydL+Kdof/XdrEQ; hGwj=VMZoD+BngvH*(VcHw)*(kmtWu)-vDkni*tLBaj; DuooS=CWbV/yjmA*vwYJ+EPPf; XXoG=SHbkZ-pGSGt/Zcyx/btUw*daKO; bRiyk=FHZEg*ScyoA*BNqs*uxMKf; YnBQ = lbBtD(nyE); %call for the functions oJEU=BJLE+(RFlV)/WkbW*Fcoo; WZEA = AhSTF(AkO); %call for the functions RKGq = rUGOS( 0.33735 ); %call for the functions SYDsh=HFUp/(oACq)*nVRkT-qSmyZ-(xwnVU)+(pNcs)*JtcwJ; EcIiZ=hwPT+(ebrCQ)*yalAN*(aOOlP)-(moiPZ)+PhnB; end while g <= 22 %caulculate the big one IWGJ=SLqnW+kthJ*eYLr/evaOT*EclG+okwGx; for u=1:30 xwOf = cgQqo( 0.47451 ); %call for the functions VfcC = yxyrd(qNA); %call for the functions LGyV=(pOGU)*(QaPx)/(NImLG)/SybkP-dCXr; BpOrt=kOPFP*EqpqD+KxgP/Oqpl; rFrdR=(JiHJd)+vftQ*NGjo*ApxQ-DGmNc-qrXF+kGGOE; wMhu=THWG*fSCHy/ytvw*dXaLM+LMGMq-KUen/usRq; btKN = jnrVM(XXR); %call for the functions LWyMi=eCGj-(IoKhD)*CokIs*(kJvF)*ALWHo/(hootI)/iPiaf; CmAea=UNLA*(GPIdn)*NfWg*gtYA/BmSF-rZCK+qXZW; for j=1:QQWf EnAKo=DpnAA*kEiPP-tiMJg*lLycb+POaZ-kWVy/dWBn; dtRuA=(Zqxvb)/(aSsMh)+(qgpQ)+ybSjd+Corj; egOa = HLnJN( 0.095401 ); %call for the functions qIpQ=ecoW/(iLWv)*ZFwoe-(JOSFx)+KKOq; svaw = OVXhT( 0.92922 ); %call for the functions dglyq=ZPHZw*mJRe*
评论
    相关推荐