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

lufeideshijie

. B7 T7 V6 ]- g/ v& ]* `下降段数据提取（附原始数据）

1 ~  V% C3 m3 @: s1 rfilename='';
: B- r# r( d8 h4 W[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
! t) y% `: `: y; B. }6 G; ZX=data_Depth;
Y=data_Load;
plot(X,Y);
7 @4 r, T$ `  r& d, ?# y; ky=Y';
dy=[-1,diff(y)]>0;
figure
) ~% W3 B& T( Wplot(dy);
%% 上升
k=diff([0 dy 0]);
ind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1

' J: x' I( [4 M. ?& r8 Y( F9 ofor ii=1:numel(ind);
    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
& K1 a5 \# L% j- u/ mend

%% 下降
k=diff([0 ~dy 0]);
+ h6 `! _' o, @1 l4 b3 g$ mind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1

for ii=1:numel(ind);
    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end

: T2 ~+ T1 b' c$ b1 ?
运算结果有1723个，我想获取20个，即plot（X,Y）中那20个下降段，并输出该段的x,y，请教大神，非常感谢！

panphp

% 读取数据
: Z. \* \2 `" P" {# Afilename='test.txt';
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);

% 平滑操作
; m9 z1 z8 C; w  m% \  k: F+ P" JN=20;
X=smooth(data_Depth,N);
Y=smooth(data_Load,N);
$ J6 t' l- Z* a$ w( u
; ~5 |& a: Y% X( B' K8 p% 找出下降段索引
  }0 e) c" q. L# wInd=(sign(diff(X))<0 & sign(diff(Y))<0);
: r1 r& r8 F9 U+ I0 k/ k! z; R4 p
% 画图
figure
plot(data_Depth,data_Load)
hold on
plot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

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