找回密码
 注册
关于网站域名变更的通知
查看: 564|回复: 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)/ B/ O# D( J& l$ j7 m/ \

    - U5 ?& B& J) l7 m3 Y3 p0 O) Hv=2000;
    ' v! f$ C0 L) Z9 |* tdx=10;+ E8 C4 t6 x) h1 ^7 U
    dt=.004;%basic model parameters2/ ^- m' V9 ]- v
    x=0:dx:2000;%x axis3! e3 C% z5 Z3 Q
    t=0:dt:2;%t axis46 [1 T4 Y) q5 s: F7 \
    seis=zeros(length(t),length(x));
    2 \/ e2 x: e7 Q# A& `% seis=event_hyp(seis,t,x,.4,700,v,1,3);  }' l9 h8 Y, j+ B( B! @- `% t* V
    flag=1;aper=3;9 T8 B6 L) J5 f9 v3 _& q
    tnot=0.4;( o/ Q4 m# k* Q5 I. Z8 P- B
    xnot=700;7 F  E- ~, t- g* U2 n: i3 M
    amat=zeros(2000,2000);" ~: Y  U( e1 {& p4 x( G
    amp=1;
    + |! B# x: h5 S' z4 a
    5 L! f% n3 U8 ^. S3 p' a% if(nargin<8)
    # t' m% V9 A% t2 f# O7 x; x%         flag=3;
      K/ `, V. O( a: A% end
    - j, H" ~' o) e+ b+ I7 u% if(nargin < 9 )
    - h4 O( B( O! ?%         aper = inf;: V2 \* z/ ], k) Y
    % end
    * Q7 Z; @% D4 O" q0 B5 R
    " ~3 l* s8 E0 B& Yv=v/2;4 O% e6 ^8 V7 m' A

    - w  I4 Q* k& W%loop over columns) b  X5 d* d' }  H
    [nsamp,nc]=size(amat);( V/ J+ `1 K( B$ h) w$ q/ H

    : l4 w% }9 i& b0 k- F! c/ |2 idt=t(2)-t(1);
    " h4 I6 r! `7 ?+ [tmin=t(1);, l+ c* C1 ~" G
    for k=1:nc
    % C) b" p$ {( d" x( n        xoff=x(k)-xnot;
    - T( u1 R  \- H; |. F        if(abs(xoff) < aper). X, H: a( B9 Z( O4 p2 v
                    tk = sqrt(tnot^2+(xoff/v)^2);
    , O8 }3 Y* O; r" P% e                a=amp;
    ( {; z, M1 c2 k. h6 F: @% C3 v                if(flag==2)
    * v0 P) j4 }4 y$ H6 \0 c                        a = tnot*a/tk;
    0 s1 p# a9 J+ S0 X                elseif(flag==3)
    9 x( \3 Z1 e' ^* v* R- _3 y7 \                        a = tnot*a/(tk^1.5);; I1 E! g  b4 F4 S6 G/ P! l
                    elseif(flag==4)
    * {. z1 J+ u0 U; v, J" k/ X! q                        a = tnot*a/(tk*tk);& Q& o4 O, S8 Q8 I
                    end; l+ @, W' [, }; l1 ]7 \3 w5 f
                    ik=(tk-tmin)/dt+1;: u) [! t, B8 j5 X5 Y. n1 n' G
                    if( between(1,nsamp,ik) )/ Y' i4 |- A/ y7 T: B+ i
                            ik1=floor(ik);: b* f1 \8 G) w9 _) }
                            ik2=ceil(ik);2 I1 k- r: W8 M! [
                            if(ik1==ik2)
    , h% {9 [$ H: H9 D3 H; D5 I                                amat(ik1,k)=amat(ik1,k)+a;) W( J9 B& U6 u5 k) M+ x% v) h
                            else
    ) Y2 m! h- H9 Q' G  U                                amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);
    - L4 R1 g( c: s                                amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);
    , r/ ]& O7 `2 C& h1 b! S+ a' L$ N                        end8 V2 J  c, J: @4 ]
                    end7 B- W3 v7 ~9 X) U, K. O' ^* }
            end
    ! t: O$ Q7 t* H/ @) s- Zend
    % X; `1 I9 F2 e  mfigure(1);imagesc(amat);colormap gray  e$ w: P, ~% c! I1 Q, Z

    7 q% ~  R* ^5 I2 V  q: V- t$ \7 E在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家
    / `, D! O% P4 `

    该用户从未签到

    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
    $ l3 K) \& f( L$ `7 y确实用错了。between函数参数类型应该是datatime类型。
    # o" [: }% G( |4 b
    请问该怎么修改呢?$ ^4 B$ e. k3 @! I- u- R

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。- ~: `( F. t/ G2 \* u
    问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。2 C/ w, z0 u6 L4 |; T! D6 U% B
    & T) U! a8 U' ^! x  a3 b
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-23 21:57 , Processed in 0.140625 second(s), 25 queries , Gzip On.

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

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

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