Page 201 - 6197
P. 201
Для знаходження коренів полінома використаємо
вбудовану в MatLab функцію roots(p), де p – вектор-
аргумент, компоненти якого коефіцієнти полінома. У нашому
випадку
>> p=[133 54 -36 10 -1];
>> y=roots(p)
y =
-0.8445
0.1379 + 0.1890i
0.1379 - 0.1890i
0.1627
Оскільки розв’язували алгебраїчне рівняння четвертого
порядку, то отримали чотири корені. Серед них вибираємо
той, який задовольняє умові задачі - y 0 . Отже,
4
*
y 0 1627, .
4
1 3y *
*
*
Знаючи y , знаходимо: y 4 0,2559 ,
4 1
2
1 y *
*
*
*
y 4 0,4186, y y 0 1627, .
2 3 4
2
Знаходимо значення Z , j 1 4, за формулою (3.79). Для
j
. У відповідності із
цього спочатку обчислимо R x *
співвідношення (3.81) маємо
0 2559, 0 4186, 0 1627, 0 1627,
5 2 5 1 .
*
R x 9,6606
0 2559, 0 4186, 0 1627, 0 1627,
За формулою (3.79) знаходимо, що
Z Z R x * 9 6606 0 1627 1,5718, , .
y
3 4 3
Тепер складемо систему рівнянь
5x 1,5718,
1
1
x 1,5718 ,
2
201