01-21-2013
Quote:
Originally Posted by
Corona688
What do you mean by polynomial interpolation?
Polynomial interpolation is an interpolation that based on three value points (two previous points and next point). So, if we see the data, there are an ID 1, 4, and 8. So, we need to find out the value of 2, 3, 5, 6, and 7 based on value 1, 4, and 8. The formula is that
(((x-x2) * (x-x3)) / ((x1-x2) * (x1-x3))) * y1 + (((x-x1) * (x-x3)) / ((x2-x1) * (x2-x3))) * y2 + (((x-x1) * (x-x2)) / ((x3-x1) * (x3-x1))) * y3
x = current ID;
x1 = the first known ID (second previous known ID); ---> 1
x2 = the second known ID (first previous known ID); ---> 4
x3 = the third known ID (next known ID); ---> 8
y1 = the first known value (the value of ID x1)
y2 = the second known value (the value of ID x2)
y3 = the third known value (the value of ID x3)
9 More Discussions You Might Find Interesting
1. UNIX for Dummies Questions & Answers
I've become obsessed with trying to get this to work.
As of yet, I am unable to figure it out. Unfortunately, I don't have Linux or UNIX available when I get home.
Anyone have tips for me on how I can pass param1 to ID via use of COUNTER and loop?
thx.
LIMIT=6
param1="999999999"
export... (0 Replies)
Discussion started by: egkumpe
0 Replies
2. Shell Programming and Scripting
Factorial calculation
Example output:
Please enter a non-negative number: 3
3! = 3 X 2 X 1
= 6
Please enter a non-negative number: 10
10! = 10 X 9 X 8 X 7 X 6 X 5 X 4 X 3 X 2 X 1
= 3628800
Please enter a non-negative number:... (1 Reply)
Discussion started by: cenco
1 Replies
3. Shell Programming and Scripting
Factorial calculation
Example output:
Please enter a non-negative number: 3
3! = 3 X 2 X 1
= 6
Please enter a non-negative number: 10
10! = 10 X 9 X 8 X 7 X 6 X 5 X 4 X 3 X 2 X 1
= 3628800
Please enter a non-negative number: -1
... (1 Reply)
Discussion started by: cenco
1 Replies
4. UNIX for Advanced & Expert Users
Hello everybody!
I have no experience with shell Programmer, but I need to compare 02 files. Txt and generate an output or a new file, after the comparisons. see:
If the column 1 of file1 is equal to column 1 of file2, and column 3 of file2 contains the column 4 of file1, output: column1... (4 Replies)
Discussion started by: He2
4 Replies
5. Shell Programming and Scripting
Hi everybody!
I have written some awk scripts that return me some results I need to process. At the moment I use openOffice to process them, but I am trying to find a more efficient solution using possibly a bash or awk script.
I have two files, file1 is in the format:
time position
... (3 Replies)
Discussion started by: Alice236
3 Replies
6. Shell Programming and Scripting
please can you help me with this script ( very very important )
what I'm trying is to write program that accepts list of user as its argument
1- If a user or more are given as arguments, the script should reset files permissions as follows:
a. Directory ~/share to 750 (if it exists).
b. All... (10 Replies)
Discussion started by: testman84
10 Replies
7. Shell Programming and Scripting
Hi ! Experts...
I just wanted to know whether it is possible in scripting...to do interpolation....
if so....have a look on my data file
I need temperature and salinity value with a bin size of 0.5 m
output looks somewhat like this
dep temp sal
0.5 25 0.077
1 25 ... (12 Replies)
Discussion started by: nex_asp
12 Replies
8. Shell Programming and Scripting
i have one function block in the beginning of my script and there are some commands inside that function which will perform some operations.
And i am invoking that function from my main script by passing some values. Sometimes it is hanging in the middle for some value.
For example:
For 1st... (3 Replies)
Discussion started by: thomasraj87
3 Replies
9. Shell Programming and Scripting
Hi I have data in the below format in two columns in excel which i will copy to notepad.
test as rec1, string
test as rec2, byteint
test as rec3, string
update date as test, datetime
name as tes2 string
I need to add trim function on all the string columns and keep the remaining... (10 Replies)
Discussion started by: pisikar
10 Replies
LEARN ABOUT DEBIAN
pdl::gsl::interp
INTERP(3pm) User Contributed Perl Documentation INTERP(3pm)
NAME
PDL::GSL::INTERP - PDL interface to Interpolation routines in GSL
DESCRIPTION
This is an interface to the interpolation package present in the GNU Scientific Library.
SYNOPSIS
use PDL;
use PDL::GSL::INTERP;
my $x = sequence(10);
my $y = exp($x);
my $spl = PDL::GSL::INTERP->init('cspline',$x,$y);
my $res = $spl->eval(4.35);
$res = $spl->deriv(4.35);
$res = $spl->deriv2(4.35);
$res = $spl->integ(2.1,7.4);
FUNCTIONS
init()
The init method initializes a new instance of INTERP. It needs as input an interpolation type and two piddles holding the x and y values to
be interpolated. The GSL routines require that x be monotonically increasing and a quicksort is performed by default to ensure that. You
can skip the quicksort by passing the option {Sort => 0}.
The available interpolation types are :
linear
polynomial
cspline (natural cubic spline)
cspline_periodic (periodic cubic spline)
akima (natural akima spline)
akima_periodic (periodic akima spline)
Please check the GSL documentation for more information.
Usage:
$blessed_ref = PDL::GSL::INTERP->init($interp_method,$x,$y,$opt);
Example:
$x = sequence(10);
$y = exp($x);
$spl = PDL::GSL::INTERP->init('cspline',$x,$y)
$spl = PDL::GSL::INTERP->init('cspline',$x,$y,{Sort => 1}) #same as above
# no sorting done on x, user is certain that x is monotonically increasing
$spl = PDL::GSL::INTERP->init('cspline',$x,$y,{Sort => 0});
eval()
The function eval returns the interpolating function at a given point. By default it will barf if you try to extrapolate, to comply
silently if the point to be evaluated is out of range pass the option {Extrapolate => 1}
Usage:
$result = $spl->eval($points,$opt);
Example:
my $res = $spl->eval($x)
$res = $spl->eval($x,{Extrapolate => 0}) #same as above
# silently comply if $x is out of range
$res = $spl->eval($x,{Extrapolate => 1})
deriv()
The deriv function returns the derivative of the interpolating function at a given point. By default it will barf if you try to
extrapolate, to comply silently if the point to be evaluated is out of range pass the option {Extrapolate => 1}
Usage:
$result = $spl->deriv($points,$opt);
Example:
my $res = $spl->deriv($x)
$res = $spl->deriv($x,{Extrapolate => 0}) #same as above
# silently comply if $x is out of range
$res = $spl->deriv($x,{Extrapolate => 1})
deriv2()
The deriv2 function returns the second derivative of the interpolating function at a given point. By default it will barf if you try to
extrapolate, to comply silently if the point to be evaluated is out of range pass the option {Extrapolate => 1}
Usage:
$result = $spl->deriv2($points,$opt);
Example:
my $res = $spl->deriv2($x)
$res = $spl->deriv2($x,{Extrapolate => 0}) #same as above
# silently comply if $x is out of range
$res = $spl->deriv2($x,{Extrapolate => 1})
integ()
The integ function returns the integral of the interpolating function between two points. By default it will barf if you try to
extrapolate, to comply silently if one of the integration limits is out of range pass the option {Extrapolate => 1}
Usage:
$result = $spl->integ($a,$b,$opt);
Example:
my $res = $spl->integ($a,$b)
$res = $spl->integ($a,$b,{Extrapolate => 0}) #same as above
# silently comply if $a or $b are out of range
$res = $spl->eval($a,$b,{Extrapolate => 1})
BUGS
Feedback is welcome.
SEE ALSO
PDL
The GSL documentation is online at
http://sources.redhat.com/gsl/ref/gsl-ref_toc.html
AUTHOR
This file copyright (C) 2003 Andres Jordan <andresj@physics.rutgers.edu> All rights reserved. There is no warranty. You are allowed to
redistribute this software/documentation under certain conditions. For details, see the file COPYING in the PDL distribution. If this file
is separated from the PDL distribution, the copyright notice should be included in the file.
The GSL interpolation module was written by Gerard Jungman.
perl v5.14.2 2012-05-30 INTERP(3pm)