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

lufeideshijie

下降段数据提取（附原始数据）

; `0 d( A. t2 M* g* Ofilename='';
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
X=data_Depth;
3 J$ |2 V/ q7 f5 M( b/ eY=data_Load;
plot(X,Y);
y=Y';
% C. l4 Y: O6 {9 P, a9 g& n! Kdy=[-1,diff(y)]>0;
& L6 h4 R6 ?# }figure
8 p! T3 C* I: k# I) L+ S' Oplot(dy);
%% 上升
! c, X: _  e( r! Fk=diff([0 dy 0]);
ind=find(k==1);
0 U  z+ @% `8 O' R& A; Cnum=find(k==-1)-ind;   %寻找连续1

, u9 ~8 S1 ?" E" w( z1 tfor ii=1:numel(ind);
% \( V+ k% |6 i( T0 ?) {    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end

% `8 _9 c1 ]4 ], v* _( R5 k: N7 W%% 下降
0 W+ a/ ]% p7 n& h) K. \k=diff([0 ~dy 0]);
ind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1
) \' f" h# H5 n5 w; _% B* I
6 n  \5 i; K# X8 y5 X! U7 jfor ii=1:numel(ind);
& t3 K, p5 x! D- Y    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end


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

panphp

% 读取数据
filename='test.txt';
[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);
8 |$ ?- u4 Y$ P4 h: oY=smooth(data_Load,N);

; M( _' @) _3 r' P9 }# s% 找出下降段索引
  T+ ?* N# Q& hInd=(sign(diff(X))<0 & sign(diff(Y))<0);
- T) J' f: E+ _$ r8 {& ?* ?
% 画图
figure
plot(data_Depth,data_Load)
hold on
4 L0 r% ^+ {& a4 |* Eplot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

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