f7_11.m

matlab时频分析配套源码第6章的程序

%读入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 ');