[求助]burg进行AR模型谱估计

[sfang9902]

这样写对吗？x1=sin(2*pi*f1*t)+0.1*sqrt(0.5)*randn(size(t))
整个程序：
clear all
fs=40;
nfft=160;
order=30;
t=0:1/fs:1;
f1=7.25;
x1=sin(2*pi*f1*t)+0.1*sqrt(0.5)*randn(size(t));
[p,f]=pburg(x1,order,nfft,fs);
p1=10*log10(p);
subplot(211);
plot(f,p1)
grid on;
xlabel('Frequency(Hz)');
ylabel('相对功率谱密度(dB/Hz)');
title('Burg PSD Estimate');
可是我得到的图形并不理想，而且感觉噪声的产生好像有问题，请信号处理的高手们帮帮忙，指点一下，本人将不胜感激！