Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ms_labs

.pdf
Скачиваний:
25
Добавлен:
10.12.2018
Размер:
4.28 Mб
Скачать

а другий індекс j – номер невідомого, при якому знаходиться даний коефіцієнт. Наприклад, коефіцієнт a23 знаходиться в другому рівнянні системи при невідомому x3 .

Числа b1, b2 , , bm називаються вільними членами системи. Коротко систему (1) можна записати так:

n

 

 

aij x j bi ,

( i 1, 2, , m ).

 

j1

 

 

Розв’язком лінійної системи (1) називається будь-яка сукупність чисел

1, 2 , , n , яка, будучи підставлена замість невідомих

x1, x2 , , xn в

рівняння даної системи, перетворює усі рівняння системи в тотожності. Система лінійних рівнянь (1) називається сумісною, якщо вона має

розв’язок. Якщо система лінійних рівнянь немає розв’язку, то вона називається

несумісною (або суперечливою).

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

Дві системи лінійних рівнянь з одним і тим же числом невідомих називаються еквівалентними, якщо вони або обидві несумісні, або обидві сумісні і мають одні і ті ж розв’язки.

Наведені нижче три типи перетворень називаються елементарними перетвореннями системи лінійних рівнянь. Це такі:

1)перестановка двох рівнянь системи;

2)множення обох частин рівняння системи на будь-яке відмінне від нуля

число;

3)додавання до обох частин рівняння одного рівняння відповідних частин другого рівняння, помножених на будь-яке число.

Виконання елементарних перетворень рівносильне вираженню однієї змінної через інші.

Для знаходження розв’язків лінійної системи рівнянь існує багато методів. Найбільш поширеними серед них є метод Крамера, метод Гауса та матричний метод.

Мета роботи – знаходження засобами Ms Excel розв’язків системи лінійних рівнянь.

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

41

лінійних алгебраїчних рівнянь. Крім того. Завданням може бути самостійне складання системи рівнянь яка має розв’язок. Для цього потрібно скласти самому визначник, який не дорівнює нулю та довільно задати вільні члени. Якщо система не розв’язуватиметься очевидно треба буде змінити значення вільних членів або коефіцієнтів при змінних.

Метод Крамера.

Нехай дано систему лінійних рівнянь, в якій число рівнянь рівне числу невідомих:

a11x1 a12x2

a1n xn b1,

 

 

 

 

 

 

b2,

 

a21x1 a22x2 a2n xn

(2)

 

 

 

 

 

 

 

 

 

 

 

 

a

x a

x a x

b .

 

 

n1 1

n2 2

nn

n

n

 

Подамо складові даної лінійної системи в такий спосіб

a11

a12

a1n

 

 

 

 

 

A a21

a22

a2n

,

 

 

 

 

 

an2

 

 

an1

ann

 

x1

 

 

 

 

 

 

x x2

 

,

 

 

 

 

xn

 

 

 

 

 

b1

 

 

 

 

 

b b2

 

,

 

 

 

 

 

bn

 

 

де A – матриця системи, x – стовпчик невідомих, b – стовпчик вільних членів.

Припустимо, що визначник системи не дорівнює нулю, тобто:

a11

a12

a1n

 

 

 

 

 

d a21

a22

a2n

0 .

 

 

 

 

 

an2

 

 

an1

ann

 

Якщо тепер замінити послідовно у визначнику d стовпчики коефіцієнтів при невідомих x j ( j 1, 2, , n ) стовпчиком вільних членів bi , то отримаємо

відповідно такі визначники:

42

 

 

b1

a12

a1n

 

 

a11

b1

a1n

 

 

 

 

 

 

 

 

 

 

d

1

b2

a22

a2n ;

d

2

a21

b2

a2n ;

 

 

 

 

 

 

 

 

 

 

 

 

an2

 

 

 

 

 

 

 

 

bn

ann

 

 

an1

bn ann

 

 

 

a11

a12

b1

a1n

 

 

 

 

a11

a12 a1,n1

