January 2019

老师好，我得到了被试的时间序列，想转换为功率谱密度曲线，请问具体怎么操作？谢谢！
我通过matlab这样做不知道对不对：
x=dlmread('*timecourse.txt')；%一个被试的时间序列
y=fft(x); %进行傅里叶转换
N=1060; %数据共有1060 timepoints
fs=1./0.392; %采样频率计算为：1./TR，我数据的TR为0.392s
f=(0:N-1)*fs/N; %计算数据的频域区间f
mag=abs(y); %计算数据的频谱幅值mag
%%%因为从f 向量内发现[ 3:333 ]这个范围频率为0.01-0.08Hz（我只分析这个频率范围），因此取此区间的数值%%%
f1=f (3:333);
mag1=mag(3:333);
power_mag1=mag1.^2; %求数据的功率值
plot (f1,power_mag1) %展示功率谱曲线

此外，我做出来的这张图，横坐标是Hz，纵坐标是什么啊？