Мы поможем в написании ваших работ!



ЗНАЕТЕ ЛИ ВЫ?

Нахождение параметров линейной функции

Поиск

 

Предположим, что зависимость между x и y линейная, т.е. приближающую функцию можно записать в виде y=ax+b.

Нужно найти такие значения a и b, для которых функция

(3.6)

минимальна.

Условия функции (3.6) запишутся так:

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

.

Суммы вычисляются по табличным данным. Для удобства вычисления можно составить расчетную таблицу:

 

 

N x Y xy x2
  x1 y1 x1y1 x12
… n x2 ... xn y2 … yn x2y2 … xnyn x22 … xn2
Σ

 

 

Нахождение параметров квадратичной функции

 

Если известно, что приближающей функцией является квадратичная функция y=ax2+bx+c, то ее коэффициенты a, b, c найдем из условия минимума функции

.

Условия минимума:

Получаем для нахождения неизвестных a, b, c систему трех уравнений, которую решаем методом Гаусса.

 

 

Расчетная таблица

 

N x y x2 x3 x4 x2y xy
  x1 y1 x12 x13 x14 x12y1 x1y1
… n x2 ... xn y2 … yn x22 … xn2 x23 … xn3 x24 … xn4 x22y2 … xn2yn x2y2 … xnyn
Σ

 

Задания. Аппроксимировать функцию, заданную таблично, некоторой функцией по методу указанному преподавателем.

 

ВАРИАНТЫ

 

№ 1 x 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8
y 3,030 3,142 3,251 3,858 3,463 3,563 3,665 3,772
№ 2 x 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8
y 3,314 3,278 3,262 3,268 3,292 3,332 3,397 3,486
№ 3 x 2,0 2,1 2,2 2,3 2,4 2,5 2,6 2,7
y 1,045 1,162 1,204 1,172 1,070 0,898 0,656 0,344
№ 4 x 0,3 0,6 0,9 1,2 1,5 1,8 2,1 2,4
y 6,715 6,735 6,750 6,741 6,647 6,649 6,645 6,636
№ 5 x 1,2 1,4 1,6 1,8 2,0 2,2 2,4 2,6
y 2,325 2,515 2,638 2,790 2,696 2,626 2,491 2,291
№ 6 x 2,1 2,3 2,5 2,7 2,9 3,1 3,3 3,5
y 1,752 1,762 1,777 1,797 1,821 1,850 1,884 1,924
№ 7 x 1,5 1,6 1,7 1,8 1,9 2,0 2,1 2,2
y 1,924 1,710 1,525 1,370 1,264 1,88 1,142 1,127
№ 8 x 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2
y 1,035 1,144 1,248 1,336 1,409 1,467 1,510 1,538
№ 9 x 2,2 2,3 2,4 2,5 2,6 2,7 2,8 2,9
y 5,785 5,685 5,605 5,545, 5,505 5,485 5,490 5,506
№ 10 x 0,5 0,7 0,9 1,1 1,3 1,5 1,7 1,9
y 4,052 4,092 4,152 4,234 4,336 4,458 4,599 4,761
№ 11 x 0,3 0,5 0,7 0,9 1,1 1,3 1,5 1,7
y 0,344 0,364 0,374 0,372 0,350 0,328 0,296 0,256
№ 12 x 0,4 0,6 0,8 1,0 1,2 1,4 1,6 1,8
y 0,205 0,235 0,249 0,245 0,225 0,190 0,140 0,076
№ 13 x 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8
y 1,044 1,161 1,203 1,172 1,070 0,896 0,654 0,342
№ 14 x 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8
y 0,525 0,625 0,678 0,681 0,640 0,552 0,432 0,362
№ 15 x 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8
y 4,230 4,253 4,256 4,240 4,205 4,150 4,075 3,980
№ 16 x 2,0 2,1 2,2 2,3 2,4 2,5 2,6 2,7
y 5,022 5,143 5,195 5,175 5,085 4,925 4,705 4,406
№ 17 x 3,0 3,1 3,2 3,3 3,4 3,5 3,6 3,7
y 1,125 1,175 1,212 1,237 1,251 1,255 1,246 1,225
№ 18 x 2,0 2,1 2,2 2,3 2,4 2,5 2,6 2,7
y 1,220 1,253 1,256 1,232 1,175 1,091 0,985 0,850
№ 19 x 0,5 0,7 0,9 1,1 1,3 1,5 1,7 1,9
y 3,150 3,171 3,181 3,179 3,165 3,140 3,105 3,059
№ 20 x 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9
y 4,018 4,025 4,035 4,048 4,063 4,080 4,099 4,120
№ 21 x 3,1 3,2 3,3 3,4 3,5 3,6 3,7 3,8
y 2,527 2,635 2,655 2,563 2,361 2,048 1,638 1,118
№ 22 x 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
y 4,030 4,142 4,251 4,358 4,468 4,561 4,465 4,762
№ 23 x 1,0 1,1 1,2 1,3 1,4 1,5 1,6 1,7
y 1,314 1,278 1,262 1,266 1,252 1,332 1,397 1,486

 

 

