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

lufeideshijie

: s3 [0 z. z3 G/ d2 A8 r4 |下降段数据提取（附原始数据）

filename='';
: a! c2 m. S% n[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
1 Y& v- T3 U# W/ ?$ ^X=data_Depth;
Y=data_Load;
# k! h  I% ^3 X2 t4 |plot(X,Y);
y=Y';
dy=[-1,diff(y)]>0;
figure
4 }  h+ v4 P0 q" Tplot(dy);
%% 上升
9 w8 [7 F  z  Sk=diff([0 dy 0]);
ind=find(k==1);
3 y' u' A% w$ z& z, j7 |num=find(k==-1)-ind;   %寻找连续1

; I+ P' h4 L5 ?! N3 M/ Jfor ii=1:numel(ind);
# a$ ?- a; w. ^- D8 ]    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
. v7 t- K9 h6 q7 uend

%% 下降
* ^" K% g1 G7 G5 tk=diff([0 ~dy 0]);
1 X  d, v: z7 N- aind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1

" C  K+ |, P. R+ J/ \- j4 }& bfor ii=1:numel(ind);
3 Q8 K. C  j1 w& ~( h  R& O    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end
1 \7 h* t; M6 D( e; ~5 ~0 n* ?2 S
: i: j! W  c7 a: {
运算结果有1723个，我想获取20个，即plot（X,Y）中那20个下降段，并输出该段的x,y，请教大神，非常感谢！

panphp

% 读取数据
filename='test.txt';
6 h+ c2 N! `. l4 }& b# z4 K+ E[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);

# w7 t1 N* A% @4 g- H% 平滑操作
N=20;
, b( _! g0 `# K8 f5 TX=smooth(data_Depth,N);
Y=smooth(data_Load,N);
2 ~2 G/ N% b+ l: y
. V. V7 Q4 _- ]- V6 B. b% 找出下降段索引
: k" D5 t4 W: Y- f$ r' v$ }2 ~, bInd=(sign(diff(X))<0 & sign(diff(Y))<0);
% I: Y4 b  W! i0 B
% 画图
8 O+ m3 x' L4 h0 u8 |6 `figure
plot(data_Depth,data_Load)
hold on
: n/ y; a) Z1 B1 fplot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

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