t_alpha.m

我认为很不错的语音处理的ｍａｔｌａｂ源代码

            alpha(set_t)=1./2+1./2.*t./(4+t.^2).^(1./2).*(1+1./2./(1+1./4.*t.^2));

        case 5;
            alpha(set_t)=1./2+1./pi.*atan(1./5.*t.*5.^(1./2))+1./pi.*t.*5.^(1./2)./(5+t.^2).*(1+2./3./(1+1./5.*t.^2));

        case 6;
            alpha(set_t)=1./2+1./2.*t./(6+t.^2).^(1./2).*(1+1./2./(1+1./6.*t.^2)+3./8./(1+1./6.*t.^2).^2);

        case 7;
            alpha(set_t)=1./2+1./pi.*atan(1./7.*t.*7.^(1./2))+1./pi.*t.*7.^(1./2)./(7+t.^2).*(1+2./3./(1+1./7.*t.^2)+8./15./(1+1./7.*t.^2).^2);

        case 8;
            alpha(set_t)=1./2+1./2.*t./(8+t.^2).^(1./2).*(1+1./2./(1+1./8.*t.^2)+3./8./(1+1./8.*t.^2).^2+5./16./(1+1./8.*t.^2).^3);

        case 9;
            alpha(set_t)=1./2+1./pi.*atan(1./9.*t.*9.^(1./2))+1./pi.*t.*9.^(1./2)./(9+t.^2).*(1+2./3./(1+1./9.*t.^2)+8./15./(1+1./9.*t.^2).^2+16./35./(1+1./9.*t.^2).^3);

        case 10;
            alpha(set_t)=1./2+1./2.*t./(10+t.^2).^(1./2).*(1+1./2./(1+1./10.*t.^2)+3./8./(1+1./10.*t.^2).^2+5./16./(1+1./10.*t.^2).^3+35./128./(1+1./10.*t.^2).^4);

        case 11;
            alpha(set_t)=1./2+1./pi.*atan(1./11.*t.*11.^(1./2))+1./pi.*t.*11.^(1./2)./(11+t.^2).*(1+2./3./(1+1./11.*t.^2)+8./15./(1+1./11.*t.^2).^2+16./35./(1+1./11.*t.^2).^3+128./315./(1+1./11.*t.^2).^4);

        case 12;
            alpha(set_t)=1./2+1./2.*t./(12+t.^2).^(1./2).*(1+1./2./(1+1./12.*t.^2)+3./8./(1+1./12.*t.^2).^2+5./16./(1+1./12.*t.^2).^3+35./128./(1+1./12.*t.^2).^4+63./256./(1+1./12.*t.^2).^5);

        case 13;
            alpha(set_t)=1./2+1./pi.*atan(1./13.*t.*13.^(1./2))+1./pi.*t.*13.^(1./2)./(13+t.^2).*(1+2./3./(1+1./13.*t.^2)+8./15./(1+1./13.*t.^2).^2+16./35./(1+1./13.*t.^2).^3+128./315./(1+1./13.*t.^2).^4+256./693./(1+1./13.*t.^2).^5);

        case 14;
            alpha(set_t)=1./2+1./2.*t./(14+t.^2).^(1./2).*(1+1./2./(1+1./14.*t.^2)+3./8./(1+1./14.*t.^2).^2+5./16./(1+1./14.*t.^2).^3+35./128./(1+1./14.*t.^2).^4+63./256./(1+1./14.*t.^2).^5+231./1024./(1+1./14.*t.^2).^6);

        case 15;
            alpha(set_t)=1./2+1./pi.*atan(1./15.*t.*15.^(1./2))+1./pi.*t.*15.^(1./2)./(15+t.^2).*(1+2./3./(1+1./15.*t.^2)+8./15./(1+1./15.*t.^2).^2+16./35./(1+1./15.*t.^2).^3+128./315./(1+1./15.*t.^2).^4+256./693./(1+1./15.*t.^2).^5+1024./3003./(1+1./15.*t.^2).^6);

        case 16;
            alpha(set_t)=1./2+1./2.*t./(16+t.^2).^(1./2).*(1+1./2./(1+1./16.*t.^2)+3./8./(1+1./16.*t.^2).^2+5./16./(1+1./16.*t.^2).^3+35./128./(1+1./16.*t.^2).^4+63./256./(1+1./16.*t.^2).^5+231./1024./(1+1./16.*t.^2).^6+429./2048./(1+1./16.*t.^2).^7);

        case 17;
            alpha(set_t)=1./2+1./pi.*atan(1./17.*t.*17.^(1./2))+1./pi.*t.*17.^(1./2)./(17+t.^2).*(1+2./3./(1+1./17.*t.^2)+8./15./(1+1./17.*t.^2).^2+16./35./(1+1./17.*t.^2).^3+128./315./(1+1./17.*t.^2).^4+256./693./(1+1./17.*t.^2).^5+1024./3003./(1+1./17.*t.^2).^6+2048./6435./(1+1./17.*t.^2).^7);

        case 18;
            alpha(set_t)=1./2+1./2.*t./(18+t.^2).^(1./2).*(1+1./2./(1+1./18.*t.^2)+3./8./(1+1./18.*t.^2).^2+5./16./(1+1./18.*t.^2).^3+35./128./(1+1./18.*t.^2).^4+63./256./(1+1./18.*t.^2).^5+231./1024./(1+1./18.*t.^2).^6+429./2048./(1+1./18.*t.^2).^7+6435./32768./(1+1./18.*t.^2).^8);

        case 19;
            alpha(set_t)=1./2+1./pi.*atan(1./19.*t.*19.^(1./2))+1./pi.*t.*19.^(1./2)./(19+t.^2).*(1+2./3./(1+1./19.*t.^2)+8./15./(1+1./19.*t.^2).^2+16./35./(1+1./19.*t.^2).^3+128./315./(1+1./19.*t.^2).^4+256./693./(1+1./19.*t.^2).^5+1024./3003./(1+1./19.*t.^2).^6+2048./6435./(1+1./19.*t.^2).^7+32768./109395./(1+1./19.*t.^2).^8);

        case 20;
            alpha(set_t)=1./2+1./2.*t./(20+t.^2).^(1./2).*(1+1./2./(1+1./20.*t.^2)+3./8./(1+1./20.*t.^2).^2+5./16./(1+1./20.*t.^2).^3+35./128./(1+1./20.*t.^2).^4+63./256./(1+1./20.*t.^2).^5+231./1024./(1+1./20.*t.^2).^6+429./2048./(1+1./20.*t.^2).^7+6435./32768./(1+1./20.*t.^2).^8+12155./65536./(1+1./20.*t.^2).^9);

        case 21;
            alpha(set_t)=1./2+1./pi.*atan(1./21.*t.*21.^(1./2))+1./pi.*t.*21.^(1./2)./(21+t.^2).*(1+2./3./(1+1./21.*t.^2)+8./15./(1+1./21.*t.^2).^2+16./35./(1+1./21.*t.^2).^3+128./315./(1+1./21.*t.^2).^4+256./693./(1+1./21.*t.^2).^5+1024./3003./(1+1./21.*t.^2).^6+2048./6435./(1+1./21.*t.^2).^7+32768./109395./(1+1./21.*t.^2).^8+65536./230945./(1+1./21.*t.^2).^9);

        case 22;
            alpha(set_t)=1./2+1./2.*t./(22+t.^2).^(1./2).*(1+1./2./(1+1./22.*t.^2)+3./8./(1+1./22.*t.^2).^2+5./16./(1+1./22.*t.^2).^3+35./128./(1+1./22.*t.^2).^4+63./256./(1+1./22.*t.^2).^5+231./1024./(1+1./22.*t.^2).^6+429./2048./(1+1./22.*t.^2).^7+6435./32768./(1+1./22.*t.^2).^8+12155./65536./(1+1./22.*t.^2).^9+46189./262144./(1+1./22.*t.^2).^10);

        case 23;
            alpha(set_t)=1./2+1./pi.*atan(1./23.*t.*23.^(1./2))+1./pi.*t.*23.^(1./2)./(23+t.^2).*(1+2./3./(1+1./23.*t.^2)+8./15./(1+1./23.*t.^2).^2+16./35./(1+1./23.*t.^2).^3+128./315./(1+1./23.*t.^2).^4+256./693./(1+1./23.*t.^2).^5+1024./3003./(1+1./23.*t.^2).^6+2048./6435./(1+1./23.*t.^2).^7+32768./109395./(1+1./23.*t.^2).^8+65536./230945./(1+1./23.*t.^2).^9+262144./969969./(1+1./23.*t.^2).^10);

        case 24;
            alpha(set_t)=1./2+1./2.*t./(24+t.^2).^(1./2).*(1+1./2./(1+1./24.*t.^2)+3./8./(1+1./24.*t.^2).^2+5./16./(1+1./24.*t.^2).^3+35./128./(1+1./24.*t.^2).^4+63./256./(1+1./24.*t.^2).^5+231./1024./(1+1./24.*t.^2).^6+429./2048./(1+1./24.*t.^2).^7+6435./32768./(1+1./24.*t.^2).^8+12155./65536./(1+1./24.*t.^2).^9+46189./262144./(1+1./24.*t.^2).^10+88179./524288./(1+1./24.*t.^2).^11);

        case 25;
            alpha(set_t)=1./2+1./pi.*atan(1./25.*t.*25.^(1./2))+1./pi.*t.*25.^(1./2)./(25+t.^2).*(1+2./3./(1+1./25.*t.^2)+8./15./(1+1./25.*t.^2).^2+16./35./(1+1./25.*t.^2).^3+128./315./(1+1./25.*t.^2).^4+256./693./(1+1./25.*t.^2).^5+1024./3003./(1+1./25.*t.^2).^6+2048./6435./(1+1./25.*t.^2).^7+32768./109395./(1+1./25.*t.^2).^8+65536./230945./(1+1./25.*t.^2).^9+262144./969969./(1+1./25.*t.^2).^10+524288./2028117./(1+1./25.*t.^2).^11);

        case 26;
            alpha(set_t)=1./2+1./2.*t./(26+t.^2).^(1./2).*(1+1./2./(1+1./26.*t.^2)+3./8./(1+1./26.*t.^2).^2+5./16./(1+1./26.*t.^2).^3+35./128./(1+1./26.*t.^2).^4+63./256./(1+1./26.*t.^2).^5+231./1024./(1+1./26.*t.^2).^6+429./2048./(1+1./26.*t.^2).^7+6435./32768./(1+1./26.*t.^2).^8+12155./65536./(1+1./26.*t.^2).^9+46189./262144./(1+1./26.*t.^2).^10+88179./524288./(1+1./26.*t.^2).^11+676039./4194304./(1+1./26.*t.^2).^12);

        case 27;
            alpha(set_t)=1./2+1./pi.*atan(1./27.*t.*27.^(1./2))+1./pi.*t.*27.^(1./2)./(27+t.^2).*(1+2./3./(1+1./27.*t.^2)+8./15./(1+1./27.*t.^2).^2+16./35./(1+1./27.*t.^2).^3+128./315./(1+1./27.*t.^2).^4+256./693./(1+1./27.*t.^2).^5+1024./3003./(1+1./27.*t.^2).^6+2048./6435./(1+1./27.*t.^2).^7+32768./109395./(1+1./27.*t.^2).^8+65536./230945./(1+1./27.*t.^2).^9+262144./969969./(1+1./27.*t.^2).^10+524288./2028117./(1+1./27.*t.^2).^11+4194304./16900975./(1+1./27.*t.^2).^12);

        case 28;
            alpha(set_t)=1./2+1./2.*t./(28+t.^2).^(1./2).*(1+1./2./(1+1./28.*t.^2)+3./8./(1+1./28.*t.^2).^2+5./16./(1+1./28.*t.^2).^3+35./128./(1+1./28.*t.^2).^4+63./256./(1+1./28.*t.^2).^5+231./1024./(1+1./28.*t.^2).^6+429./2048./(1+1./28.*t.^2).^7+6435./32768./(1+1./28.*t.^2).^8+12155./65536./(1+1./28.*t.^2).^9+46189./262144./(1+1./28.*t.^2).^10+88179./524288./(1+1./28.*t.^2).^11+676039./4194304./(1+1./28.*t.^2).^12+1300075./8388608./(1+1./28.*t.^2).^13);