MATLAB爱好者论坛-LabFans.com

MATLAB爱好者论坛-LabFans.com (https://www.labfans.com/bbs/index.php)
-   MATLAB论坛 (https://www.labfans.com/bbs/forumdisplay.php?f=6)
-   -   本人初学者,急需知道下面的程序哪里有问题,多谢 (https://www.labfans.com/bbs/showthread.php?t=4011)

sfang9902 2008-07-26 16:31

本人初学者,急需知道下面的程序哪里有问题,多谢
 
这是我设计中的一个关键例子:(望各位高手帮帮忙)
设信号为:x(t)=sin(2π20t),在信号中加入方差为0.1的白噪声序列,取采样频率fs=200Hz,采样点数N=200,用Burg算法分析信号。(应该出现谱分裂现象,可是我弄的程序不出现,调了很多次,也没出来,请各信号处理高手指点,小女子将万分感激)
我编的程序如下:
clear all
fs=200;
nfft=200;
order=20;
t=0:1/fs:1;
f1=20;

x=sin(2*pi*f1*t)+sqrt(0.1)*randn(size(t));
[p,f]=pburg(x,order,nfft,fs);
p1=10*log10(p);

plot(f,p1)
grid on;
xlabel('Frequency(Hz)');
ylabel('相对功率谱密度(dB/Hz)');
title('Burg PSD Estimate');

s060872 2008-09-04 14:39

回复: 本人初学者,急需知道下面的程序哪里有问题,多谢
 
看来你真是初学者,
clear all
fs=200;
nfft=200;
order=20;
t=0:1/fs:1;
f1=20;

x=sin(2*pi*f1*t)+sqrt(0.1)*randn(size(t));
[p,f]=pburg(x,order,nfft,fs);
p1=10*log10(p);
figure,
plot(f,p1)
grid on;
xlabel('Frequency(Hz)');
ylabel('相对功率谱密度(dB/Hz)');
title('Burg PSD Estimate');

s060872 2008-09-04 14:40

回复: 本人初学者,急需知道下面的程序哪里有问题,多谢
 
你必须加个figure语句在plot的上面!!!

sfang9902 2008-09-18 15:38

回复: 本人初学者,急需知道下面的程序哪里有问题,多谢
 
加与不加没有区别,主要是要出现谱分裂。


所有时间均为北京时间。现在的时间是 10:40

Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.