<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Analysis of the degree of polynomials in the interpolation problem</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Alicja Winnicka Institute of Mathematics Silesian University of Technology Kaszubska 23</institution>
          ,
          <addr-line>44-100 Gliwice</addr-line>
          ,
          <country country="PL">Poland</country>
        </aff>
      </contrib-group>
      <fpage>27</fpage>
      <lpage>30</lpage>
      <abstract>
        <p>-Interpolation is a tool based on finding a specific function that passes through certain points. This issue is important from the point of view of applications in data processing, and in particular finding a curve showing some initial conditions. In this work, interpolation with multitude of different degrees is analyzed due to their advantages and disadvantages.</p>
      </abstract>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>I. INTRODUCTION</title>
      <p>Mathematics gave the backbone to many fields of science,
in particular to computer science. A large number of more
and more complex algorithms or applications is based on the
determination of specific elements. To this end, a numerical
approach is used that the computer is able to perform.
However, this approach means approximation, since it is not always
possible to accurately determine any values for continuous
functions. The reason is an infinite set of values.</p>
      <p>The approximation of results is a very important aspect,
above all in biometry. In [1], the authors used curve
approximation to modify the set of points representing the signature.
Again in [2], the idea of using interpolation was used for
the purpose of verification of security protocols. Another
area where approximation is used is feature extraction [3]
and classification [4]. In addition, this type of mathematical
operations are used in the process of creating rules for fuzzy
systems [5]. In this paper, we analyze different degrees of
polynomials in the interpolation process.</p>
    </sec>
    <sec id="sec-2">
      <title>II. POLYNOMIAL INTERPOLATION</title>
      <p>Polynomial interpolation is a method in numerical analysis,
which purpose is to find the polynomial approximated to the
original function. The interpolation works, when we have a
set of points x0, x1,..., xn, which are called nodes and their
values y0, y1,...,yn. The searched polynomial must satisfied
a few requirements: its degree must be smaller or equal to
(n 1) and points (xi; yi) must belong to it. It means, that
w(x0) = y0, w(x1) = y1,..., w(xn) = yn and it will be an
approximation of searched function.</p>
      <p>Interpolation is usually used in the situation, when we have
a certain number of measurements, for example in physics
experiment, and we have to find an approximated function for
them.</p>
      <p>c 2019 for this paper by its authors. Use permitted under Creative
Commons License Attribution 4.0 International (CC BY 4.0).</p>
      <p>
        Let us assume, that polynomial w(x) has a following
formula
w(x) = anxn + an 1xn 1 + ::: + a2x2 + a1x + a0
(
        <xref ref-type="bibr" rid="ref1">1</xref>
        )
Rememebring that w(xi) = yi we can write the system of
equations
&gt;8anx0n + an 1x0n 1 + ::: + a1x0 + a0 = y0
&gt;
&gt;&gt;&gt;&gt;anx1n + an 1x1n 1 + ::: + a1x1 + a0 = y1
&lt;
&gt;&gt;&gt;anxnn 1 + an 1xnn 11 + ::: + a1xn 1 + a0 = yn 1
&gt;
&gt;
&gt;:anxnn + an 1xnn 1 + ::: + a1xn + a0 = yn
which can be calculated as the system of matrices
6264 xx::10nn: xx10nn::: 11 xx10nn::: 22 ::::::::: xx::10: :11::3577 2646a an::n: 13775 = 2664 yy::01: 3775 (
        <xref ref-type="bibr" rid="ref3">3</xref>
        )
xnn xnn 1 xnn 2 ::: xn 1 a0 yn
We have to find coefficients a0; a1,...,an, but the condition
number of this matrix may be large. Additionaly number of
equations may cause difficulty in calculating it fast enough.
In due to it, the simpler method was created – writing the
polynomial as Lagrange’s polynomial.
      </p>
      <p>A. Lagrange’s Polynomial</p>
      <p>
        The Lagrange’s polynomial is a easier way to find searched
polynomial of the function. Instead of calculating a system
of equation or matrices, we must find the polynomial using
specific formula
(
        <xref ref-type="bibr" rid="ref2">2</xref>
        )
Values yi can be replaced with values of function our f (xi)
and then polynomial will have following formula
w(x) =
n
X yi
i=0
n
Y
      </p>
      <p>x
