Code

Comparison of different power spectrum estimates

Senthilkumar December 28, 2011 Coded in Scilab

This scilab program is used to compare the frequency resolution of different power spectrum estimation methods

//Caption: Determination of Frequency Resolution of  the 
//[1]. Barlett [2]. Welch [3].Blackman Tukey Power Spectrum Estimation
//Methods
clear;
clc;
close;
Q = 10;  //Quality factor
N = 1000; //Length of the sample sequence
//Bartlett Method
F_Bartlett = Q/(1.11*N);
disp(F_Bartlett,'Frequency Resolution of Bartlett Power Spectrum Estimation')
//Welch Method
F_Welch = Q/(1.39*N);
disp(F_Welch,'Frequency Resolution of Welch Power Spectrum Estimation')
//Blackmann-Tukey Method
F_Blackmann_Tukey  = Q/(2.34*N);
disp(F_Blackmann_Tukey,'Frequency Resolution of Blackmann Tukey Power Spectrum Estimation') 
//Result
//Frequency Resolution of Bartlett Power Spectrum Estimation
//     0.0090090
//Frequency Resolution of Welch Power Spectrum Estimation
//     0.0071942  
//Frequency Resolution of Blackmann Tukey Power Spectrum Estimation 
//     0.0042735