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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x

3 Y" G7 J' c! ^6 ]下降段数据提取(附原始数据)
3 G4 a+ C9 l9 i% T4 g7 t  F& ~/ `6 ?2 `' \
filename='';
0 d. r! X! o1 \[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
. l. Q5 Q2 g' `X=data_Depth;  r+ A' U8 e! S) K" W+ O
Y=data_Load;0 o4 W8 T! T$ q; D8 p  J
plot(X,Y);
$ P) i0 ~# s. W6 f& N1 Z4 x# yy=Y';6 k. M, n! k/ C
dy=[-1,diff(y)]>0;
; p4 d! m8 N8 ~; U2 \8 Sfigure3 Y$ n3 V7 r  D1 ~
plot(dy);. s& n, h7 H0 n
%% 上升' c7 h# Q& P! j4 h/ {6 ?
k=diff([0 dy 0]);
, c' p1 {# J& B( Y. mind=find(k==1);
" i3 w- h1 w; ?2 O, j- i$ i% anum=find(k==-1)-ind;   %寻找连续1
! {9 R# {/ P( o: |' w* Q4 N
2 i3 O/ n' J$ l+ Lfor ii=1:numel(ind);
! ^5 s6 b9 e. J3 s1 b( s& v+ i% k    up{ii}=y(ind(ii):ind(ii)+num(ii)-1);" I2 @& w0 E  d$ d0 Z1 K& e
end
9 j9 I/ y+ y, _- \5 X$ x( D5 U/ h0 a+ w/ O* |+ f0 q
%% 下降+ n1 q/ ~/ v# Y5 S) Y: J( Z
k=diff([0 ~dy 0]);
6 Z/ o5 t5 K" m- oind=find(k==1);
8 N' W! H' K/ ~& gnum=find(k==-1)-ind;   %寻找连续1
# @  y3 {5 t' d5 I* S9 V. j8 I) i2 |/ @2 `9 i/ u4 i9 o( D% n. N, Q
for ii=1:numel(ind);
. u, F1 k4 G* g: g    down{ii}=y(ind(ii):ind(ii)+num(ii)-1);
$ H# j8 U: Q' G) fend* h- t. a  ~8 @  N/ Q. K

$ X' a6 }+ L7 E: w" c. C
$ {7 s/ ]7 M: ~  Q' W& u运算结果有1723个,我想获取20个,即plot(X,Y)中那20个下降段,并输出该段的x,y,请教大神,非常感谢!

该用户从未签到

2#
发表于 2020-6-10 11:34 | 只看该作者
% 读取数据
- m5 s. f( H: t# `8 z( W: T, |filename='test.txt';
7 V; W0 o0 o. X& ]' m; E[data_Depth,data_Load,data_Time,data_DV,data_LV]=textread(filename,'%n%n%n%n%n%*[^\n]','headerlines',3525);
! s) V0 |& V' y5 w
4 x! ~  D( a* L! [3 Z) S% 平滑操作
# T5 j/ ?+ d4 _* g& ~" jN=20;
1 Q0 W$ x" J& d% l& WX=smooth(data_Depth,N);8 G: H& }, D- ^3 y- I* L
Y=smooth(data_Load,N);
/ L0 l, N/ C" I# j* s5 x4 z' r; C9 k. b, y2 [
% 找出下降段索引
& b  W1 i. {1 ~2 ^4 }Ind=(sign(diff(X))<0 & sign(diff(Y))<0);7 e+ i' A0 b1 d/ w) R* r0 ]
# }2 |2 ]7 _' D/ A2 {" x
% 画图# Y" U6 ]/ l+ h/ I- y
figure
$ G/ v: k, L+ D1 z1 {& nplot(data_Depth,data_Load)3 o- b& C( l- `7 t7 S$ B* z5 |
hold on2 Q) Z0 ^% W2 ?
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 14:20 , Processed in 0.156250 second(s), 24 queries , Gzip On.

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

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

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