找回密码
 注册
关于网站域名变更的通知
查看: 502|回复: 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)
    & {: r- k! Z" N, y$ C3 g  l2 U  j9 v. w, I3 O
    v=2000;0 j3 Z- D) U8 }) x$ g) J' I
    dx=10;, i0 d$ p2 `; S3 f& K6 `6 q
    dt=.004;%basic model parameters28 Q; u& N! E! R7 U2 ]
    x=0:dx:2000;%x axis3
    + X+ d9 E7 C/ t/ t( vt=0:dt:2;%t axis4& T! L' F4 t4 f5 v) _3 F
    seis=zeros(length(t),length(x));
    , K0 Y7 ]. I2 |4 A+ [8 W% seis=event_hyp(seis,t,x,.4,700,v,1,3);
    3 R8 r8 u0 Q' [, O! Y! k: wflag=1;aper=3;
    . O4 S4 E! V% R: P+ c, K2 H1 b2 Otnot=0.4;' s# E; c2 \, z; C/ |7 y
    xnot=700;; p4 n# K$ M3 |* U
    amat=zeros(2000,2000);
    1 K6 @# w* F/ F+ G5 b3 k9 mamp=1;
    9 A1 x$ E6 w: ?1 t/ m/ c2 d2 A# v6 B0 ~3 a% F
    % if(nargin<8)
    ' L8 f" w( D2 L; |) \$ _# V$ u6 h2 T%         flag=3;
    ; F" o) G8 Q( j4 l9 y* v( o3 a% end8 y3 X% v$ s( e/ ?5 z
    % if(nargin < 9 )
    ) o& [6 Y3 O  n$ l1 k. G%         aper = inf;' Q- @' b3 z" ?8 t$ S( C# M
    % end' }% x2 ]( \0 f

    3 u- W9 d- t4 T0 _* p( yv=v/2;; ?( C0 J- @, m+ z* j: K

    / ~6 C+ b# W/ }( G! l%loop over columns7 O* ~5 \- n" a  g1 ^6 o- n
    [nsamp,nc]=size(amat);
    8 w+ _% N7 d3 j9 y( E9 p6 L, m: T" v! \$ U- W( ~; U$ ~% c
    dt=t(2)-t(1);
    ' L" U+ `2 X: Z$ o* Atmin=t(1);2 O# N7 K* W$ J% W, ^
    for k=1:nc
    " G) x$ P3 k5 f) k! l        xoff=x(k)-xnot;
    ! E! `9 }! {+ D: \% a% C        if(abs(xoff) < aper)
    9 Y5 `$ P  l3 o% k: C$ i                tk = sqrt(tnot^2+(xoff/v)^2);$ {% w0 ?$ [5 H4 h( g7 P$ \8 @+ {
                    a=amp;; [! C; C4 t! K  [' B! Q5 v
                    if(flag==2)
    5 F( I9 o7 y# ]" o: U! d3 \( z; _' s3 y                        a = tnot*a/tk;
    7 P7 F5 x0 y4 N7 J! S                elseif(flag==3)3 F' g8 V  l  K, N# f( V: C; x' ^) ?
                            a = tnot*a/(tk^1.5);( ?: L" c( [1 o8 d$ b
                    elseif(flag==4). l5 J7 k0 Q2 O& G8 j8 T
                            a = tnot*a/(tk*tk);
    ( S- D/ g% K  T                end0 G. O9 ]3 a9 s5 v/ S
                    ik=(tk-tmin)/dt+1;
    ( t2 x1 M" N- E! D2 R! b                if( between(1,nsamp,ik) )( T! _( {6 q0 ?" d
                            ik1=floor(ik);7 B8 O3 E- Q6 I
                            ik2=ceil(ik);0 a1 _/ b4 M; }" z) c& A/ F
                            if(ik1==ik2)
    * u) j, k# r" `$ i; f7 `                                amat(ik1,k)=amat(ik1,k)+a;
    5 j! X; l! n4 S+ }                        else
    ( Z  ], ?3 E2 ]8 G9 E  \                                amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);  V$ T6 n3 Z( p+ q7 v( B. y5 t
                                    amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);
      \# A4 d$ b/ i, g9 A# \  G3 J                        end
    $ a# x! h7 h- @& L7 y4 |0 L" D                end
    . R2 @& D0 Y2 Z/ m6 }: e+ `0 R' r4 \: `        end/ a( h" k. A# B$ Y2 i4 Z* W! n' z
    end$ ~$ Y# n# [4 a3 [4 q3 v/ j0 }
    figure(1);imagesc(amat);colormap gray
    . W& f5 [* V# k3 }
    $ ?" ?1 F) _: E2 D在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家! m& m6 Q' s. j

    该用户从未签到

    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
    ) ^0 L; y0 o1 n6 x0 J0 p# e确实用错了。between函数参数类型应该是datatime类型。

    9 v5 d! }5 `! h+ `! N/ y2 y请问该怎么修改呢?
    $ \7 i2 E% Q* Y( ^0 _3 m

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。2 M7 Q0 l" t: \, k
    问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。( x0 h8 A* D% G8 e, y
    1 u3 P" M) }- f" O/ a
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-8-24 06:25 , Processed in 0.140625 second(s), 27 queries , Gzip On.

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

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

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