tmp613.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 37 行

%MATLAB PROGRAM 6-13
%Transfer Function Estimate
%Create input signal
N=1024;
Nfft=256;
window=hanning(256);
noverlap=128;
dflag='none';
Fs=1000;
n=0:N-1;
t=n/Fs;
randn('state',0)
xn=sin(2*pi*50*t)+randn(1,N);
%Create a system: filter
h=ones(1,10)/10;
%Compute output signal
yn=filter(h,1,xn);
%Estimate Transfer Function by 'etf'
[HEST,f]=tfe(xn,yn,Nfft,Fs,window,noverlap,dflag);
%Compute frequency response by 'freqz'
H=freqz(h,1,f,Fs);
%Compare the results
subplot(211)
plot(f,abs(H));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Frequency Response')
axis([0 500 0 1]);
grid

subplot(212)
plot(f,abs(HEST));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Transfer Function Estimate');
axis([0 500 0 1])
grid