Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
![]() |
|
![]() |
#1 |
初级会员
注册日期: 2008-06-16
年龄: 44
帖子: 11
声望力: 17 ![]() |
![]()
这是我设计中的一个关键例子:(望各位高手帮帮忙)
设信号为: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'); |
![]() |
![]() |
![]() |
#2 |
初级会员
注册日期: 2008-09-04
年龄: 41
帖子: 13
声望力: 17 ![]() |
![]()
看来你真是初学者,
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'); |
![]() |
![]() |
![]() |
#3 |
初级会员
注册日期: 2008-09-04
年龄: 41
帖子: 13
声望力: 17 ![]() |
![]()
你必须加个figure语句在plot的上面!!!
|
![]() |
![]() |
![]() |
#4 |
初级会员
注册日期: 2008-06-16
年龄: 44
帖子: 11
声望力: 17 ![]() |
![]()
加与不加没有区别,主要是要出现谱分裂。
|
![]() |
![]() |