Lacunary Interpolation Using Quartic B-Spline

Similar documents
Deficient Quartic Spline Interpolation

Convergence of C 2 Deficient Quartic Spline Interpolation

INTERPOLATION BY QUARTIC SPLINE

Spline Solution Of Some Linear Boundary Value Problems

Solve Non-Linear Parabolic Partial Differential Equation by Spline Collocation Method

Approximations for Linear Tenth-order Boundary Value Problems through Polynomial and Non-polynomial Cubic Spline Techniques

Interpolation by Spline Functions

Cubic Spline Interpolation with New Conditions on M 0 and M n

Polynomials tend to oscillate (wiggle) a lot, even when our true function does not.

A Quintic Spline method for fourth order singularly perturbed Boundary Value Problem of Reaction-Diffusion type

Natural Quartic Spline

Generalised Mean Averaging Interpolation by Discrete Cubic Splines

1, 2

On the deviation of a parametric cubic spline interpolant from its data polygon

A Transformation Based on the Cubic Parabola y = x 3

Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon

Positivity Preserving Interpolation of Positive Data by Rational Quadratic Trigonometric Spline

Linear Interpolating Splines

MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA

Table for Third-Degree Spline Interpolation Using Equi-Spaced Knots. By W. D. Hoskins

Adaptive osculatory rational interpolation for image processing

Discrete Cubic Interpolatory Splines

specified or may be difficult to handle, we often have a tabulated data

SPLINE COLLOCATION METHOD FOR SOLUTION OF HIGHER ORDER LINEAR BOUNDARY VALUE PROBLEMS

Interpolation & Polynomial Approximation. Cubic Spline Interpolation II

Interpolation and Splines

Research Article Data Visualization Using Rational Trigonometric Spline

A Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings

Knot Insertion and Reparametrization of Interval B-spline Curves

99 International Journal of Engineering, Science and Mathematics

A New Algorithm for Developing Block Methods for Solving Fourth Order Ordinary Differential Equations

Video 11.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar

Approximation of a Fuzzy Function by Using Radial Basis Functions Interpolation

A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions

A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions

Spline Collocation Method for Solving Three Dimensional Parabolic Equations

Cubic B-spline Solution of Two-point Boundary Value Problem Using HSKSOR Iteration

Lecture 9: Introduction to Spline Curves

Lecture VIII. Global Approximation Methods: I

8 Piecewise Polynomial Interpolation

SPLINE APPROXIMATION VIA THE CONTROL POLYGON

February 2017 (1/20) 2 Piecewise Polynomial Interpolation 2.2 (Natural) Cubic Splines. MA378/531 Numerical Analysis II ( NA2 )

Rational Bezier Surface

Lecture 2.2 Cubic Splines

ASSESSMENT OF COMPLEX VARIABLE BASIS FUNCTIONS IN THE APPROXIMATION OF IDEAL FLUID FLOW PROBLEMS

Almost Curvature Continuous Fitting of B-Spline Surfaces

Computational Physics PHYS 420

Asymptotic Error Analysis

Polynomial Approximation and Interpolation Chapter 4

Special Quartics with Triple Points

Interactive Graphics. Lecture 9: Introduction to Spline Curves. Interactive Graphics Lecture 9: Slide 1

Parameterization for curve interpolation

Computation of interpolatory splines via triadic subdivision

Important Properties of B-spline Basis Functions

American International Journal of Research in Science, Technology, Engineering & Mathematics

Piecewise Polynomial Interpolation, cont d

Wavelet-Galerkin Solutions of One and Two Dimensional Partial Differential Equations

Second Triangular Hermite Spline Curves and Its Application

An introduction to NURBS

International Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research)

RATIONAL CURVES ON SMOOTH CUBIC HYPERSURFACES. Contents 1. Introduction 1 2. The proof of Theorem References 9

Approximation of 3D-Parametric Functions by Bicubic B-spline Functions

