14-02-2013, 04:53 PM
DESIGN AND ANALYSIS OF SPECTRUM ESTIMATORS
DESIGN AND ANALYSIS OF SPECTRUM ESTIMATORS.docx (Size: 29.36 KB / Downloads: 16)
PROGRAM
close all;
clear all;
randn('state',0);
fs=150;
t=(0:fs)/fs;
a=0.9;
b=0.5;
f1=50;
f2=20;
x=a*sin(2*pi*f1*t)+b*sin(2*pi*f2*t)+randn(size(t));
figure(1)
hs=spectrum.periodogram('rectangular');
psd(hs,x,'fs',fs,'nfft',1024);
xlabel('FREQUENCY(HZ)');
ylabel('POWER (dB)');
title('POWER SPECTRUM ESTIMATION OF PERIODOGRAM');
%welch
randn('state',0);
fs=1500;
t=(0:fs)/fs;
a=0.9;
b=0.5;
f1=200;
f2=600;
x=a*sin(2*pi*f1*t)+b*sin(2*pi*f2*t)+randn(size(t));
figure(2)
hs=spectrum.welch('rectangular',150,50);
psd(hs,x,'fs',fs,'nfft',1024);
xlabel('FREQUENCY(HZ)');
ylabel('POWER (dB)');
title('POWER SPECTRUM ESTIMATION OF WELCH');
%BARTLETT
figure(3)
hs=spectrum.welch('bartlett');
psd(hs,x,'fs',fs,'nfft',1024);
xlabel('FREQUENCY(HZ)');
ylabel('POWER (dB)');
title('POWER SPECTRUM ESTIMATION OF BARTLETT');
%BLACKMANN TUKEY
figure(4)
hs=spectrum.welch('blackman');
psd(hs,x,'fs',fs,'nfft',1024);
xlabel('FREQUENCY(HZ)');ylabel('POWER (dB)');
title('POWER SPECTRUM ESTIMATION OF BLACKMANNTUKEY');