请各位高手帮忙看一下粗体字是什么意思？这是一个函数？

xiaohaha007

% for the given period T, get the index for the constants
period = [0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.075, 0.09, 0.1, 0.12, 0.15, 0.17, 0.2, 0.24, 0.3, 0.36, 0.4, 0.46, 0.5, 0.6, 0.75, 0.85, 1, 1.5, 2, 3, 4, 5];
   
if (length(find(period == T)) == 0)
    index_low = sum(period<T);
    T_low = period(index_low);
    T_hi = period(index_low+1);
   
    [sa_low, sigma_low] = AS_1997_horiz(M, r_rup, T_low, is_soil, fault_type, HW, arb);
    [sa_hi, sigma_hi] = AS_1997_horiz(M, r_rup, T_hi, is_soil, fault_type, HW, arb);
  x = [T_low T_hi];
    Y_sa = [log(sa_low) log(sa_hi)];
    Y_sigma = [sigma_low sigma_hi];
    sa = exp(interp1(x,Y_sa,T));
    sigma = interp1(x,Y_sigma,T);
else
    index = find(period == T);

Happy99

试下并看下其输出
period = [0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.075, 0.09, 0.1, 0.12, 0.15, 0.17, 0.2, 0.24, 0.3, 0.36, 0.4, 0.46, 0.5, 0.6, 0.75, 0.85, 1, 1.5, 2, 3, 4, 5];
T=0.05; period<T, pause
T=0.05; sum(period<T)

xiaohaha007

回复 2 # Happy99 的帖子

呵呵，谢谢你的指点，我试试。

zhouyang664

楼主是网上找的程序吧？看看下面的有没有帮助：
http://www.stanford.edu/~bakerjw/Epsilon/AS_2008_nga.m

xiaohaha007

回复 4 # zhouyang664 的帖子

对，就是你说的那个网站，导师让先看看文章，再看程序，文章大概知道，可惜程序看不太懂，弄了好长时间没什么进展，希望你能给点建议，谢谢你了。

happy

这一部分代码比较简单，主要是检索满足条件的值
要看懂这段代码具体起到一个什么作用，那就需要了解AS_1997_horiz到底是做什么的，以及该函数中的参量到底是什么物理意义，这个应该在文章中有所体现

zhouyang664

我没看过这个文章和代码，我只是好奇这个程序，才搜索了一下，

xiaohaha007

回复 6 # happy 的帖子

好的，谢谢你的帮助！

xiaohaha007

回复 7 # zhouyang664 的帖子

哦，还是要谢你对我提出问题的关注，呵呵.