Local Approximation by Splines with Displacement of Nodes

CHAPTER 5 NUMERICAL INTEGRATION METHODS OVER N- DIMENSIONAL REGIONS USING GENERALIZED GAUSSIAN QUADRATURE

Constrained modification of the cubic trigonometric Bézier curve with two shape parameters

Effectiveness of Element Free Galerkin Method over FEM

An explicit feature control approach in structural topology optimization

Normals of subdivision surfaces and their control polyhedra

MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix.

Hyperbola for Curvilinear Interpolation

A Note On The Sparing Number Of The Sieve Graphs Of Certain Graphs

Fourier Transform of Bernstein-Bézier Polynomials

Algorithmic Aspects of Acyclic Edge Colorings

Quasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters

Operators to calculate the derivative of digital signals

DEPARTMENT - Mathematics. Coding: N Number. A Algebra. G&M Geometry and Measure. S Statistics. P - Probability. R&P Ratio and Proportion

Optimization Problems Under One-sided (max, min)-linear Equality Constraints

lecture 10: B-Splines

Element Quality Metrics for Higher-Order Bernstein Bézier Elements

(Refer Slide Time: 00:02:24 min)

(i, j)-almost Continuity and (i, j)-weakly Continuity in Fuzzy Bitopological Spaces

A NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME

Using Modified Euler Method (MEM) for the Solution of some First Order Differential Equations with Initial Value Problems (IVPs).

APPM/MATH Problem Set 4 Solutions

Sci.Int.(Lahore),29(6), ,2017 ISSN ;CODEN: SINTE

Component Connectivity of Generalized Petersen Graphs

AM205: lecture 2. 1 These have been shifted to MD 323 for the rest of the semester.

MEI Desmos Tasks for AS Pure

Transform Extreme Point Multi-Objective Linear Programming Problem to Extreme Point Single Objective Linear Programming Problem by Using Harmonic Mean

08 - Designing Approximating Curves

Math 225 Scientific Computing II Outline of Lectures

Gracefulness of a New Class from Copies of kc 4 P 2n and P 2 * nc 3

A numerical grid and grid less (Mesh less) techniques for the solution of 2D Laplace equation

EECS 556 Image Processing W 09. Interpolation. Interpolation techniques B splines

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24

MATLAB. Advanced Mathematics and Mechanics Applications Using. Third Edition. David Halpern University of Alabama CHAPMAN & HALL/CRC

Constructing Cubic Splines on the Sphere

F A S C I C U L I M A T H E M A T I C I

LINEAR CODES WITH NON-UNIFORM ERROR CORRECTION CAPABILITY

Assessing the Quality of the Natural Cubic Spline Approximation

Transcription:

General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan Hama Faraj Jwamer and 2 Bzhar Jamal 1,2 University of Sulaimani, College of Science, Department of Mathematics, Kurdistan Region, Sulaimani, Iraq karwan.jwamer@univsul.edu.iq Abstract In this paper, a quadratic B-Spline has been constructed that is an approximate solution to a function with very limited given lacunary data and approximate boundary conditions. Then the error bound for the B-spline is found. Further this construction is used to solve boundary value problem. Keywords: B-spline, error bounds, boundary value problems, approximate solution. 2000 MSC No: 65D05, 65D07. 1. Introduction Boundary value problem (BVP) is a differential equation together with some conditions imposed at different points. The BVPs have been applied in mathematics, engineering and various fields of sciences. The rapid increasing of their applications has led to formulating and upgrading several existing methods and new approaches [1]. Quartic B-spline is a piecewise polynomial of degree four satisfying third order parametric continuity [2]. In this paper, quartic B-spline is manipulated to approximate the solution of a BVP with lacunary data given on it and the boundary conditions are approximation, means that it satisfies conditions (1)- (2), see below. By presuming the B-spline to be the solution for this problem, an undetermined system of linear equations of order with being the number of uniform subintervals is built. Adding three approximate boundary conditions into this system gives a square system of that have a unique solution for this problem. This method can make use of the problem s equation to construct an error equation. Minimization of the error equation would give the value of the variable that produces the best approximation of the solution. 2. Lacunary Interpolations on Boundary Value Problem using B-spline of Degree four. In this section, we try to find an approximate solution of a function having very limited lacunary data on it, the function has the following interpolation conditions;

