【17.09.04/05】再见FFT，你好SF

果然要每天汇报实在是太困难了，尤其像现在开始正式上课之后，我争取2-3天一总结，至少保证一周一汇报吧\_(:зゝ∠)\_

# 每日研究总结
这两天主要干了3件小事：
 - 修正了之前程序里的两个重大bug
 - 确认了FFT方法不能用于M81光变曲线的分析
 - 思考了Structure Function表示的含义

##自己搬的砖
我为了要实现“按照指定的index顺序给列表排序”这个功能，一开始写了个函数：
```python
def rearrange(target,index):
    result = [ ]
    for i in index:
        result.append(target[i])
    return result
```

但发现这样输出列表的顺序并不是按照我指定的index列表中的顺序来的，直到我把改成了这样：
```python
def rearrange(target,index):
    result = [ ]
    n = len(index)
    for i in range(n):
        j = index.index(i)
        result.append(target[j])
    return result
```
我猜测是不是for循环中`for i in list`时i不一定是按照顺序从list中取的，只有`for i in range(n)`才是按顺序递增？（_待求证_）

----------
然后是在考察光变曲线的FFT时，除了之前提到的等距采样（见[17.09.02](http://blog.leanote.com/post/fangzheng42/research_170902)）外，我居然忘了采样定理这回事。一开始做出来的FFT频率只对应到几赫兹，惨不忍睹😅。
但是转念一想，如果我要提高采样频率的话，就是要在插值后的函数中比原始数据点多取好多点，对于像这样10年的观测光变曲线来说似乎是非常不合适的，因此只好放弃啦）
（所以说你导师永远是你导师\_(:зゝ∠)\_，一开始就没让用FFT，应该听话的(\*/ω╲\*)）

----------
关于Structure Function嘛，虽然我不太清楚它到底有什么统计意义，不过之前导师说有周期性的话一定会在图上有显示，于是我就将信将疑地试了试做周期为2正弦信号地SF图，似乎对应的周期时SF图上会出现一个个dip，嗯……当然这里的SF还没有经过rebin和统计平均，当然这是这周要做的事啦
![sf_sin](https://leanote.com/api/file/getImage?fileId=59af505bab64411cbe000b50)

##别人盖的楼
这两天和室友讨论到了一个问题：我们经常看到说AGN的光变时标很短，有时候只有几分或者几小时，说明发出辐射的源的大小也只有 几光分或者几光时。
以前都是直接接受这种观点，但是确实仔细一想，为什么光变时标能表征大小呢？
我隐约有点感觉是和响应有关，但又说不太清的感觉，还是要想想……
#不务正业
终于拿学邮搞到了GitHub Student Pack，愉快地在digital ocean上搞了个服务器，又可以过上离线下载，远程登电脑，云计算，以及搞各种事的生活啦！

方正42

导航

最近发表

友情链接