找回密码
 注册
查看: 471|回复: 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)
    ! ^, }; Q0 g# V" E& W/ T" B/ T" D1 H2 d
    v=2000;
    % ]6 D! {, H4 v) ?/ wdx=10;; d3 j* h0 [' M6 n# U
    dt=.004;%basic model parameters2* }) {9 \, S9 p  z  I! B8 Y. q
    x=0:dx:2000;%x axis36 |  d! Q: `2 E
    t=0:dt:2;%t axis4
    $ a+ x4 P  o' i8 M# o% Eseis=zeros(length(t),length(x));5 C3 A$ T* S& C6 Z) O% @* A( ~/ ?4 z
    % seis=event_hyp(seis,t,x,.4,700,v,1,3);$ b& ^( P! q$ [( [7 ^
    flag=1;aper=3;2 B5 u  w0 h$ R& J( I) ^, x$ U
    tnot=0.4;" v* U+ ]1 u" P0 i/ T$ Q
    xnot=700;
    9 ^5 @( Y3 y; [% t& k% ]; {* hamat=zeros(2000,2000);
    ( T1 n5 [3 N& k4 N: gamp=1;
    9 N% X/ k( L; ?' B& _
    9 y9 Q* C! J3 ]& F8 b% if(nargin<8)6 N& _$ ^" i$ d6 {
    %         flag=3;
    # I; G0 f. ~2 ]0 q% end
    0 v5 E+ N  P9 v( f! D0 ~$ N/ E% if(nargin < 9 )7 r  X) b. u+ }
    %         aper = inf;  ^1 z& h; T: p6 }0 k
    % end
    ! M  R9 t% l  ^. V3 }
    * N! J" @" ^, e: {/ c% h3 H4 Ov=v/2;
    1 N3 ]% [5 o3 O/ I' U  ]+ J) X* A& \# P. R( g3 w
    %loop over columns. ^2 `9 f6 c0 r; k5 O  m7 p! Y4 ^2 \8 Z
    [nsamp,nc]=size(amat);+ i" H# ~1 \" j& N$ |2 H& V8 F5 @/ _

    / D* a/ |$ I, K' Zdt=t(2)-t(1);! Z9 S8 T. z" K8 G2 V4 a9 C) C0 |9 i
    tmin=t(1);
    2 E) M0 G- a% X5 y' O) Ufor k=1:nc
    3 T, z6 J& n3 p        xoff=x(k)-xnot;
    ) ^1 i7 y4 d- ~) l9 a5 p        if(abs(xoff) < aper)- K  f- j( j* ~* n6 u
                    tk = sqrt(tnot^2+(xoff/v)^2);
    7 z) x, U  ^/ D6 N' ^' k                a=amp;7 h' k2 z8 h! _8 R8 b" M5 d* u
                    if(flag==2)$ D1 O. o6 _8 _9 i* X# U% s+ C: d
                            a = tnot*a/tk;7 B6 ]. r! q' o" G5 k; E! J4 P
                    elseif(flag==3)5 j1 {2 g$ ~- x; q. J( s
                            a = tnot*a/(tk^1.5);6 Z$ O+ U9 y% h: K
                    elseif(flag==4)
    0 V; u! l3 `+ k, t# a  \5 y& O                        a = tnot*a/(tk*tk);+ I5 R* J" n% W% h! L, m
                    end
    % m2 B7 ~' y4 `" U, U                ik=(tk-tmin)/dt+1;
    7 G5 e$ i! f$ t' `+ b# y4 o                if( between(1,nsamp,ik) )
    " w+ e3 }  ]: A3 I# G8 i                        ik1=floor(ik);7 k+ j2 T* Y; ^. g0 j& s6 b: X1 ^
                            ik2=ceil(ik);1 X1 |( @( w' W! T. o
                            if(ik1==ik2)
    & ^2 a- q) ]& S5 N                                amat(ik1,k)=amat(ik1,k)+a;
    . J9 B  M: x$ p5 K5 _                        else
    + ?. Y2 C) M. w7 J2 P3 u! q  P                                amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);- y+ d4 A* N6 f% x& E" I# v
                                    amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);7 J' z3 w5 |1 p
                            end0 n) |6 F- w6 [* J  p
                    end
    " n: A+ V# j/ w7 J4 c! D        end5 ]! n/ k& G. U' \0 d
    end
    $ r. ^8 C5 b8 g# J8 Mfigure(1);imagesc(amat);colormap gray0 q. r/ e% A4 p2 |' y2 h
    # ]' R9 H9 F5 ?5 v
    在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家) J3 c6 e) W# h: n2 w

    该用户从未签到

    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
    % A6 u; l: A  |# l确实用错了。between函数参数类型应该是datatime类型。

    ) ]  K7 h6 T3 }  z请问该怎么修改呢?
    1 W4 Z$ s- M* U6 m- {0 m. c

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。
    / |) V& B% X" J0 I, |问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。! x# E% ?& l  L' x% K
    & i) O+ F- ^$ S& S5 f" [
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-5-23 22:42 , Processed in 0.078125 second(s), 24 queries , Gzip On.

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

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

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