Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ

НАЦИОНАЛЬНЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ «ХПИ»

Кафедра «Вычислительной техники и програмирования»

Расчётно–графическое задание

по курсу «Теория алгоритмов и вычислительные методы»

Харьков – 2005

Исходные данные:

Вариант №	y0	y1	y2	y3	y4	y5	h	x0
64	-0.02	0.604	0.292	-0.512	-1.284	-2.04	0.5	0.3

Задача 1

Исходные данные вводятся в ЭВМ как абсолютно точные числа и представляются в ней в виде чисел с плавающей точкой с относительной погрешностью в одну миллионную. Введенные данные x₀ и y₀ служат основой формирования двух векторов x=(x₀, x₁, …, x_n) и y=(y₀, y₁, …, y_n) по рекуррентным формулам:

Вычислить скалярное произведение с := (x, y) по алгоритму:

с := 0; i := 0;

while i < n + 1 do c := c + x_i · y_i;

и оценить аналитически и численно инструментальную абсолютную и относительную погрешности.

Решение

Поскольку данные представляются в ЭВМ в виде чисел с плавающей точкой с относительной погрешностью, то

x₀ = x₀(1+δ)

y₀ = y₀(1+δ)

C₀ = x₀y₀(1+δ)

П ри i = 1

При i = 2

x₂ = x₀³(1+δ)⁵

y₂ = y₀(1+δ)³

C₂ = x₀y₀(1+δ)⁵ + x₀²(1+δ)⁷ + x₀³y₀(1+δ)¹⁰

При i = 3

x₃ = x₀⁴(1+δ)⁷

y₃ = (1+δ)⁵

C₃ = x₀y₀(1+δ)⁶ + x₀²(1+δ)⁸ + x₀³y₀(1+δ)¹¹ + x₀⁴(1+δ)¹⁴

При i = 4

x₄ = x₀⁵(1+δ)⁹

y₄ = y₀(1+δ)⁷

C₄ = x₀y₀(1+δ)⁷ + x₀²(1+δ)⁹ + x₀³y₀(1+δ)¹² + x₀⁴(1+δ)¹⁵ + x₀⁵y₀(1+δ)¹⁸

Выявим закономерность изменения C_i:

При расчете C_n без учета погрешности исходных данных и погрешности вычисления, получим

Обозначим эту сумму как S₁.

Тогда абсолютная погрешность S₂

а относительная погрешность

Оценим инструментально относительную и абсолютные погрешности при n = 10

S₁ = 0.0923071

S₂ = 1.45914·10^-6

S₃ = 1.58075·10^-5

Задача 2

Для функции g(x), заданной своими значениями в шести точках, составить таблицу всех повторных разностей. Преобразовать функцию g(x) с помощью линейного преобразования x = a + b * k в функцию G(k) с целочисленным аргументом k. В качестве проверки правильности заполнения таблицы вычислить аналитически конечную разность Δⁿg(x) = ΔⁿG(k) для n = 5.

Решение

Составим таблицу всех повторных разностей:

k	x	y	Δy	Δ2y	Δ3y	Δ4y	Δ5y
0	0.3	0.02	-1.576	0.044	-0.136	0.66	-0.54
1	1.1	-1.556	-1.532	-0.092	0.524	0.12	—
2	1.9	-3.088	-1.624	0.432	0.644	—	—
3	2.7	-4.712	-1.192	1.076	—	—	—
4	3.5	-5.904	-0.116	—	—	—	—
5	4.3	-6.02	—	—	—	—	—

Найдем формулу перехода от x к k:

В ыполним проверку, вычислив аналитически конечную разность

Δⁿg(x)= ΔⁿG(k) для n = 5:

Конечные разности, вычисленные аналитически и таблично Δⁿg(x) = ΔⁿG(k) для n = 5 совпали, следовательно, таблица повторных разностей составлена верно.

Задача 3

Таблично заданную функцию G(k) с целочисленным аргументом представить в виде разложения по факториальным многочленам (z⁽ⁿ⁾ = z · (z-1) · (z-2) · … · (z - n + 1)) и преобразовать его в степенные многочлены G(z) и G(x).

Решение

Представим функцию G(k) в виде разложения по факториальным многочленам:

Преобразуем функцию G(k) в степенной многочлен G(z):

Выполним проверку при k = 1:

0.604=0.604

Так как результаты совпали, значит степенной многочлен G(z) представлен правильно.

Преобразуем функцию G(k) в степенной многочлен G(x). Зная, что получим:

Проверим вычисления при x = 0.8:

