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

lufeideshijie

& j3 G5 r( l! H- M: e下降段数据提取（附原始数据）
* y0 ]# f( A+ L3 a; s
' I8 W" q, y4 W7 b' v: g9 [filename='';
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
2 D+ T: O4 e7 @; e! {" u3 Y7 F8 uX=data_Depth;
Y=data_Load;
plot(X,Y);
y=Y';
dy=[-1,diff(y)]>0;
figure
plot(dy);
0 Z8 k! o. ]$ Y; O& d; n) t%% 上升
k=diff([0 dy 0]);
ind=find(k==1);
2 T8 C, J& r2 {* ~9 _num=find(k==-1)-ind;   %寻找连续1
7 I. K1 h! E% R1 K
; V% d+ v) P+ m* u0 @# n8 ofor ii=1:numel(ind);
    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
end

%% 下降
k=diff([0 ~dy 0]);
ind=find(k==1);
num=find(k==-1)-ind;   %寻找连续1

! T( t( I: w& Z* z+ }for ii=1:numel(ind);
6 ?. E# u: S6 c- j; C    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
- v. G' _1 R* w5 T* r0 e) }- C$ \end
& Y% k9 a0 f* o4 t% y, p5 I

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

panphp

% 读取数据
filename='test.txt';
/ q3 e5 K! x, X3 p2 r[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);

) y: o% h# w5 J% X/ v5 G0 ?" l% 平滑操作
N=20;
X=smooth(data_Depth,N);
Y=smooth(data_Load,N);
! l" K) s' g$ H) }+ D/ Q
% 找出下降段索引
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);
; H) L0 X* ^4 E  _2 k/ N
' f/ t* d: A6 [" ]% 画图
figure
plot(data_Depth,data_Load)
hold on
plot(data_Depth(Ind),data_Load(Ind),'r.')

Colbie

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