找回密码
 注册
关于网站域名变更的通知
查看: 562|回复: 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)% N  l# f7 {* y' R, K$ h6 o

    + q/ k, X& O7 B* j+ _6 Av=2000;" E+ M2 l- ~4 z4 i. D* j
    dx=10;- z, o  e9 l8 a, o2 E+ W
    dt=.004;%basic model parameters2
    6 i8 L7 X4 B9 U2 mx=0:dx:2000;%x axis3! G0 ]# Z: j% J" Q1 R) [6 U
    t=0:dt:2;%t axis4
    ) `0 n5 x# W; B; `  }5 I) ]seis=zeros(length(t),length(x));, y$ y: }! d! n7 J( b% u6 F5 s
    % seis=event_hyp(seis,t,x,.4,700,v,1,3);
    : {  H' G. V( ~+ z- M7 W. gflag=1;aper=3;- `9 x: d8 g. h* |/ @8 C* r
    tnot=0.4;% P( B. F5 ^- y5 [
    xnot=700;: m+ V) o% ]( V
    amat=zeros(2000,2000);, G& K! }0 r+ M& L* D. r; t% J7 z. H
    amp=1;3 {( Q; P9 v; B- D! Z8 h' s5 |
    + U8 f7 K% @* _4 X' U' Y6 t
    % if(nargin<8); G, }1 ~5 R0 [9 d
    %         flag=3;
    / T1 f! e# g" o$ f5 p2 ?% end
    . B; p8 L8 @' h4 Y4 {- m2 V, O% if(nargin < 9 )( T. j! _3 Q8 o2 k! Q( s8 G
    %         aper = inf;
    . L. _! r* h6 M  b; U/ [% end
    ! Y: ?) ^! A! h1 I
    & X4 L, v9 v- u; V( M0 \v=v/2;/ u) [# p1 d3 o( W8 ?7 a
    . x% q0 l6 k0 n7 ?7 v1 k2 v& _) s
    %loop over columns8 n% \  ~9 K$ [3 T9 A
    [nsamp,nc]=size(amat);
    ) r% ^6 {9 r( X* j* m* g4 ~2 F# c! d- d, g
    dt=t(2)-t(1);
    # `+ Z  o3 G: Vtmin=t(1);
    4 Q$ m' T2 X/ g, mfor k=1:nc
    % c3 ?) y5 y5 o4 V# `' n2 a, _: a        xoff=x(k)-xnot;) z3 E, g  b0 O- d4 x) Z4 u
            if(abs(xoff) < aper); K$ f! c: Z5 J- O+ g0 i( Z$ x
                    tk = sqrt(tnot^2+(xoff/v)^2);& H& Q% @8 p6 `% v2 w
                    a=amp;
    / @! L( g& m# c/ z                if(flag==2)) v: d( E2 B4 ]0 J2 i
                            a = tnot*a/tk;
    8 A1 k1 s/ Z$ p1 P  r1 |                elseif(flag==3)! g: f  q  C! E5 Z( J4 ], U
                            a = tnot*a/(tk^1.5);0 s* A2 Q' _; O; r
                    elseif(flag==4)5 W1 H/ w+ }5 \; ^; d+ s" a
                            a = tnot*a/(tk*tk);
    0 g7 H$ L. @& o8 {                end  m% g8 N: s2 d2 k# c/ w
                    ik=(tk-tmin)/dt+1;1 _! g/ s) ?# x$ f( e
                    if( between(1,nsamp,ik) )
    9 {# n5 U9 X. u' \& [" |) r+ n                        ik1=floor(ik);0 X6 ^$ Q/ p( n0 l8 Q* `
                            ik2=ceil(ik);
    ; x- h) n, r1 E. F                        if(ik1==ik2)( u2 a) B: |; N
                                    amat(ik1,k)=amat(ik1,k)+a;/ x( t/ G4 F; n& C" c* H/ i
                            else$ O* Z4 \. M3 R+ ^6 U3 B- [; {
                                    amat(ik1,k)=amat(ik1,k)+a*(ik-ik2)/(ik1-ik2);
    + h7 @/ w" i, a4 U( Q* J8 ?                                amat(ik2,k)=amat(ik2,k)+a*(ik-ik1)/(ik2-ik1);& G: E& ^' o& L$ _& X- |& s0 c
                            end( t9 _4 E  c, r- d) y/ b! l
                    end' C% P( l8 e9 \: U6 r& E
            end. ~8 V7 j7 c$ I2 H4 X* L
    end! h, `) B' h. n) l7 K; l" C; {+ A  ^
    figure(1);imagesc(amat);colormap gray
    2 G1 b9 @9 u3 K- m. b
    . |1 O2 P9 h( I" m( t在这段代码里面的between用法是不是有些不正确,因为是学习代码所以不太懂,问问大家8 @( c# v+ s( Y3 `+ @/ \

    该用户从未签到

    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
    & ~5 e0 r, z" v. H! o9 q8 s确实用错了。between函数参数类型应该是datatime类型。
    4 O' }# u) \" H( @, S  B
    请问该怎么修改呢?
    7 k) B- N+ q. }3 {7 y5 v/ X

    该用户从未签到

    4#
    发表于 2022-11-9 13:24 | 只看该作者
    if (between(1,nsamp,ik),这个位置就是if的一个判断逻辑嘛。
    " I* l8 o* \% l. B0 ]问题是你要想知道怎么改,你首先得知道这个位置的判断逻辑是什么。- s8 l% P* Z/ Q2 |; Y  {7 h# {/ f

    - y7 ?; O- u3 C8 `
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

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

    EDA365公众号

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

    GMT+8, 2025-11-23 14:05 , Processed in 0.140625 second(s), 24 queries , Gzip On.

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

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

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