Содержимое функции func продемонстрировано на рис.1.26.
Рис. 1.26. Функция func(x)
Рекомендуется по возможности построить график. Потом уже описать функцию в М-файле. Затем использовать функцию fzero для нахождения корня.
Если график (рис.1.27) покажет, что уравнение имеет несколько корней, то следует для каждого из корней определить интервал изоляции и использовать его в вызове функции fzero (для каждого из корней по отдельности), то есть к функции придется обратиться несколько раз.
При вызове функции fzero возвращаются только вещественные числа. Комплексные (мнимые) корни она не вычисляет.
function Start_exFzero1 function y=exFzero_1(x) y=((x-1).^2).^(1/3)-(x.^2).^(1/3); end x=[-2:0.1:2]; y=((x-1).^2).^(1/3)-(x.^2).^(1/3); plot(x,y,'r-'); grid on; fzero(‘exFzero_1’,[0 1] end | |
В Command Window >> ans= 0.5000 |
Рис.1.27. График трансцендентного уравнения и решение уравнения функцией fzero.
С помощью функции polyval можно получить графическое (рис.1.28) решение уравнения.
function exPolyval p=[1 0.4 0.6 -1.6]; x=[-1:0.1:3]; y=polyval(p,x); plot(x,y,'k-'); grid on; end |