j=0^j6=i xi
xj
xj
w(x) =
n
X f (xi)
i=0
n
Y</p>
      <p>
        x
j=0^j6=i xi
xj
xj
(
        <xref ref-type="bibr" rid="ref4">4</xref>
        )
(
        <xref ref-type="bibr" rid="ref5">5</xref>
        )
      </p>
      <p>All of approximation methods are determined by an error
dependent on the features of the method. For the polynomial
interpolation method, when we have function f : [a; b] ! R
and its approximated polynomial calculated on set od nodes
x0, x1,...,xn, the error have the following formula
f (x) w(x) = Mn (x x0)(x</p>
      <p>n!
where Mn is determined by the formula:
x1):::(x
xn)
and f (n)(x) means n-th derivative of f (x).</p>
      <p>Mn = supjf (n)(x)j</p>
    </sec>
    <sec id="sec-3">
      <title>III. EXPERIMENTS</title>
      <p>To find the difference for various n and xi, we compared
polynomials for three functions:
f (x) = sin(x) 2 [0; ],
g(x) = x5 + 2x2 3 2 [ 1; 1],
h(x) = 1+0:108x2 2 [ 20; 20].</p>
      <p>
        For each of them we used another interval and nodes to
show additional features of Lagrange’s interpolation. For the
trigonometric function f (x) = sin(x) we found the second,
third, fourth and fifth degree polynomial. The effect was
showed on Fig.1. We can see on the graph, that the second
degree polynomial, linear function, is not visible. This is due
to our interval and nodes (x0 = 0 and x1 = ) and their
values, which are 0. Because it is linear function, it found the
polynomial w(x) = 0. For a quite accurate approximation of
(
        <xref ref-type="bibr" rid="ref6">6</xref>
        )
(
        <xref ref-type="bibr" rid="ref7">7</xref>
        )
sin(x) five nodes was enough. The found polynomial has a
following formula
and covers the original function on the graph.
      </p>
      <p>Second analyzed function is g(x), which is third degree
polynomial. In this way we can say, that for n = 4 it will
find precisely the same polynomial as g(x). However it finds
a little different polynomial showed on Eq.9.
w(x) = 3 6:710 16z +2x2 +3:610 15x3 +8:810 16x4 +x5
(9)</p>
      <p>We can see, that coefficients for x3 and x4 are very close to
be 0. In this case if indeed they would equal zero, we would
get the function g(x). Similar to previous function f (x), the
third degree polynomial w(x) covers original function on the
graph.</p>
      <p>The above results suggest us, that for the bigger number of
the nodes we will get polynomial, which is more approximated
to original function. However with bigger n, the calculating
last longer and becomes less effective.</p>
      <p>The third function h(x) = 1+0:108x2 is moe complicated,
because it shows us the problem called Runge’s phenomenon.
It is happening for the functions, where our nodes are
equidistant. Then with increasing the degree of polynomial, quality
of interpolation increases, but after a while rapidly deacreases.
It is the most visible at the end of the interval, in this case
[ 20; 20]. The Fig. 2 shows the effect of eqidistant nodes.
We can see, that each one graph for increased degree is less
effective atthe ends of interval. It is the most visible for tenth
degree (red color), where difference between polynomial and
original function h(x) is huge.</p>
      <p>To prevent increasing of this error, we should find another
nodes, which will be concetrated at the ends of interval. For
this function we used two set of nodes to find tenth degree
polynomial:
xi =
20; 16; 12; 8; 4; 0; 4; 8; 12; 16; 20
(10)
and
x1 = 20; 19; 17; 15; 10; 5; 0; 5; 10; 15; 17; 19; 20
(11)
where the first one shows us the Runge’s phenomenon and
the other is to reduce difference between function and its
polynomial. Comparision of both polynomials is shown on
Fig.3.</p>
      <p>Increasing density of nodes at the ends of the interval causes
aproaching the polynomial to the original function, but it also
caused less effectivness for values about [ 10; 5] and [5; 10].
To reduce these errors we should find more nodes in this
interval and calculate a polynomial with the bigger degree.
However, the calculations will last longer because of bigger set
of data, which will decrease effectiveness of the interpolation.</p>
    </sec>
    <sec id="sec-4">
      <title>IV. CONCLUSION</title>
      <p>In this paper we compared the degree of polynomials in
