Re: Interpolation
- To: mathgroup at smc.vnet.net
- Subject: [mg32309] Re: Interpolation
- From: "Allan Hayes" <hay at haystack.demon.co.uk>
- Date: Sat, 12 Jan 2002 05:18:28 -0500 (EST)
- References: <a1mbl8$2ot$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Yas, Here is an idea that might be refined: dat=Table[{t,Sin[0.75+Exp[Pi*t]]},{t,0.,1.,.001}]; ListPlot[dat] Cases[Partition[dat,3, 1],{a_,b_,c_}/;a[[2]]\[LessEqual]b[[2]]&&b[[2]]\[GreaterEqual]c[[2]]:>b] Of course the value {0., 0.983986} is not found. -- Allan --------------------- Allan Hayes Mathematica Training and Consulting Leicester UK www.haystack.demon.co.uk hay at haystack.demon.co.uk Voice: +44 (0)116 271 4198 Fax: +44 (0)870 164 0565 "Yas" <y.tesiram at pgrad.unimelb.edu.au> wrote in message news:a1mbl8$2ot$1 at smc.vnet.net... > G'day mathgroup, > I have some general questions on Interpolation as implemented in > Mathematica. Essentially, I have a series of profiles for which I have > no general function. The profile is generated by repetitive > multiplication of rotation matrices and so it is difficult (and also > impractical) to deduce an analytic function. Instead to answer the > questions I have, I have decided to use Interpolation -- a pitiful > state, nevertheless a shade better than extrapolation. The first step > involves creating the Interpolating function, > > thing1 = Interpolation[data]; > > Then I find the first derivative, > > thing2 = D[thing1[x], {x, 0, lastdatapoint}] > > Next I want to find the points where thing2 = 0, but I have run into > problems with Mathematica complaining about inverse functions etc, > although plots of thing2 versus x look fine. > > My primary question is, > > 1. The profiles that I have, have maxima and minima whose co-ordinates I > want to find, hence the differentiation step. Interpolation does a good > job reproducing the data points and in finding a derivative that can be > plotted but not when asked to Solve for thing2 = 0. How do I go about > finding the accuracy of the Interpolating Function to test whether the > value near thing2 = 0 is well behaved? > > And a secondary question is, > > 2. Is there another efficient method of estimating the co-ordinates of > the minima and maxima by computer of data sets? > > The length of each of these data sets is 1000 points and for that reason > I have not pasted it into the email. As a demonstrative example, the > Table of values generated by, > > Table[Sin[0.75 + Exp[Pi*t]], {t, 0, 1, 1/1000}]; > ListPlot[%] > > is closely related. > > Any help or comments will be appreciated. > > Thanks > Yas > >