b1

 

 

 

 

 

a22

 

 

 

 

 

 

 

a22 a2,n 1

 

 

d

3

 

a21

b2 a2n

;

d

n

 

a21

b2

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an2

 

 

 

 

 

 

 

an2 an,n1

 

 

 

 

 

an1

bn ann

 

 

 

 

an1

bn

 

Теорема Крамера. Система n лінійних рівнянь з n невідомими,

визначник якої відмінний від нуля,

завжди сумісна і має єдиний розв’язок,

який обчислюється за формулами:

x

d1

;

x

 

 

d2

; …

x

n 1

 

dn 1

;

x

 

 

dn

.

(3)

2

 

n

 

1

d

 

 

 

d

 

 

 

d

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формули (3) називаються формулами Крамера.

Приклад. Розв’язати засобами Ms Excel систему лінійних алгебраїчних рівнянь методом Крамера. Система має такий вигляд:

 

x1

 

x2

2x3

1,

 

2x1

 

x2

2x3

4,

 

 

4x

x

2

4x

3

2.

 

1

 

 

 

 

1.

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

Використовуючи матричну операцію «МОПРЕД», що знаходиться в категорії функцій «Математичні», виділяємо масив (визначника) B4:D6 і натискаємо ОК. В комірці D8 отримуємо значення d 6 , що й зображено на скриншоті справа від визначника системи.

 

1

1

2

d

2

1

2

 

4

1

4

 

 

 

 

43

2. Аналогічно обчислюємо визначники d1, d2 , d3 .

 

1

1

2

 

d1

4

1

2

 

2

1

4

 

 

 

 

 

1

2

 

1

d 2

2

4

2

 

4

2

4

 

 

 

 

1 1 1 d3 2 1 4

41 2

3.Остаточно, використовуючи формули Крамера (3) знаходимо розв’язки даної лінійної системи алгебраїчних рівнянь:

x

6

1,

x

 

 

12

2 ,

x

 

12

2 .

 

2

 

3

1

6

 

 

6

 

 

6

 

 

 

 

 

 

 

 

 

Таким чином, практично безпосереднє отримання значення визначника уможливлює розв’язок системи n лінійних алгебраїчних рівнянь з n невідомими.

Матричний метод.

I. Система n лінійних рівнянь з n невідомими.

Систему (2) можна записати у вигляді матричного рівняння:

A X B

(3)

В розгорнутому вигляді систему (3) можна подати так

44

a11

a21

an1

a12 a1n

a22 a2n

an2 ann

 

x1

 

b1

 

 

 

 

 

 

 

 

x2

 

b2

.

 

 

 

 

 

 

xn

 

 

 

 

 

 

bn

 

Розглянемо метод розв’язку системи (2) в загальному вигляді, який часом називають методом оберненої матриці. Вважаємо, що квадратна матриця Ann є не виродженою, тобто її визначник A 0 . В цьому випадку існує

обернена матриця A 1 . Помноживши зліва обидві частини матричної рівності

(3) на обернену матрицю A 1 , отримаємо:

A 1 A X A 1 B E X A 1 B E X X ,

де E – одинична матриця.

Звідси розв’язок системи методом оберненої матриці є таким

 

X A 1 B .

 

(4)

Таким чином, для розв’язку системи (2),

тобто знаходження вектора

X

необхідно знайти обернену матрицю коефіцієнтів A 1 і помножити її справа на вектор вільних членів B .

Приклад 1. Розв’язати систему

3x 2y 7,4x 5y 40.

Розв’язок.

1. Ввести матрицю A (в даному випадку розміром 2 2) в діапазон А1:В2

 

 

3

2

 

 

A

4

5

.

 

 

 

Вектор B 7

40 введіть в діапазон С1:С2

2. Знайдіть обернену матрицю A 1 . Для цього:

виділіть блок комірок під обернену матрицю. Наприклад, виділіть блок

А3:В4

(вказівником мишки при натиснутій лівій клавіші);

