找回密码
 注册
关于网站域名变更的通知
查看: 560|回复: 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)
    4 Q! z- X* r: k8 [) W3 y! z- m* t1 U, Y2 T+ n# n6 F
    v=2000;
    ; {  r8 f& e: a0 {# \8 \  g3 vdx=10;# U9 s1 `! S2 F* l; ~: A
    dt=.004;%basic model parameters2
    6 g' v% t  C: C8 Ox=0:dx:2000;%x axis3
    % O) \! d5 J2 X1 w8 q; v! ut=0:dt:2;%t axis45 K( _: ~! W& P+ W
    seis=zeros(length(t),length(x));
    - N* V5 b* O# _$ x1 h! i% seis=event_hyp(seis,t,x,.4,700,v,1,3);2 V. Q# y) `% }1 X& e
    flag=1;aper=3;
    1 f! L8 C5 ^( @tnot=0.4;7 M* H  Q$ k8 P% h
    xnot=700;
    " m9 S9 L6 \! G$ Z9 G& i" }/ r# Kamat=zeros(2000,2000);
    5 M  O$ H2 _7 U- o& n+ Q: Y7 Gamp=1;& v& g/ \- H$ u8 Z" I! b
    1 W% \, w" n: S
    % if(nargin<8)
    : }3 r' C$ S8 W5 f% l%         flag=3;4 [# m5 Z0 Z2 L' r5 U+ R
    % end2 l  X: @" |& r7 R; C! t
    % if(nargin < 9 )# U: c  m9 k$ |
    %         aper = inf;
    . D0 S8 T* W- b; m% end
    " L+ K  F' |! @! T$ |
      x6 q' q7 n4 b& Pv=v/2;8 ]) J( n" y* b( p0 m) c
    / J0 l- Y, G" u7 B0 `: o% V% f: m
    %loop over columns
    * h# s* B% a6 Y2 M$ x1 x[nsamp,nc]=size(amat);
    & Y: x4 Q) `- ?( d1 E# L5 j0 N7 D  n. v6 s
    dt=t(2)-t(1);$ [6 y! O: h) d) q. o6 i& `0 J: Z5 U
    tmin=t(1);5 ^. z/ L. ]% K( W8 s! c
    for k=1:nc
    2 U+ N0 c1 u) Q: c% |        xoff=x(k)-xnot;3 l. R7 v% I. J/ [) ]/ g
            if(abs(xoff) < aper)+ S# J" y; i2 I0 f( P
                    tk = sqrt(tnot^2+(xoff/v)^2);
    & N/ Y6 m' \  C) a. F: i3 H& J                a=amp;  t3 r2 J8 C7 n) q! F8 s
                    if(flag==2)
    4 B. h. Y3 O1 F                        a = tnot*a/tk;
    2 r! V7 o# H& w                elseif(flag==3)* |5 z. Q/ @0 V1 G; e6 o: P
                            a = tnot*a/(tk^1.5);) F; L4 \3 O8 l# M$ b9 _
                    elseif(flag==4)
    " x% |; S$ E" F- M( c9 r/ u                        a = tnot*a/(tk*tk);6 i) H3 u) x, R* P# z" h2 ]# k7 \
                    end
    & r- u) M3 I$ a! R0 [$ Q& r                ik=(tk-tmin)/dt+1;, B& r0 m, @9 X+ T# f; P/ _
                    if( between(1,nsamp,ik) )
    5 x5 s" e* V' H; ?# R                        ik1=floor(ik);
    # n" P1 f2 \, N( R0 c5 L                        ik2=ceil(ik);, e- m4 e5 S+ _1 k' p
                            if(ik1==ik2)
    8 |) f# T1 Q+ u$ @; r/ n                                amat(ik1,k)=amat(ik1,k)+a;& \9 |  ~7 I( @0 V7 f1 [& C
                            else+ s* V4 D8 i6 r# t7 i* T* `/ O
                                    amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);0 w" Y, X% z  w- i) v0 X
                                    amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);
    $ ?0 T7 x* w  x& ~6 ?5 M! N                        end
    5 X9 V9 |3 m2 z# z) |0 L                end
    " \& f6 W. v5 \/ b& k4 h5 u9 _        end1 M. D6 L! x% r
    end4 w; i+ z8 o2 A8 v4 }; R
    figure(1);imagesc(amat);colormap gray
    . z7 ~* i3 G3 f4 n( L  c# H3 M+ l  Z7 D
    在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家
    $ T0 ~' _5 ?5 L9 n

    该用户从未签到

    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:279 H0 Z: ^7 T1 N3 _- }0 i8 l
    确实用错了。between函数参数类型应该是datatime类型。

    0 m4 P  z# y1 u) x1 i: J! y% b请问该怎么修改呢?
    - @$ z# u, R9 c) A' A) S

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。
    & D! O+ C+ B2 c" y9 ~: u8 \问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。2 _2 s5 V: @2 `8 @7 G7 l2 _# c

    9 I# e; E7 h/ u- T1 G4 l" W5 y; u
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-23 13:49 , Processed in 0.156250 second(s), 27 queries , Gzip On.

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

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

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