clc; close all; clear;
t = linspace(0,pi,200)';
d = sin(2*t);
% 微分变号信息
IndMin = find(diff(sign(diff(d)))>0)+1; % 获得局部最小值的位置
IndMax = find(diff(sign(diff(d)))<0)+1; % 获得局部最大值的位置
figure; hold on; box on;
plot(t,d);
plot(t(IndMax), d(IndMax), 'k*')
plot(t(IndMin), d(IndMin), 'r^')
title('计算离散节点的波峰波谷信息', 'FontWeight', 'Bold');