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，纵坐标是什么啊？

第一期 功能磁共振数据分析培训（基础班）

老师您好，对于gretna这个工具箱，学生想请教您几个问题：
1. gretna 是否只能去除头动过大的体积，而不去除被试；
2. 对于网络分析部分的输入矩阵是得到的功能连接矩阵吗，还是Fisher‘z变化后的相关矩阵；
3.对于Modular Interaction这一网络属性中Community Index选项的输入是nodal index这个1-90的一列数吗？
望老师能不吝赐教！

杭州师范大学静息态功能磁共振影像免费培训（第四期）