натисніть на панелі інструментів Стандартна кнопку Вставка функції;

в отриманому діалоговому вікні Майстер функцій в робочому полі

Категорія

виберіть пункт Математичні, а в робочому полі Функція — ім’я

45

функції МОБР. Клацніть кнопкою ОК;

діалогове вікно МОБР, яке з’явилося мишкою відсуньте в сторону від

початкової

матриці і введіть діапазон вихідної матриці А1:В2 в

робоче поле

Масив (вказівником мишки при натиснутій лівій

клавіші). Натисніть сполучення клавіш СТRL+SНIFТ+ЕNТЕR;

якщо обернена матриця не з’явилась в діапазоні А3:В4, то треба

клацнути

вкзівником мишки в рядку формул и повторити

натискання

сполучення

клавіш

СТRL+SНIFТ+ЕNТЕR.

В результаті в діапазоні А3:В4 з’явиться обернена матриця

 

 

0.217391

0.086957

 

 

 

0.173913

0.13043

 

 

 

 

3. Множенням оберненої матриці A 1 на вектор B знаходять вектор X . Для цього:

виділяють блок комірок під результуючу матрицю (під вектор X ). Її розмірність буде m p , в даному прикладі 2 1. Наприклад, виділіть блок комірок С3:С4 (вказівником мишки при натиснутій лівій клавіші);

натискають на панелі інструментів Стандартна кнопку Вставка функції;

в отриманому діалоговому вікні Майстер функцій в робочому полі

Категорія виберіть пункт Математичні, а в робочому полі Функція

— ім’я функції МуМнож. Клацнути кнопкою ОК;

Діалогове вікно МуМнож мишкою відсуньте в сторону від вихідних матриць і введіть діапазон оберненої матриці A 1 , тобто діапазон А3:В4 в робоче поле Массив1 (вказівником мишки при натиснутій лівій клавіші), а діапазон матриці B тобто С1:С2 введіть в робоче поле Массив2. Після чого натиснути сполучення клавіш

СТRL+SНIFТ+ЕNТЕR.

якщо вектор X не з’явиться в діапазоні С3:С4, треба клацнути вказівником мишки в рядку формул и повторити натискання

СТRL+SНIFТ+ЕNТЕR.

Врезультаті в діапазоні С3:С4 з’явиться вектор X . Причому x 5 буде знаходитися в комірці С3, а y 4 – в комірці С4.

46

Можна здійснити перевірку знайденого розв’язку. Для цього знайдений вектор X необхідно підставити в вихідне матричне рівняння A X B. Перевірка здійснюється в такий спосіб.

1. Виділити блок комірок під результуючу матрицю (під вектор B ). Її розмірність буде m p , в даному прикладі 2 1. Наприклад, виділіть блок

комірок D1:D2 (вказівником мишки при натиснутій лівій клавіші).

2.Натисніть на панелі інструментів Стандартна кнопку Вставка функции.

3.В діалоговому вікні Майстер функцій на робочому полі Категорія виберіть Математичні, а на робочому полі Функція – ім’я функції МУМНОЖ. Клацніть кнопкою ОК.

4.Діалогове вікно МУМНОЖ мишкою відсуньте в бік від існуючих матриць і введіть діапазон вихідної матриці A А1:В2 в робочому полі Масив1 (вказівником мишки при натиснутій лівій кнопці), а діапазон матриці

X СЗ:С4 – на робочому полі Масив2. Після цього натиснути сполучення клавіш СТRL+SНIFТ+ЕNТЕR.

5. Якщо вектор B не з’явився в діапазоні D1:D2, то треба клацнути вказівником мишки в рядку формул і повторити натискання

СТRL+SНIFТ+ЕNТЕR.

В результате в диапазоне D1:D2 появится вектор B , і, якщо система

розв’язана правильно, то вектор,

який

появиться буде рівний вихідному

B 7

40

.

 

 

 

 

 

Приклад 2. Розв’язати систему

 

 

 

 

 

 

 

5x1 3x2

4x3

