ZangYF Group's CVs

将时间序列曲线转换为功率谱密度曲线

老师好,我得到了被试的时间序列,想转换为功率谱密度曲线,请问具体怎么操作?谢谢!
我通过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,纵坐标是什么啊?