找回密码
 注册
关于网站域名变更的通知
查看: 498|回复: 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)
    3 M( ]- L$ n: W8 C* I  s0 ?
    8 O; p$ l6 g8 N3 }$ c; y, zv=2000;3 l# ?% Y: r+ s$ J& z, t) ^' s
    dx=10;
    ' p1 g& w( L' P9 y7 J1 t; {dt=.004;%basic model parameters2. H2 g; g8 P# b
    x=0:dx:2000;%x axis3
    3 e" p9 [$ H, V4 T6 q# ~t=0:dt:2;%t axis4! [2 N/ S# a/ \" J& W2 @
    seis=zeros(length(t),length(x));6 l6 ]* D( U6 K9 I; G8 z
    % seis=event_hyp(seis,t,x,.4,700,v,1,3);; W, ?+ M+ b4 L
    flag=1;aper=3;: v; @  F( v& X# [
    tnot=0.4;
    - {. t- t& ~7 S5 Q7 x  Fxnot=700;9 T+ W; Z2 `) U& b
    amat=zeros(2000,2000);
    & p3 s3 o$ |1 i  k. n# W  Damp=1;% B5 H: E5 Y$ C0 f# R

    ! m% [7 g# S4 g2 J% if(nargin<8)# Q( p3 ?: W0 q: |4 e5 z
    %         flag=3;
    4 W: ^5 Q2 k( S9 o% end
    ) R9 ?9 Y# V1 c2 [% if(nargin < 9 )1 F* Z" D7 ?; G5 X
    %         aper = inf;- v/ }( w$ @* Q; i/ z3 h
    % end
    ) A/ P2 W1 v5 v1 g% n1 v
    * G& ?% q# r# I$ L/ [* xv=v/2;. L/ C2 ~" K2 O( L
    - s* T1 m6 U1 T' J
    %loop over columns. d0 {# {, o0 s$ Z
    [nsamp,nc]=size(amat);
    7 j$ O  R* c) k+ e8 B. t& C5 G2 X6 Y) c, y, k2 N2 c
    dt=t(2)-t(1);
    : {; c& h9 G- Y# Htmin=t(1);6 h, z7 {+ B* @9 L- `' q- o4 J# }
    for k=1:nc1 X2 \3 G2 W" Q0 Z2 I" [
            xoff=x(k)-xnot;/ r5 H3 m# F' ?# j4 ~
            if(abs(xoff) < aper)
    # `- h1 J. P" D& \, T- ~                tk = sqrt(tnot^2+(xoff/v)^2);& J7 O9 n( @8 s5 I. D
                    a=amp;
    - A4 i+ z( g' G  W) l                if(flag==2)3 w7 S' ^- J+ v# |: S: b
                            a = tnot*a/tk;
    $ m& H/ V$ [: X+ A. E4 A                elseif(flag==3)
    3 s8 \! w% ?0 {" \                        a = tnot*a/(tk^1.5);. E7 ^6 T9 ]- c2 M1 A
                    elseif(flag==4)
    8 {0 \: p7 y7 z/ C1 X" T) M# G2 Y                        a = tnot*a/(tk*tk);9 q6 d. P7 U1 U) f1 K" x4 ^
                    end
    8 \1 e" h1 b  Z8 f" _1 x                ik=(tk-tmin)/dt+1;
      P) W4 g- R, {( c' Y- q& O                if( between(1,nsamp,ik) )4 D# d. I1 C1 k0 |0 q" M# v
                            ik1=floor(ik);
    ( e- w0 \$ |  V' y5 h  I                        ik2=ceil(ik);
    1 I2 H& V6 B6 ?" y$ N) N8 t                        if(ik1==ik2)
    * s- ~1 H& Q& N+ {3 H& Q                                amat(ik1,k)=amat(ik1,k)+a;
    6 P* }/ {1 [! H                        else
    ; n+ ~+ W/ u) A9 D* ~                                amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);
    0 V, x/ x; T2 d1 ~                                amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);
    - \; H& v. B% q1 i                        end
    9 q3 a0 U- K" f. t) \& j                end* z1 l1 B8 ?; Q# C8 |
            end! c9 O4 G! J7 ?# Y+ [7 X
    end
    * ^2 b5 P0 V2 ^; Ffigure(1);imagesc(amat);colormap gray
    , l' S; [. {( t1 N5 q8 q+ ~6 F
      g! @! d, x1 d) j+ L在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家3 a; G6 r3 N1 }# n* U4 |

    该用户从未签到

    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* x# p/ E. N$ W7 I
    确实用错了。between函数参数类型应该是datatime类型。

    " y. h: b6 X, M& b' n( r( B4 H请问该怎么修改呢?
    ; V1 G1 O. F& l4 [& _+ Q. {5 ?

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。1 H+ X0 e1 O. F8 r7 V
    问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。
    6 K; e4 w, m  f) h  T+ y3 ]5 U- W+ O4 j8 B1 Q7 F* [
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-8-12 04:59 , Processed in 0.109375 second(s), 27 queries , Gzip On.

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

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

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