找地震动记录数据中的极大值和极小值

clc,clear,close all;
fs=1/0.01 ;%%%%%%%时间间隔的倒数
tspan=52.78 ;%%%%%%地震波记录总持时
t=0:1/fs:tspan;%起始数值，时间间隔，总持时
x=importdata(‘TCU045.txt’);%读入txt.格式的地震波>>>>>>>>>>输入
figure
plot(x,‘k-’)
[Pks, Locs] = findpeaks(x);% 极大值Pks, 对应下标Locs
[Pks2, Locs2] = findpeaks(-x);
Pks2 = -Pks2;% 极小值Pks2, 对应下标Locs2
hold on
plot(Locs,Pks,‘r*’);%正极大值
hold on
plot(Locs2,Pks2,‘bp’);负极大值

绘图结果：

2021.4.7.17:19