[Date Index]
[Thread Index]
[Author Index]
Re: Fitting step function
*To*: mathgroup at smc.vnet.net
*Subject*: [mg122100] Re: Fitting step function
*From*: Jacopo Bertolotti <jacopo.bertolotti at gmail.com>
*Date*: Fri, 14 Oct 2011 05:52:42 -0400 (EDT)
*Delivered-to*: l-mathgroup@mail-archive0.wolfram.com
*References*: <201110120744.DAA16264@smc.vnet.net>
Maybe I totally misunderstood your point but, if you really want to know
where the phase jumps (likely phase wraps) happen, wouldn't be a lot
easier to just look when the difference between two adjacent points is
bigger than Pi or something like that?
Here is an example code (sorry, it is not easily readable but I hope it
makes my point clear):
phasedata = Table[Arg[Exp[I x^2]], {x, 0, 10, 0.01}]; (*any phase data.
Here I am just faking some.*)
jumps = Drop[ Flatten@Position[Round[Differences[phasedata]/(2 Pi)],
Except[0]], 1]; (*find those points whose difference with the previous
one is approximately a multiple of 2 Pi. There is a spurious 0 that is
eliminated with the Drop command. *)
jumps = Table[{jumps[[i]], 0}, {i, 1, Dimensions[jumps][[1]]}]; (* give
the coordinates a y value just for the sake of plotting *)
ListPlot[phasedata, Joined -> True, Epilog -> {PointSize[Medium], Red,
Map[Point, jumps]}] (*It works!!!*)
Regards
Jacopo
On 10/12/2011 09:44 AM, premiumxy wrote:
> Hello,
>
> I have the following problem: I have a set of data about some phase information of a signal. At certain points there is a clear step of 180 degrees.
>
> I would like to fit this data to find the x-coordinates of these steps. I have the feeling that "normal" step-functions in Mathematica (such as UnitStep, HeavisideTheta, Sign...) are not varied in the fit-routines (such as NonlinearModelFit). I guess it this is caused by the functions not being continuous.
>
> I would appreciate any hints for this quite simple (?)problem. The only option that I can currently see is to create a "step-like" continuous function such as Tanh[A x] (with A>>1)
>
> Thanks!
>
>
Prev by Date:
**Re: Plot function with two arguments**
Next by Date:
**Re: Fitting step function**
Previous by thread:
**Re: Fitting step function**
Next by thread:
**Re: Fitting step function**
| |