Методы аппроксимации. Методы аппроксимации функций
МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ
ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ
ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«ВОРОНЕЖСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ»
(ФГБОУ ВПО «ВГТУ», ВГТУ)
Факультет радиотехники и электроники
Кафедра высшей математики и физико-математического моделирования
КУРСОВАЯ РАБОТА
по дисциплине: Математика
Тема: «Методы аппроксимации функций»
Разработал студент группы КП-121
И.С. Кононученко
Руководитель Кострюков С.А
ЗАДАНИЕ на курсовую работу
Тема: «Методы аппроксимации функций».
Студент группы КП-121 Кононученко Илья Сергеевич
1. Методы аппроксимации функций.
1.1. Непрерывная аппроксимация.
2. Точечная аппроксимация.
3. Интерполяционный полином Лагранжа.
4. Интерполяционный полином Ньютона.
5. Погрешность глобальной интерполяции.
6. Метод наименьших квадратов.
7. Подбор эмпирических формул.
8. Кусочно-постоянная интерполяция
9. Кусочно-линейная интерполяция.
2. Практическая часть.
2.1. Построить интерполяционный многочлен для функции f(x)=lnx- по узлам х=2; 4; 6; 8; 10; 12. Вычислить приближенное значение логарифма от 5,75. Получить оценку погрешности остаточного члена.
2.2. Функцию f(x), заданную таблицей, аппроксимировать линейной зависимостью ??(х)=Ах2+Вх+С. Найти х, для которого f(x)=10.
1. Методы аппроксимации функций
1.1 Непрерывная аппроксимация
1.2 Точечная аппроксимация
4 Интерполяционный полином Ньютона
8 Кусочно-постоянная интерполяция
9 Кусочно-линейная интерполяция
Практическая часть
2.1 Построить интерполяционный многочлен для функции f(x)=lnx-по узлам х=2; 4; 6; 8; 10; 12. Вычислить приближенное значение логарифма от 5,75. Получить оценку погрешности остаточного члена
2.2 Функцию f(x), заданную таблицей, аппроксимировать линейной зависимостью ?(х)=Ах+В, квадратичной зависимостью ?(х)=Ах2+Вх+С. Найти х, для которого f(x)=10
Список литературы
1.МЕТОДЫ АППРОКСИМАЦИИ ФУНКЦИЙ
1.1Непрерывная аппроксимация
Если исходная функция f(x) задана аналитическим выражением, то при построении аппроксимирующей функции возможно требовать минимальности отклонения одной функции от другой на некотором непрерывном множестве точек, например, на отрезке. Такой вид аппроксимации называется непрерывным или интегральным.
Теоретически для наилучшего приближения целесообразно требовать, чтобы во всех точках некоторого отрезка отклонения аппроксимирующей функции от функции f(x) было по абсолютной величине меньше заданной величины:
В этом случае говорят, что функция равномерно приближает функцию f(x) с точностью e на интервале. Практическое получение равномерного приближения представляет большие трудности, и поэтому этот способ применяется главным образом в теоретических исследованиях.
Наиболее употребительным является так называемое среднеквадратичное приближение, для которого наименьшее значение имеет величина
Потребовав обращения в нуль частных производных от М по параметрам, определяющим функцию, получают уравнения, позволяющие найти наилучшие значения этих параметров.
2 Точечная аппроксимация
Аппроксимация, при которой приближение строится на заданном дискретном множестве точек, называется точечной.
Для получения точечного среднеквадратичного приближения функции y=f(x), заданной таблично, аппроксимирующую функцию строят из условия минимума величины
где yi - значения функции f(x) в точках xi.
Основная сфера применения среднеквадратичного приближения - обработка экспериментальных данных (построение эмпирических формул).
Другим видом точечной аппроксимации является интерполирование, при котором аппроксимирующая функция принимает в заданных точках xi, те же значения yi , что и функция f(x), т.е. .
Рисунок 1
В этом случае, близость интерполирующей функции к заданной функции состоит в том, что их значения совпадают на заданной системе точек.
На рис. 1 показаны качественные графики интерполяционной функции (сплошная линия) и результаты среднеквадратичного приближения (пунктирная линия). Точками отмечены табличные значения функции f(x).
3 Интерполяционный полином Лагранжа
Лагранж предложил строить интерполяционный полином в виде разложения
где li(x) - базисные функции.
Для того, чтобы полином удовлетворял условиям Лагранжа, т.е. был бы интерполяционным, базисные функции li(x) должны обладать следующими свойствами:
) быть полином степени n
) удовлетворять условию
Лагранж показал, что функции, обладающие указанными свойствами, должны иметь следующий вид
С учетом этого выражения интерполяционный полином Лагранжа может быть записан в виде
В отличие от интерполяционного полинома в канонической форме для вычисления значений полинома Лагранжа не требуется предварительно определять коэффициенты полинома путем решения системы уравнений. Однако для каждого значения аргумента x полином Лагранжа приходится пересчитывать вновь, коэффициенты же канонического полинома вычисляются только один раз. Поэтому практическое применение полинома Лагранжа оправдано только в том случае, когда интерполяционная функция вычисляется в сравнительно небольшом количестве точек x.
Интерполяционный полином Лагранжа оказывается очень удобным для приближенного вычисления определенных интегралов. Если, например, некоторую функцию заменить интерполяционным полином Лагранжа, то определенный интеграл от нее может быть вычислен следующим образом
Значения интегралов от не зависят от f(x) и могут быть легко вычислены аналитически.
1.4 Интерполяционный полином Ньютона
Рассмотрим еще одну форму записи интерполяционного полинома
Требования совпадения значений полинома с заданными значения функции в узловых точках Ni(xi)=yi, i=0,1,…,n приводит к системе линейных уравнений с треугольной матрицей для неизвестных коэффициентов:
решить которую не составляет труда.
Интерполяционный полином называется полиномом Ньютона. Интересная особенность полинома Ньютона состоит в том, что каждая частичная сумма его первых (m+1) слагаемых представляет собой интерполяционный полином степени m, построенный по первым (m+1) табличным данным.
5 Погрешность глобальной интерполяции
Ошибка приближения функции f(x) интерполяционным полиномом n-й степени Ln(x) в точке x определяется разностью
Можно показать, что погрешность Rn(x) определяется следующим выражением
Здесь - производная (n+1) порядка функции f(x) в некоторой точке, а функция определена как
Если максимальное значение производной f (n+1)(x) равно
то для погрешности интерполяции следует оценка
Конкретная величина погрешности в точке x зависит, очевидно, от значения функции в этой точке. Качественный характер зависимости показан на рис. 2.
Рисунок 2
Вследствие описанного поведения погрешности, глобальная интерполяция в некоторых случаях может давать совершенно неудовлетворительный результат. Из рисунка видно, что погрешность интерполяции тем выше, чем ближе точка x лежит к концам отрезка. За пределами отрезка интерполяции (т.е. при экстраполяции) быстро растет, поэтому погрешность возрастает существенно.
1.6 Метод наименьших квадратов
Пусть для исходных данных xi, fi, i=1,…,N (нумерацию лучше начинать с единицы), выбран вид эмпирической зависимости: y=?(a0,a1,…,am) с неизвестными коэффициентами a0,a1,…,am . Запишем сумму квадратов отклонений между вычисленными по эмпирической формуле и заданными опытными данными:
S(a0,a1,…,am)=(?(x1,a0,a1,…,am)-fi)2
Параметры a0,a1,…,am будем находить из условия минимума функции S(a0,a1,…,am). В этом состоит метод наименьших квадратов (МНК).
Известно, что в точке минимума все частные производные от S по равны нулю:
Рассмотрим применение МНК для частного случая, широко используемого на практике. В качестве эмпирической функции рассмотрим полином
?(x)=a0+a1x+a2x2+…+amxm
Формула (1) для определения суммы квадратов отклонений примет вид:
S(a0,a1,…,am)=(a0+a1x+a2x2+…+amxm-fi)2 (2)
Вычислим производные
Приравнивая эти выражения к нулю и собирая коэффициенты при неизвестных a0,a1,…,am , получим следующую систему линейных уравнений
Данная система уравнений называется нормальной. Решая эту систему линейных уравнений, получаем коэффициенты.
В случае полинома первого порядка m=1, т.е. , система нормальных уравнений примет вид
При m=2 имеем:
Как правило, выбирают несколько эмпирических зависимостей. По МНК находят коэффициенты этих зависимостей и среди них находят наилучшую по минимальной сумме отклонений.
1.7 Подбор эмпирических формул
При интерполировании функций мы использовали условие равенства значений интерполяционного полинома и данной функции в узлах интерполяции. Если же исходные данные получены в результате опытных измерений, то требование точного совпадения не нужно, так как данные не получены точно. В этих случаях можно требовать лишь приближенного выполнения условий интерполяции. Это условие означает, что интерполирующая функция F(x) проходит не точно через заданные точки, а в некоторой их окрестности, так, например, как это показано на рис.
аппроксимация полином интерполяция формула
Рисунок 3
Тогда говорят о подборе эмпирических формул. Построение эмпирической формулы состоит из двух этапов подбора вида этой формулы, содержащей неизвестные параметры a0,a1,…,am, и определение наилучших в некотором смысле этих параметров. Вид формулы иногда известен из физических соображений (для упругой среды связь между напряжением и деформацией) или выбираются из геометрических соображений: экспериментальные точки наносятся на график и примерно угадывается общий вид зависимости путем сравнения полученной кривой с графиками известных функций. Успех здесь в значительной степени определяется опытом и интуицией исследователя.
Для практики важен случай аппроксимации функции многочленами, т.е. F(x)=a0+a1x+a2x2+…+amxm .
После того, как выбран вид эмпирической зависимости степень близости к эмпирическим данным определяется, используя минимум суммы квадратов отклонений вычисленных и экспериментальных данных.
1.8 Кусочно-постоянная интерполяция
На каждом отрезке интерполяционный многочлен равен константе, а именно левому или правому значению функции.
Для левой кусочно-линейной интерполяции
F(x)= fi-1, если xi-1 ?x Для правой кусочно-линейной интерполяции F(x)= fi-1, если xi-1 Легко понять, что условия интерполяция выполняются. Построенная функция является разрывной, что ограничивает ее применение. Для левой кусочно-линейной интерполяции имеем графическое представление Рисунок 4 1.9 Кусочно-линейная интерполяция На каждом интервале функция является линейной Fi(x)=kix+li. Значения коэффициентов находятся из выполнения условий интерполяции в концах отрезка: Fi(xi-1)=fi-1, Fi(xi-1)=fi . Получаем систему уравнений: kixi-1+ li= fi-1, kixi+ li= fi , откуда находим ki=li= fi- kixi . Следовательно, функцию F(x) можно записать в виде: F(x)= x+ fi- kixi , если, т.е. Или F(x)=ki ·(x-xi-1)+fi-1, ki = (fi - fi-1) / (xi - xi-1), xi-1 ? x ? xi, i=1,2,...,N-1 При использовании линейной интерполяции сначала нужно определить интервал, в который попадает значение x, а затем подставить его в формулу. Итоговая функция будет непрерывной, но производная будет разрывной в каждом узле интерполяции. Погрешность такой интерполяции будет меньше, чем в случае кусочно-постоянной интерполяции. Иллюстрация кусочно-линейной интерполяции приведена на рисунке Рисунок 5 2. ПРАКТИЧЕСКАЯ ЧАСТЬ 2.1 Построим интерполяционный многочлен для функции f(x)=lnx- по узлам х=2; 4; 6; 8; 10; 12. Формула для вычисления данного многочлена выглядит следующим образом: где n- количество узлов. Рассчитаем значения базисных полиномов. Формула для расчета базисных полиномов: Запишем значения узлов функции: Вычислим значения функций f(x) в соответствующих узлах: f(x0)==0.6931471805599453-1.5=-0.8068528194400547(x1)= =1.386294361119891-1.25=0.136294361119891(x2)= =1.791759469228055-1.1666666666666667=0.625092802561388(x3)= =2,079441541679835-1.125=0.954441541679835(x4)= =2.302585092994045-1.1=1.202585092994045(x5)= =2.484906649788-1.083333333333333=1.401573316454667 Рассчитаем значения соответствующих базисных полиномов: Запишем формулу вычисления многочлена f(x)=lnx- по полученным данным: L(x)=f(x0)·l0(x)+ f(x1)·l1(x)+ f(x2)·l2(x)+ f(x3)·l3(x)+ f(x4)·l4(x)+ f(x5)·l5(x). Подставим в формулу полученные значения: L(x)=((- 0.8068528194400547) ·(x-4)(x-6)(x-8)(x-10)(x-12)+ +0.136294361119891·5(x-2)(x-6)(x-8)(x-10)(x-12)- 0.625092802561388·10· · (x-2)(x-4)(x-8)(x-10)(x-12)+ 0.954441541679835·10(x-2)(x-4)(x-6)(x-10)(x-12)-1.202585092994045·5(x-2)(x-4)(x-6)(x-8)(x-12)+ 1.401573316454667· ·(x-2)(x-4)(x-6)(x-8)(x-10)=0,000443792912875·x5-0.001895922201567·x4+ 032520620421826·x3-0.289410042490318·x2+1.50294940468648·x-2.886362165898854 Рисунок 6 L(x)= 0.000443792912875·x5-0.001895922201567·x4+ 032520620421826·x3-0.289410042490318·x2+ 50294940468648·x-2.886362165898854 Из рисунка видно, что графики функций совпадают. Вычислим приближенное значение логарифма от 5,75 с точностью до 0,001. Воспользуемся разложением Пользуясь формулой посчитаем приближенное значение логарифма: Получим оценку погрешности остаточного члена: Формула нахождения остаточного члена в других точках: Rn(x)=f(x)-Ln(x). Подставим значения и вычислим остаточный член: Rn(x)= -0.234721044665224-(-0.149875603361276)= 0.0122 Для абсолютной погрешности интерполяционной формулы Лагранжа можно получить следующую оценку: 0122374?9.9512361 Из оценки следует, что выбирая достаточно большое число точек разбиения можно получить результат с необходимой точностью. Функцию f(x), заданную таблицей аппроксимируем линейной зависимостью ?(х)=Ах+В, квадратичной зависимостью ?(х)=Ах2+Вх+С.
x10151720f(x)371117Решение:
Для решения этой задачи воспользуемся методом наименьших квадратов. Система нормальных уравнений для линейной зависимости (x)=Ax+B: Учитывая, что n=4: ; Решаем систему линейных уравнений: Следовательно, линейная зависимость будет иметь вид: Рассмотрим квадратичную зависимость?(х)=Ах2+Вх+С. Система нормальных уравнений имеет вид: Найдем не подсчитанные суммы: Следовательно, квадратичная зависимость будет иметь вид: Рисунок 7 Функция, заданная таблицей. Линейная зависимость Квадратичная зависимость По графику найдем значение х, для которого f(x)=10. Список литературы 1. Кириллова С.Ю. Вычислительная математика/Кириллова С.Ю. Изд-во Владим. гос. ун-та, 2009. -102с. 2. Справочное пособие по приближенным методам решения задач высшей математики/ Л.И. Бородич, А.И. Герасимович, Н.П. Кеда и др.; под ред. Л.И. Бородич.- М.: Высшая школа, 1986. -189с. 3. Тюканов, А.С. Основы численных методов: учеб. пособие для студентов. Изд-во РГПУ им. А.И. Герцена, 2007. -226с. Нужна помощь по изучению какой-либы темы?
Наши специалисты проконсультируют или окажут репетиторские услуги по интересующей вас тематике. АППРОКСИМАЦИЯ ФУНКЦИЙ Информатика, кибернетика и программирование Как упростить вычисление известной функции fx или же ее характеристик если fx слишком сложная Ответы на эти вопросы даются теорией аппроксимации функций основная задача которой состоит в нахождении функции y=x близкой т. Обоснование способов нахождения удачного вида функциональной зависимости и подбора параметров составляет задачу теории аппроксимации функций. В зависимости от способа подбора параметров получают различные методы аппроксимации; наибольшее распространение среди них получили интерполяция и среднеквадратичное... В окружающем нас мире все взаимосвязано, поэтому одной из наиболее часто встречающихся задач является установление характера зависимости между различными величинами, что позволяет по значению одной величины определить значение другой. Математической моделью зависимости одной величины от другой является понятие функции
y=f(x)
.
В практике расчетов, связанных с обработкой экспериментальных данных, вычислением
f(x),
разработкой вычислительных методов, встречаются следующие две ситуации:
1. Как установить вид функции
y=f(x),
если она неизвестна? Предполагается при этом, что задана таблица ее значений, которая получена либо из экспериментальных измерений, либо из сложных расчетов.
2. Как упростить вычисление известной функции
f(x)
или же ее характеристик, если
f(x)
слишком сложная?
Ответы на эти вопросы даются теорией аппроксимации функций,
основная задача
которой состоит в нахождении функции
y=
(x)
, близкой (т.е. аппроксимирующей) в некотором нормированном пространстве к исходной функции
y=f(x).
Функцию
(x)
при этом выбирают такой, чтобы она была максимально удобной для последующих расчетов.
Основной подход
к решению этой задачи заключается в том, что
(x)
выбирается зависящей от нескольких свободных параметров, т.е. , значения которых подбираются из некоторого условия близости
f(x)
и
(x)
.
Обоснование способов нахождения удачного вида функциональной зависимости и подбора параметров составляет задачу
теории аппроксимации функций
.
В зависимости от способа подбора параметров получают различные
методы аппроксимации;
наибольшее распространение среди них получили
интерполяция
и
среднеквадратичное приближение,
частным случаем которого является
метод наименьших квадратов.
Наиболее простой, хорошо изученной и нашедшей широкое применение в настоящее время является
линейная аппроксимация
,
при которой выбирают функцию, линейно зависящую от параметров, т.е. в виде обобщенного многочлена:
. (3.1)
Здесь известная система линейно независимых (базисных) функций. В качестве в принципе могут быть выбраны любые элементарные функции, например: тригонометрические, экспоненты, логарифмические или комбинации таких функций. Важно, чтобы система базисных функций была
полной
, т.е. обеспечивающей аппроксимацию
f(x)
многочленом (3.1) с заданной точностью при.
Приведем хорошо известные и часто используемые системы. При интерполяции обычно используется система линейно независимых функций. Для среднеквадратичной аппроксимации удобнее в качестве брать ортогональные на интервале [ 1, 1] многочлены Лежандра:
Заметим, что если функция задана на отрезке [
a, b
]
,
то при использовании этой системы необходимо предварительно осуществить преобразование координат, приводящее интервал к интервалу.
Для аппроксимации периодических функций используют ортогональную на [
a, b
] систему тригонометрических функций. В этом случае обобщенный многочлен (3.1) записывается в виде.
Интерполяция
является одним из способов аппроксимации функций. Суть ее состоит в следующем. В области значений
x
a, b
], где функции
f
и
должны быть близки, выбирают упорядоченную систему точек (узлов) (обозначим), число которых равно количеству искомых параметров. Далее параметры подбирают такими, чтобы функция совпадала с
f(x)
в этих узлах, (рис.3.1), для чего решают полученную систему из
n
алгебраических в общем случае нелинейных уравнений.
В случае линейной аппроксимации (3.1) система для нахождения коэффициентов линейна и имеет следующий вид:
(3.2)
Система базисных функций, используемых для интерполяции, должна быть
чебышевской
, т.е. такой, чтобы определитель матрицы системы (3.2) был отличен от нуля и, следовательно, задача интерполяции имела единственное решение.
Для большинства практически важных приложений при интерполяции наиболее удобны обычные алгебраические многочлены, ибо они легко обрабатываются.
Интерполяционным многочленом
называют алгебраический многочлен степени
n-1
, совпадающий с аппроксимируемой функцией в выбранных
n
точках.
Общий вид алгебраического многочлена
. (3.3)
Матрица системы (3.2) в этом случае имеет вид
, (3.4)
и ее определитель (это определитель Вандермонда) отличен от нуля, если точки
x
i
разные. Поэтому задача (3.2) имеет единственное решение, т.е. для заданной системы различных точек существует единственный интерполяционный многочлен.
Погрешность аппроксимации
функции
f(x)
интерполяционным многочленом степени, построенным по
n
точкам, можно оценить, если известна ее производная порядка
n,
по формуле
(3.5)
Из (3.5) следует, что при
h
0 порядок погрешности
p
(см. подразд. 1.4) при интерполяции алгебраическим многочленом равен количеству выбранных узлов
p=n.
Величина
может быть сделана малой как за счет увеличения
n,
так и уменьшения
h.
В практических расчетах используют, как правило, многочлены невысокого порядка
(n
6
)
, в связи с тем что с ростом
n
резко возрастает погрешность вычисления самого многочлена из-за погрешностей округления.
Один и тот же многочлен можно записать по-разному, например. Поэтому в зависимости от решаемых задач применяют различные виды представления интерполяционного многочлена и соответственно способы интерполяции.
Наряду с общим представлением (3.3) наиболее часто в приложениях используют интерполяционные многочлены в форме Лагранжа и Ньютона. Их особенность в том, что не надо находить параметры, так как многочлены в этой форме прямо записаны через значения таблицы.
. (3.6)
Здесь текущая точка, в которой надо вычислить значение многочлена, разделенные разности порядка
k,
которые вычисляются по следующим рекуррентным формулам:
Схема расчета многочлена Ньютона представлена на рис. 3.2.
Вычисления по интерполяционной формуле (3.6) для
n
> 3 используют редко. Обычно при интерполяции по заданной таблице из
m>3
точек применяют квадратичную
n=3
или линейную
n=2
интерполяцию. В этом случае для приближенного вычисления значения функции
f
в точке
x
находят в таблице ближайший к этой точке
i
-узел из общей таблицы, строят интерполяционный многочлен Ньютона первой или второй степени по формулам
(3.7)
и за значение
f(x)
принимают (линейная интерполяция
) или (квадратичная интерполяция
). Схема расчета для линейной и квадратичной интерполяции приведена на рис. 3.3.
(3.8)
Многочлены выбраны так, что во всех узлах, кроме
k
-го, они обращаются в ноль, в
k
-м узле они равны единице:
Поэтому из выражения (3.8) видно, что.
Схема расчета интерполяционного многочлена Лагранжа представлена на рис. 3.4.
Следует отметить, что ввиду громоздкости многочлены Ньютона и Лагранжа уступают по эффективности расчета многочлену общего вида (3.3), если предварительно найдены коэффициенты.
Поэтому когда требуется производить много вычислений многочлена, построенного по одной таблице, оказывается выгодно вначале один раз найти коэффициенты и затем использовать формулу (3.3). Коэффициенты находят прямым решением системы (3.2) c матрицей (3.4), затем вычисляют его значения по экономно программируемой формуле (алгоритм Горнера)
. (3.9)
Схема расчета интерполяционного многочлена общего вида по формуле (3.9) с прямым решением системы (3.2) приведена на рис.3.5.
Находить коэффициенты многочлена (3.3) можно, не решая прямо систему (3.2), а используя разложение коэффициентов Лагранжа (3.8):
.
(3.10)
Рекуррентные формулы для нахождения коэффициентов:
(3.11)
получаются из вида многочленов, если использовать очевидное представление
Схема алгоритма вычисления коэффициентов многочлена общего вида по формулам (3.10), (3.11) представлена на рис. 3.6.
Суть среднеквадратичной аппроксимации заключается в том, что параметры функции подбираются такими, чтобы обеспечить минимум квадрата расстояния между функциями
f(x)
и в пространстве (см. подразд. 1.3), т.е. из условия
. (3.12)
В случае линейной аппроксимации (3.1) задача (3.12) сводится к решению СЛАУ для нахождения необходимых коэффициентов:
(3.13)
Здесь скалярные произведения в
L
2.
Матрица системы (3.13) симметричная, и ее следует решать методом квадратного корня.
Особенно просто эта задача решается, если выбрана
ортогональная система функций
, т.е. такая, что
Тогда матрица СЛАУ (3.13) диагональная и параметры находятся по формуле
В этом случае представление (3.1) называется
обобщенным рядом Фурье
, а называются
коэффициентами Фурье
.
МНК является частным случаем среднеквадратичной аппроксимации. При использовании МНК в области значений
x
, представляющей некоторый интервал [
a, b
]
,
где функции
f
и
должны быть близки, выбирают систему различных точек (узлов)
x
1
, ...,
x
m
, число которых обычно больше, чем количество искомых параметров, Далее, потребовав, чтобы сумма квадратов невязок во всех узлах была минимальна (рис. 3.7):
находим из этого условия параметры. В общем случае эта задача сложная и требует применения численных методов оптимизации. Однако в случае линейной аппроксимации (3.1), составляя условия минимума суммы квадратов невязок во всех точках (в точке минимума все частные производные должны быть равны нулю):
(3.14)
получаем систему
n
линейных уравнений относительно
n
неизвестных следующего вида:
Или. (3.15)
Здесь векторы-таблицы функ-ций. Элементы матрицы
G
и вектора в (3.15) определяются выражениями
скалярные произведения векторов.
Система (3.15) имеет симметричную матрицу
G
и решается методом квадратного корня.
При аппроксимации по МНК обычно применяют такие функции, которые используют особенности решаемой задачи и удобны для последующей обработки.
Схема расчета коэффициентов многочлена вида (3.3) по методу наименьших квадратов представлена на рис. 3.8.
Приведем пример аппроксимации по МНК. Предположим, что известна таблица значений
f(x)
: {
x
1
=
1
,
y
1
=
0.5
,
x
2
=
2
,
y
2
=1.2
,
x
3
=3
,
y
3
=0.8}, т.е.
m
=3. Требуется найти параметры аппроксимирующей функции вида, (n
=2).
Составляем сумму квадратов невязок
Условия минимума (3.1
4
):
Приводя подобные члены, получим окончательно систему двух уравнений с симметричной матрицей относительно неизвестных
c
1
и
c
2
:
Решая ее, находим.
На рис. 3.9 приведена таблица функции
f(x)
и полученная по МНК функция
(x
).
Порядок расчета по МНК следующий. Вначале по исходной таблице формируется матрица
G
и рассчитываются коэффициенты (см. рис. 3.8) (в качестве
j
k
(x
) здесь берется функция
x
k
-
1
). Затем, используя полученные коэффициенты, рассчитывается значение функции в искомой точке (см. рис. 3.5, б).
Во всех вариантах (табл. 3.1.) требуется аппроксимировать заданную исходную функцию
f(x)
многочленом на интервале [
a, b
]
.
Задано количество неизвестных параметров
n
, вид аппроксимации и
m
количество точек, в которых задана функция. Таблица исходной функции
y
i
=
f
(x
i
)
вычисляется в точках
Используя полученную таблицу требуется вычислить значения функций и погрешность в точках, построить графики и проанализировать качество полученной аппроксимации.
Таблица 3
.
1
N вар.
Функция f(x)
Вид аппроксимации
МНК
Ньютона PN
Лагранжа PL
Общего вида POG
Общего вида POL
Линейная PNL
МНК
Квадратичная PNS
Ньютона PN
Лагранжа PL
Общего вида POG
Общего вида POL
МНК
Квадратичная PNS
Линейная PNL
1. Как ставится задача линейной аппроксимации функций?
2. Что такое интерполяция, ее геометрическая интерпретация?
3. Напишите интерполяционный многочлен Ньютона 2-го порядка.
4. Напишите интерполяционный многочлен Лагранжа 2-го порядка.
5. Как осуществляется аппроксимация по методу наименьших квадратов и его геометрическая интерпретация?
6. Выведите СЛАУ относительно коэффициентов функции по МНК для таблицы.
PAGE 30
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны. Размещено на http://www.allbest.ru/
Введение
Аппроксимация (от латинского "approximate" -"приближаться")- приближенное выражение каких-либо математических объектов (например, чисел или функций) через другие более простые, более удобные в пользовании или просто более известные. В научных исследованиях аппроксимация применяется для описания, анализа, обобщения и дальнейшего использования эмпирических результатов. Как известно, между величинами может существовать точная (функциональная) связь, когда одному значению аргумента соответствует одно определенное значение, и менее точная (корреляционная) связь, когда одному конкретному значению аргумента соответствует приближенное значение или некоторое множество значений функции, в той или иной степени близких друг к другу. При ведении научных исследований, обработке результатов наблюдения или эксперимента обычно приходиться сталкиваться со вторым вариантом. При изучении количественных зависимостей различных показателей, значения которых определяются эмпирически, как правило, имеется некоторая их вариабельность. Частично она задается неоднородностью самих изучаемых объектов неживой и, особенно, живой природы, частично обуславливается погрешностью наблюдения и количественной обработке материалов. Последнюю составляющую не всегда удается исключить полностью, можно лишь минимизировать ее тщательным выбором адекватного метода исследования и аккуратностью работы. Поэтому при выполнении любой научно-исследовательской работы возникает проблема выявления подлинного характера зависимости изучаемых показателей, этой или иной степени замаскированных неучтенностью вариабельности значений. Для этого и применяется аппроксимация - приближенное описание корреляционной зависимости переменных подходящим уравнением функциональной зависимости, передающим основную тенденцию зависимости (или ее "тренд"). При выборе аппроксимации следует исходить из конкретной задачи исследования. Обычно, чем более простое уравнение используется для аппроксимации, тем более приблизительно получаемое описание зависимости. аппроксимация алгебраический интерполяционный полином Поэтому важно считывать, насколько существенны и чем обусловлены отклонения конкретных значений от получаемого тренда. При описании зависимости эмпирически определенных значений можно добиться и гораздо большей точности, используя какое-либо более сложное, много параметрическое уравнение. Однако нет никакого смысла стремиться с максимальной точностью передать случайные отклонения величин в конкретных рядах эмпирических данных. Гораздо важнее уловить общую закономерность, которая в данном случае наиболее логично и с приемлемой точностью выражается именно двухпараметрическим уравнением степенной функции. Таким образом, выбирая метод аппроксимации, исследователь всегда идет на компромисс: решает, в какой степени в данном случае целесообразно и уместно "пожертвовать" деталями и, соответственно, насколько обобщенно следует выразить зависимость сопоставляемых переменных. Наряду с выявлением закономерностей, замаскированных случайными отклонениями эмпирических данных от общей закономерности, аппроксимация позволяет также решать много других важных задач: формализовать найденную зависимость; найти неизвестные значения зависимой переменной путем интерполяции или, если это допустимо, экстраполяции. 1. Теоретическое описание задачи
Получить аналитическое описание графически заданных зависимостей концентрации дырок р-типа от температуры в образцах кремния с примесью бора (график 1 и 2) методами Лагранжа, Ньютона, Форсайта и сравнить точности каждого из методов при решении данной задачи. Исходные данные для выполнения курсовой работы:
Рис.1. Зависимость концентрации дырок р от температуры в образцах кремния с примесью бора. 2. Теоретическое описание методов решения
Аппроксимацией
(приближением) функции f(X) называется нахождение такой функции g(X) (аппроксимирующей функции
)
, которая была бы близка заданной. Критерии близости функций f(X) и g(X) могут быть различные. В том случае, когда приближение строится на дискретном наборе точек, аппроксимацию называют точечной
или дискретной.
В том случае, когда аппроксимация проводится на непрерывном множестве точек (отрезке), аппроксимация называется непрерывной
или интегральной
. Примером такой аппроксимации может служить разложение функции в ряд Тейлора, то есть замена некоторой функции степенным многочленом. Наиболее часто встречающим видом точечной аппроксимации является интерполяция
(в широком смысле). Пусть задан дискретный набор точек X i (i=0,1,…,n), называемых узлами интерполяции
, причем среди этих точек нет совпадающих, а также значения функции Y(X) в этих точках. Требуется построить функцию g(X),
проходящую через все заданные узлы. Таким образом, критерием близости функции является g(X i)=y i . В качестве функции g(X) обычно выбирается полином, который называют интерполяционным полиномом
. В том случае, когда полином един для всей области интерполяции, говорят, что интерполяция глобальная
.
В тех случаях, когда между различными узлами полиномы различны, говорят о кусочной
или локальной интерполяции
.
Найдя интерполяционный полином, мы можем вычислить значения функции f(X) между узлами (провестиинтерполяцию в узком смысле слова
), а также определить значение функции f(X) даже за пределами заданного интервала (провести экстраполяцию
). Следует иметь в виду, что точность экстраполяции обычно очень невелика особенно при удалении от заданного интервала. Аппроксимация функций с помощью алгебраических интерполяционных полиномов.
Задача аппроксимации функции с помощью алгебраического интерполяционного полинома формулируется следующим образом. Пусть аналитическое выражение функции Y=f(X) неизвестно, заданы только ее значения Y 1 ...Y N в точках X 1 ...X N некоторого отрезка . Необходимо найти полином степени n для которого выполняются условия: Так как в точках X j значение функции Y j и значение полинома P n (X j) должны совпадать между собой, то неизвестные коэффициенты полинома можно найти путем решения системы уравнений (1.2) Интерполяционная формула Лагранжа.
Одну из простейших формул интерполяции позволяет построить метод Лагранжа. По условию находим полином P N -1 (X) степени (N-1), который в N точках совпадает с N значениями функции f(X). Если найти систему полиномов {ц j (X)}, каждый из которых в точке Xj равен единице, а в остальных точках равен нулю, то интерполяционный полином можно представить в виде Это следует из того, что Последовательность функций {ц j (X)} такого типа называется фундаментальной системой полиномов. По предположению полином ц j (X) в точках X k при k?j обращается в нуль, поэтому его можно представить в виде где С J - некоторая постоянная. Учитывая, что ц j (X j)=1, получим Отсюда следует, что интерполяционный полином Лагранжа
имеет вид Формула Лагранжа при N?4 становится громоздкой при практическом использовании, так как в нее входит произведение П(х). Рассмотрим случай выбора узлов интерполяции, когда формула значительно упрощается. Предположим, что функция Y=f(X) задана на отрезке [-1,1]. Далее полученные результаты обобщим на случай произвольного отрезка . Сначала введем полиномы Чебышева 1-го и 2-го рода. По определению полиномы Чебышева 1-го и 2-го рода задаются с помощью формул: Если узлами интерполяции являются нули полинома T N (X), то интерполяционный полином Лагранжа можно представить в виде Предположим, что узлами интерполяции являются нули полинома Чебышева 2-го рода U N (X). В этом случае интерполяционную формулу Лагранжа можно представить в виде Интерполяционная формула Ньютона.
На практике для аппроксимации функции часто используется интерполяционный полином Ньютона. Этот полином вводится с помощью разделенных разностей различных порядков, найденных по значениям функцииY1,…,Y N в точках X 1 ,…,X N . По определению разделенные разности первого порядка равны Разности 2-го порядка определяются с помощью разностей 1-го порядка: Разделенные разности n-го порядка можно представить в виде Отсюда следует, что разделенная разность является симметричной функцией относительно узлов X j , т.е. не зависят от порядка расположения входящих в нее переменных X j . Построим интерполяционный полином Ньютона. Пусть X-произвольная точка отрезка . Рассмотрим разность первого порядка Из этого выражения находим значение функции в точке: Разность второго порядка имеет вид Подставив это выражение в (1.15) получим Разность 3-го порядка позволяет представить (1.19) в виде Продолжая процесс подстановки, получим выражение которое можно представить в следующей форме Полином P N -1 (X) является интерполяционным, так как имеют место равенства Этот полином называется интерполяционным полиномом Ньютона
, а R N -1 (X) - остаточным членом формулы Ньютона. Так как по значениям функции в некоторых точках можно единственный интерполяционный полином, то полином Ньютона путем перегруппировки его членов можно преобразовать в интерполяционный полином Лагранжа, для которого каждое из слагаемых суммы (1.18) зависит от всех узлов интерполяции, произвольный m-й член полинома Ньютона зависит только от m первых узлов. Поэтому для полинома Ньютона добавление новых узлов интерполяции приводит лишь к появлению новых слагаемых полинома, без изменения первоначальных. На практике часто используется интерполяционный полином Ньютона, представленный в виде который называется формулой Ньютона интерполирования назад. Метод наименьших квадратов с помощью ортогональных полиномов,
интерполяционная формула Форсайта
Рассмотренные выше методы позволяют аппроксимировать функции, заданные экспериментальными данными, с помощью интерполяционных многочленов. На практике интерполяционные формулы применяются в тех случаях, когда ошибки в данных можно не учитывать и число N точек X j является малым. При больших N эти формулы становятся громоздкими, а также возникают трудности, связанные с неустойчивостью интерполяционного процесса на концах отрезка . В реальных задачах ошибки в экспериментальных данных необходимо учитывать. Если зависимости между параметрами являются достаточно гладкими, то даже при больших N часто нет необходимости выбирать для аппроксимации функций можно использовать метод наименьших квадратов (VYR)/ Предположим, что функция Y=f(X) задана на отрезке экспериментальными значениями где е j - некоррелируемые случайные величины, имеющие нулевое математическое ожидание и дисперсию у 2 . При аппроксимации функции Y=f(X) алгебраическим полиномом (1.28) с помощью МНК по экспериментальным данным необходимо оценить коэффициенты a k полинома таким образом, чтобы сумма квадратов была минимальной. Алгебраический полином (1.28) является частным случаем общей линейной модели Оценка коэффициентов общей линейной модели сводится к решению системы нормальных уравнений X T Xa = X T Y, которая для приближения (1.28) имеет вид Оценки коэффициентов. Решение системы (1.31) существует, если определитель Дn+1 ? 0. В рамках современной математики задача аппроксимации с помощью полинома (1.28) формулируется как задача оценки коэффициентов модели, представляющей собой комбинацию функций некоторой подсистемы Lm системы базисных функций L={1,X,X 2 ,…,X n ,…}. Влияние плохой обусловленности значительно уменьшается, если вместо модели (1.28) рассматривается модель представляющая собой линейную комбинацию элементов подсистемы Lm системы ортогональных полиномов L={1,ц 1 (X), ц 2 (X),…, ц n (X),…}. Система полиномов ортогональна на в следующем смысле В настоящее время разработано несколько подходов к построению ортогональных полиномов. Одной из наиболее простых является система полиномов Чебышева. На практике удобной является так же система ортогональных полиномов Форсайта
где и и выбираются из условия ортогональности. Если аппроксимирующая функция имеет вид (1.34), то Эту систему уравнений можно представить в матричной форме с - (n+1)-мерный вектор-столбец неизвестных параметров модели (1.34). Из условия ортогональности матрица системы нормальных уравнений является диагональной: Из линейной алгебры известно, что матрица, обратная к диагональной, также является диагональной., причем ее элементы равны обратным величинам диагональных элементов исходной матрицы. Поэтому учитывая, что решение нормальной системы уравнений можно найти по формуле получим оценки коэффициентов модели (1.34) Оценки коэффициентов не коррелированны между собой и имеют дисперсии где - дисперсия случайных ошибок эксперимента. При решении практических задач степень аппроксимирующего полинома обычно неизвестна. Если функция Y=f(X) аппроксимируется с помощью полинома (1.28), то выбор его степени часто осуществляется следующим образом. Начиная с некоторого малого числа n 0 , выбирается возрастающая последовательность целых чисел n 1 ,n 2 ,n 3 ,…,n p ,… и для этих степеней путем решения системы (X T X)a=X T y вычисляются коэффициенты полинома. Для каждого значения n с помощью найденных оценок вычисляются остаточные дисперсии При увеличении n остаточная дисперсия сначала обычно убывает, а позже наступает момент, когда она начинает возрастать. Поэтому степень полинома n выбирается равной значению n, при котором остаточная дисперсия является минимальной. Система нормальных уравнений взвешенного метода наименьших квадратов имеет вид (1.29), где Задачу аппроксимации функции взвешенным методом наименьших квадратов можно также решить, используя ортогональные полиномы. При этом коэффициенты модели (1.34) выбираются из условия минимума функции Если полиномы образуют ортогональную систему функций, то матрица нормальной системы является диагональной, а решение системы имеет вид Ортогональные полиномы можно найти с помощью рекуррентных формул 3. Расчёт и графики по каждому методу аппроксимации
Результат работы программы расчета по методу Лагранжа:
Для кривой 1: Для кривой 2: Рис.2. Зависимость концентрации дырок р от температуры в образцах кремния с примесью бора, полученная по результатам метода Лагранжа. Результат работы программы расчета по методу Ньютона:
Для кривой 1:Для кривой 2: Рис.3. Зависимость концентрации дырок р от температуры в образцах кремния с примесью бора, полученная по результатам метода Ньютона. Результат работы программы расчета по методу Форсайта:
Степень алгебраического полинома М=2
Для кривой 1: Для кривой 2: Рис.4. Зависимость концентрации дырок р от температуры в образцах кремния с примесью бора, полученная по результатам метода Форсайта. Степень алгебраического полинома М=3
Для кривой 1: Для кривой 2: Рис.5. Зависимость концентрации дырок р от температуры в образцах кремния с примесью бора, полученная по результатам метода Форсайта. Степень алгебраического полинома М=4
Для кривой 1: Для кривой 2: Рис.6. Зависимость концентрации дырок р от температуры в образцах кремния с примесью бора, полученная по результатам метода Форсайта. 4. Сводный график
Рис.7. Сводный график. 1 - заданный график для «1» 2 - заданный график для «2» 3 - аппроксимированный график по методу Лагранжа для «1» 4 - аппроксимированный график по методу Лагранжа для «2» 5 - аппроксимированный график по методу Ньютона для «1» 6 - аппроксимированный график по методу Ньютона для «2» 7 - аппроксимированный график по методу Форсайта со степенью полинома 2 для «1» 8 - аппроксимированный график по методу Форсайта со степенью полинома 2 для «2» 9 - аппроксимированный график по методу Форсайта со степенью полинома 3 для «1» 10 - аппроксимированный график по методу Форсайта со степенью полинома 3 для «2» 11 - аппроксимированный график по методу Форсайта со степенью полинома 4 для «1» 12 - аппроксимированный график по методу Форсайта со степенью полинома 4 для «2» 5. Анализ точности
Аппроксимация по методу Лагранжа
Погрешность составляет 0,13 % Погрешность составляет 0,05 % Аппроксимация по методу Ньютона
Погрешность составляет 0,15 % Погрешность составляет 0,06 % Аппроксимация по методу Форсайта со степенью 2
Погрешность составляет 0,71% Погрешность составляет 0,62 % Аппроксимация по методу Форсайта со степенью 3
Погрешность составляет 1,21% Погрешность составляет 0,5 % Аппроксимация по методу Форсайта со степенью 4
Погрешность составляет 0,12 % Погрешность составляет 0,05% Заключение
В своей практике инженер часто сталкивается с необходимостью аналитически описать экспериментально полученные зависимости, представленные графически или таблично. Для этого используются методы аппроксимации, соответствующее алгоритмическое и программное обеспечение. В данной курсовой работе мы на практике знакомились с различными методами аппроксимации, соответствующим алгоритмическим и программным обеспечением. Проведя все расчеты можно сделать вывод, что из рассмотренных методов аппроксимации, метод наименьших квадратов со степенью 4 является самым точным. Список литературы
1. Бронштейн и Семендяев Справочник по математике для ВТУЗов. - М., 1986г. 2. Шафрин Ю.А. Информационные технологии.- М., Лаборатория базовых знаний, 2000г. 3. Конева Н.Е. Информационные технологии в электронике. Методические указания к лабораторному практикуму для студентов специальности 210105 - Электронные приборы и устройства. МГОУ, 2009г. 4. Норенков И.П. Системы автоматизированного проектирования. Учебное пособие, Высшая школа, Москва, 1986г. 5. Вллах И., Сингхал К. Машинные методы анализа и проектирования электронных схем. - М., Радио и связь, 1988г. 6. Чау Л.О., Пен-Ман Лиин. Машинный анализ электронных схем. - М.,Мир. 7. Шур Т. Решение инженерных задач на ЭВМ, - М., 1982г. 8. Чахмахсазян Е.А., Бармаков Ю.Н., Голденберг А.Э. Машинный анализ интегральных схем. - М., Советское радио, 1974г. 9. Бахвалов Н.С., Лапин А.Р., Чижонков Е.В. Численные методы в задачах и упражнениях. - М., Высшая школа, 2000г. 10. Бахвалов Н.С., Кобельков Г.М., Поспелов В.В. Сборник задач по методам вычислений. - М.,Издательство МГУ,1989г. 11. Конева Н.Е. Информационные технологии в электронике. Методические указания к курсовой работе для студентов специальности 210105 - Электронные приборы и устройства. МГОУ, 2009г. Размещено на Allbest.ru Проведение аппроксимации данных с помощью Excel, расчет площадей (отдельно для выпуклой и вогнутой кривых периферического, серединного и корневого сечения) и целевой функции V с целью нахождения полного объема бетонной строительной конструкции. контрольная работа , добавлен 26.01.2010 Определение уравнений динамики и передаточных функций элементов системы автоматического управления. Дискретизация последовательного корректирующего звена методом аппроксимации операции интегрирования. Анализ устойчивости автоматической системы управления. курсовая работа , добавлен 27.02.2014 Балансировка ротора машин и балансировка гибких роторов как задача оценивания дисбалансов. Условие допустимости одной статической балансировки. Оценивание методом наименьших квадратов. Целевая функция метода наименьших квадратов и численные эксперименты. дипломная работа , добавлен 18.07.2011 Материальный баланс и расход абсорбента. Определение коэффициента диффузии ацетона в воде. Поверхность массопередачи, формула для её расчета. Определение геометрических параметров абсорбера с помощью уравнения массопередач и через высоту единиц переноса. курсовая работа , добавлен 05.11.2012 Исследование влияния скорости печати на качество оттисков по совмещению красок при многокрасочной флексографской печати. Математическое моделирование как приближённое описание реальных объектов с помощью математических выражений, его главные этапы. контрольная работа , добавлен 14.04.2011 Нагрев металла перед прокаткой. Автоматизация процесса нагрева металла. Выбор системы регулирования давления. Первичный измерительный преобразователь перепада давления. Метод наименьших квадратов. Измерение и регистрация активного сопротивления. курсовая работа , добавлен 25.06.2013 Особенности статической настройки, использование пробных заготовок с помощью рабочего калибра. Настройка по пробным заготовкам с помощью универсального измерительного инструмента. Ее проведение с учетом переменных систематических погрешностей и без них. презентация , добавлен 26.10.2013 Разработка и компоновочные схемы токарных многоцелевых станков. Привод главного движения. Обработка фасонной поверхности с помощью копира. Управление фрикционными муфтами с помощью кулачка. Регулирование подачи с помощью конуса Нортона и гидропривода. реферат , добавлен 02.07.2015 Получение тонкопленочных покрытий в вакууме, термическое и магнетронное испарение. Конструирование жидкофазного магнетрона с помощью AutoCAD. Методы исследования параметров тонких пленок. Измерение толщины тонкопленочных покрытий с помощью профилометра. дипломная работа , добавлен 15.06.2012 Усовершенствование шлифовальной операции технологического процесса обработки хвостовой части метчика с помощью методов технического творчества. Совершенствование шлифования цилиндрической поверхности с помощью мозгового штурма и метода проб и ошибок. Математической
моделью зависимости одной величины от
другой является понятие функции y=f(x)
.
Аппроксимацией
называется
получение некой функции, приближенно
описывающей какую-то функциональную
зависимость f(x),
заданную
таблицей значений, либо заданную в виде,
неудобном для вычислений. При этом эту
функцию выбирают такой, чтобы она была
максимально удобной для последующих
расчетов. Основной
подход
к
решению этой задачи заключается в том,
что функция fi(x)
выбирается
зависящей от нескольких свободных
параметров c1,
c2, …, cn,
значения
которых подбираются из некоторого
условия близости f(x)
и fi(x)
.
Обоснование способов нахождения удачного
вида функциональной зависимости и под-
бора параметров составляет задачу
теории
аппроксимации функций
.
В зависимости от способа подбора
параметров получают различные методы
аппроксимации
,
среди которых наибольшее распространение
получили интерполяция
и
среднеквадратичное
приближение
.
Наиболее
простой является линейная
аппроксимация
,
при которой
выбирают функцию линейно зависящую от
параметров, т. е. в виде обобщенного
многочлена:
.
Интерполяционным
многочленом
называют
алгебраический многочлен степени n-1
,
совпадающий с аппроксимируемой функцией
в n
выбранных
точках. Погрешность
аппроксимации
функции
f(x)
интерполяционным
многочленом степени n-1
,
построенным по n
точкам, можно
оценить, если известна ее производная
порядка n.
Суть
среднеквадратичной
аппроксимации
заключается в том, что параметры функции
подбираются такими, чтобы обеспечить
минимум квадрата расстояния между
функциямиf(x)
и
fi
(x
,
c
).
Метод
наименьших квадратов
является частным случаем среднеквадратичной
аппроксимации. При использовании метода
наименьших квадратов аналогично задаче
интерполяции в области значений x
,
представляющей некоторый интервал [a,
b
],
где функции
f(x)
и
fi(x)
должны быть
близки, выбирают систему различных
точек (узлов) x1,
..., x
m, число
которых больше, чем количество искомых
параметров. Далее, требуют чтобы сумма
квадратов невязок во всех узлах была
минимальна. Следует
отметить, что ввиду громоздкости
многочлены Ньютона и Лагранжа уступают
по эффективности расчета многочлену
общего вида. Поэтому, когда требуется
производить многократные вычисления
многочлена, построенного по одной
таблице, оказывается выгодно вначале
один раз найти коэффициенты с. Коэффициенты
находят прямым решением системы с, затем
вычисляют его значения по алгоритму
Горнера. Недостатком такого вида
аппроксимации является необходимость
решения системы линейных алгебраических
уравнений. Лагранжем
была предложена своя форма записи общего
интерполяционного алгебраического
многочлена в виде, не требующем решения
системы линейных алгебраических
уравнений. Следует отметить, что ввиду
громоздкости многочлены Ньютона и
Лагранжа уступают по эффективности
расчета многочлену общего вида. Ньютоном
была предложена форма записи общего
интерполяционного алгебраического
многочлена в виде, не требующем решения
системы линейных алгебраических
уравнений. Следует отметить, что ввиду
громоздкости многочлены Ньютона и
Лагранжа уступают по эффективности
расчета многочлену общего вида.Репетиторство
Отправь заявку
с указанием темы прямо сейчас, чтобы узнать о возможности получения консультации.ТЕМА 3. АППРОКСИМАЦИЯ ФУНКЦИЙ
3.1. Зачем нужна аппроксимация функций?
3.2. Что такое интерполяция?
3.3. Какие бывают многочлены и способы интерполяции?
Интерполяционный многочлен Ньютона (PN)
Линейная (PNL) и квадратичная(PNS) интерполяция
Интерполяционный многочлен Лагранжа (PL)
Интерполяция общего вида, использующая прямое решение системы (3.2) методом Гаусса (POG)
Интерполяция общего вида, использующая расчет коэффициентов многочлена (3.3) через многочлен Лагранжа (POL)
3.4. Что такое среднеквадратичная аппроксимация?
Метод наименьших квадратов (МНК)
3.5. Варианты заданий
3.6. Контрольные вопросы
А также другие работы, которые могут Вас заинтересовать
33218.
Классификация режимов работы судового электрооборудования в зависимости от продолжительности рабочего цикла
83.58 KB
Судовое электрооборудование будет работать надежно, если оно не только правильно сконструировано, но и правильно используется. СЭО используется правильно, если оно соответствует условиям работы судового механизма, устройства и т.п.
33224.
Электромагнитные измерительные приборы
13.15 KB
Магнитоэлектрический прибор измерительный прибор непосредственной оценки для измерения силы электрического тока напряжения или количества электричества в цепях постоянного тока Электри́ческая мо́щность физическая величина характеризующая скорость передачи или преобразования электрической энергии. Ток изменяющийся во времени по значению и направлению называется переменным. В практике применяют периодически изменяющийся по синусоидальному закону переменный ток.
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Подобные документы
Интерполяция общего вида
Интерполяционный многочлен Лагранжа
Интерполяционный многочлен Ньютона