f7_11.m

来自「工具箱提供了各种常用的信号处理」· M 代码 · 共 28 行

%读入EEG信号
load eeg.dat;
sig=eeg(2,1:512)';
N=length(sig);
%采样频率
fs=312.5;

%时域波形
figure(1);
plot((1:N)/fs,sig,'LineWidth',2);
xlabel('时间 t/s');
ylabel('电压 A/uV');

%计算谱图
[tfr,t,f]=tfrsp(sig);

%显示
figure(2);
contour(t/fs,f(1:N/2)*fs,tfr(1:N/2,:));
xlabel('时间 t/s');
ylabel('频率 f/Hz');

figure(3)
mesh(t/fs,f(1:N/2)*fs,tfr(1:N/2,:));
xlabel('时间 t/s');
ylabel('频率 f/Hz');
zlabel('幅值 A/ uV ');