花了两天时间,终于研究明白了!
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。
全部