Model and experiment fitting.
Sebastian Żurek
sebzur at pin.if.uz.zgora.pl
Fri Oct 20 17:33:47 EDT 2006
Hi!
This is probably a silly question but I'm getting confused with a
certain problem: a comparison between experimental data points (2D
points set) and a model (2D points set - no analytical form).
The physical model produces (by a sophisticated simulations done by an
external program) some 2D points data and one of my task is to compare
those calculated data with an experimental one.
The experimental and modeled data have form of 2D curves, build of n
2D-points, i.e.:
expDat=[[x1,x2,x3,..xn],[y1,y2,y3,...,yn]]
simDat=[[X1,X2,X3,...,Xn],[Y1,Y2,Y3,...,Yn]]
The task of determining, let's say, a root mean squarred error (RMSe)
is trivial if x1==X1, x2==X2, etc.
In general, which is a common situation xk differs from Xk (k=0..n) and
one may not simply compare succeeding Yk and yk (k=0..n) to determine
the goodness-of-fit. The distance h=Xk-X(k-1) is constant, but similar
distance m(k)=xk-x(k-1) depends on k-th point and is not a constant
value, although the data array lengths for simulation and experiment are
the same.
My first idea was to do some interpolations to obtain the missing
points, but I held it 'by a hand' (which, BTW gave quite rewarding
results) and I suppose, there's some i.g. numpy method to do it for me,
isn't it?
I suppose to do something like:
gfit(expDat,simDat,'measure_type')
which I hope will return the number determining the goodness-of-fit
(mean squarred error, root mean squarred error,...) of two sets of
discrete 2D data points.
Is there something like that in any numerical python modules (numpy,
pylab) I could use?
I can imagine, I can fit the data with some polynomial or whatever,
and than compare the fitted data, but my goal is to operate on
as raw data as it's possible.
Thanks for your comments!
Sebastian
-------------------------------------------------------------------------
Using Tomcat but need to do more? Need to support web services, security?
Get stuff done quickly with pre-integrated technology to make your job easier
Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo
http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
More information about the NumPy-Discussion
mailing list