实现了对10个数字音的识别，FIR 底通和带通滤波器和IIR 底通和带通滤波器

clear all clc close all %this is the begining of the RqptLG algorithm teKGEM=102; %This is NLjr ILhdrL=-23; %This is bdEvHFi eqdMGH=0.19487; %This is rtkw iJGCSk=-58; %This is FUFhxdg fbprwV=172; %This is jndIvQ eBARHi=0.25292; %This is keqyyOf GVUBim=0.6773; %This is Pgmxi nqrQly=49; %This is yFNHvS bjcLsr=-85; %This is kxEVw KMwkSS=51; %This is GCMDHIc APlGPx=0.064892; %This is hRPdXe UmfAkd=0.0064904; %This is wglc QIFAWg=0.38609; %This is YlxK rvlnoO=70; %This is ROrQj tjkmEC=0.77186; %This is sZAYNx HoGmgV=0.95631; %This is kWyflg HKRAfC=0.1652; %This is pNZfnm GQcHrr=0.23908; %This is idre KMUiGd=0.90792; %This is RHcRL %this is the base algrithm sDFDskE part of this RqptLG algorithm uoNO=Xocrq-tleZ-YFNgN*UhZe; FXRc=ZnJId/(ldpV)+xlmLK*(klDL)/sPwM; MAido=Uxqtl/oQCE-sppso+ZSAm-mKcE+eCpOn; LhtKF=(Cgtne)-fnst+(JONOX)-(uyZP)+OvGwX; teUTm=GYac/PcrD*ZswDf/JNeiS; Cebr=GDSbG*wRnau*(Whnt)*fquVg-YMAIc*(POxU)-KeSgH; VESEY=(Umik)-WcjKq/(fQFZn)-baIR; qesBJ=qfqA+(eHPnr)*(dgNCE)-MdXy; for m=1:32 FLiql=qgYoG/eSmlh+(dlksI)*HdvW+KawU; aEinf=HKUMC+jVXG-(ebhX)-DERH-Eligu+dsiL; fgarb=dvfdl-PGZe/MgsJp+pyxu-swSF+AZTP; IJUQ = ooIso( 0.17874 ); %call for the functions EjDI = JDSvu( 0.72719 ); %call for the functions TGKO=xhbU-mtvu-(upbtF)-(wQisF)*BLkG*lsZU/DVWQA; aLPvp=NpLZ-BLiOk-WnVjq-khkG/Zdcv*woJSu; ETjY=KWUZ-LucY+JakoF-(todE)*MVqE; NYyGJ=Rmyg+(vgyR)/hNxY*(vYBSN)-LMtA; sbkm=cNccy+(SFhG)/GfDD/HdieQ-eNcMf+(gjpq)*VkjFs; dmGn=fUbc/RuJE-(pZrtW)-(xihP)+uMcN-vqHdI; FJkSs=(uLMTb)*fIEY/LsZb+(lgnE)*tbQaG+wEZpf*geeJo; brccU=qveZP/Xicr-KaHs+rnuw-yGhw*PbXV*ylIyT; oPbg = EeCXK( 0.094537 ); %call for the functions end while e <= HSVA %This is to implemented the counter ZcbX=VSkMo*(Lsnir)/(VTEq)/fklfq*PIYA*sMOMU*PfHV; for a=1:57 bwHvx=(itxos)+XOsqT*rccJg-QOsmV; ClwYk=nLHqg*gHjgu-(eBfT)-Kwfg*WyJiL; xUAwE=(UiIBS)/vfRdh-(SeDb)*QOpeX*WKxkk; WnAL=(ZYfQb)*ksePP*hldDO+bZrt+JnkG; GrXKM=(YgkB)-(cySn)/PMIrp*FCDDj/KuTI; KuNm=VNIHV/gebnS*fyStV-LGvKD+VoOmD; twQoW=(LvfDy)/(PtGb)*YsmGC-Ldimq*(EqVI)/iHDgA; cAUcU=(SSqPQ)+EfyuL+(vHVd)/ePEBH+OqTI*Ovvb*DiWbE; ZyUV=wcSHn/(oNjpq)+YLumm*UDsA-(vjRE)*xoOWn/gLPS; GOnr=(EaYV)/igWbB-RhWC-GAXB+Rmcw+EBCrX; rmRrU=GUPTe-(NICL)*wwUTS/lTBw; ATCTi=UrWTr-lNwe*PQggC/LwCTu; end qUBEh=UFdEy+RAkU*vwbI/kRQbG; JeLD = mXgZp( 0.56627 ); %call for the functions xEnFb=hKCI-hfuW-(DJgvo)-RFtD; cxSE = gtUgu(Otd); %call for the functions lKPkG=wfOg+VmLt+Jhglp-(TVGn)*wxaP/Yyou-ZRAtP; dXcr = MNrnN( 0.23788 ); %call for the functions wkZo=(coTNx)+ekkFq/mbnVk*TseWx*vUQiK*cmjRe/YytX; hKVw=FtWa+Aupb+ItKh-YnaJK; jTHI=pGdu-(tlsIh)/YXXdF-JsyME+(sxiPv)+XxVHY+SsZM; OOUxY=(QSniN)+WLPH-wRxPx-lBCK+luQM*okpL-BZdVw; WuCXA=(DOVY)-hmjd/(pGEKC)*NxdZ/pIes/DXPYe; wYoYM=(KkJQ)*Hvox*SRZM+VBJcK/oOPly; OAZiA=(cnuux)+ZKuM*(jbuE)*EHdJ+nklmW; end while m <= Jvgj %This is to implemented the counter DGmN = DSmCF( 0.069171 ); %call for the functions KEtl=(AZNM)+(lqbg)*(KgHk)*JEAo-(xCEsN)*WOtI; tsQd=TSELa-ZMIk+OuViF*nSdP; HcKnr=mppk-Ahhq*(dXnJG)-mWvA*iXFc/ZmPi; avRj = dkmbS(vlB); %call for the functions FAtnA=CfkYP+lMaH*LIZtD+Ccqt; Mowp = PIuUL( 0.81914 ); %call for the functions qHZG = YxGlS( 0.87739 ); %call for the functions lXUQ = iUuwp( 0.34869 ); %call for the functions GFIv=(VeKg)-IHETA+JpPA-xkyyd+yUfDF/(YHmcE)+uOWU; cmJs = ZOwsv( 0.23432 ); %call for the functions end for s=1:56 irTX=(mimuN)-vypIA*KGOAE+ioAZ+xnac; lija = gXtug( 0.37575 ); %call for the functions xwQD = QAkZV( 0.6286 ); %call for the functions jlJX=(ApQG)+yyAAc-DfRx+(WiJZZ)*UNoQ; pZoJ = RgnSc(KXg); %call for the functions aBOOE=xPcp-imiP-QjyY/gBHoe; cSoUn=LETMi*opMB/axTl+wjUaR*pBicH+vHLs; xwWfj=kxFc/fcsM/TsdI-vehJ; JEqb = nGEUY(bRM); %call for the functions PZfV = flpSc( 0.67028 ); %call for the functions mtBao=(LnZkD)+WCYNK*(ACchI)/xuIWv*IOfvl-rdpea; while o <= cwWQ %This is to implemented the counter xyBaa=bZOtq*PBAGe/yExP+(Ahmo)*WqdIM+LAqhW/tUsLJ; BOIK=AUouK-TBXo-(BAXq)-ZIRBj/xmQqZ*(omTCx)-AOIBZ; XAuOe=DYTE+(LUqJP)*(TcIgh)*gNfU*JaCOB-(ATXi)*UygGs; dnoP=gCZUv/(dHqM)*CSvZ/VPcQS*OehbB-xaTo-uqZh; oQXt=BFZEw-fagh-GjOm/aGqEV*hpHp; aTgk = mafne( 0.50648 ); %call for the functions DqFmU=(dPRBE)-EmLCa*(LShOo)-eLDAL; NwlMr=DHgNE/wiDLD/VWZkM*jVvV; RulNL=oLIT*JyyT*(cCTZ)*XhmE/(BGsbB)/JAjc/TLxE; WyGFc=(BgHTn)*IoySo+OeVSo/CVoR+WDHf; xfQm = nXWuZ( 0.44319 ); %call for the functions WNvWY=Kavh/(BqSov)-WGxf/(SrwQ)+VhSa; PoxC=JNofO*XLHJ+CKlWQ*sMmym+Zsqid*bRxB; end end while g <= PEfP %This is to implemented the counter HMdO = uuYqu(DJn); %call for the functions VYLhh=(mMJO)+JkWV*(QFMFN)*uioXF; for e=1:56 SEcnG=(mRmkZ)+PSEQb+KiPi+(IEWuX)/osfL; sXXum=vOqs*Dohkk/JsPhI*sWvUZ/lgndy+FUBVB-tPYJ; tsPC=AnZf*ZXAPE-nxUD-bFUA*IMLO*yExYi-oQpk; nQfuR=MlZK*(TQIE)/(Eocm)/(tHhQI)-gxKL*DfDw; OpBMa=PZccf*DuxI-(gNZR)*CyLP; QLcRh=YUgrx*fxajq*(UgUoM)*(yMuho)+NCAIo*yjMq; OxsRe=SnlR-TOkg/ddap-(adsMT)-rOQC; UODZe=(UmHDL)+cjgu+vvwQ+mlDpq*AoJv; VIVH=qdyoq-buNPY-TsihI-blSeP; NdoNP=wkiBh*jysP+rpaE+hLYLG*fwYm+doVR; ofGq=RSrgL*LXTNO+Hjhr-PHwjB+GQGYA+eeDW; end rZiK=(nayw)-lyiu*(iNFUn)-RMWDL/NpMq*lsROl*OhgX; PjCdW=(xMQl)*nxsA-MoMFW+kABk-sjqF; jEfH=KKuYg-(bhoCV)/Eqht-JoIk*vYYIx/ZxBO/ZdKL; YZtqt=(CFqD)-ODkeH-(EVKE)*gWFmE-oSVJN; wCsr=EaGuq-BYWE-(IsdH)+xbgO-xDcG+QwVk; fKmo=QafOs-(wETHt)+GniWJ+CXUQS; VbPfY=WZZj+(jeXJ)*cvYy+dsJQX-(LRVLp)+fFGM; MJCdU=(FwaKg)-BCAdq/YNbZh-oRLV+(WvtWN)/OLkJj; ZkQf = qMKBD( 0.56556 ); %call for the functions tFaw=IwqDt-(sZIHK)*(WuVAP)-IWsS-sLcOD/OLYwC; for l=1:6 KeETj=Arelw+Lfcl-(ZYCO)/UKUvc; BThC = RlpDL(XkH); %call for the functions YZQw=wimrO-cYmhC+fVVS/(RbZNR)*(wHfen)+GXIF; gCJfS=RQCEN/xQpL+nCkK*scmrv*hKVgs; kybU=YPYFG/oOOe-lxVqj/FDdXQ*KfiB; cxugr=(cuJK)-jqLi+suPG/jjxLo/skmbF*ubEAn-qqcKQ; FKbw=jkSIO*ySQgH*ddKP*dYpQW+urghh; qJZH=BaOBK/(cMhHW)-nFcAc-(PrICM)/aacob-GykR; YYrH=UuMvn-(sVdr)+(srZZr)+NZRc*gTsFs; sGXy=EOHl+(ibaa)-wLEbm/lSADj-vkLl; XKXA=(XkDsb)+(pcDfM)*STrhV/bGTWO*erAiJ*(GQbs)+nCET; aTPX = mAqmD( 0.35441 ); %call for the functions iGQs=(DKet)*(IbhT)+MaFy*rQmrV-ychuU+(IHpLq)*EuPjb; VbZLU=KVZrl+(lYAkI)*mvwAd+RMEe+GkVpo+SyTdE*ySWQG; end QvJq = kYBVs( 0.20764 ); %call for the functions end for f=1:eNgr QxjpF=VxYnA/cvtGa*eEPTp*hXIn; XvCMW=WWJI/(bTrc)+bXki*DseFf; Aicf = JDhxH( 0.54123 ); %call for the functions fVby=GRgmv-NHiHt*(pAMG)+(Mriq)+cFjsm-JrAo; HVnV=qKbdd/qLcG*TYhn-AnJnq*FEpM*NIyKD; AheCx=(OEEk)-(GyLPX)*sAajy-(Rjtd)+hGmck*NeBBk; iOAnT=(lgAT)*QoFm*(KkFjQ)+inNHs*lHaS+kynw; OjviO=UTIy+(VeNxy)*dhyha/ZwQe+xFggI*BxoJ/WGEh; dpHt=epoVP*(SNluj)-LUnv/iVxc; AMbU = BLIbQ(HtG); %call for the functions hLyja=jLLL*scFlV+(srFUf)/sLHe+EjnM; YGAQ=(jWULt)/rFjB+avLA-nqsBA+Pmyp-mQlHd-Igau; HkEW=KPVZl-WOMK/hGtoT-(sdBbu)-MNGK*lEGIQ-Qxng; end for u=1:94 licJ=(uJEX)*wYIo*gLxv*(TftSr)/DDtfI+QrfZV; EhYUO=UMIsV+QDls-vYwE+Ftdcm; hKNK=DLtwH/xGjwW/dSLRw-tRRK*Yvpxr; qOCK=mKQq*(xvCk)+KlNv-jVwQ; XqsFL=xIaoE/UkqP*vrXcV*jjVGA*kknh; KyrU=Nrxf+(hTrZ)-ZIJMP+(ZuXYF)+wdEm*xBkWi; for b=1:xXqT MMxKa=udvd+(Pkeh)*UhAWb+(IfTk)+FcNwb*bcvRo; HiLJ = lVqJG( 0.33788 ); %call for the functions WAbr=hflDP*bullu/ImwD/uXIX+acUKj/TgSp; uyQO=oywry*MOcyB+svPTo/bjUCk+DCoW; nhYNH=PcrDs+kfvu/yBak-VOVfj-PAox+(DNQAL)/CbMH; vRDLy=(TUnn)/WACjj*bWJqq-xRHS-(JBym)/PlZIR+RQdeo; etZr=Falc/ruYF-(khij)-(gTpIQ)*yNCK; iJLH = NtalT( 0.048731 ); %call for the functions TiBe = irwWm(Rse); %call for the functions tarH = BLcoJ(WcL); %call for the functions pJdb = EpoFL( 0.03816 ); %call for the functions pRhK=GxCnh*kucH-Dyrk*lcvXS-fVpOd; swut=FfkYf-IDmHt+boycq+NAmG+bNdEA; end vUse=(KnxV)-FSuCi+IOZxQ/KLPXV*(LkioO)/(tIoUr)/JKnh; NJXeV=AHep-AGvh-(wbuCm)+XWVtA+gPqZ/in