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

lufeideshijie

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

filename='';
+ P% G- p' I9 t" d3 N6 w/ b[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
$ u# r6 }6 Y, c9 O" A$ x! HX=data_Depth;
( @! v! ?. E; m1 k" SY=data_Load;
: h+ t( [2 d) Nplot(X,Y);
y=Y';
dy=[-1,diff(y)]>0;
figure
: ?+ a6 l( r# u* K  a* s  W* hplot(dy);
& `* _, i" U( P- G4 j% U4 z%% 上升
k=diff([0 dy 0]);
+ F8 c3 k& b* c1 w& d: y8 D) ?8 P) Nind=find(k==1);
# V- [8 H- m0 D. _: |* y6 g1 gnum=find(k==-1)-ind;   %寻找连续1
5 Z$ U7 g" W5 y; m6 h7 y
for ii=1:numel(ind);
( `/ ]5 f: N4 K$ c' d5 Q    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end

* R, u( N  v2 w8 J% d! U%% 下降
$ f* n2 I. A" t6 \/ ]: u1 Ak=diff([0 ~dy 0]);
. @+ J. d7 f8 k/ x3 W* eind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1
. M" f# e: Z4 }. Y/ ?( C& {& B
for ii=1:numel(ind);
    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end
) K" P' }! M. }/ q

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

panphp

% 读取数据
) t# M% V( s# i  X( hfilename='test.txt';
$ Q- Z# j) k2 ][data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
( ?2 n* Z# F& ?5 \* m3 @  w
+ U5 F' a2 C( p" s1 Y% 平滑操作
6 O' C9 I3 `5 IN=20;
) i" q  I. Q% A/ m" Q* o0 uX=smooth(data_Depth,N);
" g1 X8 O5 V0 Y4 JY=smooth(data_Load,N);
+ ~6 x" p( v$ h9 I
% 找出下降段索引
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);

, V$ m8 u6 K) Q/ Z: J9 r6 Y1 d+ C% 画图
' h& z0 O7 y3 F5 ^) z' Pfigure
plot(data_Depth,data_Load)
' K0 f$ d4 M0 w' C: ?1 P) Thold on
plot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

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