Path: news.mathworks.com!newsfeed-00.mathworks.com!solaris.cc.vt.edu!news.vt.edu!news.glorb.com!postnews.google.com!r33g2000yqn.googlegroups.com!not-for-mail
From: Adam Chapman <adamchapman1985@hotmail.co.uk>
Newsgroups: comp.soft-sys.matlab
Subject: Re: An Interesting Real-World Problem...
Date: Thu, 2 Jul 2009 11:55:40 -0700 (PDT)
Organization: http://groups.google.com
Lines: 35
Message-ID: <13c05c8e-4af1-4053-b21a-a45a0137acf8@r33g2000yqn.googlegroups.com>
References: <h2iuip$nm2$1@fred.mathworks.com>
NNTP-Posting-Host: 87.113.152.34
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable
X-Trace: posting.google.com 1246560940 27586 127.0.0.1 (2 Jul 2009 18:55:40 GMT)
X-Complaints-To: groups-abuse@google.com
NNTP-Posting-Date: Thu, 2 Jul 2009 18:55:40 +0000 (UTC)
Complaints-To: groups-abuse@google.com
Injection-Info: r33g2000yqn.googlegroups.com; posting-host=87.113.152.34; 
	posting-account=QATl1QoAAADZYVAWORqO_aubcdzudO6s
User-Agent: G2/1.0
X-HTTP-UserAgent: Mozilla/4.0 (compatible; MSIE 7.0; Windows NT 6.0; 
	Trident/4.0; SLCC1; .NET CLR 2.0.50727; Media Center PC 5.0; .NET CLR 
	1.1.4322; InfoPath.2; .NET CLR 3.5.30729; .NET CLR 3.0.30618; 
	OfficeLiveConnector.1.3; OfficeLivePatch.0.0),gzip(gfe),gzip(gfe)
Xref: news.mathworks.com comp.soft-sys.matlab:552495


On Jul 2, 7:34 pm, "Geoffrey Akien" <geoff.ak...@gmail.com> wrote:
> I've been using MATLAB for around 6 months now, and only now have I needed to do something that is actually mathematically-related!  I'm a chemist, so high-level maths is not my strong point.
>
> I have some real-world data that I'd like to pick out the dramatic change in gradient of:
>
> http://www.flickr.com/photos/40058045@N07/3681613485/
> (wish I could include this in the post but never mind)
>
> Because the data is quite noisy (and not random) this is quite awkward.  In that picture there are about 5000 points.  You cannot smooth too much (in practice >300) without losing the information about the point at which the gradient changes (the smoothed data starts to cut the corner).  Even when the data is smoothed, the presence of non-random noise means that the "flat" part of the curve is not actually flat, and there are several stationary points (calculated thanks to STATIONARY from FEX) along that flat part.
>
> I've also tried fitting a least-squares linear regression using POLYFIT for each of the data points with ~500 points either side of it.  This was with the aim of picking out the part of the graph with the lowest R2 value - this I assumed would represent the fast gradient change.  However, this still does not pick out the right point, selecting a point about 0.1 hours before it should be, and I can't think why this is the case.
>
> Any other good suggestions?
>
> Thanks
>
> Geoff

Is the gradient of the sharp drop always the same? Im thinking a
macauley step function might be useful