Исчисления методами Лагранжа Рунге Кутта Ньютона и Гаусса
Рефераты >> Математика >> Исчисления методами Лагранжа Рунге Кутта Ньютона и Гаусса

2.2.Решение

Решение будем производить методом Лагранжа.Oцениваем шаг

h=xi+1 -xi

В этой таблице h=const.Для интерполяции функции с произвольно задаными узлами выбираем интерполяционный многочлен Лагранжа:

;

Выражения,называемые коэффициентами Лагранжа:

Далее построим матрицу Лагранжа:

Обозначим произведение строк через ,а произведение элементов главной диагонали через ,тогда :

Вычислим её:

отсюда:

Пn+1=4,00384 .10-9

D0=7,68488.10-6 D5=1.1475.10-8

D1=-1.84275.10-7 D6= -1.16944.10-8

D2= 4.2525.10-8 D7=2.3625.10-8

D3=2.92313 10-9 D8= -8.91.10-8

D4= -7.0875.10-9 D9=7.86713.10-7

Далее по формуле:

,

имеем

В результате проделанной работы мы произвели интерполяцию функции заданной таблицей 2.1 и получили значение функции в точке х=0,38 y=0,683860.

О справедливости полученного результата мы можем судить из того ,что точка х=0,38 находиться точками х=0,30 и х=0,40 и искомое значение должно находиться между соответствующими значениями этих точек. Полученное значение y=0,683860 находиться в пределах между y(0.30)=0.670320 и y(0.40)=0.740818.

Следовательно решение верно.

3.Задача 3

3.1.Постановка задачи

Решить систему линейных уравнений:

9.3x1+(1.62+a)x2+6.1x3+1.9x4=-12.65+b;

4.92x1+7.45x2+(9.7-a)x3+2.46x4=10.21;

4.77x1+(6.21+a)x2+9.04x3+2.28x4=13.45;

3.21x1+(2.65-a)x2+3.69x3+6.99x4=-10.35.

методом Гаусса. Все расчёты ведите с тремя значащими цифрами после запятой.

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

3.2.Решение

Перепишем систему линейных уравнений в виде:

9.3x1+(1.62+0.8)x2+6.1x3+1.9x4=-12.65+3.6;

4.92x1+7.45x2+(9.7-0.8)x3+2.46x4=10.21;

4.77x1+(6.21+0.8)x2+9.04x3+2.28x4=13.45;

3.21x1+(2.65-0.8)x2+3.69x3+6.99x4=-10.35.

9.3x1+2.42x2+6.1x3+1.9x4=-9.05;

4.92x1+7.45x2+8.9x3+2.46x4=10.21;

4.77x1+7.01x2+9.04x3+2.28x4=13.45;

3.21x1+1.85x2+3.69x3+6.99x4=-10.35.

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

а15,а25,а35,а45---свободные члены

---суммирующий (контрольный) коэффициент

Прямой ход. Заполнение таблицы:

1.Запишем аij в четырёх строках и пяти столбцах раздела 1 таблицы(i=1,2,3,4,j=1,2,3,4,5)

2.Стимулирующие аi6 запишем в столбце å (столбец контроля)

3.Вычисляем b1j=a1j/a11 (j=1,2,3,….6) и запишем в пятой строке раздела 1

4.Вычисляем и проверяем совпала ли она с b16 c вычисления ведутся с постоянным количеством знаков после запятой). В противном случае проверяем действия пункта 3.

5.Вычисляем b1ij(1)=aij-ai1.b1j(i=2,3,4 , j=2,3,….6) и записываем их в в первые три строки раздела 2.

6.Проверка. Сумма элементов каждой строки и должен совпасть с указанной в п.4 точностью, иначе надо проверить п.5.

7.Вычисляем и записываем в четвёртой строке раздела 2

8.Проверка как в п.4.

9.Вычисляем и записываем в первые две строки раздела 3.

10.Проверка как в п.4.

11.Вычисляем (j=3,4,5,6) и записываем в третьей строке раздела 3.

12.Проверка как в п.4.

13. Вычисляем и записываем в первую строку раздела 4.

i

ai1

ai2

ai3

ai4

ai5

åai6

1

1

2

3

4

9.3

4.92

4.77

3.21

1.0

2.42

7.45

7.01

1.85

0.2602

6.1

8.9

9.04

3.69

0.6559

1.9

2.46

2.28

6.99

0.2043

-9.05

10.21

13.45

-10.35

-0.9731

10.67

33.94

36.55

5.39

1.1473

2

2

3

4

6.1698

5.7688

1.0148

1.0

5.6730

5.9114

1.5846

0.9195

1.4548

1.3055

6.3342

0.2358

14.9977

18.0918

-7.2263

2.4308

28.2953

31.0775

1.7073

4.5861

3

3

4

0.6069

0.6515

1

-0.0547

6.0949

-0.0901

4.0690

-9.6931

6.7045

4.6212

-2.9467

7.6144

4

5

4

1

1

1

6.1536

1

-14.0611

-2.2850

6,4986

-3.0059

-3.9866

-7.9075

-1.2850

7,4986

-2.0059

-2.9866


Страница: