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

[sfang9902]

请教一下产生噪声的小程序，
如果有信号x(t)=sin(2*π*7.25t+50),即信号幅值为1，初相角50，并在信号中加入最大幅值为0.1的随机噪声（信噪比SNR=20），选取对模型AR阶数order=30用burg算法对信号进行分析。
怎样产生噪声呢？谢谢。（如果有高手对这个小程序很了解的话，请指教），多谢！

[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');
可是我得到的图形并不理想，而且感觉噪声的产生好像有问题，请信号处理的高手们帮帮忙，指点一下，本人将不胜感激！

[zzxinlovewl]

你对信噪比snr的定义理解不对，它不是直接幅度的比值，而是一个dB的式子。

[脚踏实地]

考虑调整阶数

[侯明月]

我有点不明白：为什么要在白噪声前还要另外再乘以 根号0.5？？
信噪比已经满足20db了的~~
我刚才试了一下，有根号0.5的图，很平缓，没有的，峰和谷都很明显。。。你试试吧~~