matlab 求函数极值的代码

x=linspace(0,100);%x的取值范围
y=x.*sin(x);%函数
a=find(diff(sign(diff(y)))<0)+1;%判断是极大值
subplot(2,1,1);plot(x,y,x(a),y(a),'r*')%画出函数图
fprintf('下面是极大值x y组\n')
fprintf('x=%-13f',x(a))%输出
fprintf('\n')
fprintf('y=%-13f',y(a))
fprintf('\n')
fprintf('下面是极小值x y组\n')
b=find(diff(sign(diff(y)))>0)+1;%判断是极小值
subplot(2,1,2);plot(x,y,x(b),y(b),'r*')%
fprintf('x=%-13f',x(b))%输出
fprintf('\n')
fprintf('y=%-13f',y(b)