MATLAB函数工具箱提供了tfe函数用来实现基于经典谱估计的系统辨识。调用方式：（1）Txy=tfe(x,y,NFFT,Fs,window)使用welch平均周期图法，根据输入变量x和输出变量y来估计系统的传递函数。参数NFFT用来指定FFT运算所采用的点数。如果x和y都是实信号、NFFT为偶数，则Txy的长度为NFFT/2+1；如果x和y都是实信号、NFFT为奇数，则Txy的长度为(NFFT+1)/2；如果x或y是复信号，则Txy的长度为NFFT；参数window用来指定所采用的窗函数。窗函数的长度必须与向量x的长度一样大。（2）[Txy,f]=tfe(x,y,NFFT,Fs,window,noverlap)返回和传递函数的估计同样大