![]() |
RTXI 1.3
|
00001 // 00002 // file = sam_spec.cpp 00003 // 00004 00005 #include <stdlib.h> 00006 #include <iostream> 00007 #include "complex.h" 00008 #include "sam_spec.h" 00009 #include "fft.h" 00010 00011 SampleSpectrum::SampleSpectrum( complex* time_signal, 00012 double samp_intvl, 00013 int num_samps ) 00014 :PsdEstimate( num_samps, 00015 samp_intvl) 00016 { 00017 complex *sample_spectrum = new complex[num_samps]; 00018 00019 fft(time_signal,sample_spectrum,num_samps); 00020 for(int samp_idx=0; samp_idx<num_samps; samp_idx++) 00021 { 00022 Psd_Est[samp_idx] = mag_sqrd(sample_spectrum[samp_idx]); 00023 } 00024 return; 00025 } 00026