2 0.777373 -3.32063 Search
3 0.822627 -4.06625 search
4 0.768 -3.17712 search
5 0.832 -4.23184 search
6 0.754745 -2.98039 search
7 0.845255 -4.47271 search
8 0.736 -2.71444 search
9 0.864 -4.82695 search
10 0.70949 -2.36229 search
11 0.89051 -5.35561 search
12 0.672 -1.9106 search
13 0.928 -6.16014 search
14 0.618981 -1.35979 search
15 0.981019 -7.41582 search
16 0.544 -0.743367 search
17 1.056 -9.43876 search
18 0.437961 -0.157497 search
19 1.16204 -12.8248 search
20 0.288 0.215057 search
Looking for a zero in the interval [0.288, 1.162]
21 0.308 0.190464 interpolation
22 0.460543 -0.256865 interpolation
23 0.37295 0.0607713 interpolation
24 0.397562 -0.010605 interpolation
Zero found in the interval: [0.288, 1.162].
x =
0.3976
>> [x fval]=fzero('Func',0.8,optimset('TolX',10^-5,'disp','iter'))
Func-count x f(x) Procedure
1 0.8 -3.6824 initial
2 0.777373 -3.32063 search
3 0.822627 -4.06625 search
4 0.768 -3.17712 search
5 0.832 -4.23184 search
6 0.754745 -2.98039 search
7 0.845255 -4.47271 search
8 0.736 -2.71444 search
9 0.864 -4.82695 search
10 0.70949 -2.36229 search
11 0.89051 -5.35561 search
12 0.672 -1.9106 search
13 0.928 -6.16014 search
14 0.618981 -1.35979 search
15 0.981019 -7.41582 search
16 0.544 -0.743367 search
17 1.056 -9.43876 search
18 0.437961 -0.157497 search
19 1.16204 -12.8248 search
20 0.288 0.215057 search
Looking for a zero in the interval [0.288, 1.162]
21 0.302415 0.198003 interpolation
22 0.467237 -0.288821 interpolation
23 0.369452 0.0698667 interpolation
24 0.398884 -0.0148167 interpolation
25 0.393735 0.00136079 interpolation
26 0.394168 2.27096e-005 interpolation
27 0.394188 -3.91718e-005 interpolation
Zero found in the interval: [0.288, 1.162].
x =
0.3942
fval =
2.2710e-005
>> X = fzero('Func',3, optimset('disp','final'))
Zero found in the interval: [0.28471, 4.92].
X =
0.3942
Для нахождения корней полинома в пакете MATLAB предусмотрена соответствующая функция roots( ), возвращающая вектор-столбец, компоненты которого являются корнями полинома (действительными или комплексными).
Обращение к функции к данной функции имеет следующий вид:
r = roots(c)
Здесь
c
вектор-строка, содержащая значения
коэффициентов полинома
.
Пример 2.4.
Решение уравнения
с использованием функцииroots(
).
>> c=[1 -11 1 1 0.1];
>> roots(c)
ans =
10.8998
0.3942
-0.1470 + 0.0409i
-0.1470 - 0.0409i