刚才问题的补充
这样写对吗?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');
可是我得到的图形并不理想,而且感觉噪声的产生好像有问题,请信号处理的高手们帮帮忙,指点一下,本人将不胜感激!
|