EDA365电子论坛网
标题:
有关下降段数据提取问题请教大神。
[打印本页]
作者:
lufeideshijie
时间:
2020-6-10 10:53
标题:
有关下降段数据提取问题请教大神。
8 K8 a. \) H% v9 M6 R- o" }
下降段数据提取（附原始数据）
$ n n$ X7 o) r' U$ N7 O# \9 d4 W
4 k1 ?' R! O3 @+ V4 b% H7 ?8 Z
filename='';
7 [& O2 u/ h6 K$ q
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
4 {0 C, }) s! i$ t0 i
X=data_Depth;
. d, D9 W1 X# K
Y=data_Load;
( Q E: t. C* i) @0 r
plot(X,Y);
1 I8 j, c& w9 A8 C
y=Y';
/ u+ m% Y) e3 w7 T
dy=[-1,diff(y)]>0;
P3 l5 t/ u* S, d# r: f
figure
2 Y5 a4 F0 P q+ j& O. N
plot(dy);
2 Q- \7 N" o2 K9 }
%% 上升
/ z6 D8 X$ z/ R- |( g
k=diff([0 dy 0]);
9 s1 ^' R& s$ q ~( V
ind=find(k==1);
1 _' \. m- Z) C" T/ S8 t0 ~
num=find(k==-1)-ind; %寻找连续1
2 x% p- m( C% r$ ]- }- d* q
4 Z) A2 ~3 U* u/ n7 ] Z s
for ii=1:numel(ind);
- h. S# R! I8 w+ o9 z
up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
7 K T* x" Y/ Z8 N- g
end
' k4 d7 r" D2 ?+ ^8 D) U
8 t+ y% J- C; i
%% 下降
3 _7 i( s1 L# n! R7 J
k=diff([0 ~dy 0]);
! b0 \8 o5 ^( l1 Z' z
ind=find(k==1);
! Y+ [$ Z& i8 N! u m+ \
num=find(k==-1)-ind; %寻找连续1
! @+ u& f' j- E2 @+ M4 n8 Z/ M0 A
, c& A/ b# E J4 b! \! W" v
for ii=1:numel(ind);
: H6 l3 }$ r5 f* Y- S+ H
down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
" C1 X- I0 A# }2 s
end
) x$ s$ u' u. \( F9 i
* E8 U( P% ?# T0 ] Y9 x( i: f
4 A: l/ I9 o' w; L
运算结果有1723个，我想获取20个，即plot（X,Y）中那20个下降段，并输出该段的x,y，请教大神，非常感谢！
作者:
panphp
时间:
2020-6-10 11:34
% 读取数据
8 |% j) a B) ?
filename='test.txt';
3 ?$ Q3 R4 e. F4 q9 a$ Q9 f
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
/ C, u$ k" D, p
; X3 C! o# [: I" v/ ]
% 平滑操作
4 N9 ?9 p o+ Q" _: |) w
N=20;
" V) o2 z) ^; ~# E6 y( S
X=smooth(data_Depth,N);
& w1 n. {" m1 \
Y=smooth(data_Load,N);
0 J+ x) K* G: b0 ?. G8 `2 _" S
9 ], l2 k4 p/ p; R
% 找出下降段索引
- x% b O3 T) ?7 q8 J: o0 e
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);
8 x, ^" M- V# E) {
1 [) ?0 P4 L, O! x/ ?- l
% 画图
- `6 g: P/ g0 g0 ?0 _5 x Q
figure
2 d( W) n8 N7 W2 L9 s+ Y: ?
plot(data_Depth,data_Load)
- L( b2 t$ @) z/ [
hold on
1 v/ u5 }7 J; \2 x
plot(data_Depth(Ind),data_Load(Ind),'r.')
作者:
Colbie
时间:
2020-6-10 17:27
你提取的数据有误吧，我看down中不只有单个的数值，还有矩阵！
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2