Burg算法功率谱估计,采样点数nfft不会影响MAX值和位置.

s=data22(6,:)

s=s(50:150)

nfft=2048

fs=500 EEG

[px1,f1]=pburg(s,6,nfft,fs) figure

plot(f1,px1) title( CI——Burg算法功率谱估计 )

xlabel( Frequency(Hz) ) ylabel( Power )



IndMax=find(diff(sign(diff(px1)))<0)+1

MAX=px1(IndMax)

loc=f1(IndMax)



figure plot(f1,px1) title( CI——Burg算法功率谱估计 )

xlabel( Frequency(Hz) ) ylabel( Power )

hold on plot(loc(1),MAX(1), * )

hold on plot(loc(2),MAX(2), * )

hold on plot(loc(3),MAX(3), * )

-　Burg algorithm for power spectrum estimation, sampling points nfft MAX value and does not affect the position of s = data22 (6, :) . S = s (50:150)　nfft = 2048　fs = 500　EEG [px1, f1] = pburg (s, 6, nfft, fs)　figure　plot (f1, px1)　title (' CI- Burg algorithm for power spectral estimation' )　xlabel (' Frequency (Hz)' )　ylabel (' Power'　)　IndMax = find (diff (sign (diff (px1))) < 0)+1　MAX = px1 (IndMax)　loc = f1 (IndMax)　figure　plot (f1, px1)　title (' CI- Burg algorithm for power spectral estimation ' )　xlabel ('　Frequency (Hz) ' )　ylabel ('　Power ' )　hold on　plot (loc (1), MAX (1),' * ' )　hold on　plot ( loc (2), MAX (2), ' *' )　hold on　plot (loc (3), MAX (3), ' *' )