Question

脑电波信号处理MatLab代码模板，小白看不懂，跪求高手解读，谢谢……

请问以下这段代码里面，load数据以后，怎么就忽然冒出了A和B和C，为什么不能直接进行进行FFT变换。A\B\C分别都是什么呢？然后i和j的值是怎么取的呢？是依据什么让它们分别等于64和100呢…………小白第一次接触matlab，跪求大神指点迷津………
基于MATLAB的代码程序
% 脑电信号的特征分析
function NaoDian_Analysis
clc;
close all;
data0 = load('Competition_test');
A = data0.X;
B = A(:,:,1);
% size(B)
C = zeros(6400,1);
for i = 1:100
for j = 1:64
C(i*j) = B(i,j,1); % 提取脑电信号
end
end
plot(C);title('原始脑电信号时域图');xlabel('time');ylabel('幅值');
% 求频谱
fs=1000; %自己设置采样频率
N=3000; %自己设置采样点数
NFFT = 2^nextpow2(N);%转化为2的基数倍
f=fs/2*linspace(0,1,NFFT/2); %求出FFT转化频率
E_change=fft(C,NFFT)/N; %进行FFT变换
figure
plot(f,2*abs(E_change(1:NFFT/2)));title('脑电信号频域图');xlabel('Frequre');ylabel('频谱值');

% 求功率谱
window=boxcar(length(C)); %矩形窗
nfft=1024;
[Pxx,f]=periodogram(C,window,nfft,fs);%直接法
figure
plot(f,10*log10(Pxx));title('脑电信号功率域图');xlabel('Frequre');ylabel('功率谱值');

Answer 1

我解释下：

data0 = load('Competition_test');%data0是个struct结构，包含X域
A = data0.X;%把X域赋给A
B = A(:,:,1);%取A的第一页给B
% size(B)%如果前面的“%”去掉，会看到B的大小，可能是100*64
C = zeros(6400,1);%生成全0矩阵，
for i = 1:100
for j = 1:64
C(i*j) = B(i,j,1); % 提取脑电信号
end
end