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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
3 F0 H3 Z# F6 u8 v+ z
下降段数据提取(附原始数据)
: z% }0 p; O! X2 Q. A4 O3 ^
6 b. r  |; l2 z& S6 @filename='';: K6 L/ {: k, U+ L
[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);& Z( q8 t6 K& c' G+ g7 V) |
X=data_Depth;  w% b. y" B: G# r1 K- c
Y=data_Load;1 C$ R: o6 G, p
plot(X,Y);8 o2 A* z- D4 Y8 F; V6 U5 B! i
y=Y';; q; E+ [# y) H
dy=[-1,diff(y)]>0;* T& }! w8 X( O. e4 p3 C
figure# d; o% S2 b; F1 G8 d) g
plot(dy);
) d$ |6 U% C! L" }4 f, F' O%% 上升& @2 y) j5 o5 Z) \% \7 A6 w% i
k=diff([0 dy 0]);
; `. M& G$ h4 Q5 ]ind=find(k==1);
6 Y7 t0 H5 g# w& @+ w, ]+ wnum=find(k==-1)-ind;   %寻找连续1
+ b& z' Z( U) t: u/ f! X6 h3 {+ M2 A  x# Z% D' B* W$ R1 a: W( W
for ii=1:numel(ind);
, _3 x0 U! v7 \% m2 r# B. X    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);
. H" \) _- K9 ]( Fend
" f$ t4 r8 @/ j( M& o4 f0 `! D, _3 Z5 Q! ^2 I# S# i% t' @
%% 下降% |8 L( E  w! h( {; X; z
k=diff([0 ~dy 0]);( L: U4 O& g4 x6 K! [
ind=find(k==1);9 B5 |" Y1 x" `1 S; O$ b, c) ]
num=find(k==-1)-ind;   %寻找连续1' ?: d" U& E: J" C1 j& j
) z( k9 u; C1 M- ^% \; p% \
for ii=1:numel(ind);
9 A1 P% u& U0 |! i4 h    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
$ {0 N8 s5 W/ w1 vend
% t" Z* P$ u4 }! e7 T6 ^1 ^
+ w: D$ ]3 w6 F: |
* ?8 W+ b. ~  x/ Q8 D: ^5 R/ t运算结果有1723个,我想获取20个,即plot(X,Y)中那20个下降段,并输出该段的x,y,请教大神,非常感谢!

该用户从未签到

2#
发表于 2020-6-10 11:34 | 只看该作者
% 读取数据
. w! {- _% V4 J- N: M  h  D2 Sfilename='test.txt';
; [( t! Z. a; w! c7 _[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);9 C* d" m. [! e' I$ ~2 T" S, N
/ j# b* q# w8 g! \; f" C! w4 d
% 平滑操作
. f; I) ]$ P3 J. `/ p: \( ^+ S3 _N=20;7 G9 G' Q& \  g8 ~3 b3 U
X=smooth(data_Depth,N);7 h5 \  j- z+ E/ H
Y=smooth(data_Load,N);/ ^6 n1 a; Y2 d: n. u" J
% C0 P- R, Z+ R' c9 e
% 找出下降段索引: |+ n3 Q+ l* x+ |6 n  L1 W
Ind=(sign(diff(X))<0 & sign(diff(Y))<0);
4 Z& v$ |6 a8 R4 O( ^1 t, ~) d* \" B  ?8 i( p  h3 z" t
% 画图
$ @( `% e5 [; k! D2 lfigure, g; D' R: D: @5 C6 a% N0 s+ s
plot(data_Depth,data_Load)8 g# [# ^; H1 I) Z( d# u
hold on' d  }( @- Q. c: H. s
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-24 08:31 , Processed in 0.109375 second(s), 23 queries , Gzip On.

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

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

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