Lab 1. Fundamentals of the theory of errors

Lab 1. Fundamentals of the theory of errors

Objective: Learn some basic concepts of the theory of errors.

Task 1: Dan number . Find the sum of a number of analytically. Calculate

partial sums of and find the value of errors in values N = 10, 10², 10³, 10⁴.

option			option

Algorithm and implementation of solutions TO MS EXCEL

Algorithm

Find the sum of a number of analytically.

Using the function. , calculate the partial sums of the indicated values N.

For each N calculate the absolute error.

Let а – exact value, а * – the exact value of the approximate value of a certain value. Absolute error of the approximate value а * It is the quantity .

Determine the number of correct digits in .

 

The number of significant digits and * called all the numbers in its records, starting with the first non-zero left.

Significant digits of a * is called valid if the absolute error does not exceed the number of discharge units corresponding to this figure.

Approximate number can be represented as a finite decimal fraction .

Then, if the figure in the image of a * true, then the inequality , , more often .

Analytical solution

Dan number

.

Implementation in MS Excel

3.1 To calculate the sum of the series to create the n auxiliary column: Edit - Fill - Progression...:

3.1 Carry out the calculation:

View the worksheet with the result of the calculation

Worksheet with formulas

Notes: The curly brackets indicate that the corresponding formula is derived array, ie, using a combination of Ctrl + Shift + Enter...

Control questions:

 

1. What is called absolute and relative accuracy of the approximate number?

2. What is called the significant digit number?

3. What is called the true number of characters?

4. Formulate valuation rules limiting errors when performing operations on the approximate numbers.

 

Lab number 2

Subject: Solution of nonlinear equations. bisection method.

Task:

1) Separate the roots of the equation graphically and analytically.

2) Clarify the roots (all!) Of the bisection method with a precision specify the number of partitions of the interval.

Questions of self-control

1) How to separate the roots of the equation?

2) What should be the value of the step at the roots of the department?

3) What conditions must be met for the application of the method of bisection of the segment?

 

option	the equation		option	the equation

 

A sample implementation of laboratory work number 2

(The solution of nonlinear equations. Bisection method.)

Formulation of the problem.

Find the root of the nonlinear equation bisection method with a precision .

The solution of the problem. Separate the root of the equation on the interval graphical method. For this function tabeliruem in this segment

We have , , , ,

We select the segment , containing isolated root, to clarify that the bisection method applicable under the scheme , , , , . believing , , as well as the condition to stop the interval bisection , up a table

						корень	погрешность	Усл.ост.
1,00000000	3,00000000	2,00000000	1,29583687	-1,17168626	0,15888308		1,00000000	нет
2,00000000	3,00000000	2,50000000	0,15888308	-1,17168626	-0,48776781		0,50000000	нет
2,00000000	2,50000000	2,25000000	0,15888308	-0,48776781	-0,15924305		0,25000000	нет
2,00000000	2,25000000	2,12500000	0,15888308	-0,15924305	0,00119806		0,12500000	нет
2,12500000	2,25000000	2,18750000	0,00119806	-0,15924305	-0,07868831		0,06250000	нет
2,12500000	2,18750000	2,15625000	0,00119806	-0,07868831	-0,03866032		0,03125000	нет
2,12500000	2,15625000	2,14062500	0,00119806	-0,03866032	-0,01870977		0,01562500	нет
2,12500000	2,14062500	2,13281250	0,00119806	-0,01870977	-0,00875050		0,00781250	нет
2,12500000	2,13281250	2,12890625	0,00119806	-0,00875050	-0,00377488		0,00390625	нет
2,12500000	2,12890625	2,12695313	0,00119806	-0,00377488	-0,00128807		0,00195313	нет
2,12500000	2,12695313	2,12597656	0,00119806	-0,00128807	-0,00004492		0,00097656	нет
2,12500000	2,12597656	2,12548828	0,00119806	-0,00004492	0,00057659		0,00048828	нет
2,12548828	2,12597656	2,12573242	0,00057659	-0,00004492	0,00026584		0,00024414	нет
2,12573242	2,12597656	2,12585449	0,00026584	-0,00004492	0,00011046		0,00012207	нет
2,12585449	2,12597656	2,12591553	0,00011046	-0,00004492	0,00003277	2,12591553	0,00006104	да
2,12591553	2,12597656	2,12594604	0,00003277	-0,00004492	-0,00000608	2,12594604	0,00003052	да
2,12591553	2,12594604	2,12593079	0,00003277	-0,00000608	0,00001335	2,12593079	0,00001526	да

 

Approximate solution , error , number of iterations .

Therefore, the approximate value of the root is equal to .

We write the approximate value of the root only true significant digits in the narrow sense.

We have , , . Rounded до . We get , , .

Let us find the number of correct digits for . We have , , . As , we obtain an approximate value of the root of the number of correct digits .

Answer:

