Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

ABOUT PROPERTIES OF ORTHOGONAL FINITE FUNCTIONS AND ABOUT THEIR APPLICATION IN ALGORITHMS OF NUMERICAL METHODS

Leontev V.L. 1
1 Ulyanovsk State University
The article is devoted to development of fundamental direction of approximation theory. This direction concern to construction of systems of grid basis finite functions with compact supports. The properties of orthogonal functions on every concrete grid are investigated in the book [1]. Some orthogonal finite functions [1] are the generalization of the theory of B-splines and the addition of the theory of orthogonal finite wavelets [2]. Analogies do not exist for some orthogonal finite functions [1] in 2D and 3D cases. The possibilities of formation of new properties of orthogonal finite functions and the possibilities of application of such functions in algorithms of numerical methods are considered. In first part of the article it is supposed the existence of the system of basis orthogonal grid finite functions, which are orthogonal also to first derivatives of all another functions of this basis on every grid. The theorem about impossibility of application of such functions in algorithms of numerical methods is proved. The possibility of construction of such systems of basis finite functions is investigated in second part of the article. The modification of grid basis finite functions is made for this purpose. It is proved, that the modifications of such type are not give grid basis finite functions with new properties.
orthogonal grid finite functions
basis functions
numerical methods
1. Leontev V.L. Ortogonal’nye finitnye funkcii i chislennye metody / V.L.Leontev. – Uljanovsk: Izd-vo Uljanovskogo gos. un-ta, 2003. 178 р.
2. Leontev V.L. Metody konechnykh ehlementov, osnovannye na ispolzovanii obobshhennykh funkcijj Kuranta v teorii uprugikh kolebanijj // Problemy dinamiki i prochnosti ehlektro- i ehnergo-mashin: tezisy dokl. Vseros. nauchnogo seminara (S.Peterburg, 18-20 maja 1993 g.). S.-Peterburg: izd-vo Instituta problem mashinovedenija RAN, 1993. рp. 21–22.
3. Leontev V.L. Ob odnom obobshhenii funkcijj Kuranta // Teorija funkcijj i priblizhenijj: trudy 7-jj Saratovskojj zimnejj shkoly (30 janvarja – 4 fevralja 1994 g.). Chast 3. Saratov: Izd-vo SGU, 1995. рp. 36-40.
4. Leontev V.L., Lukashanec N.Ch. Setochnye bazisy ortogonalnykh finitnykh funkcijj // Zhurnal vychislitel’nojj matematiki i matematicheskojj fiziki. 1999. T. 39, no. 7. рр. 1158–1168.
5. Leontev V.L. Ob ortogonal’nykh finitnykh funkcijakh i o chislennykh metodakh, svjazannykh s ikh primeneniem // Obozrenie prikladnojj i promyshlennojj matematiki. 2002. T. 9, vyp. 3. pp. 497–504.
6. Daubechies I. Comm. Pure and Appl. Math., 1988, 41, pp. 909–996.

Известны [1] сеточные базисы ортогональных финитных функций (ОФФ) (первый такой базис построен В.Л. Леонтьевым в 1993 году [2], другие подобные базисы созданы и исследованы затем, например, в [3, 4, 5]), состоящие из финитных непрерывных функций, и такие, что скалярные произведения функционального пространства L2 любых двух различных функций каждого такого базиса на каждой конкретной сетке равны нулю. Ранее (в 1988 году) были построены финитные ортогональные вейвлеты И. Добеши [6], но они в отличие от базисных функций [1] не могут быть записаны в аналитической форме, строятся приближенно с заданной точностью с помощью итерационного алгоритма кратноразрешающего анализа, не обладают симметрией (четностью, нечетностью), имеют существенно более сложную структуру и низкую гладкость. Базисные функции [6] построены с целью их использования в спектральном анализе сигналов, где успешно применяются, но не адаптированы к их использованию в алгоритмах численно-аналитических методов (проекционно-сеточных (ПСМ) и вариационно-сеточных (ВСМ)). Функции [1] созданы специально для их применения в алгоритмах численно-аналитических методов, связанных со смешанными вариационными принципами (Рейсснера, Ху-Васидзу) или с обобщающими эти вариационные принципы смешанными проекционными условиями. Показано [1], что использование ОФФ устраняет основной недостаток смешанных ВСМ или ПСМ – в несколько раз более высокое число узловых неизвестных по сравнению с ВСМ, основанными на вариационном принципе Лагранжа. Ортогональность таких базисных финитных функций позволяет исключить большую часть узловых неизвестных величин до начала решения глобальной сеточной системы линейных алгебраических уравнений (СЛАУ) на ЭВМ (число исключаемых величин в 2-4 раза превышает число остающихся в СЛАУ узловых величин), тем самым по числу арифметических операций, необходимых для получения решения СЛАУ, ВСМ, основанные на вариационном принципе Рейсснера, оказываются близкими к ВСМ, основанным на вариационном принципе Лагранжа. При этом качество решений для всех величин (не только для основных неизвестных функций, например, перемещений, но и для их производных, например, деформаций и напряжений) в смешанном ВСМ или ПСМ является существенно более высоким и по точности, и по гладкости, чем в ВСМ, связанном с вариационным принципом Лагранжа.

