matlab 求函数的极值点和拐点

发布网友 发布时间：2022-04-22 05:24

我来回答

共2个回答

热心网友 时间：2022-06-01 20:06

展开1全部花了两天时间，终于研究明白了！

function hh
global dy1 dy2
y='x^2*sin(x^2-x-2)'
dy1=diff(y)
dy2=diff(y,2)
subplot(3,1,1)
ezplot(y,[-2 2])
subplot(3,1,2)
ezplot(dy1,[-2 2]),hold on,plot(-2:2,zeros(length(-2:2)))
subplot(3,1,3)
ezplot(dy2,[-2 2]),hold on,plot(-2:2,zeros(length(-2:2)))
x01=fsolve(@myfun1,[-1.5 -0.7 0 1.6])
x02=fsolve(@myfun2,[-1.9 -1.3 -0.5 1.3])
function f1=myfun1(x)
global dy1
f1=subs(dy1);%very inportamt!!!!!;
function f2=myfun2(x)
global dy2
f2=subs(dy2);%very inportamt!!!!!;

结果：
y =

x^2*sin(x^2-x-2)

dy1 =

2*x*sin(x^2-x-2)+x^2*cos(x^2-x-2)*(2*x-1)

dy2 =

2*sin(x^2-x-2)+4*x*cos(x^2-x-2)*(2*x-1)-x^2*sin(x^2-x-2)*(2*x-1)^2+2*x^2*cos(x^2-x-2)

Optimization terminated: first-order optimality is less than options.TolFun.

x01 =

-1.5326 -0.7315 0 1.5951

Optimization terminated: first-order optimality is less than options.TolFun.

x02 =

-1.9240 -1.2650 -0.4742 1.2404

参考资料：http://hi.baidu.com/chemical%5Fengineering/album/item/d6f7367ace36fefe2f73b3e1.html

热心网友 时间：2022-06-01 20:06

先求其导数，令导数的式子=0 得到的解便是极值点，如果超出了[-2,2]的范围，还要把-2与2比较。