|
原帖由 tingxin81 于 2008-1-27 23:21 发表
根据‘songzy41’提供的资源,编辑M文件如下:
clear
clc
close all hidden
format long
waveFile='sunday.wav';
[speech, sf, nbits]=wavread(waveFile);
index1=4000;
nfft=1000;
index2=index1+nfft-1; ...
你的第一个问题虽然已经解决,但是为什么要试验呢,MATLAB HELP文件中已经明确说明了rceps的算法:
Algorithm
-
rceps is an M-file implementation of algorithm 7.2 in [2], that is,
y = real(ifft(log(abs(fft(x)))));
-
Appropriate windowing in the cepstral domain forms the reconstructed minimum phase signal:
w = [1;2*ones(n/2-1,1);ones(1-rem(n,2),1);zeros(n/2-1,1)];
ym = real(ifft(exp(fft(w.*y))));
而且MATLAB中还有另外一个求倒频谱的函数cceps()
[ 本帖最后由 uestclzl 于 2009-12-21 03:35 编辑 ] |
评分
-
1
查看全部评分
-
|