Главная » Просмотр файлов » Н.Н. Калиткин - Численные методы

Н.Н. Калиткин - Численные методы (1133437), страница 14

Файл №1133437 Н.Н. Калиткин - Численные методы (Н.Н. Калиткин - Численные методы) 14 страницаН.Н. Калиткин - Численные методы (1133437) страница 142019-05-12СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 14)

На каждой итерации вес известен по предыдущей итерации, поэтому коэффициенты а~ь", 6!о легко находятся пз условия минимума квадратичной формы. Практика показывает, что коэффициенты наилучшего приближения слабо зависят от выбора веса, поэтому обычно итерации сходятся быстро.. СРЕДНЕКВАДРАТИЧНОЕ ПРИБЛИЖЕНИЕ В указанных диапазонах изменения аргумента погрешность первой формулы ие превышает 0,44ю а погрешность второй формулы — 2,44де. Таким образом, точность этих аппроксимаций вполне достаточна для многих 'практических приложений. в) Положим у(х)=агс!йх при О=.х (со. Эта функция монотонна, причем у(х) =х прн х-РО и у(+ оо)=п/2. Легко построить дробь <р(х)=х '(1+ — х), удовлетворяющую тем же условиям.

Она дает грубую аппроксимацию арктан~енса; локальная погрешность в точке х=! составляет 304~4. Несложное видоизменение этой формулы ага!Их-х,г'1/ 1+( — х) дает вчетверо лучшую точность. г) Тангенс в первой четверти можно грубо аппроксимировать формулой х ~( — — х), передающей поведение вблизи нуля а наличие полюса прн х=л/2. д) В задачах рассеяний часто встречается одна из специальных функций— интегральная экспонента: Р е-' 1 жч хь Е! (х)= ~ — бт=!п — — — С+ 7 ( — 1)" 4 —. =3 т = эз ьь!' х ь=! Ряд, в который она разлагается, сходится при любых положительных значе.

ниях аргумента. Но только при х - ! сходимость достаточно быстрая, и ряд пригоден для вычисления функции. Если учесть асимптотику Е!(х)=е'х)х при х-ьоо, то рациональную аппроксимацию при хзи! целесообразно искать в следующем виде: л зд)г л Е1(х) — — у ааль 7 Ьчхч, а„=Ь„=1, х е=э е=о (б 1) где не полииомиальная часть асимптотики выделена отдельным множителем. Оказывается, уже л=з, т. е. шесть свободных коэффициентов обеспечивают точность 10 444.

Отметим, что рациональными функциями при небольшом числе коэффициентов можно удовлетворительно аппроксимировать функции с разрывами производной вроде у(х) =(х~, которые плохо поддаются аппроксимации другими способами. Первый ряд абсолютно сходится, но при х)1 сходимость очень медленная; второй ряд сходится асимптотически при больших значениях х. Заменяя первые члены каждого ряда дробями, получим Ф (х) бх при хс!, )' и (3+ха) Ф (х) ~ 1— 2х „4 е при х~!. )' и (1+ 2хз) (гл, И АППРОКСИМАЦИЯ ФУНКЦИЙ 9 3, Равномерное приближение 1. Наилучшие приближения.

Поскольку чебышевская норма сильнее нормы Б, то принято считать, что равномерная аппроксимация лучше аппроксимации в среднем. Поэтому поиску равномерных и особенно нпллучших равномерных приближений, определяемых условием Л(у, гр)=пип, где Л(у, гр)= гпах (у(х) — гр(х)!, (52) а к<э где минимум ищется на множестве функций гр (х), посвящено много работ. В частности, получены следующие результаты (доказательства большинства из них приведены в учебнике И.

С. Березина и Н. П. Жидкова (4]). а) Гели выбрала линейная аппроксимация (37) с чебышевской системой функций фь(х), то равномерное наилучшее приближение единственно"). Дока. зательство существования наилучшего приближения для этого случая было приведено и й 2, и, 1.

б) Чтобы обобщенный многочлеи гр(к) по чебышевской системе функций яра (л), ! -- к - и, был наилучшим равномерным прибчижением к у'(х) на [а, Н, необходимо и достаточно, чтобы иа этом отрезке нашлась не менее и+1 таких тачек, в которых погрешность б(х)=у(х) — гр(х) попеременно принимает значения + Л и — Л (р, гр), Следовательно, погрешность имеет на (а, б) ие менее п нулей, как и у миогочленов наилучшего среднеквадратичного приближения, Впервые этот результат бьш получен П. Л, Чебышевым в 1859 г. для алгебраических мнагочленов.

