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

метки: Корень, Численный, Расчет, Метод, Отрезка, Инженерный, Функция, Формула

Решение алгебраических, нелинейных и трансцендентных уравнений

Пусть задана непрерывная функция f(х) и требуется найти все или некоторые корни уравнения

f(x)=0. (1)

Эта задача распадается на несколько задач. Во-первых, надо исследовать количество, характер и расположение корней. Во-вторых, найти приближенные значения корней. В-третьих, выбрать из них интересующие нас корни и вычислить их с требуемой точностью.

Первая и вторая задачи решаются аналитическими и графическими методами.

Когда ищутся только действительные корни уравнения, то полезно составить таблицу значений f(x).

Если в двух соседних узлах таблицы функция имеет разные знаки, то между этими узлами лежит нечетное число корней уравнения (по меньшей мере, один).

Если эти узлы близки, то, скорее всего, корень между ними только один. Но выявить по таблице корни чётной кратности ^** Кратными корнями называются совпадающие по значению корни. Например, известно, что уравнение y=x² имеет корень x=0, но правильнее утверждать, что это уравнение имеет 2 кратных корня x₁ =0 и x₂ =0. сложно. По таблице можно построить график функции у=f(х) и графически найти точки его пересечения с осью абсцисс. Этот способ более нагляден и дает неплохие приближенные значения корней. Во многих задачах техники такая точность уже достаточна. В технике еще популярны графические методы решения уравнений (номография).

Построение графика позволяет выявить даже корни чётной кратности.

Иногда удается заменить уравнение (1) эквивалентным ему уравнением (х)=(х), в котором функции y ₁ =(х) и y₂ =(х) имеют несложные графики. Например, уравнение хsinх—1=0 удобно преобразовать к виду sinx=l/x. Абсциссы точек пересечения этих графиков будут корнями исходного уравнения.

Приближенные значения корней уточняют различными итерационными методами. Рассмотрим наиболее эффективные из них.

Метод половинного деления (дихотомия)

Пусть мы нашли такие точки a и b что f(a)f(b)0, т. е. на отрезке [a,b] лежит не менее одного корня уравнения. Найдем середину отрезка x _c =(a+b)/2 и вычислим f(x_c ).

Из двух половин отрезка выберем ту, для которой f(x_c )f(a или b)0, т.е. отрезок на котором функция меняет знак. Затем новый отрезок опять делим пополам и выберем ту половину, на концах которой функция имеет разные знаки, и т. д. (рис. 1).

Приближённые методы решения алгебраического уравнения

