СОДЕРЖАНИЕ
ВВЕДЕНИЕ
1. Задача 1
1.1 Постановка задачи
1.2 Решение
2. Задача 2
2. 1. Постановка задачи
2. 2. Решение
3. Задача 3
3. 1. Постановка задачи
3. 2. Решение
4. Задача 4
4. 1. Постановка задачи
4. 2. Решение
СПИСОК ЛИТЕРАТУРЫ
ВВЕДЕНИЕ
Основой автоматизации умственного труда человека является широкое внедрение вычислительной техники во все сферы деятельности человека. Применение ЭВМ ускорило процесс математизации науки и техники. Расширяется круг профессий,для которых математическая грамотность и наличие практических навыков применения ЭВМ становятся необходимыми.
Решение технической или научной задачи включает её математическое описание на языке уравнений, функций. Очень часто математическая формулировка задачи может оказаться непереводимой на язык ЭВМ,так как ЭВМ выполняет только арифметические действия.
Численный метод решения задачи –это определённая последовательность операций над числами, язык которого числа и арифметические действия. Численные методы легко реализуются на ЭВМ,что делает эти методы мощным и универсальным инструментом. Процесс решения инженерной задачи на ЭВМ сложный и длительный. Он включает в себя этапы, требующие от разработчика профессиональной подготовки и грамотности. Для снижения трудоёмкости, на всех типах ЭВМ создан мощный аппарат технологической поддержки работы пользователя ЭВМ.
1. Задача 1
1. 1. Постановка задачи
Необходимо графически определить один корень уравнения. Уточнить корень уравнения с точностью Е=0,001 методом Ньютона. Дано нелинейное уравнение :
tg(ax+b)=x2
где a=0,5 и b=0,2
1. 2. Решение
Для того,чтобы определить корень,преобразуем уравнение к виду :
tg(0. 5x+0. 2)=x2
Построим графикидвух функций :
y1= tg(0. 5x+0. 2) и y2=x2;
Кривые на рис. 1 описаны следующим образом:
1) y1= tg(0. 5x+0. 2) функция периодическая,её значения сведём в таблицу 1. 1
Таблица 1. 1.
x | -3. 1 | -3 | -2 | -1 | 0 | 1 | 2 | 2. 1 | 2. 2 |
y | -4. 45 | -2. 57 | -1. 02 | -0,3 | 0,2 | 0,84 | 2. 57 | 3. 0 | 3. 6 |
2) y2=x2 – парабола
y2=0 когда x=0
y2=4 при x=±2
По графику определяем,что уравнение имеет несколько корней. Для уточнения корня выберем интервал [0,1]. Уточняем корень по формуле Ньютона:
xn+1= xn-
Необходимо выбрать начальное значение x0, исходя из условия сходимости:
f(x0)f "(x0)>0
f(x)= tg(0. 5x+0. 2) – x2
Проверяем условия сходимости для x=0 :
f(0)f"(0)<0,условие не соблюдается
Проверяем условие сходимости для x=1. 0 :
f(0)f"(0)>0,условие соблюдается
берём за x0=1
и условие:
Т=
Решение запишем в виде таблицы:
n | x n | f(x n) | f '(x n) | T<E 10-1 | |
0 | 1. 000000 | -0. 158000 | -1. 151000 | 0. 137271 | Нет |
1 | 0. 862728 | -0. 013000 | -0. 976000 | 0. 013119 | Нет |
2 | 0. 849416 | -0. 000467 | -0. 958000 | 0. 000487 | Нет |
3 | 0. 848929 | -0. 000009 | -0. 958000 | 0. 000009 | Да |
4 | 0. 848920 |
В результате проделанной работы мы определили один корень уравнения вида tg(0. 5x+0. 2)=x2 графически,а затем уточнили его методом Ньютона и получили
X=0. 848929
Вывод по решению:
В результате проделанной работы мы определили один корень уравнения
Tg(0. 5x+0. 2)=x2графически, а затем уточнили его методом Ньютона и получили x=0. 848929
2. Задача 2
2. 1. Постановка задачи
Выбрать формулу интерполяции и с её помощью определить значение функции в точке x=0,38. Функция задана в виде таблицы 2. 1,Степень интерполяционного многочлена равна 3.
Таблица 2. 1
0,15 | 0,860708 |
0,25 | 0,778801 |
0,30 | 0,740818 |
0,40 | 0,670320 |
0,45 | 0,637628 |
0,55 | 0,576950 |
0,60 | 0,548812 |
0,65 | 0,522046 |
0,70 | 0,496585 |
0,75 | 0,472237 |
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,38y=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.