bfre2.m

里面囊括了基于matlab滤波器设计的各种.m文件

 
% Y = bfre2(h0,h1,g0,g1,L,l)
% BFRE ---- Branch Frequency Response Evaluation(2 Channel)
% Each part impulse response of system is h0,h1,g0 and
% g1. We produce some sinusoidal signals as input with
% the length L. There are 1024 sinusoidal signals
% altogether. The output Y is a matrix, whose columns
% are FFTs of output signals. The first column is lowest
% frequency output and the last column is highest
% frequency output. At the end, the function will produce
% a contour drawing for the input output mapping.
% The function PLAYB is used to play back the contour.
% L will be more than 200, l is between 2 and 10.

function Y = bfre2(h0,h1,g0,g1,M,N,L)

omega = linspace(0,pi,1024);
omega = omega(:);
n = linspace(0,L,L+1);
input = cos(omega*n);

for I = 1:1024
 x = input(I,:);
 Y(:,I) = sysrep(h0,h1,g0,g1,M,N,x);
end
%Y = abs(fft(Y));

%playdb(Y,l)