... с 1 =с 2 =с Найденная точка интересна тем, что она является единственной общей точкой для всех отрезков построенной последовательности Используя непрерывность функции f (x ), докажем, что она является корнем уравнения ... предположениях уравнение (1.1) имеет в точности один корень на интервале (a , b ). 2. Метод дихотомии Этот метод ещё называется методом вилки. Нам необходимо найти корень уравнения ( ...

Если требуется найти корень с точностью , то продолжаем деление пополам до тех пор, пока длина отрезка не станет меньше 2. Тогда середина последнего отрезка даст значение корня с требуемой точностью. Дихотомия проста и очень надежна: к простому корню она сходится для любых непрерывных функций f(x), в том числе недифференцируемых; при этом она устойчива к ошибкам округления. Скорость сходимости невелика: за одну итерацию точность увеличивается примерно вдвое, т. е. уточнение трех цифр требует 10 итераций (т.к. длина отрезка, на котором лежит корень, после 10 итераций равна 1/2 ¹⁰ =1/102410^-3 ).

Зато точность ответа гарантируется.

Перечислим недостатки метода.

1. Для начала расчета надо найти отрезок, на котором функция меняет знак.

2. Если в этом отрезке несколько корней, то заранее неизвестно, к какому из них сойдется процесс (хотя к одному из них сойдется).

3. Метод неприменим к корням четной кратности.

4. Для корней высокой нечетной кратности он сходится, но менее точен и хуже устойчив к ошибкам округления, возникающим при вычислении f(x).

5. Наконец, на системы уравнений дихотомия не обобщается.

Утверждение 1.

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

Чётно кратный корень это корень уравнения вида

(x+a) ^{2 n} =0, где n — целое, n[0,]. (2)

Решением этого уравнения будет корень x=-a кратности 2n. В общем виде уравнение может иметь как чётно, так и нечётно кратные корни. Можно записать общий вид уравнения имеющего (k+m) только действительных корней так:

(x+x ₁ )^{2 n 1} (x+x₂ )^{2 n 2} …(x+x_k )^{2 nk} (x+x_{k +1} )^{2 n ( k +1)+1} (x+x_{k +2} )^{2 n ( k +2)+1} …(x+x_{k + m} )^{2 n ( k + m )+1} =0, (3)

где n1,…,n(k+m) [0,] — целые числа; x ₁ x₂ … x_{k + m} .

В уравнении (3) k чётно кратных и m нечётно кратных корней. Оно раскладывается на (k+m) уравнений, из которых легко получаются корни. Если задать начальный отрезок [-x ₁ -r,-x₁ +r], где r — мало, и проверить условие смены знака функции на его границах, то обнаружим, что знак не меняется в силу чётности степени. А если аналогично проверить нечётно кратные корни, то получим обратную ситуацию.

Следствие 1.

Если корень имеет чётную кратность, то на границах бесконечно малого отрезка с центром в этом корне функция имеет одинаковые знаки.

Технология возделывания картофеля в заданных условиях

... остается еще низкой - 9,8 т/га, при использовании интенсивной технологии около 20 т/га. В России районировано 155 сортов картофеля, различающихся по срокам созревания и хозяйственному назначению. По ... как бы продолжением стебля Корневая система имеет ростковые (глазковые) или первичные корни, образующиеся в начале прорастания клубней. Другие, так называемые пристолонные клубни, появляющиеся в ...

Следствие 2.

Если корень имеет нечётную кратность, то на границах бесконечно малого отрезка с центром в этом корне функция имеет разные знаки.

Пусть на заданном отрезке [a,b] лежит 1 корень чётной кратности, тогда в силу следствия 1 на границах отрезка знак меняться не будет, что означает остановку выполнения итераций и недостижение необходимой точности. Если же на отрезке [a,b] лежит 1 чётно кратный корень и 1 нечётно кратный корень, то чётно кратный корень будет просто игнорирован методом, т.к. условие смены знака являющееся также основным условием, с помощью которого определяется корень на текущем полуотрезке, в силу следствия 1 не выполнится. Следовательно, чётно кратный корень не может быть найден с помощью данного метода.

Утверждение 2., Доказательство.

Первым вариантом постоянства знака функции на границах отрезка является отсутствие корня на нём, поэтому исключим этот случай как тривиальный, будем считать, что на отрезке хотя бы один корень существует. Вторым вариантом — существование чётного количества корней.

Если f(a)f(b)0, то продолжать итерации невозможно, т.к. условие смены знака не подтверждается. Если же, тем не менее, на первом шаге не проверять условие смены знака и разделить отрезок пополам, то может возникнуть ситуация, в которой корни распределяться по чётному количеству в каждой половине отрезка. А чётное количество корней означает чётное количество пересечений оси Ox, даже если существуют кратные корни.

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

Утверждение 3., Доказательство.

В силу утверждения 1 будем рассматривать только корни нечётной кратности. Так как функция меняет знак на концах отрезка, предположим, f(a)?0, f(b)?0. Тогда если f(x _c )?0, то для дальнейшего приближения выберем отрезок [x_c ,b], т.к. f(b)f(x_c )?0. Если же f(x_c )?0, то для дальнейшего приближения выберем отрезок [a,x_c ], т.к. f(a)f(x_c )?0.Для второго случая, когда f(a)?0, f(b)?0 аналогично доказывается существование одного из полуотрезков, на котором функция меняет знак. Из чего следует, что после каждой итерации для одного из полуотрезков условие смены знака обязательно будет выполнено. Следовательно, нет причин для остановки итерационного процесса, который завершится лишь по достижении заданной точности.

Построим блок-схему алгоритма вычисления корня уравнения вида (1) с помощью метода дихотомии. Пусть на начальном отрезке [a,b] функция меняет знак, т.е. на этом отрезке существует нечётное количество нечётно кратных корней. Пример такой функции изображён на рис. 1. Необходимо найти корень x _т с точностью е. Будем считать x_т точным значением корня, x_ч — значение корня полученное данным методом, тогда задача считается выполненной, если x_ч [x_т -е,x_т +е].

1

Дихотомия применяется тогда, когда требуется высокая надежность счета, а скорость сходимости малосущественна., Метод простых итераций

Основной принцип метода заключается в том, что уравнение (1) представляется в виде:

Сила (физическая величина)

... равным нулю. На этом основано измерение величины силы, когда она компенсируется силой, величина которой известна . Второй закон Ньютона позволяет измерять величину силы. Например, знание массы планеты и ... столетий. Эти несоответствия устранил в XVII в. Исаак Ньютон, используя для описания силы математические методы. Механика Ньютона оставалась общепринятой на протяжении почти трехсот лет.[5] ...

x=ц(x), (4)

где ц(x) можно определить многими способами, например, так:

ц(x)=x-бf(x), б=const, или

ц(x)=x+ш(x)f(x),

где ш(x) — произвольная, непрерывная, знакопостоянная функция на отрезке [a,b].

Метод простых итераций в силу (4) определяется следующей рекурсивной формулой:

1

x _n+1 =ц(x_n ), где n=0,1,2,… (5)

Здесь n имеет смысл номера итерации, x ₀ — некоторое начальное приближение. Из (5) видно, что если x_n >x_т , то этот предел и есть корень уравнения (рис. 2).

Пусть в окрестности точки x _т (x_т -Д,x_т +Д), где Д>0 функция ц(x) удовлетворяет условию Липшица:

|ц(x ₂ )-ц(x₁ )|?q|x₂ -x₁ | (6)

для любых x ₂ ,x₁ (x_т -Д,x_т +Д),

0<q<1, (7)

при этом x ₀ (x_т -Д,x_т +Д), (8)

причём, (x _т -Д,x_т +Д) [a,b].

В связи с этим допущением можно сделать несколько утверждений.

Утверждение 1.

Полученные с помощью (5) x _n (x_т -Д,x_т +Д) для любого целого n?0.

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

В силу (8)

|x ₀ -x_т |<Д,

из (5),(6) и (7) получим

|x ₁ -x_т |=|ц(x₀ )-ц(x_т )|?q|x₀ -x_т |,

|x ₂ -x_т |=|ц(x₁ )-ц(x_т )|?q|x₁ -x_т |?q² | x₀ -x_т |,

|x _n -x_т |=|ц(x_n-1 )-ц(x_т )|?q|x_n-1 -x_т |?q² | x_n-2 -x_т |?…? qⁿ | x₀ -x_т |<Д, (9)

Из последнего неравенства |x _n -x_т |<Д следует, что для любого целого n?0 верно утверждение 1.

Утверждение 2.

Последовательность {x _n } сходится при n>? к пределу x_т , являющемуся корнем уравнения.

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

В силу (5) и (6):

|x _{n + m} -x_n |=|ц(x_{n + m -1} )-ц(x_{n -1} )|?q|x_{n + m -1} -x_{n -1} |?q² |x_{n + m -2} -x_{n -2} |?…?qⁿ |x_m -x₀ | (10)

Следовательно, в силу (7) для любого целого

m?0

Признак Коши., Утверждение 3.

На интервале (x _т -Д,x_т +Д) существует только 1 корень уравнения (2).

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

Пусть существует 2 корня

x ₁ и x₂ ?x₁ , тогда в силу (4) x₁ =ц(x₁ ), x₂ =ц(x₂ ).

Из (6) следует |x ₂ -x₁ |=|ц(x₂ )-ц(x₁ )|?q|x₂ -x₁ |, т.е. |x₂ -x₁ |?q|x₂ -x₁ |.

Производительность труда и методы её определения

... определяет уровень производительности труда на данном предприятии. В первой главе курсовой работы показана экономическая сущность и значение производительности труда, представлены показатели и методы измерения производительности труда, а также приведена методика планирования роста производительности труда с ...

Из этого следует, что q=1, а это противоречит (7).

Следовательно, корень только 1.

Следует так же отметить, что если ц(x) имеет на интервале (x _т -Д,x_т +Д) непрерывную производную, то (6) выполняется, когда

|ц?(x)|?q => |ц?(x)|<1 для x(x _т -Д,x_т +Д).

(11)

Если |ц?(x _т )|<1, x₀ (x_т -Д,x_т +Д), то итерации сойдутся.

Если |ц?(x _т )|>1, то в силу непрерывности производной |ц?(x)|>1 и на некотором интервале (x_т -Д₁ ,x_т +Д₁ ), поэтому итерации не сойдутся к корню. Если же |ц?(x_т )|<1, но x₀ (x_т -Д,x_т +Д), т.е. |ц?(x₀ )|>1, то в данном случае о сходимости процесса можно судить только после дополнительного анализа функции ц(x).

Как и при поиске решения методом дихотомии будем считать задачу выполненной, если найденное некоторое значение x _ч [x_т -е,x_т +е], где е — заданная точность. Для определения того, когда можно прекратить итерации, т.е. когда достигнута заданная точность, подробнее рассмотрим неравенство (9).

По сути, нам необходимо добиться выполнения следующего неравенства:

|x _{n +1} -x_т | ? е, (12)

а в силу того, что из (9) можно получить неравенство

|x _{n +1} -x_т | ? q^{n +1} |x₀ -x_т |,

где q можно определить как для любого целого n?1, выведем условие достижения заданной точности (12).

Введём обозначения

д _{n +1} =|x_{n +1} -x_т |. Из (9) ,

а так же очевидно, что д ₀ -д_{n +1} =|x₀ -x_{n +1} |=о, тогда:

д ₀ =д_{n +1} +о

д _{n +1} =gд=g(д_{n +1} +о)

• (13)

Неравенство (13) является условием остановки процесса итераций, т.е. условием достижения заданной точности. В завершение рассмотрения данного метода остаётся только построить блок-схему его алгоритма. Будем считать, что |ц?(x _т )|<1, x₀ (x_т -Д,x_т +Д).

1

1

Метод касательных (Ньютона)

1

Метод Ньютона называют также методом касательных и методом линеаризации. Суть метода заключается в том, что в точке приближения к функции строится касательная (Рис. 3).

Следующая точка приближения — это точка пересечения полученной прямой с осью Ox. Процесс продолжается вплоть до достижения заданной точности.

Из рисунка очень легко получить итерационную формулу метода, используя геометрический смысл производной. Если f(x) имеет непрерывную производную f(x)?0, тогда получим

Аналогично получаем x ₂ , x₃ , и т.д. Таким образом, можем записать общую формулу:

Метод Ньютона можно рассматривать, как частный случай метода простых итераций, если задать

В этом случае

Условие сходимости метода простых итераций можно переписать для метода Ньютона следующим образом:

(14)

Рассуждения по поводу выбора начального приближения в методе Ньютона такие же, как и в методе простых итераций, только вместо (11) используется (14).

Для данного метода также применимо (13).

Напишем блок-схему алгоритма метода Ньютона.

1

1

Метод секущих

1

В данном методе, в отличие от метода Ньютона, проводятся не касательные, а секущие (Рис. 4).

• (15)

В качестве начального приближения необходимо задать не только x ₀ , но и x₁ . Метод секущих имеет одно преимущество перед методом Ньютона — здесь не нужно вычислять производную. Но этот метод имеет также существенные недостатки. Сходимость итераций может быть немонотонной не только вдали от корня, но и в малой окрестности корня.

В знаменателе формулы (15) стоит разность значений функции. Вдали от корня это не существенно; но вблизи корня, особенно корня высокой кратности, значения функции малы и очень близки. Возникает потеря значащих цифр, приводящая к «разболтке» счёта. Это ограничивает точность, с которой можно найти корень; для простых корней это ограничение не велико, а для кратных может быть существенным.

От «разболтки» страхуются так называемым приёмом Гаврика. Выбирают не очень малое е, ведут итерации до выполнения условия |x _{n +1} -x_n |<е и затем продолжают расчёты до тех пор, пока |x_{n +1} -x_n | убывают. Первое же возрастание обычно означает начало «разболтки»; тогда расчёт прекращают и последнюю итерацию не используют. Все ограничения по сходимости итераций для данного метода такие же, как и в методах простых итераций и Ньютона. А вот определение достижения заданной точности, как видно из описания метода, затруднительно, и, даже, возможна ситуация, когда из-за «разболтки» счёта заданная точность не будет достигнута никогда. При использовании метода секущих в явном виде определить точность трудно, поэтому используют косвенный метод. Считают, что вблизи корня |x_{n +1} -x_n ||x_т -x_{n +1} |. Конечно эта оценка весьма примерна, но при больших n (в идеале при n>?) это так и есть. Построим блок-схему алгоритма данного метода с использованием приёма Гаврика.

1

1

Численные методы вычисления определённых интегралов, Метод левых прямоугольников

1

Эти методы численного интегрирования основываются на геометрическом смысле интеграла. Как известно интеграл есть площадь криволинейной трапеции ограниченной подынтегральной функцией f(x) на отрезке [a,b]. Для вычисления интеграла I отрезок [a,b] разбивают на n отрезков длиной h. На каждом отрезке криволинейную трапецию приближают прямоугольником, так как его площадь можно легко вычислить. Затем суммируют все полученные площади, получая тем самым приближённое значение интеграла. На рис. 5 проиллюстрирован, так называемый, метод левых прямоугольников. Его название объясняется тем, что высота прямоугольника f(x) вычисляется в левой границе отрезка h. Выведем формулу для приближённого вычисления интеграла I.

Как видно из рисунка площадь первого слева прямоугольника S ₀ =f(x₀ )h. Площадь следующего S₁ =f(x₁ )h. Легко заметить, что площадь i-го прямоугольника S_i =f(x_i )h. Всего таких прямоугольников n, нумерация их ведётся от 0 до n-1. Таким образом, приближённое значение интеграла, полученное этим методом, вычисляется по формуле:

• (16)

Оценим погрешность формулы (16) на отрезке [x _i ,x_i +h]. Для этого разложим функцию f(x) по формуле Тейлора, выбирая x_i за центр разложения и предполагая наличие у функции требуемых по ходу рассуждений непрерывных производных.

f(x)=f(x _i )+(x-x_i )f(x_i )+… (17)

В формуле (17) для определённости отбросим члены, содержащие производные высших порядков, т.е. 2-ю и выше. Тогда получим:

, т.е.

(18)

Чтобы получить общую погрешность метода на отрезке [a,b] просуммируем все R _{Л i} :

(19)

Поскольку в (17) были отброшены члены, содержащие более высокие степени длины интервала, то выражение остаточного члена (19) является асимптотическим, т.е. выполняющимся при h>0 с точностью до членов более высокого порядка малости. Но для справедливости этой оценки необходимо существование непрерывной f(x); если f(x) кусочно-непрерывна, то удаётся сделать лишь мажорантную оценку

(20)

Эту формулу можно записать иначе

, (21)

где n — количество отрезков, на которые разбит отрезок [a,b].

Перед началом вычислений по данному методу необходимо определиться с h или, что, по сути, то же самое, с n. Пусть известны границы отрезка [a,b], задана подынтегральная функция f(x) и точность е, с которой должен быть вычислен интеграл. Чтобы определить h и n воспользуемся простым и очевидным неравенством:

|R _Л |?е, (22)

откуда получаем для непрерывной f(x):

(23)

либо если учесть, что получим

• (24)

Для кусочно-непрерывной f(x) с учётом (20) и (21) справедливо следующее

(25)

(26)

Метод правых прямоугольников

1

Этот метод похож на предыдущиу. Отличие в том, что высота прямоугольников вычисляется по правой границе (Рис. 6).

Выводы формул для данного метода аналогичны предыдущему. Основные отличия заключаются в нумерации. Формула метода правых прямоугольников выглядит следующим образом:

• (27)

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

Метод средних прямоугольников

1

Чтобы уменьшить погрешность методов левых и правых прямоугольников был предложен метод средних, т.е. метод в котором высота прямоугольника вычисляется в середине отрезка h (Рис. 7).

Обращаясь к рисунку легко увидеть, что площади прямоугольников вычисляются по следующим формулам:

(28)

Оценим погрешность формулы (28) на отрезке [x _i ,x_i +h]. Прежде всего, получим разложение по формуле Тейлора для данного метода. Центр разложения в данном случае будет точка . Тогда получим:

(29)

В формуле (29) отбросим члены, содержащие 3-ю и выше производные. Получаем:

I _i -I_{С р. i} =R_{С р. i}

, т.е.

• (30)

Аналогично тому, как были получены формулы (20), (21), (23), (24), (25) можно получить формулы, по которым оцениваются n и h, для данного метода. Запишем окончательные формулы для непрерывной f(x).

(31)

(32)

Для кусочно-непрерывной f(x) получаем мажорантные оценки:

(33)

(34)

Сравнивая формулы (19) и (30) погрешностей методов приходим к выводу, что погрешность метода средних во много раз ниже погрешности метода левых или правых прямоугольников, т.е. метод средних во много раз точнее и для достижения заданной точности требует меньше машинного времени. Поэтому блок-схему построим именно для метода средних.

Введём обозначение. [n] — целая часть n, полученная путём отбрасывания дробно части. Для вычисления интеграла дано е, [a,b], f(x), f(x).

Прежде всего, через е необходимо получить n. Пользоваться формулой (32) мы не будем, т.к. в ней необходимо вычислять интеграл. Воспользуемся мажорантной оценкой (34).

Для этого перед вычислением интеграла необходимо найти , но это легко осуществимо с помощью небольшого циклического процесса. Для него понадобиться знание количества разбиений , на которых будем вычислять значения второй производной. Обозначим это количество n ₁ .

1

1

Метод трапеций

1

Метод трапеций основан на том, что криволинейная трапеция приближается прямолинейной (Рис. 8).

Т.е. площади вычисляются по следующей формуле:

Таким образом, получаем общую формулу трапеций:

• (35)

Теперь оценим погрешность метода. Вывод формулы погрешности аналогичен выводу (30), поэтому приведём сразу окончательную оценку погрешности для непрерывной f(x):

; (36)

для кусочно-непрерывной f(x):

(37)

Оценивая n и h, получим:

(38)

(39)

Для кусочно-непрерывной f(x):

(40)

(41)

Как видно из полученных формул погрешность метода средних примерно вдвое меньше погрешности метода трапеций. Очевидно, что если функция определена на всём интервале, лучше пользоваться методом средних, метод трапеций используют обычно для функций определённых только в узлах сетки. Знаки главного члена погрешности у формул трапеций и средних разные. Поэтому, если есть расчёты по обеим формулам, то точное значение интервала лежит, как правило, в вилке между ними. Деление этой вилки как 2:1 даёт уточнённый результат близкий к тому, который получается при использовании более точного метода Симпсона.

Метод Симпсона

Этот метод основан на том, что функция f(x) приближается на отрезке [x _i -h, x_i +h] параболой (причём x_i отстоит от x_{i +1} на расстоянии 2h) (Рис.9).

То есть через заданные точки проводится парабола. Но известно, что уравнение параболы имеет вид с(x)=ax² +bx+c, т.е. чтобы определить коэффициенты a, b, c необходимо решить систему из трёх уравнений, а для этого необходимо знать координаты как минимум трёх точек, через которые проходит парабола с(x).

В связи с этим, в отличие от предыдущих методов, для вычисления площади отдельной криволинейной трапеции понадобиться не две, а три точки. Для вывода формулы Симпсона рассмотрим подробнее i-ю криволинейную трапецию ограниченную сверху параболой с_i (x).

1

Для упрощения расчётов сделаем следующие преобразования. Перенесём i-ю криволинейную трапецию в начало координат так, что x _i =0 (следовательно, x_i -h=-h, x_i +h=h).

Очевидно, что её площадь не изменится. Найдём её площадь (Рис. 10).

Для нахождения площади криволинейной трапеции необходимо знать форму кривой ограничивающей её. В нашем случае это парабола . Для нахождения коэффициентов a_i , b_i , c_i составим систему из трёх уравнений. Для простоты обозначим

a _i ?a, b_i ?b, c_i ?c, f(x_i -h)?f₀ , f(x_i )?f₁ , f(x_i +h)?f₁ .

Тогда

Таким образом, получаем

• (42)
• (43)

Таким образом, получаем общую формулу Симпсона:

Для упрощения этой формулы учтём, что

где a и b — границы отрезка интегрирования [a, b]. С учётом этого получим:

где n — количество отрезков длиною 2h. То есть количество отрезков h, на которые разбит отрезок [a, b] должно быть обязательно чётным. Длина отрезка [a, b] равна 2nh.

Для оценки погрешности формулы Симпсона на отрезке [x _i -h, x_i +h] разложим функцию f(x) по формуле Тейлора до пятого члена с центром разложения в точке x_i . То есть

(44)

Тогда с учётом (43) и (44) получим:

Т.е.

, (45)

Сравнив формулы (45), (36), (30) и (19) увидим, что метод Симпсона является наиболее точным из всех описанных методов численного интегрирования.

Учитывая, что для метода Симпсона , запишем формулы оценки n и h для заданного е. Для непрерывной :

Мажорантные оценки для кусочно-непрерывной :

Построим блок-схему.

1

1

1

Приближение функций

Как известно, функцию можно задать многими способами. Самые распространённые из них — это аналитический, графический и табличный. В первых двух способах функция обычно определена на бесконечном числе точек, третий же способ задаёт лишь их конечное количество, например, экспериментальные данные или результаты сложных вычислений. В связи с этим встаёт вопрос получения промежуточных значений таблично заданной функции. Для их нахождения используют два способа: интерполяция (экстраполяция) ^** Экстраполяция отличается от интерполяции тем, что этим методом пытаются вычислить следующее за последним значение или стоящее до первого заданного, а не промежуточное как при интерполяции. и аппроксимация. При использовании интерполяции считают табличные значения точными и интерполяционный многочлен должен точно проходить через заданные точки (Рис. 11).

Аппроксимирующая же функция лишь приближается как можно ближе к заданным точкам, но не обязательно через них проходит (Рис. 13).

Это связано с тем, что табличные данные считаются приближенными, и точное прохождение функции через эти точки будет только увеличивать ошибку.

Интерполяция

Итак, как было сказано выше, задачей интерполяции является поиск такого многочлена, график которого проходит через заданные точки.

Пусть функция y=f(x) задана с помощью таблицы (табл. 1).

Таблица 1

x	x 0	x 1	x 2	…	x n
y	y 0	y 1	y 2	…	y n

Необходимо получить многочлен P _n (x) такой, чтобы выполнялось условие:

P _n (x_i )=y_i . (46)

Для этого зададимся конкретным видом многочлена. Пусть P _n (x) имеет следующий вид:

P _n (x_i )=a₀ +a₁ x+a₂ x² +…+a_n xⁿ . (47)

Для того, чтобы определить коэффициенты a ₀ , a₁ ,… необходимо решить систему из n уравнений с n неизвестными:

(48)

Полином с коэффициентами, полученными путём решения системы (48) называют интерполяционным полиномом Лагранжа и обозначают L _n (x).

Решение системы (48) весьма трудоёмко, поэтому интерполяционный полином Лагранжа представляют в виде линейной комбинации многочленов степени n:

1

• (49)

Необходимо, чтобы каждый многочлен l _i (x) обращался в нуль во всех узлах интерполяции, за исключением i-го, в котором он равен 1 (Рис. 12).

Если l_i (x) удовлетворяет таким условиям, то в i-ом узле интерполяции многочлен L_n (x) примет значение y_i , что удовлетворяет условию поставленной задачи. Таким условиям удовлетворяет многочлен вида:

• (50)

Таким образом, интерполяционный многочлен Лагранжа можно представить следующей общей формулой:

(51)

С использованием формулы (51) построим блок-схему алгоритма метода интерполяции функции полиномом Лагранжа.

1

1

Аппроксимация

1

В задачу аппроксимации входит нахождение такой функции y=f(x), что расстояния между заданными точками y _i и значениями f(x_i ) были минимальными (Рис. 13).

Обозначим отклонение:

е _i =y_i -f(x_i )

В качестве оценки общего отклонения кривой f(x) от табличных данных (Табл. 1) можно было бы взять сумму отклонений е _i , но отклонения могут быть разными по знаку и, не смотря на большие е_i их сумма может быть близка к нулю. Очевидно, что необходимо брать сумму абсолютных значений отклонений, но на практике неудобно пользоваться этой функцией, поэтому в качестве критерия оценки отклонения кривой берут сумму квадратов отклонений:

• (52)

Для определения функции f(x) необходимо, во-первых, задать её общий вид, например, f(x)=ax+b, во-вторых, подставив f(x) в (52) и минимизировав у, найти коэффициенты (a и b).

Такой метод определения коэффициентов для функции f(x) называется методом наименьших квадратов. Наиболее часто встречающиеся виды функции f(x) для метода наименьших квадратов приведены в таблице 2. Формула y=f(x) называется эмпирической формулой или уравнением регрессии y на x.

Таблица 2

Общий вид функции	Аналитическая формула	Вид регрессии
y=f(x,a,b)	y=ax+b	линейная
y=f(x,a,m)	y=ax m	степенная
y=f(x,a,m)	y=ae mx	показательная
y=f(x,a,b)		дробно-линейная
y=f(x,a,b)		логарифмическая
y=f(x,a,b)		гиперболическая
y=f(x,a,b)		дробно-рациональная

Рассмотрим подробнее метод наименьших квадратов на примере линейной регрессии, т.е. общий вид функции такой: f(x)=ax+b. Требуется найти методом наименьших квадратов коэффициенты a и b. Для определения минимального у необходимо приравнять нулю частные производные этой функции по параметрам a и b. Для случая линейной регрессии формула (52) приводится к следующему виду:

• (53)

Условие экстремума запишется так:

• (54)

Отсюда получаем систему из двух уравнений с двумя неизвестными.

(55)

(56)

(57)

Итак, коэффициенты линейной регрессии вычисляются по формулам (57), но для определения параметров других видов регрессии выводить аналогичные формулы необязательно. Можно воспользоваться уже полученными. Для этого будем приводить все виды регрессии к линейному виду.

Степенная регрессия

Прологарифмируем степенную зависимость при условии, что x>0 и a>0.

В общем виде линейную зависимость будем записывать так:

Y=AX+B. (58)

Для линейной регрессии Y=y A=a, X=x, B=b. Для регрессии степенной:

То есть для того, чтобы воспользоваться формулами (57) необходимо вместо подставлять , вместо — . Тогда по первой формуле определяющей значение параметра A получим значение m, а по второй — значение . Следовательно, остаётся небольшое преобразование для получения a: .

Перепишем формулы (57) в общем виде с учётом (58).

• (59)

Показательная регрессия

Сделаем следующее преобразование:

• Тогда (a>0, y>0).

Дробно-линейная регрессия.

(y?0).

Логарифмическая регрессия

(x>0).

Гиперболическая регрессия

(x?0).

Дробно-рациональная регрессия

(x?0, y?0).

Построим блок-схему алгоритма вычисления коэффициентов a и b для линейной регрессии.

1

Список использованной литературы

[Электронный ресурс]//URL: https://inzhpro.ru/kursovaya/chislennyie-metodyi-v-injenernyih-raschetah/

1. Калиткин Н.Н. Численные методы. — М.: Наука, 1978. — 512с.

2. Самарский А.А. Введение в численные методы. — М.: Наука, 1982. — 272с.

3. Хемминг Р.В. Численные методы/ Пер. с англ. М.: Наука, 1968. — 400с.

4. Численные методы анализа Б. П. Демидович, И. А. Марон и Э. З. Шувалова. — М.: Наука, 1967. — 368с.

5. Баранов А.В., Рябчук Э.В. Численные методы в инженерных задачах. — Волгоград.: Политехник, 1988. — 128с.

6. Фильчаков П.В. Справочник по высшей математике. — Киев: Наукова думка, 1975. — 500с.