exn519b.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 14 行

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第五章例5-1-9b程序
% 求曲线的极值（用导数为零方法）
% 西安电子科技大学出版社出版  陈怀琛编著 2007年7月
%
clear,clf,
x=linspace(0,pi,10000);
y=2*sin(2*x).^2+5/2*x.*cos(x/2).^2; 
plot(x,y),grid on,hold on
h=pi/(10000-1);Dy=diff(y)/h;	% 先用数值方法求导数
plot(x(1:end-1),Dy, '-.r')		% 绘制导数曲线
kh=input('过零点检测门限取h的kh倍,kh= ')
nm=find(abs(Dy)<kh*h)+1, ym=y(nm)   % 求导数过零点， 
nm=find(abs(Dy)<8*h), ym=y(nm),xm=x(nm)
set(gcf,'color','w')