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

lufeideshijie

/ @0 B% E  E( T( _" Z% S/ b下降段数据提取（附原始数据）

8 M  q: C) {) T0 x# _9 Q) ]6 Yfilename='';
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
# S' N4 K. f( j' aX=data_Depth;
1 D# a3 {) q/ E' B2 i  Y5 SY=data_Load;
7 u! }  W# X  A- \4 Pplot(X,Y);
7 p$ ~$ l7 I& e1 Xy=Y';
. M( ~" y/ I7 N7 `* {, }) Ndy=[-1,diff(y)]>0;
figure
" @, e0 h% n. S( R/ |/ k. |6 o! Fplot(dy);
' `" |8 Z* p4 m0 i%% 上升
7 v4 Q6 x1 Y6 k1 H8 nk=diff([0 dy 0]);
ind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1
* I6 @9 p0 ^( i) ~
/ Q4 c0 ~% B- x/ g* t, V& n8 k6 B5 ~for ii=1:numel(ind);
# ^% M$ U. M' z" Z' J  B& Y" `    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end
, S- R3 e5 ], `; R, b
%% 下降
+ B7 j% t8 p5 p. Q+ e# \7 P6 nk=diff([0 ~dy 0]);
; W: X8 c4 p9 Q; l' Gind=find(k==1);
6 B& ^5 n) `3 _0 i# ~& [3 Wnum=find(k==-1)-ind;   %寻找连续1
) {- e- T) l1 e. F4 n) s
for ii=1:numel(ind);
( ?* _* e( u) _9 _' O    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
3 i, M( \& `0 z" lend

1 e& n) [2 s8 j6 c. ]
运算结果有1723个，我想获取20个，即plot（X,Y）中那20个下降段，并输出该段的x,y，请教大神，非常感谢！

panphp

% 读取数据
# |" W; k8 [0 g8 E+ {2 Mfilename='test.txt';
, @! q5 f% Z, M2 e: r[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);

% 平滑操作
( \( ], @2 Z) _" x3 `( X/ CN=20;
X=smooth(data_Depth,N);
Y=smooth(data_Load,N);

- s) e0 m& b! I# M0 X% 找出下降段索引
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);

1 p$ u" Z7 }' p' M$ ^' k* C% 画图
figure
plot(data_Depth,data_Load)
6 r% P; A$ m1 k- X( uhold on
plot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

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