%-------------------------------------------------------------------------
% exa100802_tfe.m,  for example 10.8.2,
% to test tfe.m  ;
%-------------------------------------------------------------------------
clear all;

% 求出的可作為已知的系統(tǒng)的幅頻響應；
Fs=1;N=256;
A=[1 1  0 0];
f=[0 .6 .7 1];
weigh=[1 10];	
b=remez(42,f,A,weigh);
[h,w]=freqz(b,1,N,1);
h=20*log10(abs(h));
subplot(221);plot(w,h);grid;

r=rand(4096,1);
x=filter(b,1,r);
M=N*2;

% 由函數(shù) tfe 辨識出的系統(tǒng)的幅頻響應；
[H,w]=tfe(r,x,M,Fs,hamming(M),0,'mean');
H=20*log10(abs(H));
subplot(222);
plot(w,H);grid;