Рис. 1.25. Решение уравнения при различных интервалах изоляции корня

Содержимое функции 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

Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: