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

关于这段代码between的用法

[复制链接]
  • TA的每日心情
    开心
    2022-1-21 15:22
  • 签到天数: 1 天

    [LV.1]初来乍到

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

    EDA365欢迎您登录!

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

    x
    % function amat=event_hyp(seis,t,x,0.4,700,v,1,3)
    ) v! S5 D: j: m7 F, S
    ) K. L. g& N: G! Xv=2000;
    % N$ Y* [5 ?6 c0 b4 }dx=10;
    0 K% _- H% n- T/ O5 e8 bdt=.004;%basic model parameters2  R5 e& ^6 ]; u
    x=0:dx:2000;%x axis3' v0 L& R* E- _4 Q# s2 I) i
    t=0:dt:2;%t axis49 N3 W' U4 A. h
    seis=zeros(length(t),length(x));
    4 }( y3 i* E. S4 b- K% seis=event_hyp(seis,t,x,.4,700,v,1,3);9 w5 t+ N* }- j' h
    flag=1;aper=3;. D( p" p* q- V$ R2 y
    tnot=0.4;. D) [, V/ f4 j& b
    xnot=700;
    7 j5 c6 v# N5 Y5 @( N  ]  v! \  pamat=zeros(2000,2000);
    * y: h) N- Z4 L5 l4 F% G7 {amp=1;; M. s4 i  z5 P3 p
    ' ?/ E2 h/ l4 l& b  W% ^: X+ g
    % if(nargin<8)4 [2 j* V( d& J( L
    %         flag=3;) H+ I. y* v! M$ H2 C2 E
    % end* p7 O( U& I6 t. X. O5 E' G! _( @! Y
    % if(nargin < 9 )
    ; `) w! S0 ~  x  \2 J$ X%         aper = inf;
    # E) j3 u9 Y% ]( z, t' N( I3 |. [! P" v% end5 X' M' a( b/ P, G$ W: A
    . _& K9 X" b, f5 g
    v=v/2;
    + ^& [5 o5 l& p8 `5 s" h2 D5 F
    + E9 D4 p. o3 d%loop over columns
    ' s2 r- C9 L2 y; T- e[nsamp,nc]=size(amat);2 n$ s6 l7 }) L4 a* d! }9 L

    & _) B  j8 n0 I7 r7 A7 K& {dt=t(2)-t(1);
    ' z4 ^7 R7 |( {. atmin=t(1);
    " E6 r( Y  q: E# D; J( [3 Mfor k=1:nc
    + Q7 P! A8 p( a# o, F) Z  R/ z        xoff=x(k)-xnot;$ M% B# T9 W, ~! `& o
            if(abs(xoff) < aper)$ b3 ?5 E6 I, J( m3 G
                    tk = sqrt(tnot^2+(xoff/v)^2);
    + B7 k7 _& _. [+ R% F  Y4 |                a=amp;
    : j/ K8 D3 J; |' G9 q0 h* A$ k                if(flag==2)
    % ]4 k: C3 {4 b7 Z; V8 o, B                        a = tnot*a/tk;8 t' a; W4 s1 p& s+ X
                    elseif(flag==3)
    * q% P2 _+ G" L4 A: N+ y1 @* n                        a = tnot*a/(tk^1.5);5 v0 T2 r7 {% ~1 O) m
                    elseif(flag==4)1 J4 `; g+ ]/ D. n* ]
                            a = tnot*a/(tk*tk);
    ' N4 g0 J0 s- f; @* T  Y                end4 H# |7 a( F, Q6 W# U) S; W
                    ik=(tk-tmin)/dt+1;
    $ ?+ s/ c7 \1 _9 P6 W                if( between(1,nsamp,ik) )3 x' v5 n8 T8 U( _: G% w
                            ik1=floor(ik);2 c! `5 a! C1 u5 c! h
                            ik2=ceil(ik);- M% X6 V8 H3 W, x- l
                            if(ik1==ik2)
    $ N5 ^% @6 c- G                                amat(ik1,k)=amat(ik1,k)+a;
    2 `2 o# W& L& K4 |! W                        else
    ( c, N9 g8 q! c- l                                amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);: X7 ?7 a" i' H0 j; s- y
                                    amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);$ P3 T* |$ Y9 ^$ ]! K: M6 `
                            end
    - _: m" P( `4 X) a3 m8 _  i! x; D- B                end6 _( P! A% a8 E! ?% r
            end# f  v9 c4 j7 s1 H  w2 \+ ]" j2 S
    end2 L4 F" Z- m* T6 ?! ^# @
    figure(1);imagesc(amat);colormap gray
    9 O4 D2 D: H$ R: @& j9 [
    ! o# u! P/ u2 A/ L( V- J  H在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家
    8 N9 E' [+ A  _. x. r0 o$ P( i

    该用户从未签到

    2#
    发表于 2022-11-9 11:27 | 只看该作者
    确实用错了。between函数参数类型应该是datatime类型。

    点评

    请问该怎么修改呢?  详情 回复 发表于 2022-11-9 13:17
  • TA的每日心情
    开心
    2022-1-21 15:22
  • 签到天数: 1 天

    [LV.1]初来乍到

    3#
     楼主| 发表于 2022-11-9 13:17 | 只看该作者
    li205212021 发表于 2022-11-9 11:278 L& D& S! U( l5 j" Y
    确实用错了。between函数参数类型应该是datatime类型。
    " H: ^2 ?3 m* H( l2 k8 p0 a# z
    请问该怎么修改呢?* H9 i9 }4 |& H5 l/ ]% A

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。. ]$ K3 X; z' Q$ j, D" ~
    问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。* k, N9 W* m! i8 @, r1 `' B  p

    : k7 u8 u: S$ H* K0 P6 D/ I. {
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-23 16:37 , Processed in 0.140625 second(s), 24 queries , Gzip On.

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

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

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