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

Тип работы:
Реферат
Предмет:
Физико-математические науки


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 519. 651: 519. 63
О СВОЙСТВАХ ОРТОГОНАЛЬНЫХ ФИНИТНЫХ ФУНКЦИЙ И ОБ ИХ ИСПОЛЬЗОВАНИИ В АЛГОРИТМАХ ЧИСЛЕННЫХ МЕТОДОВ Леонтьев В. Л.
Ульяновский государственный университет, Ульяновск, e-mail: alex_lion@rambler. ru
Статья посвящена развитию фундаментального направления теории аппроксимации. Это направление относится к построению систем сеточных базисных функций с компактными носителями. Свойство ортогональности таких функций на каждой конкретной сетке исследовано в книге [1]. Некоторые ортогональные финитные функции [1] являются обобщением теории В-сплайнов и дополнением теории ортогональных финитных вейвлетов [2]. Некоторые ортогональные финитные функции [1] в двумерном и трехмерном случаях не имеют аналогов. Здесь рассматриваются возможности формирования новых свойств ортогональных финитных функций и возможности использования функций с такими свойствами в алгоритмах численных методов. В первой части статьи предполагается существование системы базисных сеточных финитных функций ортогональных на каждой конкретной сетке, функции которой ортогональны также первым производным всех других функций этого базиса на каждой сетке. Доказана теорема о невозможности использования такой системы функций в алгоритмах проекционных численных методов. Во второй части статьи исследуется возможность построения таких систем базисных финитных функций в общем случае. Проводится модификация сеточных базисных финитных функций с этой целью. Доказывается, что не существует модификаций указанного типа, приводящих к системе ортогональных сеточных базисных финитных функций, каждая из которых ортогональна также первым производным всех других функций этого базиса на каждой сетке.
Ключевые слова: ортогональные сеточные финитные функции, базисные функции, численные методы
ABOUT PROPERTIES OF ORTHOGONAL FINITE FUNCTIONS AND ABOUT THEIR APPLICATION IN ALGORITHMS OF NUMERICAL METHODS Leontev V.L.
Ulyanovsk State University, Ulyanovsk, e-mail: alex_lion@rambler. ru
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.
Keywords: orthogonal grid finite functions, basis functions, numerical methods
Известны [1] сеточные базисы ортогональных финитных функций (ОФФ) (первый такой базис построен В. Л. Леонтьевым в 1993 году [2], другие подобные базисы созданы и исследованы затем, например, в [3, 4, 5]), состоящие из финитных непрерывных функций, и такие, что скалярные произведения функционального пространства L2 любых двух различных функций каждого такого базиса на каждой конкретной сетке равны нулю. Ранее (в 1988 году) были построены финитные ортогональные вейвлеты И. Добеши [6], но они в отличие от базисных функций [1] не могут быть записаны в аналитической форме, строятся приближенно с заданной точностью с помощью итерационного алгоритма кратноразрешающего анализа, не обладают симметрией (четностью, нечетностью). имеют существенно более сложную структуру и низкую гладкость. Базисные функ-
ции [6] построены с целью их использования в спектральном анализе сигналов, где успешно применяются, но не адаптированы к их использованию в алгоритмах численно-аналитических методов (проекционносеточных (ПСМ) и вариационно-сеточных (ВСМ)). Функции [1] созданы специально для их применения в алгоритмах численно-аналитических методов, связанных со смешанными вариационными принципами (Рейсснера, Ху-Васидзу) или с обобщающими эти вариационные принципы смешанными проекционными условиями. Показано [1], что использование ОФФ устраняет основной недостаток смешанных ВСМ или ПСМ — в несколько раз более высокое число узловых неизвестных по сравнению с ВСМ, основанными на вариационном принципе Лагранжа. Ортогональность таких базисных финитных функций позволяет исключить большую часть узловых неизвестных
величин до начала решения глобальной сеточной системы линейных алгебраических уравнений (СЛАУ) на ЭВМ (число исключаемых величин в 2−4 раза превышает число остающихся в СЛАУ узловых величин), тем самым по числу арифметических операций, необходимых для получения решения СЛАУ, ВСМ, основанные на вариационном принципе Рейсснера, оказываются близкими к ВСМ, основанным на вариационном принципе Лагранжа. При этом качество решений для всех величин (не только для основных неизвестных функций, например, перемещений, но и для их производных, например, деформаций и напряжений) в смешанном ВСМ или ПСМ является существенно более высоким и по точности, и по гладкости, чем в ВСМ, связанном с вариационным принципом Лагранжа.
1. В последние годы в связи с развитием алгоритмов численных методов возник вопрос: возможно ли построение базисных ОФФ, обладающих дополнительными свойствами. В [1] и более ранних работах В. Л. Леонтьева отмечалось, что ряд таких ОФФ имеет неиспользованные для придания функциям свойства ортогональности свободные параметры, которые можно использовать для появления у этих функций дополнительных свойств. Однако, оказывается, что эти дополнительные свойства не могут планироваться произвольно. Здесь показывается, что невозможно создание финитных функций Лагранжева базиса (в таком базисе каждая отдельная функция базиса определяет значение глобальной аппроксимации функции или ее производной только в одном узле), которые, во-первых, обладают свойствами: взаимно-ортогональными являются сами базисные функции, а также каждая функция является ортогональной первым производным всех других функций этого базиса на каждой конкретной сетке, и, во-вторых, предназначены для полноценного использования их в проекционных условиях.
Пусть
7 + і& lt- х1 + 2
& lt- … (
— координаты узлов сетки, расположенных на действительной оси.
Теорема. Пусть {ф (х)} - сеточный Ла-гранжев «базис» в Ь2 (& amp-) такой, что
ф (х) = 5. ,
тогда если функции ф. (х) обладают следующими свойствами ортогональности
I
& lt-рі(х)<-р](х)сІх=Аду-
Ф7 (х)гіх=Вду,
СІХ
(і)
где А, В — некоторые числовые константы, то они не могут быть использованы для аппроксимации неизвестных функций при построении приближенных решений в проекционных условиях, связанных с невязками дифференциальных уравнений, содержащих неизвестные функции и их первые производные.
Доказательство.
Пусть /1(х), /2(х) — непрерывные функции, первая из которых строго возрастает, а вторая строго убывает на интервале (х. -2, х^), причем
/і(х) = Дх). Пусть также
(2)
Н-1
-1
— аппроксимации функций У1(х), /2(х) на интервале (х х) линейными комбинациями функций ф-(х). Здесь, А (х) — строго возрастающая функция- /2 (х) — строго убывающая функция, причем
где Аи), № - постоянные коэффициенты, равные значениям /1(х), /2(х) в узлах х
Скалярные произведения первых производных этих аппроксимаций и ф (х)
|ЇС (х)1 *у+1 ф,(х)^= [ (¦/+1 2 Г^фДхЛ /У
J хі-1 А ^ СІХ — /
Ф](х)& amp- = В?
& lt-0
¦*/+1 і 1
СІХ
Ф ]{х)(Ьс= |
х]- V / х}-
равны в силу (1) и (2). Это означает, что в проекциях первых производных двух раз-
^Ф,(х)
СІХ
Гг
(& lt-)
Ф7 (х)с!х = врг
(г)
личных функций — строго возрастающей и строго убывающей, отличных от нуля
и
и имеющих различные знаки, — информация о первых производных двух различных в окрестности точки х. функций при проецировании аппроксимаций этих первых производных теряется. Таким образом, использование подобных функций в качестве базисных в проекционных условиях методов типа Бубнова-Галеркина применительно к невязкам дифференциальных уравнений невозможно. Следовательно, расщепление системы сеточных уравнений на совокупность несвязанных уравнений с помощью такого ПСМ невозможно.
2. Рассмотрим общий случай. Пусть fx) — материнская локальная (финитная) функция, обладающая на отрезке [-1, 1] свойствами
f-1) = /1) = 0- f0) = A,
где A не равно нулю и равно, например, единице, и обладающая также свойством
fx) + fix — 1) = f0) V х? [0, 1], (3)
необходимым для наличия у последовательности сеточных наборов функций, порожденных этой материнской функцией, полного набора свойств системы базисных функций. За пределами отрезка [-1, 1] функция fx) тождественно равна нулю.
Предположим, что сеточные функции, порожденные материнской функции fx), обладают не только известным свойством ортогональности [1] самих сеточных функций на каждой конкретной сетке (внутригрупповой взаимной ортогональностью финитных базисных функций)
|/Ы/Ь-1)Л=о,
(4)
о
Г df{x) {dxfbс-1)
dx-0.
(5)
мую функцией (-аф (х)) и поэтому свойство (3) выполняется для модифицированных сеточных функций автоматически в случае, если до модификации сеточные функции им обладают. Последнее здесь предполагается, поскольку рассматривается модификация базисных функций с целью изучения вопроса возможности возникновения у них дополнительного свойства, а именно, свойства ортогональности производных сеточных функций, с одной стороны, и самих сеточных функций, с другой стороны. Далее предполагается также, что модифицирующие функции непрерывны и имеют в точках х = 0,± 1 значения, равные нулю.
Условия (4), (5) с учетом такой модификации записываются в виде
/ (/ (х) + «ф (х))(/ (х -1) — аф (х))& lt-&- = 0-
0 ил ил
Пусть
?(х) = fx) + аФ (х)
на [0, 1], тогда последние два условия с учетом (3) записываются в виде 1
|у (х) (/(0)-? = 0-
I
dx
но и дополнительным свойством внутригрупповой ортогональности производных сеточных функций самим сеточным функциям
Выполнение первого из двух последних условий обеспечивается при определенном выборе числовых значений коэффициента а, которые, как показано в [1], существуют. Рассмотрим второе условие, которое переписывается в виде 1
|& lt-%Ы/(0)-0,5ч/2Ы)=0
или
Пусть на отрезке [-1, 1] материнская функция модифицируется, во-первых, некоторой финитной функцией аф (х) (а — постоянный коэффициент), имеющей конечный носитель, не выходящий за пределы отрезка [0, 1], а, во-вторых, финитной функцией (-аф (-х)) с конечным носителем, не выходящим за пределы отрезка [-1, 0]. Первая модифицирующая функция симметрична относительно линии х = 0,5, вторая — относительно линии х = -0,5. В этом случае функция _Дх — 1) на отрезке [0, 1] имеет соответствующую модификацию, определяе-
[& gt-(*)/(0)-0,5чг2Ы]'=0,
откуда получается
?(1)/(0)-0,5|/2(1)-|/(0)/(0) + 0,5|/2(0) = - _/2 (0)+0,5/2 (о) = -0,5/2 (о) ^ 0,
поскольку
?(0) = Х0) = А ф 0- ?(1) = /1) = 0,
а следовательно, второе условие не выполняется при любом выборе модифицирующей функции ф (х) и значений параметра а.
Заключение
Известные примеры численно-аналитических «методов», в рамках которых последние три года строятся и используются «базисные» системы финитных функций, «обладающих» свойствами ортогональности как функций, так и их первых производных, являются неверными в своей основе. Такое заключение основано здесь не только на представленных формальных математических исследованиях, но и на том, что в итоге применения в рамках подобного «метода» процедуры проецирования, связанной с сеткой, каждое уравнение во внутреннем узле сетки оказывается не имеющим никакой связи ни с уравнениями в других узлах, ни с граничными условиями. Возникает естественный вопрос: почему подобное расщепление системы проекционных уравнений возможно для нефинитных функций, например, как известно, для тригонометрических? Это возможно потому, что, несмотря на отделение проекционных уравнений друг от друга в смысле наборов содержащихся в них коэффициентов из состава аппроксимирующих точное решение линейных комбинаций, каждое возникающее после применения процедуры проецирования уравнение, соответствующее отдельной базисной тригонометрической функции, учитывает внешнее воздействие на моделируемую систему на всей области краевой задачи в силу нефинитности базисной функции и при этом любое из таких уравнений может учитывать в общем случае граничные условия.
Список литературы
1. Леонтьев В. Л. Ортогональные финитные функции и численные методы / В. Л. Леонтьев. — Ульяновск: Изд-во Ульяновского гос. ун-та, 2003. — 178 с.
2. Леонтьев В. Л. Методы конечных элементов, основанные на использовании обобщенных функций Куранта в теории упругих колебаний // Проблемы динамики и прочности электро- и энерго-машин: тезисы докл. Всерос. научного семинара (С.- Петербург, 18−20 мая 1993 г.). — СПб.: Изд-во Института проблем машиноведения РАН, 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.
References
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, osnovan-nye na ispolzovanii obobshhennykh funkcijj Kuranta v teorii uprugikh kolebanijj // Problemy dinamiki i prochnosti ehlek-tro- 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 orto-gonalnykh finitnykh funkcijj // Zhurnal vychislitel’nojj matema-tiki 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.
Рецензенты:
Вельмисов П. А., д.ф. -м.н., профессор, заведующий кафедрой «Высшая математика» Ульяновского государственного технического университета, г. Ульяновск-
Журавлев В. М., д.ф. -м.н., профессор, профессор кафедры «Теоретическая физика» Ульяновского государственного университета, г. Ульяновск-
Попов Ф. А., д.т.н., профессор, зам. директора по ИТ, Бийский технологический институт, г. Бийск.
Работа поступила в редакцию 26. 10. 2012.

ПоказатьСвернуть
Заполнить форму текущей работой