Скачиваний:
4
Добавлен:
27.04.2021
Размер:
293.4 Кб
Скачать

няння визначаємо з нижнього рядка таблиці, у даному випадку х = 1,4400003 з точністю 0,0002586. При потребі збільшити точність продовжуємо заповнювати рядки таблиці. Зауважимо, що кількість знаків після коми може також залежати від ширини стовпця та встановлених шрифтів відображення даних у комірках.

A

B

С

1Ітераційний метод розв'язування рівняння 10sin x 3x 5 = 0

2Послідовні наближення розв’язку

3

Попереднє значення

Наступне значення

Похибка розв’язку

 

хn

хn+1

хn+1- хn

4

1

1,1382366

0,1382366

5

1,13823662

1,2523303

0,1140937

6

1,25233028

1,3323244

0,0799942

7

1,33232443

1,3816575

0,0493331

8

1,38165753

1,4095304

0,0278729

9

1,40953039

1,4244611

0,0149307

10

1,42446108

1,4322233

0,0077622

11

1,43222331

1,4361948

0,0039715

12

1,43619483

1,4382101

0,0020152

13

1,43821007

1,4392283

0,0010182

14

1,43922831

1,4397417

0,0005134

15

1,4397417

1,4400003

0,0002586

Рис. 6. Отримані в комірках числові значення

З рис. 3 видно, що рівняння має два розв’язки. Ми отримали один з них. Змінюючи початкове наближення – значення у комірці А4 – ми отримуватимемо у таблиці лише наближення знайденого розв’язку. Другий розв’язок таким чином знайти не вдалося.

Якщо не вдалося отримати інші розв’язки, або якщо ітераційний процес не збігається, то це не означає, що інших розв’язків, або розв’язку взагалі немає. Можливо, у цьому випадку невдало зроблене перетворення рівняння, або вибране невідповідне початкове наближення. Існує можливість замість рівняння x = f (x) для отримання розв’язку задачі використати еквівалентне йому рівняння

x = f 1 (x) ,

де f 1 (x) – функція, обернена до функції f (x) . Суть такої заміни стає зрозумілою, якщо зважити, що розв’язку рівняння x = f (x) відповідають точки перетину графіків функцій y = x та y = f (x) . Але функції y = f (x) та

11

y = f 1 (x) є симетричними відносно прямої

y = x і тому перетинаються з

прямою

y = x

у тих

самих точках. Тому розв’язки

рівнянь x = f (x) та

x = f 1 (x) – однакові. Проте якщо для одного з них ітераційний метод не

збігається, то для другого може збігатися. Крім того, можна отримати ті

розв’язки, які не були знайдені.

 

 

 

 

 

 

 

 

 

Повернемося до нашого прикладу і продовжимо розв’язування рівняння.

Якщо

f (x) =

10sin

x 5 ,

то

обернена

функція

має

вигляд

 

 

 

3

2

 

 

 

 

 

 

 

 

 

 

 

3x +5

 

 

 

 

 

 

 

 

 

 

f 1 (x) = arcsin

10

 

. На рис.7 показані графіки цих функцій.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, що розв’язки рівнянь

x =

10sin

x 5

та

x =

 

3x +5 2

3

 

arcsin

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

однакові,

оскільки обидві криві на рис.3.7 перетинають пряму y = x

у тих

самих точках.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

1,6

 

 

 

 

 

 

 

 

 

 

 

 

 

1,4

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

y =10sin

x 5

 

 

 

 

 

 

 

 

0,8

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,6

 

 

 

 

 

 

3x +5 2

 

 

 

 

 

0,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y = arcsin

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

0

0,2

0,4

0,6

0,8

1

1,2

1,4

1,6

 

 

Рис. 7. Обернені функції симетричні відносно прямої y = x

Розв’яжемо рівняння, записане другому вигляді, а саме

x= arcsin 3x +5 2 .

10

Відповідно до такого вигляду рівняння сплануємо внесення формул у комірки робочого аркуша (рис. 8).

12

