• 정보과학회지
• /
• 제1권1호
• /
• pp.72-74
• /
• 1983

Traditionally, polynomials have been used to approximte functions with prescribed values at a number of points(called the knots) on a given interal on the real line. The method of splines recently developed is more flexible. It approximates a function in a piece-wise fashion, by means of a different polynomial in each subinterval. The cubic spline gas ets origins in beam theory. It possessed continuous first and second deriatives at the knots and is characterised by a minimum curvature property which es rdlated to the physical feature of minimum potential energy of the supported beam. Translated into mathematical terms, this means that between successive knots the approximation yields a third-order polynomial sith its first derivatives continuous at the knots. The minimum curvature property holds good for each subinterval as well as for the whole region of approximation This means that the integral of the square of the second derivative over the entire interval, and also over each subinterval, es to be minimized. Thus, the task of determining the spline lffers itself as a textbook problem in discrete computer programming, since the integral of ghe square of the second derivative can be obviously recognized as the criterion function whicg gas to be minimized. Starting with the initial value of the function and assuming an initial solpe of the curve, the minimum norm property of the curvature makes sequential decision of the slope at successive knots (points) feasible. It is the aim of this paper to derive the cubic spline by the methods of computer programming and show that the results which is computed the all the alues in each subinterval of the spline approximations.

• 한국근적외분광분석학회:학술대회논문집
• /
• 한국근적외분광분석학회 2001년도 NIR-2001
• /
• pp.1041-1041
• /
• 2001

Removal of variability in spectra data before the application of chemometric modeling will generally result in simpler (and presumably more robust) models. Particularly for sparsely sampled data, such as typically encountered in diode array instruments, the use of Savitzky-Golay (S-G) derivatives offers an effective method to remove effects of shifting baselines and sloping or curving apparent baselines often observed with scattering samples. The application of these convolution functions is equivalent to fitting a selected polynomial to a number of points in the spectrum, usually 5 to 25 points. The value of the polynomial evaluated at its mid-point, or its derivative, is taken as the (smoothed) spectrum or its derivative at the mid-point of the wavelength window. The process is continued for successive windows along the spectrum. The original paper, published in 1964 [1] presented these convolution functions as integers to be used as multipliers for the spectral values at equal intervals in the window, with a normalization integer to divide the sum of the products, to determine the result for each point. Steinier et al. [2] published corrections to errors in the original presentation [1], and a vector formulation for obtaining the coefficients. The actual selection of the degree of polynomial and number of points in the window determines whether closely situated bands and shoulders are resolved in the derivatives. Furthermore, the actual noise reduction in the derivatives may be estimated from the square root of the sums of the coefficients, divided by the NORM value. A simple technique to evaluate the actual convolution factors employed in the calculation by the software will be presented. It has been found that some software packages do not properly account for the sampling interval of the spectral data (Equation Ⅶ in [1]). While this is not a problem in the construction and implementation of chemometric models, it may be noticed in comparing models at differing spectral resolutions. Also, the effects on parameters of PLS models of choosing various polynomials and numbers of points in the window will be presented.

• Structural Engineering and Mechanics
• /
• 제2권3호
• /
• pp.245-256
• /
• 1994

A new axisymmetric crack model is proposed on the basis of p-version of the finite element method limited to theory of small scale yielding. To this end, axisymmetric stress element is formulated by integrals of Legendre polynomial which has hierarchical nature and orthogonality relationship. The virtual crack extension method has been adopted to calculate the stress intensity factors for 3-D axisymmetric cracked bodies where the potential energy change as a function of position along the crack front is calculated. The sensitivity with respect to the aspect ratio and Poisson locking has been tested to ascertain the robustness of p-version axisymmetric element. Also, the limit value that is an exact solution obtained by FEM when degree of freedom is infinite can be estimated using the extrapolation equation based on error prediction in energy norm. Numerical examples of thick-walled cylinder, axisymmetric crack in a round bar and internal part-thorough cracked pipes are tested with high precision.

• 대한수학회지
• /
• 제43권3호
• /
• pp.579-591
• /
• 2006

