exa100802_tfe.m

来自「功率谱估计的一些matlab例程」· M 代码 · 共 28 行

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

% 求出的可作为已知的系统的幅频响应；
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;

% 由函数 tfe 辨识出的系统的幅频响应；
[H,w]=tfe(r,x,M,Fs,hamming(M),0,'mean');
H=20*log10(abs(H));
subplot(222);
plot(w,H);grid;