130 Karwan Hama Faraj Jwamer et al. With the following lacunary boundary conditions (1) The B-spline is non-zero at five knots, we can find the value of at the nodal points by differentiating it with respect to the value (x) and its first three derivatives at the nodal points can be tabulated as in table 1. Table 1: Values of (x) and its first three derivatives at the nodal points 0 1 11 11 1 0 0 0 0 0 0 0 From the definition of B-Spline, the solution of Eq.1 using Quartic B-spline is approximated by; (2) then and Now without loss of generality, we can re-write Eq. 3 as below, where and all other s are zero, By shifting the B spline to the right side by step, mathematically meaning; for all, then Eq. 3 can be re-written as follows

Lacunary Interpolation Using Quartic B-Spline 131 Doing the same steps on first and second derivatives of the B-spline in Eqs.4 and 5 respectively, gives and From the lacunary conditions and by substituting the values of following equations are formulated; at the knots using table 1, the.. (9) From Eq.9, there are unknowns to be founded and equations, so we need to write three more equations which are the lacunary boundary conditions, to have; (10) Eqs. 9 and 10 forms a square system of follows system represented in the matrix form as

132 Karwan Hama Faraj Jwamer et al. 1 11 11 1 0 0 0 0... 0 ) -4-12 12 4 0 0 0 0... 0 ) 12-12 -12 12 0 0 0 0... 0 ) 0 1 11 11 1 0 0 0... 0 ) 0 0 1 11 11 1 0 0 0 ) 0 0 0 1 11 11 1 0... 0 = ) (11).............................. 0 0... 0 0 0 0 1 11 11 1 ) 0 0... 0 0 0 0-4 -12 12 4 ) Above Matrix is invertible, to show this we can use the principal of mathematical induction, this yields that above system has a unique solution. With this, the construction of a fourth degree B-spline is completed which is an approximate solution of the problem given by Eqs.1 and 2. 3. Error bound In this section, an error bound of the fourth degree B-spline is formulated Denote Eq.11 by, and Let be an approximate solution of, and be an approximate solution of.

Lacunary Interpolation Using Quartic B-Spline 133 Lemma 1 If is an matrix as defined in Eq.11, then Proof: It easy to show that by mathematical induction the matrix A is invertible, then It is known that, or =, then ( ) Lemma 2 Let be another approximate solution of using exact boundary conditions, then, Proof: For the approximate solution of, another matrix system could be obtained as follows;, and, Then ( ) ( ) ( ) ( ) Now using properties of norm, the following yields; Lemma 3 The following inequalities are true for i), ii) Proof For iii) and iv), and from Table 1, for i), ii), iii), iv). Lemma 4 Let be a quartic B-spline approximate solution of with lacunary and approximate boundary conditions, and let Ŝ(x) be another approximate solution of with lacunary boundary conditions, then for, the following are true

134 Karwan Hama Faraj Jwamer et al. i), ii), iii), iv). Proof Using Lemma 2 and Lemma 3, and for, we have i) For, For ii) For, iii) For ; iv) Finally for ;, we have the following expansions, that will be used in subsequent section. ( ) (12) where Theorem 1 Let be a fourth degree B-Spline approximation solution of, then the following inequalities are true; i) ( ) ii) ( ( )) iii) Proof: Using lemma 4, prove of (iv) is as follows; Assigning as and as,and using Eq. 12, then the last term of above equation can be obtained as follows; ( ) ( )

