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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
1 |6 G5 X& @8 v# r
下降段数据提取(附原始数据); ?- E7 J$ B5 \( a4 R  _
% g) z% K% y$ P) f8 X5 u
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);8 a" O- @( u7 m
y=Y';0 a0 X/ u: a6 [, c% _' [
dy=[-1,diff(y)]>0;5 ^' S, N4 d; l2 q: R: T
figure
: ?+ a6 l( r# u* K  a* s  W* hplot(dy);
& `* _, i" U( P- G4 j% U4 z%% 上升( E: I# x8 D' Q' x( ^8 _
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 y4 v" S  C, V5 l
for ii=1:numel(ind);
( `/ ]5 f: N4 K$ c' d5 Q    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);. t: m( R  f; J: i
end- d, P& G. D% n( L! p3 \$ V

* 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);) ~+ Q, V5 u/ T* v, d8 _
num=find(k==-1)-ind;   %寻找连续1
. M" f# e: Z4 }. Y/ ?( C& {& B% t( |  h" M! Q
for ii=1:numel(ind);1 ]& n, E7 M( y2 V( o
    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);8 s- A& u, p1 g7 H' M4 f6 A
end
) K" P' }! M. }/ q4 H( I0 w% q) G% w

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

该用户从未签到

2#
发表于 2020-6-10 11:34 | 只看该作者
% 读取数据
) 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  t1 B% M/ }, O4 X$ D
% 找出下降段索引8 G0 m) S* ]. t! `
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);1 Q' g7 e5 I  F# ]% q' J; B

, V$ m8 u6 K) Q/ Z: J9 r6 Y1 d+ C% 画图
' h& z0 O7 y3 F5 ^) z' Pfigure  f( p# o- C# {1 V
plot(data_Depth,data_Load)
' K0 f$ d4 M0 w' C: ?1 P) Thold on4 z0 f3 O# v. P( a9 M
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-11-5 13:54 , Processed in 0.156250 second(s), 23 queries , Gzip On.

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

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

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