polynomial interpolation. This is accurate method for function,
which are similar to polynomials and usually the better result
is given for the bigger number of nodes. Howevere with
increasing number of nodes and degree of polynomial, the
time for calculations increases too.</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>D.</given-names>
            <surname>Połap</surname>
          </string-name>
          and
          <string-name>
            <surname>M.</surname>
          </string-name>
          <article-title>Woz´niak, “Flexible neural network architecture for handwritten signatures recognition</article-title>
          ,”
          <source>International Journal of Electronics and Telecommunications</source>
          , vol.
          <volume>62</volume>
          , no.
          <issue>2</issue>
          , pp.
          <fpage>197</fpage>
          -
          <lpage>202</lpage>
          ,
          <year>2016</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>M.</given-names>
            <surname>Rocchetto</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Viganò</surname>
          </string-name>
          , and
          <string-name>
            <given-names>M.</given-names>
            <surname>Volpe</surname>
          </string-name>
          , “
          <article-title>An interpolation-based method for the verification of security protocols</article-title>
          ,
          <source>” Journal of Computer Security</source>
          , vol.
          <volume>25</volume>
          , no.
          <issue>6</issue>
          , pp.
          <fpage>463</fpage>
          -
          <lpage>510</lpage>
          ,
          <year>2017</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>P.</given-names>
            <surname>Upchurch</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Gardner</surname>
          </string-name>
          , G. Pleiss,
          <string-name>
            <given-names>R.</given-names>
            <surname>Pless</surname>
          </string-name>
          ,
          <string-name>
            <given-names>N.</given-names>
            <surname>Snavely</surname>
          </string-name>
          , K. Bala, and
          <string-name>
            <given-names>K.</given-names>
            <surname>Weinberger</surname>
          </string-name>
          , “
          <article-title>Deep feature interpolation for image content changes</article-title>
          ,”
          <source>in Proceedings of the IEEE conference on computer vision and pattern recognition</source>
          ,
          <year>2017</year>
          , pp.
          <fpage>7064</fpage>
          -
          <lpage>7073</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>S.</given-names>
            <surname>Niklaus</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Mai</surname>
          </string-name>
          , and
          <string-name>
            <given-names>F.</given-names>
            <surname>Liu</surname>
          </string-name>
          , “
          <article-title>Video frame interpolation via adaptive convolution</article-title>
          ,”
          <source>in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition</source>
          ,
          <year>2017</year>
          , pp.
          <fpage>670</fpage>
          -
          <lpage>679</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>L.</given-names>
            <surname>Yang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Chao</surname>
          </string-name>
          , and
          <string-name>
            <given-names>Q.</given-names>
            <surname>Shen</surname>
          </string-name>
          , “
          <article-title>Generalized adaptive fuzzy rule interpolation</article-title>
          ,
          <source>” IEEE Transactions on Fuzzy Systems</source>
          , vol.
          <volume>25</volume>
          , no.
          <issue>4</issue>
          , pp.
          <fpage>839</fpage>
          -
          <lpage>853</lpage>
          ,
          <year>2016</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>R.</given-names>
            <surname>Damaševicˇius</surname>
          </string-name>
          , C. Napoli, T. Sidekerskiene˙, and
          <string-name>
            <surname>M.</surname>
          </string-name>
          <article-title>Woz´niak, “Imf mode demixing in emd for jitter analysis</article-title>
          ,
          <source>” Journal of Computational Science</source>
          , vol.
          <volume>22</volume>
          , pp.
          <fpage>240</fpage>
          -
          <lpage>252</lpage>
          ,
          <year>2017</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>M.</given-names>
            <surname>Wozniak</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Napoli</surname>
          </string-name>
          , E. Tramontana, G. Capizzi,
          <string-name>
            <given-names>G. L.</given-names>
            <surname>Sciuto</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R. K.</given-names>
            <surname>Nowicki</surname>
          </string-name>
          , and
          <string-name>
            <given-names>J. T.</given-names>
            <surname>Starczewski</surname>
          </string-name>
          , “
          <article-title>A multiscale image compressor with rbfnn and discrete wavelet decomposition</article-title>
          ,” in
          <source>2015 International Joint Conference on Neural Networks (IJCNN)</source>
          . IEEE,
          <year>2015</year>
          , pp.
          <fpage>1</fpage>
          -
          <lpage>7</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>