生成归一化频率的频谱图

根据提供的引用内容，我们可以使用MATLAB中的pwelch函数来生成归一化频率的频谱图。pwelch函数可以计算信号的功率谱密度（PSD），并且可以自动进行零填充和窗口处理。

下面是一个示例代码，其中a是一个包含信号的矩阵，Nf是FFT的点数，p是AR模型的阶数：

% 计算p阶AR模型的功率谱
Par = sigma(p) ./ fftshift(abs(fft([1, a(p, :)], Nf)).^2);

% 计算m阶AR模型的输入方差
for m = 2:p
    sigma(m) = sigma(m-1) * (1 - abs(k(m))^2);
end

% 生成归一化频率的频谱图
fs = 1000; % 采样率
win = hann(Nf)'; % 汉宁窗
[Pxx, f] = pwelch(a, win, [], Nf, fs); % 计算PSD
plot(f/fs, 10*log10(Pxx)); % 绘制频谱图
xlabel('Normalized Frequency');
ylabel('Power/Frequency (dB/Hz)');

其中，pwelch函数的第一个参数是信号矩阵，第二个参数是窗口函数，第三个参数是重叠长度，第四个参数是FFT的点数，第五个参数是采样率。