sinesum

SINESUM Estimates a signal as a sum of pure sine waves of known frequency

function [output m p]= sinesum(signal, freq, Fs)

SINESUM Estimates a signal as a sum of pure sine waves of known frequency
'signal' - Input signal you wish to estimate
'freq' - Vector of known frequencies with significant contributions
If no vector is specified, vector is assumed to be a
fundamental (primary freq of signal) plus harmonics (number of
harmonics specified by the constant 'harmonics'
'Fs' - Sampling rate of input signal
'output' - Estimate of signal using pure sine functions
'm' - Vector of frequency component amplitudes
'p' - Vector of frequency component phases (reconstructed
signal must be reversed in time as of 4/24/2012)

This effectively eliminates "noise" at frequencies other than those specified in the vector 'freq'. The resulting signal, 'output' is a sum of pure sinusoids. This function is useful for extracting information about the strength of specific frequency components as well as the amplitude of a noisy signal.

Created on: April 24, 2012
Author: Srikanth Iyer
Email: [e-mail]