In this paper we obtain a sufficient and necessary condition for an analytic function f on the unit ball B with Hadamard gaps, that is, for $f(z)\;=\;{\sum}^{\infty}_{k=1}\;P_{nk}(z)$ (the homogeneous polynomial expansion of f) satisfying $n_{k+1}/n_{k}{\ge}{\lambda}>1$ for all $k\;{\in}\;N$, to belong to the weighted Bergman space $$A^p_{\alpha}(B)\;=\;\{f{\mid}{\int}_{B}{\mid}f(z){\mid}^{p}(1-{\mid}z{\mid}^2)^{\alpha}dV(z) < {\infty},\;f{\in}H(B)\}$$. We find a growth estimate for the integral mean $$\({\int}_{{\partial}B}{\mid}f(r{\zeta}){\mid}^pd{\sigma}({\zeta})\)^{1/p}$$, and an estimate for the point evaluations in this class of functions. Similar results on the mixed norm space $H_{p,q,{\alpha}$(B) and weighted Bergman space on polydisc $A^p_{^{\to}_{\alpha}}(U^n)$ are also given.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제23권3호
• /
• pp.211-236
• /
• 2019

This paper is devoted to the study and investigation of the Runge-Kutta discontinuous Galerkin method for a system of differential equations consisting of two hyperbolic conservation laws. The numerical coupling flux which is used at a given interface (x = 0) is the upwind flux. Moreover, in the linear case, we derive optimal convergence rates in the $L_2$-norm, showing an error estimate of order ${\mathcal{O}}(h^{k+1})$ in domains where the exact solution is smooth; here h is the mesh width and k is the degree of the (orthogonal Legendre) polynomial functions spanning the finite element subspace. The underlying temporal discretization scheme in time is the third-order total variation diminishing Runge-Kutta scheme. We justify the advantages of the Runge-Kutta discontinuous Galerkin method in a series of numerical examples.

• 대한수학회논문집
• /
• 제9권2호
• /
• pp.467-489
• /
• 1994

Let $\Omega$ be a closed and bounded polygonal domain in R$^2$, or a closed line segment in R$^1$ with boundary $\Gamma$, such that there exists an invertible mapping T : $\Omega$ \longrightarrow $\Omega$ with the following correspondence: x$\in$$\Omega$ ↔ x = T(x) $\in$$\Omega$, (1.1) and (1.2) t $\in$ U$\sub$p/($\Omega$) ↔ t = to T$\^$-1/ $\in$ U$\sub$p/($\Omega$), where $\Omega$ denotes the corresponding reference elements I = [-1,1] and I ${\times}$ I in R$^1$ and R$^2$ respectively, (1.3) U$\sub$p/($\Omega$) = {t : t is a polynomial of degree $\leq$ p in each variable on $\Omega$}, and (1.4) U$\sub$p/($\Omega$) = {t : t = to T $\in$ U$\sub$p/($\Omega$)}.(omitted)

• 대한토목학회논문집
• /
• 제12권4_1호
• /
• pp.67-76
• /
• 1992

축대칭(軸對稱) 선형강성(線形彈性) 응력해석을 위해 p-version 유한요소법에 기초한 계층적(階層的) 정식화 과정이 제안되었다. 이 방식은 적분형 르장드르 다항식을 사용하여 절점좌표값을 갖지 않는 절점을 추가하여 형상함수의 조합형태로 변위함수(變位)를 근사시키는 방법이다. 형상함수(形狀函數)가 계층적 성질을 갖기 때문에 강성도(剛性度)행렬과 하중벡터도 계층적이 된다. 본 연구에서 제안된 요소(要所)의 장점(長點)은 다음과 같다. 첫째, 개선된 수치연산의 효율성이며 둘째, 요소간에 서로 다른 차수(次數)의 형상함수를 사용할 수 있고 셋째, p-세분화를 할 때 저차(低次)일 때 계산된 값을 그대로 사용할 수 있다. 수치예제를 통해 제안된 요소의 정확도(正確度), 효율성(效率性), 모델링의 간편성(簡便性), 적용성(適用性) 및 변위와 응력 그리고 에너지 Norm등을 사용하여 그 우월성을 입증하고 있다. 몇 가지 예제의 해석결과는 이미 발표된 논문과 아울러 해석적 방법에 의한 결과와 비교되었다.

• 한국컴퓨터정보학회논문지
• /
• 제24권10호
• /
• pp.57-64
• /
• 2019

이 논문에서는 인터벌 그래프 컬러링 역문제 중 다항시간 안에 풀이 가능한 경우에 대해 연구한다. 인터벌 그래프의 컬러링 역문제는 주어진 인터벌 그래프를 K개의 서로 다른 색깔로 색칠할 수 없는 경우를 가정하며, 다음과 같이 정의된다. 주어진 인터벌 그래프가 K개의 색깔을 이용해서 모두 칠해질 수 있도록 인터벌 그래프와 연관되어 있는 인터벌 시스템을 최소한으로 수정하는 문제이다. 인터벌 시스템에서 두 인터벌이 부분적으로라도 서로 겹쳐있는 구간이 있을 경우 두 인터벌에 해당하는 노드들이 엣지로 연결되어 있음을 의미하고, 따라서 이 경우에는 해당 노드들을 같은 색깔을 이용해 칠할 수 없다. 따라서 겹쳐져 있는 인터벌들을 이동시켜 해당 그래프의 chromatic number를 바꿀 수 있다. 본 논문에서는 인터벌의 길이가 모두 1 또는 2이며, 인터벌의 이동이 본래 위치 대비 오른쪽으로만 가능하다는 제한이 있는 경우에 대해 집중 탐구한다. 이 문제를 해결하는 다항시간 알고리즘으로 sorting과 선입선출 방식을 사용하는 2단계 알고리즘을 제안한다.