1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
|
function [yf, yt, ydiff] = tolstoy(x, h, L)
% compute the convolution between two signals using the time-partitioned
% overlap-save algorithm.
% x -- signal (vector)
% h -- impulse response (vector)
% L -- fft size (IR frames are L/2, signal stepping is L/2 with 50% overlap)
% yf -- result using frequency domain convolution
% yt -- result using time domain convolution
% ydiff -- error between yf and yt
%
% @author gregor heinrich :: arbylon
% @date Nov./Dec. 2006
% TODO: time-varying impulse responses.
% TODO: trade-off calculation for frame sizes.
close all
% standard values
if nargin==0
L = 2^5;
x=ones(20,1);
h=ones(3,1);
%breaks truncation rule
L = 2^5;
x=ones(112,1);
h=ones(38,1);
end
x = x(:);
h = h(:);
fprintf('signal length: %i, IR length: %i\n', length(x), length(h));
fprintf('1. convolution in time domain.\n');
tic
yt = conv(x,h);
t1 = toc;
%plot(yt);
%title('time-domain convolution y_t');
%figure
K = L / 2;
fprintf('2. convolution in frequency domain.\n');
fprintf('signal frame length: %i / stepping: %i, IR frame length: %i\n', L, K, K);
tic
yf = tols(x, h, L);
t2 = toc;
fprintf('comparison of calculation time: time-domain: %gs, freq-domain: %gs.\n', t1, t2);
plot(yf);
title('frequency-domain convolution y_f');
%figure
len = min([length(yt),length(yf)]);
ydiff = yf(1:len)-yt(1:len);
%plot(ydiff);
%title('error plot y_t - y_f');
error = sum(ydiff.^2);
fprintf('total squared error: %g.\n', error);
fprintf('done.\n');
end % test
|