Вносимо і копіюємо формули у комірки таблиці з модифікацією відносних адрес комірок, виконуючи дії аналогічні тим, які проводилися при заповненні таблиці (рис. 5). Отримуємо результати обчислень за внесеними у комірки формулами (рис. 9).

A

B

С

1Ітераційний метод розв'язування рівняння 10sin x 3x 5 = 0

2Послідовні наближення розв’язку

3

Попереднє значення

 

Наступне значення хn+1

Похибка розв’язку

 

хn

 

 

 

хn+1- хn

4

1

=(ASIN((3*A4+5)/10))^2

=ABS(A4-B4)

5

=B4

=(ASIN((3*A5+5)/10))^2

=ABS(A5-B5)

6

=B5

=(ASIN((3*A6+5)/10))^2

=ABS(A6-B6)

 

 

 

Рис. 8. Розташування і вигляд формул у комірках робочого аркуша

 

 

 

 

 

 

 

 

A

 

B

 

С

1Ітераційний метод розв'язування рівняння 10sin x 3x 5 =0

2Послідовні наближення розв’язку

3

Попереднє значення

Наступне значення

Похибка розв’язку

 

хn

хn+1

хn+1- хn

4

1

0,85987642

0,1401236

5

0,85987642

0,73991738

0,119959

6

0,73991738

0,6506887

0,0892287

 

14

0,50453622

0,50321513

0,0013211

15

0,50321513

0,50247451

0,0007406

Рис. 9. Результати обчислень

Ітераційний процес збіжний, оскільки у стовпчику С значення похибки розв’язку зменшується. Отримано наближення другого розв’язку рівняння

х = 0,50247451 з точністю 0,0007406.

Записуємо відповідь задачі: х1 ≈ 1,4400003, х2 ≈ 0,50247451.

4.2. Проведення циклічних обчислень

В електронній таблиці передбачена можливість багаторазових переобчислень, чим можна скористатись для розв’язання рівняння x = f (x) за ітера-

ційною формулою xn+1 = f (xn ) .

Щоб задати режим ітераційних обчислень, виконаємо команду головного меню Сервис | Параметры… і на вкладці Вычисления діалогового вікна Параметры встановимо прапорець для пункту итерации. Крім того, доці-

13

льно в розділі Вычисления поставити перемикач з позиції автоматически у позицію вручную для ручного запуску ітераційних обчислень натисканням функціональної клавіші F9.

Внесемо у комірку А4 початкове значення x0 . У комірку В4 внесемо формулу обчислення f (A4) . Після цього змінимо дані у комірці А4, записавши у неї формулу обчислення f (B4) . Таким чином, у комірці А4 обчислення проводяться за значенням у комірці В4, а у комірці В4 для обчислень використаємо значення з комірки А4. Незважаючи на такі перехресні звертання, зациклювання не відбудеться, оскільки гранична кількість ітерацій та точність обчислень, до досягнення якої відбуваються переобчислення, в електронній таблиці задаються на вкладці Вычисления, і після запуску циклічного обчислювального процесу він обов’язково завершиться. Якщо перемикач обчис-

лень знаходиться у позиції

вручную,

після внесення формул у комірки ро-

бочого аркуша натискаємо функціональну клавішу F9, щоб запустити ітера-

ційний процес вручну.

 

 

 

 

 

 

 

 

Для прикладу ще раз отримаємо розв’язок рівняння записаного як

 

x =

10sin x 5

та x

 

 

3x +5

2

 

3

 

= arcsin

10

 

.

 

 

 

 

 

 

 

 

 

Виконаємо такі дії.

 

 

 

 

 

 

 

 

1. Задамо ітераційний режим обчислень.

 

Виконуємо команду

Сервис | Параметры… , Вычисления і встанови-

мо прапорець для пункту итерации.

 

 

 

2. Приймемо

х = 1

як початкове значення і заповнимо комірки таблиці

такими даними.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

B

 

1

 

 

 

1

 

 

 

 

=(10*SIN(КОРЕНЬ(A1))-5)/3

 

 

2

 

 

 

1

 

 

 

 

=(ASIN((3*A2+5)/10))^2

 

 

Отримаємо наступне наближення розв’язку.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

B

 