Subject: Solution of nonlinear equations. method of chords and tangents method

Task: 1) Separate the roots of the equation graphically and programmatically.

2) Clarify the roots of this method with the accuracy .

3) Draw a diagram of the method to each root of the equation.

 

option	The equation		option	The equation

Questions of self-control

1) What is the idea of the method of bisection of the segment? Geometric illustration.

2) How do you calculate the approximate root of the equation and what is its uncertainty?

3) As a result of the error depends on the choice of the approximate solution?

Лабораторная работа № 4

Тема: Решение нелинейных уравнений. Метод итерации.

Задание: 1) Отделить корни уравнения графически и программно.

2) Уточнить один из корней уравнения методом итерации с точностью , указать число итераций.

3) Нарисовать схему применения метода итерации к данному корню уравнения.

Вопросы самоконтроля.

1) Как отделяются корни уравнения?

2) Какой должна быть величина шага при отделении корней?

3) Какие условия должны быть выполнены для применения метода итерации?

4) Какова идея метода итерации? Геометрическая иллюстрация.

5) Какое условие должно выполняться для сходимости итерационной последовательности?

6) Как находится равносильное уравнение, применяемое для итерационного процесса? Критерий выбора равносильного уравнения.

7) Как определяется погрешность метода итерации при заданной точности?

8) Какие положительные и отрицательные стороны метода итерации (сравнить с методом деления отрезка пополам)?

Вариант	Уравнение		Вариант	Уравнение

 

Теоретические сведения

Дана система n линейных уравнений с n неизвестными

 

(1)

 

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

Наиболее распространенным является метод Гаусса, согласно которому путем последовательного исключения неизвестных система (1) приводится к треугольному виду:

 

(2)

 

Приведение матрицы системы к виду (2) называется прямым ходом. Вычисление неизвестных – обратным ходом. Необходимость округлять в промежуточных вычислениях приводит к тому, что возникает очень большая погрешность округления, искажающая результат. Существует несколько видов вычислительных схем метода Гаусса, в различной степени уменьшающих погрешность округления. Наиболее эффективной является “схема с выбором главного элемента”.

 

Прямой ход

Исключение неизвестных в прямом ходе осуществляется по этапам. На каждом i-м этапе (i=1, 2, …, n) среди коэффициентов при неизвестном выбираем наибольший по абсолютной величине – “главный элемент”. Строка, его содержащая, называется главной. Затем главную строку прибавляем ко всем остальным строкам, предварительно умножив ее на специально подобранные числа так, чтобы коэффициенты при во всех строках, кроме главной, обратились бы в нуль.

Главная строка i-го этапа в дальнейших преобразованиях не участвует, поэтому для приведения системы (1) к виду (2) нужно проделать исключение неизвестных.

 

Обратный ход

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

 

Прямой ход

В первом столбце поставим номера разделов.

Раздел первый

 

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

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

.

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

5. Элементы всех строк первого столбца, кроме главного, делим на него и результат с противоположным знаком записываем в столбец , так, , , .

 

Раздел второй

6. Исключаем элементы первого столбца, кроме главного. Для этого главную строку, умноженную на , прибавляем к i-й неглавной строке, результат записываем в раздел второй. Так, первая строка второго раздела получается из первой неглавной строки первого раздела:

 

 

=

 

.

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

8. Аналогично вычисляем вторую и третью строки второго раздела

+

 

=

 

.

+

 

=

 

.

 

9. Проводим контроль, как в п. 7.

 

10. Среди коэффициентов при неизвестном выбираем главный (-1,150), подчеркиваем главную строку и записываем ее в качестве второй в таблицу для обратного хода.

Далее, делим на главный элемент все числа, стоящие с ним в одном столбце (0,300 и -0,300), результаты записываем с противоположным знаком в качестве . Затем производим исключение коэффициентов при , как в пп. 6 – 8.

 

11. Третий раздел будет иметь две строки. Главный элемент выбираем среди коэффициентов при и поступаем, как в п. 6. Получаем строку последнего четвертого раздела. Записываем ее в таблицу обратного хода, сделав предварительно контроль, как в п. 7.

 

Обратный ход

1. Последняя строка таблицы соответствует уравнению . Отсюда находим .

2. Контроль: вместо свободного члена берем , находим . Проверяем выполнение условия :

.

 

3. Найденное значение подставляем в следующее уравнение, соответствующее третьей главной строке: , отсюда

 

.

 

Контроль: ,

 

.

4. Найденные значения и подставим в следующее уравнение, получим . Также по известным получим .

 

5. Из первого уравнения получим .

Проверка. Значения неизвестных подставим в каждое из уравнений исходной системы и подсчитаем невязки (см. разд. 3):

 

 

 

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

Ответ: Решением данной системы будет вектор (в ответе оставляем только 1 запасную цифру).

 

1. Матричные функции Excel

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