в) Для функции р(х), имеющей р непрерывных производных„ причем уан(х) удовлетворяет у~повию Липшица с константой (р, д. Джексоном в 1911 г. получены некоторые оценки скорости сходимости наилучших равномерных приближений. При аппроксимации алгебраическим многочленом и-й степени на отрезке — 1 ~ х -= 1; Лл (Сее)ЙЯЯ)р/Ь 2л (Я+1) ппчг) =О(1)лг Я) (53) а при аппроксимации периодической функции с периодом 2п тригонометрическим многочленом такой нсе степени: Л„),(С,Уя)Я =О(1ря (54) где Се — универсальная константа (Са(137), С.

Н. Бернштейн доказал, что пз сходнмости приближений со скоростью О (1)па' Яте), е ) О, следует наличие у функции ограниченной р-1-1-й производной, поэтому оценки Джексона почти иеулучшаемы. Таким образом, эти приближения для достаточно гладких функций быстро сходятся при л-я-со, а для липшиц-непрерывных, ио не гладких функций следует полагать и = О, т. е. для иих приближения сходятся яяедлеяяно. Для произвольной функции, непрерывной на конечном отрезке а ек х«=. б, равномерные приближеиня алгебраическими и тригонометрическими чногочленами также сходятся (теорема, доказанная К. Вейерштрассом в 1885 г„но скорость сходимости, как показал С.

Н. Бернштейн в 1938 г., может быть скгшь угодно малан, 1(яяеиио, как бы медленно ни убывали члены монотонной последовательности б„-я-о, б„=..б„ы"-О, всегда найдется такая непрерывная функция р(х), для которой Л(у, )чя„(х))=б„. Соответствующая оценка джексона для ашебраической аппроксимации произвольной функции, непрерывной РАВНОМЕРНОЕ ПРИНЛИЖЕН1ГЕ прн — 1 ~ х~ 1 (н тем самым равномерно-непрерывной), есть Ьл ~ (')зСо+ 2) ы (2)л), (55) а для тригонометрической аппроксимации непрерывной функции с пернодом 2л: а ~ (,2 Со+ 2) оэ (2л,'л), /1 (56) где ы — модуль непрерывности функцвн.

Наилучшее равномерное прнблнженне рациональной функцией (отношеннеч многочленов) имеет таков нсе порядок точнсстн, как в оценках (55) — (56), где под л надо подразумевать полное число свободных коэффициентов, которое на единицу меньше суммарной степени числителя н знаменателя. г) Многочлены наилучшего равномерного прнблнження не обеспечивают хорошей сходнмостн (а иногда н просто сходнмостн) производных гр' (х) к у' (х).

Если нужна сходнмость производных, то приходится строить другие много- члены, которые имеют меньшую скорость сходвмости Например, многочлены С. Н. Бернштейна в В„(х)= ~~ С (1 — х)" ехеу! — ), 0(а~1, %т А /й1 а (,л)' е=-о (57) равномерно сходятся к любой непрерывной функции у(х), но не быстрее чем О (1)л), сколь бы гладкой функння нн была; зато если существует непрерыв- ная пронэводная у'Я' (х), то производные многочленов С. Н. Берншгейна В'„М (х) равномерно сходятся к ней на указанноч отрезке прн л-ьсо. Сходные оценки сугцествуют н для наилучших аппроксимаций алгебраическими многочленамн на отрезке — 1(х~ 1. Из неравенства (58) следует, что при небольших и погрешность многочлеиов наилучшего среднеквадратичного приближения даже в )! ))с несильна превосходит погрешность многочленов наилучшего равномерного приближения (например, при и ~ 12 не более чем в 7 раз).

Из оценок (53) — (56) следует, что для функций с непрерывными старшими производными, не слишком большими по абсолютной величине, наилучшие равномерные .приближения обеспечивают высокую точность уже при небольших и 5 —:1О. Значит, для таких функций наилучшие среднеквадратичные приближения будут обеспечивать в )! !!с почти ту же точность, что и наилучшие равномерные приближения. Только для недостаточно гладких функций д) Наибольший практический интерес представляет соотношение между точностями, достигаемыми при наилучшей равномерной и наилучшей среднеквадратичной аппроксимациях. Пусть для произвольной функции у (х) с периодом 2и тригонометрический многочлен наилучшего равномерного приближения есть )7„(х). Гхоказана (см.

монографию В. Л. Гончарова 191, стр. 186), что тригонометрический многочлен наилучшего среднеквадратичного приближения той же степени Я„(х) имеет погрешность не более: ))у (х) — Я„(х) ))с ( (4,5+! п и) )!у (х) — )сг„(х)))с. (58) лппгоксимлция Фгнкцип [гл. и среднеквадратичные приближения не сходятся нли плохо сходятся в 1! 1!с, но в этом случае и наилучшие равномерные приближения сходятся настолько медленно, что практически их трудно использовать.

Описанные в 9 2 алгоритмы нахождения наилучших среднеквадратичных приближений намного проще, чем известные алгоритмы нахождения наилучших равномерных приближений, По всем указанным причинам на практике много удобнее искать наилучшие среднеквадратичные, а пе равномерные приближения; как отмечалось в 9 2, для улучшения их сходимости следует явно выделять в простой форме основные особенности функции и ее младших производных и аппроксимировать оставшуюся достаточно гладкую часть. К нахождению равномерных приближений прибегают в основном при разработке алгоритмов для стандартных программ вычисления функций, когда добиваются очень высокой точности при минимальном числе членов суммы. 2. Нахождение равномерного приближения. Для функции, заданной на отрезке [а, 61, не найдено способа определения коэффициентов наилучшего равномерного приближении за конечное число действий.

Рассмотрим простой итерационный процесс нахождения коэффициентов. Чебышевскую норму можно рассматривать как предел ~! !!ь Р при р-«со и единичном весе. В пространстве (р задачу нахождения наилучшего приближения зу — ~рЦ =гп(п удобно решать и итерированием веса: ь ~ ргп (к) [р (х) — ~р<" и (х)11 г(х = пп'и, й рпп (х) = ! У(х) — ф'~ (х) ',г-'; (59) для начала итерационного процесса можно положить рпп (х) = !. Если ~р (х) является обобщенным многочленом, то на каждой итерации задача на минимум опять сводится к решению системы линейных (относительно коэффициентов а„) уравнений. Для решения полной задачи зу — сг!!с=гп!и надо выбрать последовательность р — «оо, для каждого фиксированного р провести итерации (59) до сходимости, а затем в коэффициентах а~~) произвести предельный переход при р-«оо (т. е. оценить, начиная с какого р, коэффициенты перестают меняться в пределах заданной точности при дальнейшем увеличении р).

Двойной предельный переход требует больших численных расчетов. Поэтому целесообразно объединить предельные переходы з-«со и р-«со, Для этого на первой итерации по з положим р = 2, на второй возьмем р = 4, на третьей — р = 6 и т. д.

ЗАДАЧИ Вместо (59) получим следующую задачу: ь ~ реп (х) (у (х) — гр'"и (х) )з г(х =- пни, р"1 (х) = (у (х) — гр(<1 (х))з<, з = О, 1, 2, (ЕО) Здесь начальное условие для итераций рго) (х) ==1 получается естественно при 8=0. Этот итерационный процесс не исследован теоретически и мало опробован н практических расчетах, но поскольку обычно коэффициенты аппроксимации слабо зависят от выбора веса, то следует ожидать быстрой сходимости процесса.

ЗАДАЧИ 1. Доказать, что разделенная разность и-го порядка выражается через узловые значения функции следующим образом: л а у(хз «1, ..., хз)= ~ у(х») П (х» — хй-'. »-о (=а тза» 2. Вывести. оценку (11). а. написать оценки погрешности типа (11) для трех случаев интерполяцн- онного лшогочлена Эрмита 7.й степени: <У< (х; х„хь ..., х,), <т< (х; х<, х„хг, х„«,, «, «э, «,) и <У< («; «<, х<, «э, «<, «т, «ь «ь «т); суавнить нк по(задки точ- ности и численные коэффициенты.

Характеристики

Тип файла
PDF-файл
Размер
21,94 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Зачем заказывать выполнение своего задания, если оно уже было выполнено много много раз? Его можно просто купить или даже скачать бесплатно на СтудИзбе. Найдите нужный учебный материал у нас!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6417
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее