开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > 语音合成与识别 > 查看源码
levension.asv
资源名称：speech.rar　[点击查看]
上传用户：ay_070428
上传日期：2014-12-04
资源大小：11427k
文件大小：2k
源码类别：

语音合成与识别

开发平台：

Matlab

levension.asv：源码内容
  1. s1=wavread('ljb');
  2. s2=s1(13000:15000);
  3. s=s2(800:1055);
  4. d=hanning(256);
  5. w=s.*d;
  6. t=fft(w,256);
  7. h=fftshift(t);
  8. nn=log(abs(h));
  9. gg=nn(128:255);
  10. p=4;
  11. for j=1:p
  12.     r(j)=w(j:256)'*w(1:256-j+1);
  13. end
  14. [A E K]=levinson(r,p);
  15. a1=fft(A,256);
  16. hh=log(p*E)-log(abs(a1));
  17. hh1=hh(1:128);
  18. p=8;
  19. for j=1:p
  20.     r(j)=w(j:256)'*w(1:256-j+1);
  21. end
  22. [A E K]=levinson(r,p);
  23. a1=fft(A,256);
  24. hh=log(p*E)-log(abs(a1));
  25. hh2=hh(1:128);
  26. f=[0:127]*8000/256;
  28. p=12;
  29. for j=1:p
  30.     r(j)=w(j:256)'*w(1:256-j+1);
  31. end
  32. [A E K]=levinson(r,p);
  33. a1=fft(A,256);
  34. hh=log(p*E)-log(abs(a1));
  35. hh3=hh(1:128);
  36. p=16;
  37. for j=1:p
  38.     r(j)=w(j:256)'*w(1:256-j+1);
  39. end
  40. [A E K]=levinson(r,p);
  41. a1=fft(A,256);
  42. hh=log(p*E)-log(abs(a1));
  43. hh4=hh(1:128);
  44. p=20;
  45. for j=1:p
  46.     r(j)=w(j:256)'*w(1:256-j+1);
  47. end
  48. [A E K]=levinson(r,p);
  49. a1=fft(A,256);
  50. hh=log(p*E)-log(abs(a1));
  52. hh5=hh(1:128);
  54. p=30;
  55. for j=1:p
  56.     r(j)=w(j:256)'*w(1:256-j+1);
  57. end
  58. [A E K]=levinson(r,p);
  59. a1=fft(A,256);
  60. hh=log(p*E)-log(abs(a1));
  61. hh6=hh(1:128);
  63. subplot(3,3,1)
  64. plot(f,gg);
  65. hold on
  66. plot(f,hh1)
  67. xlabel('Frequency/Hz');
  68. ylabel('Magnitude/dB');
  69. text(3000,3,'p=4');
  70. subplot(3,3,2)
  71. plot(f,gg);
  72. hold on
  73. plot(f,hh2)
  74. xlabel('Frequency/Hz');
  75. ylabel('Magnitude/dB');
  76. text(2500,3,'p=8');
  77. subplot(3,3,3)
  78. plot(f,gg);
  79. hold on
  80. plot(f,hh3)
  81. xlabel('Frequency/Hz');
  82. ylabel('Magnitude/dB');
  83. text(2500,3,'p=12');
  84. subplot(3,3,4)
  85. plot(f,gg);
  86. hold on
  87. plot(f,hh4)
  88. xlabel('Frequency/Hz');
  89. ylabel('Magnitude/dB');
  90. text(2500,3,'p=16');
  91. subplot(3,3,5)
  92. plot(f,gg);
  93. hold on
  94. plot(f,hh5)
  95. xlabel('Frequency/Hz');
  96. ylabel('Magnitude/dB');
  97. text(2500,3,'p=20');
  98. subplot(3,3,6)
  99. plot(f,gg);
  100. hold on
  101. plot(f,hh6)
  102. xlabel('Frequency/Hz');
  103. ylabel('Magnitude/dB');
  104. text(2500,3,'p=30');
  111.