1

 

 

 

1

 

 

 

 

 

1,138236616

 

 

2

 

 

 

1

 

 

 

 

 

0,859876421

 

 

3. Змінимо формули у комірках таблиці.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

B

 

1

 

 

 

=В1

 

 

 

=(10*SIN(КОРЕНЬ(A1))-5)/3

 

 

2

 

 

 

=В2

 

 

 

=(ASIN((3*A2+5)/10))^2

 

14

Відбудеться ітераційний обчислювальний процес і отримаємо в комірках нові значення.

 

A

B

1

1,440262105

1,440262173

2

0,502474514

0,502059699

Згідно з обчисленнями в таблиці значення коренів рівняння х = 1,440262173 та х = 0,502059699, із заданою точністю обчислень співпадають з розв’язками, отриманими для цих коренів у попередньому пункті.

4. Натискаємо функціональну клавішу F9. При цьому ітераційний процес продовжується, і виконається кількість ітерацій, вказана на вкладці Вычисления діалогового вікна Параметры, якщо при цьому не буде досягнута задана відносна точність обчислень. Спостерігаємо за зміною обчислених значень коренів рівняння і врешті отримаємо х = 1,440262242 та х = 0,501532457.

Підвищення точності обчислень можна досягнути, змінивши параметри ітераційного процесу. Відкриємо вікно Параметры і переходимо на вкладку Вычисления. У текстових полях цієї вкладки збільшимо граничну кількість ітерацій та зменшимо значення відносної похибки обчислень на кілька порядків. Циклічний процес проходитиме більшу кількість повторень.

4.3. Використання сервісу “Підбір параметра”.

Алгебраїчні рівняння можна розв’язувати в Excel, використовуючи сервіс, який дозволяє автоматизувати пошук розв’язку рівняння, записаного у вигляді F(x) = 0 та контролювати точність обчислень.

Проілюструємо використання сервісу добирання параметра при знаходженні розв’язків рівняння 10sin x 3x 5 = 0 .

Для цього прикладу F(x) =10sin x 3x 5 . Виконаємо такі дії.

1.

