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

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

语音合成与识别

开发平台：

Matlab

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