0.6045128 ≈ 0.604

Так как результаты совпали, то вычисления сделаны верно.

Задача 4

Вывести аналитическое выражение суммы для функции целочисленного аргумента G(z). Проверить правильность вычисления полученного выражения прямым суммированием табличных значений G(k), k = 0, 1, 2, 3, 4, 5 (m = 5).

Решение.

Для вычисления значения суммы используем функцию G(z) в виде разложения по факториальным многочленам, полученным в задаче 3:

г де

Для проверки, просуммируем значения G(k) из таблицы:

-0.02 + 0.604 + 0.292 - 0.512 - 1.284 - 2.04 = - 2.96

- 2.96 = - 2.96

Так как результаты вычисления аналитического выражения и суммы табличных значений G(k) совпали, значит аналитическое выражение для суммы выведено правильно.

Задача 5

Составить таблицу упорядоченных разделенных разностей для g(x). Проверить правильность таблицы для разделенной разности [x₀; x₁; x₂; x₃] по формуле ее аналитического представления.

Решение

Составим таблицу упорядоченных разделенных разностей для g(x):

xi	g(xi)	[xi; xi+1]	[xi; xi+1; xi+2]	[xi; xi+1; xi+2; xi+3]	[xi; xi+1; xi+2; xi+3; xi+4]	[xi; xi+1; xi+2; xi+3; xi+4;xi+5]
0.3	-0.02	1.248	-1.872	0.592	0.0533333	-0.1567999
0.8	0.604	-0.624	-0.984	0.6986666	-0.3386666	—
1.3	0.292	-1.608	0.064	-0.0213333	—	—
1.8	-0.512	-1.544	0.032	—	—	—
2.3	-1.284	-1.512	—	—	—	—
2.8	-2.04	—	—	—	—	—

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

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

Задача 6

Получить интерполяционные многочлены Лагранжа и Ньютона, проходящие через первые четыре точки таблично заданной функции G(x), и сравнить их степенные представления.

Решение

Для нахождения интерполяционного многочлена Лагранжа используем формулу

где n = 3.

Проведем проверку вычислений, подставив x=0.8 в интерполяционный многочлен Лагранжа, получим y₁=0.604

Интерполяционный многочлен Ньютона находится по формуле:

l_n(x) = g₀ + (x-x₀)[x₀;x₁] + (x-x₀)(x-x₁)[x₀;x₁;x₂] + … +

+(x-x₀)(x-x₁)∙ …∙(x-x_n-1)[x₀;x₁;x₂;…;x_n]

П одставив в формулу g_i и x_i получим:

Интерполяционные многочлены Ньютона и Лагранжа совпадают.

Проведем проверку вычислений, подставив x=0.8 в интерполяционный многочлен Ньютона, получим y₁=0.604

Задача 7.

Вывести выражения для вычисления второй производной в точке x=x₃ в виде функций:

где ∆ⁿg(0) и g(x_n) для n = 0,1,…,5 соответственно значения разностей в точке x = x₀ и ординаты g(x_n) = g_n из задачи N2. Значения производной вычисленные по выведенным формулам, сравнить с вычисленным значением производной, найденной путем дифференцирования интерполяционного многочлена G(x):

Решение

Для вычисления производной воспользуемся оператором д ифференцирования:

Выражение для вычисления производной в точке x₀ имеет вид:

Для того, чтобы преобразовать его к выражению для вычисления производной в точке x₃, применим оператор сдвига:

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

Тип файла документ

Документы такого типа открываются такими программами, как Microsoft Office Word на компьютерах Windows, Apple Pages на компьютерах Mac, Open Office - бесплатная альтернатива на различных платформах, в том числе Linux. Наиболее простым и современным решением будут Google документы, так как открываются онлайн без скачивания прямо в браузере на любой платформе. Существуют российские качественные аналоги, например от Яндекса.

Будьте внимательны на мобильных устройствах, так как там используются упрощённый функционал даже в официальном приложении от Microsoft, поэтому для просмотра скачивайте PDF-версию. А если нужно редактировать файл, то используйте оригинальный файл.

Файлы такого типа обычно разбиты на страницы, а текст может быть форматированным (жирный, курсив, выбор шрифта, таблицы и т.п.), а также в него можно добавлять изображения. Формат идеально подходит для рефератов, докладов и РПЗ курсовых проектов, которые необходимо распечатать. Кстати перед печатью также сохраняйте файл в PDF, так как принтер может начудить со шрифтами.

Список файлов ответов (шпаргалок)