Внесемо в комірку А1 значення 0 (початкове наближення розв'язку).

2.

У комірку В1 внесемо формулу =10*SIN(КОРЕНЬ(А1)) – 5 – 3*А1.

Укомірці В1 матимемо значення 0.

3.Виконаємо команду Сервис | Подбор параметра. Відкриється діало-

гове вікно Подбор параметра (рис. 10).

4.У поля вікна вводимо адресу комірки, що містить формулу (В1), значення, яке повинно отриматися у комірці з формулою (0), а також адресу комірки з початковим наближенням розв’язку (А1).

15

Рис.10. Вигляд вікна Подбор параметра з внесеними адресами комірок

5. Клацаємо мишею на кнопці ОК.

Якщо розв'язок знайдено, то у вікні Результат подбора параметра маємо відповідне повідомлення. У даному випадку

Розв'язок знайдено Підібране значення 0

Поточне значення –1,21293Е-05.

6.Клацаємо мишею на кнопці ОК. Значення x (0,50152732) матимемо в комірці А1, а відповідне йому значення формульного виразу (–1,21293E-05)

у В1.

7.Уведемо в комірку А1 інше початкове значення 1 і виконаємо дії з підбору параметра. Отримаємо значення x = 1,44003339, для якого значення формули 0,000341106.

Отримано два різних корені рівняння x1 ≈ 0,50152732 та x2 ≈ 1,44003339. Оскільки рівняння може мати неєдиний розв’язок, то підбирання параметра доцільно проводити декілька разів, починаючи з різних допустимих значень невідомої змінної.

Якщо розв’язку рівняння не існує, то у випадку спроби його визначити підбиранням параметра, у вікні Результат подбора параметра отримаємо повідомлення, що розв’язку не знайдено.

5.ЗАВДАННЯ ДЛЯ РОЗВ’ЯЗУВАННЯ НЕЛІНІЙНИХ АЛГЕБРАЇЧНИХ

РІВНЯНЬ

Знайти наближено корені рівняння ітераційним методом, а також використовуючи сервіс Excel “Підбір параметра”.

16

Варіанти рівнянь

 

 

Рівняння

варіанта

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

 

 

ln(x + 2) = x

2.

e

x

 

x

 

 

 

+0,8 = x

 

 

 

 

sin

2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

 

ln x +sin

3

x

 

 

 

 

=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

4.

 

 

ex2 / 2

= x

5.

 

(x3 +3)(ex ex ) =5

6.

 

arctg

x +0,5 = x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

cos(

x +2) +

 

x +3 =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.

 

arcsin x =

 

 

 

 

 

x

 

 

1 2

9

 

(2 +sin 4x)x =1

10

 

e( x 3) = x +2

11

 

2x + log2

 

2x 3

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

12

 

ex +1 ln

 

x +1,5

 

= 2

 

 

 

 

 

 

 

 

 

13

 

 

x2 2 =10x ln x

14

 

 

 

 

 

1

 

 

 

x +sin

 

= 2

 

 

 

 

 

 

x

 

15

 

 

x3 +10 cos x = 3

16

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

sin x ln(x +0,5) 2 = 0

 

17

 

 

ex = ex +2,5

18

 

arcsin(2x +1) = x2

19

 

 

(2x tg x)2 =1

20

 

 

ln x +1 = sin x

17

6.ЗМІСТ ЗВІТУ

1.Номер, тема і мета лабораторної роботи.

2.Короткі теоретичні відомості про методи розв’язування рівнянь.

3.Завдання свого варіанта роботи.

4.Розподіл даних і формул для обчислень у комірках ЕТ у вигляді таб-

лиці.

5.Послідовність дій при виконанні роботи на комп’ютері.

6.Результати розв’язку задачі, а також формули обчислення.

7.Короткі висновки.

7.КОНТРОЛЬНІ ЗАПИТАННЯ

1.Викладіть методику розв’язування систем лінійних алгебраїчних рівнянь.

2.Дайте поняття про матричні функції Excel, які потрібні для розв’язування систем лінійних алгебраїчних рівнянь.

3.Які ви знаєте способи розв’язування нелінійних алгебраїчних рівнянь?

4.Викладіть методику реалізації ітераційної розрахункової схеми в електронній таблиці для розв’язування нелінійних алгебраїчних рівнянь.

5.Викладіть методику проведення циклічних пере обчислень для розв’язування нелінійних алгебраїчних рівнянь.

6.Викладіть методику використання сервісу «Підбір параметра» для розв’язування нелінійних алгебраїчних рівнянь.

БІБЛІОГРАФІЧНИЙ СПИСОК

7.Гера Б.В., Риков В.А. Використання табличного процесора Excel для розв’язування інженерних задач : методичний посібник. – Дніпропет-

ровськ, 2007. – 87 с.

8.Гера Б.В., Риков В.А. Методичні вказівки з дисципліни “Обчислювальна техніка і програмування”. ч. 3. Використання табличного процесора Excel для розв’язування рівнянь. // Дніпропетр. нац. ун-т залізнич. трансп. ім. акад. В.Лазаряна; – Дніпропетровськ, 2007. – 11 с.

9.Глинський Я.М. Практикум з інформатики. – Львів: Деол, СПД Глин-

ський, 2005. – 296 с.

10.Вальдрат О.Л., Чаповська Р.Б. Робота з Microsoft Excel 2000. – Київ: ЦУЛ, Фітосоціоцентр, 2002. – 186 с.

11.Ларсен Р.У. Инженерные расчеты в Excel. – М.: Изд.дом “Вильямс”, 2002. – 544 с.

12.Обшта А.Ф., Каленюк П.І., Клочко Н.Ф., Гоблик Н.М. Excel for Windows9x: Лабораторний практикум. – Львів: Асоціація діагностичні системи, 2000. – 160 с.

13.Microsoft Office Excel 2003. Учебный курс. / В. Кузьмин. – СПб.: Пи-

тер; Киев: Издательская группа BHV, 2004. – 493 c.

18

Соседние файлы в предмете Вычислительная техника в инженерных и экономических расчетах