4. РЕШЕНИЕ СИСТЕМ ЛИНЕЙНЫХ УРАВНЕНИЙ

МЕТОДОМ ГАУССА

 

Существует несколько приемов решения систем линейных уравнений. Один из таких приемов называется методом Гаусса или методом последовательного исключения неизвестных. Этот метод используется и в электронных вычислительных машинах. Ниже описывается способ решения системы линейных уравнений по методу Гаусса с использованием микрокалькуляторов.

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

В общем виде система трех линейных уравнений с тремя неизвестными имеет вид:

(1)

 

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

Разделим все члены первого уравнения системы (1) на a11, затем все члены второго уравнения – на a21, все члены третьего уравнения – на a31. Если a21 или a31 равен нулю, то соответствующее уравнение входит в новую систему без изменений.

Затем вычтем из каждого последующего уравнения первое полученное уравнение. В результате получим систему:

(2)

Предположим, что в системе (2) хотя бы один из коэффициентов a22 или a32 отличен от нуля. Пусть таким является, например, коэффициент a22. разделим почленно второе уравнение системы (2) на a22, третье – на a32. затем из третьего полученного уравнения вычтем второе полученное уравнение. В результате этих преобразований система примет вид:

(3)

Для того, проводить поэтапный контроль за результатами преобразований уравнений, введем новые определения. Число, равное суме всех коэффициентов данного уравнения, включая его свободный член, будем называть строчной суммой. Контрольной суммой будем называть то число, которое на первом этапе совпадает с соответствующей строчной суммой, а на последующих этапах получается из контрольных сумм предыдущего этапа в результате того же преобразования, которое проводится над соответствующим уравнением. Например, если уравнение получается делением предыдущего уравнения на какое-либо число, то и новая контрольная сумма получается из соответствующей предыдущей контрольной суммы делением на это число. Легко видеть, что строчная и контрольная суммы одного и того же уравнения должны быть практически равными. В их записи различие может быть только в запасных знаках.

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

При замене правых частей системы уравнений соответствующими строчными суммами получается система уравнений такая, что если x1, x2, x3 есть решение первоначальной системы, то x1+1, x2+1, x3+1 будет решение новой системы.

В самом деле, пусть, например, S1 – строчная сумма первого уравнения, а x1, x2, x3 – решение системы (1), тогда

 

.

Таким образом, числа удовлетворяют первому уравнению новой системы (1).

Наряду с системами (1), (2), (3), будем иметь в виду еще три системы (1’), (2’), (3’). Левые части параллельных пар систем (1) и (1’), (2) и (2’), (3) и (3’) соответственно совпадают, а правые части систем (1’), (2’), (3’) получаются заменой свободных членов системы (1), (2), (3) соответствующими строчными суммами. Системы (1), (2), (3) эквивалентны между собой; системы (1’), (2’), (3’) – также эквивалентны между собой. Причем система (3’) имеет вид:

(3’)

Используя «обратный ход», из системы (3) последовательно находим x3, x2, x1; параллельно из системы (3’) последовательно находим . Числа должны отличаться на 1 (расхождение в запасных знаках во внимание не принимается).

