|
输入:EEG_file:脑电图信号文件韦尔奇:你想要的长度seg_len:每个段的长度Fs:采样率 |
|
输出:Wel_psd: PSD的Welch估计开始（1）初始化:MPI_Init;排名←MPI_Comm_rank;大小←MPI_Comm_size;overlap_len = seg_len / 2;（2）读EEG_file（3)EEG_file中的数据存储到数组welch[n]；（4)计算有多少段:n_segs（5）for node = 1 to size pardo //并行进程(6）将汉明视窗应用于韦尔奇[n］/ /n间startPos stopPos);startPos和stopPos分别是某个数据段的起始点和结束点。（7）对加窗数据进行FFT。结果存储在变量fft_x中;（8）计算fft_x振幅的平方。结果存储在变量程序中;（9)求和功率谱Pxy + =程序;(10）结束了(11）if (rank = = 0)//主进程（12）从每个过程中接收Pxy，并计算所有周期图的平均值。（13）/ /其他奴隶的过程(14）发送Pxy到主进程结束 |
|