Lacunary Interpolation Using Quartic B-Spline 135 Hence ( ( )) where Similar to the proof of (iv) the followings can be proved; iii) Assigning as and as, then the last term of the above equation can be obtained as follows; ( ) ( ) ( ( )) Hence ( ( ))where. Applying similar techniques complete the proof of the theorem ii) ( ) i) ( ) Example 1: The analytical solution is given by Table 2 compares the numerical results between present B-Spline method and the computational method used in [12]

136 Karwan Hama Faraj Jwamer et al. Table 2 Exact Solution [12 ] B-Spline Max Error 1/8-0.1239381121 2.37 8.94 1/16-0.0437870146 5.75 3.46 1/32-0.0312344196 1.47 5.59 Error of the present B-spline method is slightly better than the error obtained in [12], this confirms that the B- Spline method is a precise one for the models of BVPs stated in section 2. Example 2: with The analytical solution of above problem is; The results of maximum absolute error [ ] are tabulated in Table 3 Table 3 Exact Solution [12 ] B-spline Error 1/8 0.07372 1.29 7.53 1/16 0.08767 3.08 4.21 1/32 0.08858 7.54 3.17 The example has been solved by Yogesh and Punkja. [12], The numerical results shown in Table 3 shows encouraging results of our method. Conclusion Based on the investigation done in the present study, it can be concluded that, B-Spline is fast, flexible and precise to be used to find approximate solutions of boundary value problems with given limited lacunary interpolation condition and approximate boundary condition, below are some conclusions; 1) B- Spline produced an approximation of analytical solution of the problem with respect to the selected subinterval. 2) B-Spline is a good tool to be used to solve lacunary interpolation problems for Boundary Value Problems. 3) It is considerable that as the subintervals are increasing, and h being small, the approximation by B-Spline is more precise and has potential to give good approximation solution for boundary value problems. 4) B-Spline of degree four is precise and flexible enough to become an approximate solution of BVPs.

Lacunary Interpolation Using Quartic B-Spline 137 5) The errors in the numerical examples are quite good, slightly better than other errors obtained from other approximation tools. References [1] S. Shafie and A. A. Majid, Approximation of cubic B- spline interpolation method, shooting and finite difference methods for linear problems on solving linear two-point boundary value. 2012. [2] R. K. Saeed, K. H. Jwamer, and F. K. Hamasalh, Introduction to numerical analysis, University of Sulaimani, Sulaimani, Kurdistan Region- Iraq. 2015. [3] V. Dahiya, Exploring B-spline functions for numerical solution of mathematical problems, Int. J. of Multidisciplinary Research and Development, 2(1) 2015, : 452-458 [4] P. J. Davis, Interpolation and Approximation, Blaisdell, New York. 1961. [5] C. De-Boor, A practical guide to splines, Springer-Verlag, Berlin. 1978. [6] Y. P. Dubey, Best error bounds for splines of degree six, Int. Journal of Math. Analysis, 5(24) 2011: 1201-1209. [7] Y. P. Dubeyand A. Shukla, The deficient C 1 quartic spline interpolation, Research Inventy, Int. J. of Engineering and Sciences, 2(9) 2013, 24-30. [8] C. A. Hall and W. W. Meyer, Optimal error bounds for cubic spline interpolation, J. Approx. Theory, 16: (1976) 105-122. [9] K. H. Jwamer, Minimizing error bounds in (0,2,3) lacunary interpolation by sextic spline function, Journal of Mathematics and Statistics, 3(4) 2007,: 249-256. [10] S. S. Rana and Y. P. Dubey, Best error bounds of deficient quantic spline interpolation, Indian J. Pure Appl. Math., 28(10) 1997, 1337-1344. [11] Saeed, R.K. (1990), A study of lacunary interpolation by splines, MSC Thesis, Salahaddin University/Erbil, Iraq. [12] Yogesh Gupta et al, A Computational Method for Solving Two Point Boundary Value Problems of Order Four. Int. J. Comp. Tech. Appl, 2 (5) 2011, 1426-1431