Code covered by the BSD License  

Highlights from
Frequency, amplitude, phase and mean value of sine wave

4.5

4.5 | 2 ratings Rate this file 33 Downloads (last 30 days) File Size: 2.46 KB File ID: #19902

Frequency, amplitude, phase and mean value of sine wave

by

 

13 May 2008 (Updated )

The function sinfapm evaluates parameters of sampled sine wave

| Watch this File

File Information
Description

The function sinfap.m evaluates frequency, amplitude, phase and mean value of a uniformly sampled harmonic signal
   x(t) = a.sin(2.pi.f.t + phi) + x_m
It uses a vector version of 3-point formulae derived by application of
Z-transform (see [1]) for finding amplitude and frequency of a signal.
If more than two output parameters are to be determined, all of them are optimized in the least squares sense by the function LMFnlsq.

Calls:
frq = sinfapm(x,fs); % Get only frequaency of sine-wave
[frq,amp] = sinfapm(x,fs); % Get frequency and amplitude
[frq,amp,phi] = sinfapm(x,fs); % Get frequency, amplitude and phase
[frq,amp,phi,ave] = sinfapm(x,fs); % ditto plus mean value
The set of more than two output parameters can be found by calling
[frq,amp,phi] = sinfapm(x,fs,Name_1,Val_1,Name_2,Val_2, ...);
[frq,amp,phi,ave] = sinfapm(x,fs,Name_1,Val_1,Name_2,Val_2, ...);

Input arguments:
x % vector of samples
fs % sampling frequency [Hz]
Name_i % name of the i-th optional parameter for optimization
Val_i % value of the i-th optional parameter (see function LMFnlsq)
Output arguments:
frq % frequency of x [Hz]
amp % amplitude of x
phi % phase in radians
Examples:
[f,a,phi,ave] = sinfapm([1.3633;-.2428;-0.9705;1.8130;-1.9631],10);
% f = 4.0000
% a = 2.0000
% phi = 0.7500
% ave = -2.2806e-005
[f,a,phi] = sinfapm([.707,1,.707,0],20,'Xtol',1e-4);
% f = 2.5001
% a = 0.9999
% phi = 0.7853 % pi/4 = 0.785398...

Acknowledgements

This file inspired Sinide Parameter Identification Of A Sine Wave From A Measured Signal.

MATLAB release MATLAB 7.6 (R2008a)
Other requirements Function LMFnlsq from FEX, Id. 17534
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (8)
18 Aug 2014 Miroslav Balda

The function sinfapm is not maintained any more. Use the function sinide instead, which is much more stable. See

sinide: www.mathworks.com/matlabcentral/fileexchange/45567

13 Mar 2014 will wehner

Thank you.
Unfortunately, I believe there are still problems with the function.

I emailed you an example (miroslav@balda.cz) with estimation of real data. It uses your function and compares to estimation with nlinfit.m

Will

14 Feb 2014 Miroslav Balda

@Will
I am sorry that the function behaved wrong. A reason for it was a hight sensitivity of the function to inaccurate initial guess of the sought parameters. Since it was impossible to improve it by little modification, I decided to build a new function. It will appear in a short time under the name sinide.m.

30 Jan 2014 will wehner

t=0:.01:1;
y=2*sin(2*pi*10*t+3)+5;

[f,a,phi,ave] = sinfapm(y,100)

Returns
f=2.486
a=-0.1283
phi=1.5719
ave=5.0035

11 Oct 2009 John

When I implement this function for noisy sine vawes, the amplitudes seem to be very far from reality. Often, amplitudes have even minus signs(?). Any ideas why ?

18 Apr 2009 Pablo

txs a lot =)).
i need it that mean

24 Jul 2008 david sweeney

How do I get this file to run????

31 May 2008 Ramesh Pachar  

Contact us