2700,

 

 

 

2x1

x2

 

x3

900,

 

 

 

 

 

 

3x

2x

2x

1600.

 

 

 

1

2

 

3

 

1.

Ввести матрицю A (в даному випадку розміром 3 3) в діапазон А1:С3

 

 

 

5

3

4

 

 

 

A

2

1

1

 

 

 

 

 

 

 

 

 

 

 

3

2

2

 

 

 

 

 

 

 

 

Вектор B 2700

900 1600 ввеcти в діапазон D1:D3.

2.

Знайдіть обернену матрицю A 1 . Для цього:

47

виділити блок комірок під обернену матрицю. Наприклад, виділити

блок

А4:С6

(вказівником мишки при натиснутій лівій клавіші);

 

натисніть на панелі інструментів Стандартна кнопку Вставка функції;

в отриманому діалоговому вікні Майстер функцій в робочому полі

Категорія

виберіть пункт Математичні, а в робочому полі Функція — ім’я функції МОБР. Клацніть кнопкою ОК;

діалогове вікно МОБР, яке з’явилося мишкою відсуньте в бік від

початкової матриці і введіть діапазон вихідної матриці А1:С3 в робоче поле Масив (вказівником мишки при натиснутій лівій клавіші). Натисніть сполучення клавіш СТRL+SНIFТ+ЕNТЕR;

якщо обернена матриця не з’явилась в діапазоні А4:С6, то треба клацнути вкзівником мишки в рядку формул и повторити натискання сполучення клавіш СТRL+SНIFТ+ЕNТЕR.

В результаті в діапазоні А4:С6 з’явиться обернена матриця.

3. Множенням оберненої матриці A 1 на вектор B знаходять вектор X . Для цього:

виділяють блок комірок під результуючу матрицю (під вектор X ). Її розмірність буде m p , в даному прикладі 3 1. Наприклад, виділіть блок комірок D4:D6 (вказівником мишки при натиснутій лівій клавіші);

натисніть на панелі інструментів Стандартна кнопку Вставка функції;

в отриманому діалоговому вікні Майстер функцій в робочому полі

Категорія виберіть пункт Математичні, а в робочому полі Функція

— ім’я функції МуМнож. Клацнути кнопкою ОК;

Діалогове вікно МуМнож мишкою відсуньте в сторону від вихідних матриць і введіть діапазон оберненої матриці A 1 , тобто діапазон А4:С6 в робоче поле Массив1 (вказівником мишки при натиснутій лівій клавіші), а діапазон матриці B тобто D1:D3 введіть в робоче поле Массив2. Після чого натиснути сполучення клавіш

СТRL+SНIFТ+ЕNТЕR.

48

якщо вектор X не з’явиться в діапазоні D4:D6, треба клацнути вказівником мишки в рядку формул и повторити натискання

СТRL+SНIFТ+ЕNТЕR.

Врезультаті в діапазоні D4:D6 з’явиться вектор X . Причому x1 200

буде знаходитися в комірці D4, а x2 300 – в комірці D5, а x3 200 – в

комірці D6.

Рекомендується зробити перевірку, підставивши знайдені значення в рівняння системи. Для цього знайдений вектор X треба підставити в початкове матричне рівняння A X B.

Приклад 3. Знайти розв’язки системи

2x1 3x2 4x3 5x4 54x1 5x2 2x3 7x4 23x1 5x2 6x3 4x4 72x1 6x2 6x3 5x4 2

1. Ввести значення елементів матриць A і B рівняння в комірки Excel.

2. Обчислити обернену матрицю за допомогою матричної функції МОБР.

49

3. Перемножити обернену матрицю A 1 на матрицю B з допомогою матричної функції МУМНОЖ (Порядок множення є важливим – першою повинна бути матриця A 1 , а другою матриця B )

4. Перевірити правильність отриманої матриці розв’язків (коренів) X .

II. Система m лінійних рівнянь з n невідомими.

Система т лінійних рівнянь з п невідомими має такий вигляд:

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]