matlab 功率谱密度 汉宁窗_功率谱密度PSD求法

大致知乎查了下PSD与振幅谱的区别，因为在做引力波数据分析时，不论matched filtering 还是whitening(其实，这个就是matched filter所做的信号还原）都要用到PSD。

大概看了知乎，大概理解为是噪声振幅谱平方的期望即为PSD。

那么我用两种方法来实验：

1.采用多次振幅谱平方求均值作为PSD

对于引力波GW150914（数据可以在官网events页面下载）：

'''数据列表
data_et,    4096s数据
time_et,    4096s时间
template    32s  模板
data        32s  数据
time        32s  时间
freq        32s  频率
S_n         32s  噪声谱
df          32s  频域间隔
dt          32s  时域间隔
times(a,b,S,df)  内积
overlap(a,b,S,df)
rho_sq(h,s,df)