作为基本图形的生成,圆有很多种方法实现,这是一个很简单但有典型的算法,希望能给人一点启发．

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) 返回和传递函数的估计同样大小、一一对应的线性频率f。参数noverlap用来指定数据分段间重叠的样本数。 （3）tfe(x,y,...dflag) 参数dflag用来指定对x和y的预处理方式，其取值有 dflag=linear：去掉加窗数据中的最佳直线拟合； dflag=mean：去掉加窗数据中的均值； dflag=none：不做处理。 （4）tfe(…) 没有输出参数，在当前图形窗口里绘制出传递函数的估计结果图。 以上这些参数的默认值分别为： NFFT=256； Fs=2； noverlap=0； window=hanning（NFFT）； dflag=none。