1. В последние годы в связи с развитием алгоритмов численных методов возник вопрос: возможно ли построение базисных ОФФ, обладающих дополнительными свойствами. В [1] и более ранних работах В.Л. Леонтьева отмечалось, что ряд таких ОФФ имеет неиспользованные для придания функциям свойства ортогональности свободные параметры, которые можно использовать для появления у этих функций дополнительных свойств. Однако, оказывается, что эти дополнительные свойства не могут планироваться произвольно. Здесь показывается, что невозможно создание финитных функций Лагранжева базиса (в таком базисе каждая отдельная функция базиса определяет значение глобальной аппроксимации функции или ее производной только в одном узле), которые, во-первых, обладают свойствами: взаимно-ортогональными являются сами базисные функции, а также каждая функция является ортогональной первым производным всех других функций этого базиса на каждой конкретной сетке, и, во-вторых, предназначены для полноценного использования их в проекционных условиях.

Пусть

{… < xj–2 < xj–1 < xj < xj + 1 < xj + 2 < …}

– координаты узлов сетки, расположенных на действительной оси.

Теорема. Пусть {φi(x)} – сеточный Лагранжев «базис» в  L2 (R) такой, что

φi(xj) = δij,

тогда если функции φi(x) обладают следующими свойствами ортогональности

Eqn134.wmf

Eqn135.wmf (1)

где A, B – некоторые числовые константы, то они не могут быть использованы для аппроксимации неизвестных функций при построении приближенных решений в проекционных условиях, связанных с невязками дифференциальных уравнений, содержащих неизвестные функции и их первые производные.

Доказательство.

Пусть f1(x), f2(x) – непрерывные функции, первая из которых строго возрастает, а вторая строго убывает на интервале
(xj–2, xj+2), причем

f1(xj) = f2(xj). (2)

Пусть также

Eqn136.wmf

Eqn137.wmf

– аппроксимации функций f1(x), f2(x) на интервале (xj–1, xj+1) линейными комбинациями функций φi(x). Здесь Eqn138.wmf – строго возрастающая функция; Eqn139.wmf – строго убывающая функция, причем

Eqn140.wmf

где Eqn141.wmf, Eqn142.wmf – постоянные коэффициенты, равные значениям f1(x), f2(x) в узлах xi.

Скалярные произведения первых производных этих аппроксимаций и φj(x)

Eqn143.wmf

и

Eqn144.wmf

равны в силу (1) и (2). Это означает, что в проекциях первых производных двух различных функций – строго возрастающей и строго убывающей, отличных от нуля и имеющих различные знаки, ‒ информация о первых производных двух различных в окрестности точки xj функций при проецировании аппроксимаций этих первых производных теряется. Таким образом, использование подобных функций в качестве базисных в проекционных условиях методов типа Бубнова‒Галеркина применительно к невязкам дифференциальных уравнений невозможно. Следовательно, расщепление системы сеточных уравнений на совокупность несвязанных уравнений с помощью такого ПСМ невозможно.

2. Рассмотрим общий случай. Пусть
f(x) ‒ материнская локальная (финитная) функция, обладающая на отрезке [–1, 1] свойствами

f(–1) = f(1) = 0; f(0) = A,

где A не равно нулю и равно, например, единице, и обладающая также свойством

