找回密码
 注册
关于网站域名变更的通知
查看: 542|回复: 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)
    0 L1 R# Z* ]. }: x) [0 |& F3 ~- R  Q; Y7 y) [& ~1 p7 P. T
    v=2000;
    " q& G/ g5 K3 `9 Wdx=10;  ^; ~' M" k$ Q; J
    dt=.004;%basic model parameters2( U. r; B% `5 T0 t6 w
    x=0:dx:2000;%x axis3
    ) ^5 W8 J+ w1 s- Vt=0:dt:2;%t axis4
    ! k# U4 ]* j4 x+ G; _& s7 qseis=zeros(length(t),length(x));
    ( U3 F. K8 x1 H# P! M  Z& @% seis=event_hyp(seis,t,x,.4,700,v,1,3);& y1 K4 t. L9 \0 y4 L
    flag=1;aper=3;% y& u. F! p, n: l% c0 x
    tnot=0.4;
    ' z+ M  E9 Y5 b0 q7 ]% `xnot=700;8 ~% _7 w& L* ~- g3 O+ {( v1 s+ P8 Z
    amat=zeros(2000,2000);
    9 A5 I# D7 ~& W4 }" S' J' aamp=1;
    $ W2 H; q9 `+ M& Z, }
    $ E& m( A* ?( n- p7 L% if(nargin<8)9 ^  ?) m# b! Q' |: i; E
    %         flag=3;
    8 B; T4 K. q" G! v' d% V% end
    / k, T2 G# k0 t& V% if(nargin < 9 )
    ) h, s  _# b5 N( z# r; d%         aper = inf;
    , k6 S! f! w; T2 U( w% end
    ' T5 G! \" @- t1 x2 o, o9 u! q, }+ `8 }
    v=v/2;
    ! S; a. O' R* c4 W* u) }3 V4 x# F; L: ~9 b1 {
    %loop over columns9 d2 W; B  t2 n8 a+ O9 i
    [nsamp,nc]=size(amat);
    8 E0 e0 h/ I/ m1 G! W# W% k, Q5 ]/ _3 A: e; z7 @: l+ ^
    dt=t(2)-t(1);
    3 F. F& m, ^+ z  ~tmin=t(1);
    ! O" L% n' T$ E% [1 kfor k=1:nc
    ; |( ?) W' l+ ~        xoff=x(k)-xnot;
      H; c1 F  J+ s$ ~5 m) _' \- e        if(abs(xoff) < aper)
    9 A5 b& O# E( @                tk = sqrt(tnot^2+(xoff/v)^2);
    & q6 a* @* j4 y                a=amp;; ?' @/ y2 ^! m2 r8 L
                    if(flag==2)$ Y9 H; M/ k! t6 w
                            a = tnot*a/tk;* w* N% g0 i0 k' `) M; ?& u) }7 ]$ h
                    elseif(flag==3)
    ; f, Z: K# n+ e0 Q; k. ^& D                        a = tnot*a/(tk^1.5);2 j" \( ^& |6 n; ^7 t  v
                    elseif(flag==4)
    - J" I" p1 @3 N$ a8 D9 h" b1 X+ {                        a = tnot*a/(tk*tk);1 ^9 [, P2 y9 b, E) I. @4 t: Y
                    end, Z3 V* U8 W8 S" M$ A8 q, Y6 D2 I
                    ik=(tk-tmin)/dt+1;
    - C$ V) R- S* H: R8 N                if( between(1,nsamp,ik) )6 X+ x8 q5 |  v- j; l0 S8 E
                            ik1=floor(ik);! e* z) G; [  j6 ~' s, w
                            ik2=ceil(ik);& L- a" ?5 k1 U
                            if(ik1==ik2)% `1 F$ c9 _, }
                                    amat(ik1,k)=amat(ik1,k)+a;8 }# S' Z$ U- X6 |5 L$ |
                            else
    " f& i+ ^. y( [2 ^3 V. ~2 H                                amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);4 a  z* l: e. c
                                    amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);
    6 Y4 l- y$ G3 C# m& f4 x' q                        end
    8 }# w7 N: R: L0 n, e! [                end7 o8 Y. E; ]1 x6 f
            end
      P6 z+ n+ j# N$ P: i1 p8 y  F# v- Vend! E4 V/ S3 [* h/ A/ ^. m! r
    figure(1);imagesc(amat);colormap gray2 p+ [! h: X5 r5 N& u
    & O  T1 m& _8 r2 F
    在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家) l0 T  \8 k; S

    该用户从未签到

    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:27
    3 ^4 F$ {) U' M* U5 F' X确实用错了。between函数参数类型应该是datatime类型。
    % |# X4 X! x4 h5 u
    请问该怎么修改呢?  |  f! b5 A" T+ I+ K' f/ n+ {! D

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。5 l6 @. U# R; x, z/ N8 x  R) ~
    问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。
    8 k5 W2 Z# ]+ m6 |! H
    # k5 a0 \2 S, N4 J6 r/ H/ Y
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-10-28 00:10 , Processed in 0.125000 second(s), 27 queries , Gzip On.

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

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

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