对比滤波器的四种实现结构：直接I,直接II,直接I转置，直接II转置

% 对比滤波器的四种实现结构
% 直接I,直接II,直接I转置，直接II转置
clc; clear; close all;

fs = 1000;
fc = 35;        % Low pass -3dB cutoff freq.
w = fc/(fs/2);
[b,a] = butter(2,w);
%filtinfo(b,a,2048,fs);

x =  randn(1,100000000);
tic;
y_filter = filter(b,a,x);
toc;

plot(x); hold on; plot(y_filter);
tic;
y_1 = filt_struct_direct_I(b,a,x);
toc;

tic;
y_2 = filt_struct_direct_II(b,a,x);
toc;

tic;
y_1t = filt_struct_direct_I_transpose(b,a,x);
toc;

tic;
y_2t = filt_struct_direct_II_transpose(b,a,x);
toc;
return;
figure;
plot(y_2 - y_filter);
title('error, direct II');
figure;
plot(y_1 - y_filter);
title('error, direct I');
figure;
plot(y_1t - y_filter);
title('error, direct I transpose');
figure;
plot(y_2t - y_f