f(x) + f(x – 1) = f(0) ∀ x ∈ [0, 1], (3)

необходимым для наличия у последовательности сеточных наборов функций, порожденных этой материнской функцией, полного набора свойств системы базисных функций. За пределами отрезка [–1, 1] функция f(x) тождественно равна нулю.

Предположим, что сеточные функции, порожденные материнской функции f(x), обладают не только известным свойством ортогональности [1] самих сеточных функций на каждой конкретной сетке (внутригрупповой взаимной ортогональностью финитных базисных функций)

Eqn145.wmf (4)

но и дополнительным свойством внутригрупповой ортогональности производных сеточных функций самим сеточным функциям

Eqn146.wmf (5)

Пусть на отрезке [–1, 1] материнская функция модифицируется, во-первых, некоторой финитной функцией αφ(x) (α – постоянный коэффициент), имеющей конечный носитель, не выходящий за пределы отрезка [0, 1], а, во-вторых, финитной функцией (–αφ(–x)) с конечным носителем, не выходящим за пределы отрезка [–1, 0]. Первая модифицирующая функция симметрична относительно линии x = 0,5, вторая – относительно линии x = –0,5. В этом случае функция f(x – 1) на отрезке [0, 1] имеет соответствующую модификацию, определяемую функцией (–αφ(x)) и поэтому свойство (3) выполняется для модифицированных сеточных функций автоматически в случае, если до модификации сеточные функции им обладают. Последнее здесь предполагается, поскольку рассматривается модификация базисных функций с целью изучения вопроса возможности возникновения у них дополнительного свойства, а именно, свойства ортогональности производных сеточных функций, с одной стороны, и самих сеточных функций, с другой стороны. Далее предполагается также, что модифицирующие функции непрерывны и имеют в точках x = 0,± 1 значения, равные нулю.

Условия (4), (5) с учетом такой модификации записываются в виде

Eqn147.wmf

Eqn148.wmf

Пусть

ψ(x) = f(x) + αφ(x)

на [0, 1], тогда последние два условия c учетом (3) записываются в виде

Eqn149.wmf

Eqn150.wmf

Выполнение первого из двух последних условий обеспечивается при определенном выборе числовых значений коэффициента α, которые, как показано в [1], существуют. Рассмотрим второе условие, которое переписывается в виде

Eqn151.wmf

или

Eqn152.wmf

откуда получается

Eqn153.wmf

поскольку

ψ(0) = f(0) = A ≠ 0; ψ(1) = f(1) = 0,

а следовательно, второе условие не выполняется при любом выборе модифицирующей функции φ(x) и значений параметра α.

Заключение

Известные примеры численно-аналитических «методов», в рамках которых последние три года строятся и используются «базисные» системы финитных функций, «обладающих» свойствами ортогональности как функций, так и их первых производных, являются неверными в своей основе. Такое заключение основано здесь не только на представленных формальных математических исследованиях, но и на том, что в итоге применения в рамках подобного «метода» процедуры проецирования, связанной с сеткой, каждое уравнение во внутреннем узле сетки оказывается не имеющим никакой связи ни с уравнениями в других узлах, ни с граничными условиями. Возникает естественный вопрос: почему подобное расщепление системы проекционных уравнений возможно для нефинитных функций, например, как известно, для тригонометрических? Это возможно потому, что, несмотря на отделение проекционных уравнений друг от друга в смысле наборов содержащихся в них коэффициентов из состава аппроксимирующих точное решение линейных комбинаций, каждое возникающее после применения процедуры проецирования уравнение, соответствующее отдельной базисной тригонометрической функции, учитывает внешнее воздействие на моделируемую систему на всей области краевой задачи в силу нефинитности базисной функции и при этом любое из таких уравнений может учитывать в общем случае граничные условия.

Рецензенты:

Вельмисов П.А., д.ф.-м.н., профессор, заведующий кафедрой «Высшая математика» Ульяновского государственного технического университета, г. Ульяновск;

Журавлев В.М., д.ф.-м.н., профессор, профессор кафедры «Теоретическая физика» Ульяновского государственного университета, г. Ульяновск;

Попов Ф.А., д.т.н., профессор, зам. директора по ИТ, Бийский технологический институт, г. Бийск.

Работа поступила в редакцию 26.10.2012.