找回密码
 注册
关于网站域名变更的通知
查看: 367|回复: 2
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-6-10 10:53 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
; {& z* b3 l/ i  r8 X! G4 `
下降段数据提取(附原始数据)0 j% i* Q* {4 }# b9 _3 r! c
+ a# G/ N( {& y% x' {
filename='';' T1 @, d7 N: D7 |# Q3 x
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);4 I" E$ A+ M9 Z5 E. m
X=data_Depth;
, u* {6 j; ?( ]0 o! E9 _: ^Y=data_Load;
5 P" f" |+ i% Aplot(X,Y);( ^% o) j5 R' _- }2 l6 @. b6 g
y=Y';
' A  u" K# {: N3 I8 Y5 C; m6 Ndy=[-1,diff(y)]>0;5 c3 Z) q$ V) ]
figure
  v$ m1 d+ u+ g$ O/ N! rplot(dy);
$ r6 ?$ D" R9 ]& _/ u' R) S%% 上升0 \1 a- O1 `" s7 _* p
k=diff([0 dy 0]);& r- t* e* z, i
ind=find(k==1);& g, t' i6 O6 s6 m
num=find(k==-1)-ind;   %寻找连续1) _4 s/ V1 _% N/ Y: n

3 y8 _- ]$ ]+ ]4 \# V* y+ ^for ii=1:numel(ind);
  ?+ ~7 a9 I3 u' [: X    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);) j$ u; d/ w% H
end
5 M9 S  c4 q: X4 F3 R  z6 E
  H2 V6 ]" b$ u1 o8 W8 @: B%% 下降
4 y4 p/ X( s/ B; [' @* ?k=diff([0 ~dy 0]);
  ?3 e. Z9 V- z; P4 aind=find(k==1);
7 A) a( E- o; C+ N$ hnum=find(k==-1)-ind;   %寻找连续14 P; U8 A: Q( p# R

) ^% s  T6 ?3 m. p3 b0 Nfor ii=1:numel(ind);9 U5 U6 J3 ?0 N
    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);* Z( ?) I6 ~# ^3 s0 B% Q
end. `4 s' I& s# Z" C1 h2 X0 L: l
$ J( V4 \& R+ c  R4 L, ?/ P

+ ?6 ^$ p" j' N' o4 c$ m/ D运算结果有1723个,我想获取20个,即plot(X,Y)中那20个下降段,并输出该段的x,y,请教大神,非常感谢!

该用户从未签到

2#
发表于 2020-6-10 11:34 | 只看该作者
% 读取数据. m& H+ l7 P6 Q% ]6 G: y7 q( t. Y. E
filename='test.txt';
/ Q# r; T" ^3 ~, C6 U[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);: x* v7 H, z5 R2 e: d0 _( @3 n7 |
. _) h# b7 {* A; A0 n- h6 s
% 平滑操作- E1 R: e! W+ T9 C6 Q# V3 |& W
N=20;
) g* U' n6 s1 LX=smooth(data_Depth,N);; f' U4 E& P" `2 M' e) j8 {
Y=smooth(data_Load,N);
# u3 g" a( N- M/ \; v" Z: x! e
) S. P& b; Z: {! U% 找出下降段索引* g% i7 L' A! }( k( {
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);3 |( _& ?) \+ ?2 e% `9 ?

  \3 z2 b! a9 k8 I$ }" n% 画图
% z4 i( o. r0 ]figure9 [3 b1 \+ R4 B4 b( \; _
plot(data_Depth,data_Load)
1 ~" M3 [; K9 E  d, L% W& Uhold on- {; ~: G; T& L8 S5 h& ^
plot(data_Depth(Ind),data_Load(Ind),'r.')
  • TA的每日心情
    开心
    2019-11-20 15:05
  • 签到天数: 2 天

    [LV.1]初来乍到

    3#
    发表于 2020-6-10 17:27 | 只看该作者
    你提取的数据有误吧,我看down中不只有单个的数值,还有矩阵!
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-8-19 18:20 , Processed in 0.109375 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表