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