有一个脑电数据，是三维数组，存为一个2×701×100的矩阵，每个维度分别对应：导联（2个导联）×信号（701个采样点）×试次（100个试次，即重复刺激100次）；第1导联为FCz电极，第2导联为Cz电极；对其第1导联的数据，对其100个试次的脑电信号，分别进行fft频谱分析，提供matlab程序

好的，这是一个简单的matlab程序，可以实现你的需求：

% 加载数据
load('data.mat'); % 假设数据保存在 data.mat 文件中

% 取出第1导联的数据
data = squeeze(data(1,:,:)); % 取出第1个维度，即第1导联

% 采样频率为 1000Hz
fs = 1000;

% 对每个试次的信号进行fft
fft_data = fft(data,[],2);

% 取出频率范围在0-1000Hz的部分
fft_data = fft_data(:,1:(size(data,2)/2+1));

% 计算每个试次的功率谱密度
psd_data = abs(fft_data).^2 / (size(data,2)*fs);

% 计算平均功率谱密度
mean_psd = mean(psd_data,1);

% 绘制波形
f = linspace(0,fs/2,size(data,2)/2+1);
plot(f,mean_psd);
xlabel('Frequency (Hz)');
ylabel('Power/Frequency (dB/Hz)');
title('Power Spectral Density');

请将数据保存在 `data.mat` 文件中，然后运行该程序即可得到波形。