有关下降段数据提取问题请教大神。

lufeideshijie

1 ~. ~; d/ L  ^# [( w5 Q; a下降段数据提取（附原始数据）
+ L; K9 k( Y  n$ X& c
filename='';
5 R5 @3 U3 F" T( x1 W$ {* u( H+ @[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
! T% t+ W) q$ Z, s2 cX=data_Depth;
Y=data_Load;
plot(X,Y);
y=Y';
dy=[-1,diff(y)]>0;
figure
plot(dy);
%% 上升
k=diff([0 dy 0]);
ind=find(k==1);
7 F+ Q  p4 T4 ?( Jnum=find(k==-1)-ind;   %寻找连续1

/ T' s2 S# X' f, ?for ii=1:numel(ind);
; {* {7 t9 ?# t$ ~& s# g1 k    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
+ y3 F" O' f2 h/ K7 N" k: w4 Tend
6 l7 q6 X: f' \9 ?' Q
( @+ H% f6 o% H! w- T& \5 m5 B%% 下降
k=diff([0 ~dy 0]);
ind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1
4 w. x1 G# P( p0 f9 T7 m
for ii=1:numel(ind);
; }/ G  K4 l8 J  V    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
& s& }& V: b' s& {  B- oend
6 a3 d% U& P0 N

运算结果有1723个，我想获取20个，即plot（X,Y）中那20个下降段，并输出该段的x,y，请教大神，非常感谢！

panphp

% 读取数据
filename='test.txt';
9 X5 W5 e9 s% Q. Z6 w[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);

% 平滑操作
N=20;
X=smooth(data_Depth,N);
1 k4 L2 W. H1 \7 @& kY=smooth(data_Load,N);

% y9 @9 g2 @' |' d, l5 c% 找出下降段索引
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);

2 p; g4 u- u% h, ]+ u2 u% 画图
figure
plot(data_Depth,data_Load)
. O4 m* n0 Q) c( Whold on
plot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

你提取的数据有误吧，我看down中不只有单个的数值，还有矩阵！