МУМНОЖ(<матрица1>;<матрица2>) – возвращает произведение матриц.

МОБР(<матрица>) – возвращает матрицу, обратную к данной.

МОПРЕД(<матрица>) – вычисляет определитель исходной квадратной матрицы.

Порядок обращения к матричным функциям:

Ø выделить блок, где будет размещен результат матричной операции;

Ø в мастере функций выбрать нужную категорию и нужную функцию;

Ø убрать окно соответствующей функции (с помощью кнопки) или отбуксировать его в сторону от исходных данных;

Ø выделить исходную матрицу (бегущая пунктирная линия);

Ø одновременно нажать клавиши Shift+Ctrl+Enter.

1. Реализация метода Гаусса средствами приложения Excel

Пример 1.1. Найти решение системы линейных алгебраических уравнений

, (1.8)

используя алгоритм метода Гаусса.

 

Последовательность действий:

Возьмем чистый лист Excel, назовем его Гаусс. Введем расширенную матрицу системы, как показано на рис.1.1 в ячейки А3:D5.

Первый этап, приведение матрицы системы к треугольной.

1. Поделим элементы первой строки на а_{11 .}Для этого в ячейку А7 введем формулу А7=А3/A$3$ (*смотри сноску) и скопируем ее вправо до конца строки.

2. Умножим элементы первой строки на (-а₂₁) и прибавим ко 2-ой строке. Для этого введем формулу А8=А7(-А$4$)+А4 и скопируем ее вправо до конца строки.

3. Умножим элементы первой строки на (- а₃₁ ) и прибавим к 3-ей строке. Для этого введем формулу А9=А7(-А$5$)+А5 и скопируем ее вправо до конца строки. Таким образом исключили неизвестное х₁ из 2-го и 3-го уравнений системы (смотри 1-ый шаг рис.1.1).

4. Осталось исключить неизвестное х₂ из 3-го уравнения системы. Для этого реализуем описанный выше алгоритм для 2-ой и3-ей строк (смотри 2-ой шаг рис.1.1).

 

Рис.1.1.

На этом первый этап метода Гаусса, закончен, матрица системы приведена к треугольной.

Второй этап. Здесь последовательно найдем неизвестные, начиная с последней строки. Для этого в ячейки G12:G14 запишем формулы:

G4=D13/C13 (для вычисления x ₃);

G3=D12-C12*G4 (для вычисления x₂);

G2=D11-C11*G4-B11*G3 (для вычисления x₁).

 

Последовательность действий

1.

Возьмем новый лист (а можно и на том же). Заготовим таблицу, как показано на рис.1.2.

Рис.1.2.

2. Заготовим ячейки А7:С7, где будет сформировано решение системы (х₁, х₂, х₃). Первоначально они остаются пустыми, т.е. равными нулю. Однако для контроля правильности вводимых далее формул, удобно ввести в эти ячейки какие-либо значения, например единицы. Эти значения можно рассматривать как нулевое приближение решения системы, .

3. Введем коэффициенты системы (матрицу А) в ячейки А3:С5.

4. В столбец D введем выражения для вычисления левых частей исходной системы. Для этого в ячейке D 3 введем и скопируем вниз до конца таблицы формулу: D3=СУММПРОИЗВ (A3:C3;$A$7:$C$7).

Используемая функция СУММПРОИЗВ принадлежит категории Математические.

5. В столбец Е запишем значения правых частей системы матрицу .

6. В столбец F введем невязки в соответствии с формулой (1.9), т.е. введем формулу F3=D3-E3 и скопируем ее вниз до конца таблицы.

7. Будет не лишним проверить правильность вычислений для случая .

8. Зададим команду меню Сервис\Поиск решения. В окне Поиск решения (рис.1.3) в поле Изменяя ячейки укажем блок $А$7:$С$7, а в поле Ограничения – $F$3:$F$5 =0. Для этого надо щелкнуть на кнопке Добавить и ввести эти ограничения.

9. Щелкнем на кнопке Выполнить.

Рис. 1.3.

Полученное решение системы (1.8) х₁ =1; х₂ =-1 х₃ =2 записано в ячейках А7:С7, рис.1.2.

 

Вопросы самоконтроля.

1) Как отделяются корни уравнения?

2) Как используется функция СУММПРОИЗВ?

3) Какой должна быть величина шага при отделении корней?

4) Какие условия должны быть выполнены для применения метода половинного деления отрезка?

 

Задания к лабораторным работам № 5-7

 

Найти решение данной системы

 

 

№ варианта	Коэффициенты при неизвестных	Свободные члены
	0,11270	-2,39990	8,95146	0,75000	8,60527
9,58778	-3,45350	0,24300	1,46840	16,40216
0,86400	4,23700	-2,50200	-1,72927	-15,88846
-0,28427	-4,58674	-1,85970	0,14940	10,90588
	1,11270	-3,02270	-10,91328	1,06140	11,56420