Научный журнал

Фундаментальные исследования

ISSN 1812-7339

"Перечень" ВАК

ИФ РИНЦ = 1,674

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Леонтьев В.Л. 1

---

1 Ульяновский государственный университет

Статья посвящена развитию фундаментального направления теории аппроксимации. Это направление относится к построению систем сеточных базисных функций с компактными носителями. Свойство ортогональности таких функций на каждой конкретной сетке исследовано в книге [1]. Некоторые ортогональные финитные функции [1] являются обобщением теории В-сплайнов и дополнением теории ортогональных финитных вейвлетов [2]. Некоторые ортогональные финитные функции [1] в двумерном и трехмерном случаях не имеют аналогов. Здесь рассматриваются возможности формирования новых свойств ортогональных финитных функций и возможности использования функций c такими свойствами в алгоритмах численных методов. В первой части статьи предполагается существование системы базисных сеточных финитных функций ортогональных на каждой конкретной сетке, функции которой ортогональны также первым производным всех других функций этого базиса на каждой сетке. Доказана теорема о невозможности использования такой системы функций в алгоритмах проекционных численных методов. Во второй части статьи исследуется возможность построения таких систем базисных финитных функций в общем случае. Проводится модификация сеточных базисных финитных функций с этой целью. Доказывается, что не существует модификаций указанного типа, приводящих к системе ортогональных сеточных базисных финитных функций, каждая из которых ортогональна также первым производным всех других функций этого базиса на каждой сетке.

Статья в формате PDF

645 KB

ортогональные сеточные финитные функции

базисные функции

численные методы

1. Леонтьев В.Л. Ортогональные финитные функции и численные методы / В.Л. Леонтьев. – Ульяновск: Изд-во Ульяновского гос. ун-та, 2003. – 178 c.

2. Леонтьев В.Л. Методы конечных элементов, основанные на использовании обобщенных функций Куранта в теории упругих колебаний // Проблемы динамики и прочности электро- и энерго-машин: тезисы докл. Всерос. научного семинара (С.- Петербург, 18-20 мая 1993 г.). – CПб.: Изд-во Института проблем машиноведения РАН, 1993. – С. 21–22.

3. Леонтьев В.Л. Об одном обобщении функций Куранта // Теория функций и приближений: труды 7-й Саратовской зимней школы (30 января – 4 февраля 1994 г.). Часть 3. – Саратов: Изд-во СГУ, 1995. – C. 36–40.

4. Леонтьев В.Л., Лукашанец Н.Ч. Сеточные базисы ортогональных финитных функций // Журнал вычислительной математики и математической физики. – 1999. – Т. 39, № 7. – C. 1158–1168.

5. Леонтьев В.Л. Об ортогональных финитных функциях и о численных методах, связанных с их применением // Обозрение прикладной и промышленной математики. – 2002. – T.9, Вып. 3. – C. 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) обладают следующими свойствами ортогональности

(1)

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

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

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

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

Пусть также

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

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

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

и

равны в силу (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] самих сеточных функций на каждой конкретной сетке (внутригрупповой взаимной ортогональностью финитных базисных функций)

(4)

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

(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) с учетом такой модификации записываются в виде

Пусть

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

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

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

или

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

поскольку

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

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

Заключение

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

Рецензенты:

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

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

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

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

Библиографическая ссылка