[논문 리뷰] Least-squares Solutions of Linear Differential Equations

Daniele Mortari|arXiv (Cornell University)|2017. 02. 25.

Differential Equations and Numerical Methods참고 문헌 6인용 수 11

한 줄 요약

이 논문은 비정상계수를 가진 비동차 선형 미분방정식을 해결하기 위한 최소제곱법을 제시한다. 이 방법은 미지의 함수 g(t)를 포함하는 제약 조건을 내장한 표현식을 사용하며, g(t)는 체비셰프 다항식을 통해 근사된다. 이 접근법은 초기값 문제에서 MATLAB의 ODE45보다 다섯 자리 수준 높은 정확도를 달성하며, 잔차와 조건수 분석을 통해 경계값 문제에서 해가 없음, 유일함, 무한해임을 효과적으로 구분한다.

ABSTRACT

This study shows how to obtain least-squares solutions to initial and boundary value problems to nonhomogeneous linear differential equations with nonconstant coefficients of any order. However, without loss of generality, the approach has been applied to second order differential equations. The proposed method has two steps. The first step consists of writing a constrained expression, introduced in Ref. \cite{Mortari}, that has embedded the differential equation constraints. These expressions are given in term of a new unknown function, $g (t)$, and they satisfy the constraints, no matter what $g (t)$ is. The second step consists of expressing $g (t)$ as a linear combination of $m$ independent known basis functions, $g (t) = \mathbfξ^T \mathbf{h} (t)$. Specifically, Chebyshev orthogonal polynomials of the first kind are adopted for the basis functions. This choice requires rewriting the differential equation and the constraints in term of a new independent variable, $x\in[-1, +1]$. The procedure leads to a set of linear equations in terms of the unknown coefficients vector, $\mathbfξ,$ that is then computed by least-squares. Numerical examples are provided to quantify the solutions accuracy for initial and boundary values problems as well as for a control-type problem, where the state is defined in one point and the costate in another point.

연구 동기 및 목표

• 비반복적이고 최소제곱 기반의 방법을 개발하여 비정상계수를 가진 선형 미분방정식과 임의의 초깃값 또는 경계조건을 해결한다.
• 닫힌형 해가 존재하지 않을 경우 비동차 선형 미분방정식에 대한 일반적인 해법의 부족을 해결한다.
• 동일한 제약 조건 표현 형식을 사용하여 초깃값 문제와 경계값 문제를 통합된 프레임워크로 다룰 수 있도록 한다.
• 잔차와 조건수 분석을 통해 알려진 해, 알려지지 않은 해, 해가 없음, 또는 무한해임을 정확하게 식별할 수 있도록 한다.
• 향후 비선형 방정식, 편미분방정식, 최적의 기저 선택 등으로의 확장을 위한 기반을 마련한다.

제안 방법

• 이 방법은 특수 함수 βi(t, t̲)를 통해 미분방정식과 경계/초기 조건을 내장함으로써, 해가 반드시 만족하도록 제약 조건이 포함된 표현식을 사용한다.
• 미지의 해 y(t)는 y(t) = g(t) + Σβi(t, t̲)(y_t̲i^(d̲i) - g_t̲i^(d̲i))로 표현되며, 여기서 g(t)는 근사되어야 할 자유 함수이다.
• 함수 g(t)는 g(t) = ξᵀh(t)로 전개되며, h(t)는 첫 번째 종류의 m개의 정규수직 체비셰프 다항식으로 이루어진 벡터이다.
• 미분방정식은 g(t)에 대해 재기록되며, 그 결과 방정식은 N개의 점에서 샘플링되어 ξ에 대한 선형 시스템을 형성한다.
• 계수 벡터 ξ는 샘플링 점들에서 잔차의 최소제곱 최소화를 통해 계산된다.
• 이 방법은 일반성을 잃지 않고 2차 미분방정식에 적용되었으며, 이는 임의의 차수로 확장 가능하다.

실험 결과

연구 질문

• RQ1 비정상계수를 가진 선형 미분방정식과 임의의 제약 조건에 대해 비반복적이고 최소제곱 기반의 방법을 개발할 수 있는가?
• RQ2 미지의 함수에 관계없이 항상 만족하도록 제약 조건을 해의 형태에 직접 내장할 수 있는가?
• RQ3 체비셰프 다항식이 이러한 방정식의 해를 근사하는 데 안정적이고 정확한 기저로 기능할 수 있는가?
• RQ4 이 방법은 경계값 문제에서 해가 유일함, 없음, 또는 무한해임을 어떻게 구분하는가?
• RQ5 기존의 ODE45와 같은 표준 ODE 솔버와 비교해 최소제곱 해의 정확도와 조건수 행동은 어떠한가?

주요 결과

• 최소제곱 방법은 알려진 해가 있는 초깃값 문제에서 MATLAB의 ODE45보다 다섯 자리 수준 높은 정확도를 달성한다.
• 최소제곱 시스템의 잔차 평균과 조건수는 경계값 문제에서 해가 없음과 무한해임을 신뢰성 있게 구분할 수 있다.
• 해의 오차는 전통적 방법과 달리 통합 구간 전체에 걸쳐 약간의 균일한 분포를 보이며, 시간이 지남에 따라 오차가 누적되지 않는다.
• 이 방법은 초깃값 문제와 경계값 문제에 동일하게 적용 가능하여 통합된 해법 프레임워크를 제공한다.
• 이 방법은 반복 계산이 필요 없고 기본 해를 구하는 데도 필요 없기 때문에, 닫힌형 해가 존재하지 않는 방정식에 적합하다.
• 비선형 방정식, 편미분방정식, 가중 최소제곱 형식으로의 확장 가능성은 보여주지만, 아직 연구가 진행 중인 분야이다.