test_kernel.asv

功能是MATLAB 语音文件存为数组的形式

win_type = 'hamming';
n = 512;    
m = 0:n-1;
switch(lower(win_type))
    case 'rectangle'
        win = randsrc(1, length(m), [1]);
    case 'hamming'
        win = 0.54 - 0.46 * cos((2* pi * m)/(n - 1));
    case 'blackman'
        win = 0.42 - 0.5 * cos((2 * pi * m)/(n - 1)) + 0.08 * cos((4 * pi * m)/(n - 1));
    case 'hanning'
        win = 0.5 * (1 - cos((2 * pi * m) / (n - 1)));
    case 'barlett'
        m = floor((n - 1)/2);
        win = zeros(1, n);
        win(1:m + 1) = 2 * [0:m] /(n -1);
        win(m+2:n) = 2 - 2 * [m + 1:(n -1)] /(n - 1);
    otherwise
        fprintf(1, 'unknown window, operation is terminated!\n');
        return;
end;
for k = 1:7
    for j = 0.1:0.05:0.45
        jj =0.45
        h = zeros(1, n);
		h(1) = 1;
		h(8) = 0.45;
        h(8 + k) = j;
		hh = h .* win;
		hh = abs(ifft(log(abs(fft(hh)))));
		atten(j) = hh(8);
    end;
    [val, idx] = max(atten);
    max_atten(k, 1) = val;
    max_atten(k, 2) = 0.1 + (idx - 1) * 0.05;
end;
[val, idx] = max(max_atten(:, 1));
fprintf(1, 'max complex:%f, echo:%d, atten:%f\n', val, 7 + idx, max_atten(idx, 2));
% h = zeros(1, n);
% h(1) = 1;
% h(8) = 0.45;
% hh = h .* win;
% hh = abs(ifft(log(abs(fft(hh)))));
% hh(8)