Все расчеты выразим в виде схемы (см. таблицу 1). Для компактности схемы систему, эквивалентную данной, запишем только при выполнении «обратного хода».

I этап: запись данной системы, вычисление строчных (контрольных) сумм.

II этап: деление каждого уравнения на коэффициент при x1, вычисление строчных и контрольных сумм.

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

IV этап: деление 2-го и 3-го уравнения на коэффициент при x2, вычисление строчных и контрольных сумм.

V этап: вычитание 2-го уравнения из 3-го, вычисление строчных и контрольных сумм.

VI этап: вычисление .

VII этап: вычисление .

VIII этап: вычисление .

 

Пример. Решить систему трех линейных уравнений с тремя неизвестными с точностью до 0,001.

Составим схему расчета и контроля (см. таблицу 2).

 

 

 

ТАБЛИЦА 1

СХЕМА РЕШЕНИЙ СИСТЕМ ПО МЕТОДУ ГАУССА

 

 

Прямой ход
Система Строчные суммы Контрольные суммы Этапы
I
II
III
IV
V

, где

.

 

 

Обратный ход   Контроль Этапы
VI     VII   VIII

 

ТАБЛИЦА 2

Коэффициенты при неизвестных Свободные члены Строчные суммы Контрольные суммы
x1 x2 x3
5,434 -1,304 2,236 1,925 8,305 -0,489 3,241 5,435 4,283 48,500 50,230 30,780 59,100 62,666 36,900 59,100 62,666 35,900
  0,35425 -6,36886 -0,21023 0,59642 -4,16794 1,84135 8,92528 -38,51993 13,23301 10,87292 -48,05673 15,86414 10,87292 -48,05673 15,86414
- - -6,73211 -0,56448 -4,76436 1,24493 -47,44521 4,30773 -58,93268 4,98818 -58,93268 4,98818
- -   0,70865 -2,20544 7,05703 -7,63132 8,76568 -8,83677 8,76568 -8,83677
- - -2,91409 -14,68835 -17,60245 -17,60245

 

-2,91409 x3=-14,68835

x2+0,70865 x3=7,05703

x1+0,35425 x2+0,59642 x3=8,92528

 

Решение данной системы Решение контрольной системы Контроль
x3=-14,68835:(-2,91409)=5,04045 x3=-17,60244:(-2,91409)=6,04045 x3= x3+1
x2=7,05703-0,70865*5,04045=3,48512 x2=8,76568-0,708657*6,04045=4,48512 x2= x2+1
x1=8,92528-0,59642*5,04045-0,35425*3,48512=4,68446 x1=10,87596-0,59642*6,04045-0,35425*4,48512=5,68446 x1= x1+1

Ответ: x1=4,684, x2=3,485, x3=5,040.

 

СОДЕРЖАНИЕ

 

1. Численное интегрирование …………………………………………………. 3

1.1. Метод прямоугольников ……………………………………………...… 3

1.2. Метод трапеций …………………………………………………………. 4

1.3. Метод парабол (метод Симпсона) ……………………………………... 5

2. Приближенное решение алгебраических и трансцендентных уравнений.. 9

2.1. Уточнение корней методом половинного деления ………………….. 11

2.2. Метод итераций ………………………………………………………... 13

2.3. Метод хорд ……………………………………………………………... 13

2.4. Метод касательных (метод Ньютона) ………………………………... 16

2.5. Комбинированный метод хорд и касательных ………………………. 18

3. Аппроксимация функций ………………………………………………….. 20

3.1. Математическая постановка задачи интерполирования ……………. 20

3.2. Интерполяционный полином Лагранжа ………………………..……. 21

3.3. Линейная интерполяция …………………………………………….… 23

3.4. Метод наименьших квадратов ……………………………………...… 24

3.5. Нахождение параметров линейной функции ……………………...… 26

3.6. Нахождение параметров квадратичной функции................................ 27

4. Решение систем линейных уравнений методом Гаусса …………………. 30

 



Поделиться:


Последнее изменение этой страницы: 2016-12-16; просмотров: 551; Нарушение авторского права страницы; Мы поможем в написании вашей работы!

infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 3.145.66.104 (0.007 с.)