Yule-Walker AR法的功率谱估计：pyulear函数

（1）Pxx=pyulear（x，order）：用Yule-Walker AR法对离散时间信号x进行功率谱估计。输入参数order为 AR模型的阶数。如果x为实信号，则返回结果为“单边”功率谱；如果x为复信号，则返回结果为“双边”功率谱。

（2）Pxx=pyulear（x，order，NFFT）：参数NFFT用来指定FFT运算所采用的点数

如果x为实信号、NFFT为偶数，则Pxx的长度为（NFFT/2+1）；

如果x为实信号、NFFT为奇数，则Pxx的长度为（NFFT+l）/2；

如果x为复信号，则Pxx的长度为NFFT；

NFFT的默认值为256。

（3）［Pxx，w］=pyulear（…）：输出参数w为和估计PSD的位置一一对应的归一化角频率，单位为rad/sample，其范围如下：

如果x为实信号，则w的范围为［0，pi］；

如果x为复信号，则w的范围为［0，2*pi］。

（4）［Pxx，f］=pyulear（…，Fs）：同时返回和估计PSD的位置一一对应的线性频率f，单位为Hz，参数Fs为采样频率，单位也是Hz。当Fs为空矩阵“［］”时，则使用默认值1 Hz。输出参数f的范围如下：

如果x为实信号，则f的范围为［0，Fs/2］；

如果x为复信号，则f的范围为［0，Fs］。

（5）［Pxx，w］=pyulear（…，‘twosided’）：在［0，2*pi］区间上进行功率谱的“双边”估计；twosided也可以由onesided代替。

（6）［Pxx，f］=pyulear（…，Fs，‘twosided’）：在［0，Fs］区间上进行功率谱的“双边”估计；twosided也可以由onesided代替。

（7）pyulear（…）：没有输出参数，在当前图形窗口里绘制出PSD估计结果图，坐标分别为dB和归一化频率。

［例4-4］用Yule-Walker AR 法进行PSD 估计，结果如图4-14 所示。

Fs=500；t=0：1/Fs：1；

x=sin（2*pi*60*t）+4*sin（2*pi*110*t）+sin（2*pi*210*t）+randn（size（t））；

pyulear（x，20，［］，Fs）。

图4-14 Yule-Walker AR法功率谱估计结果图