Skip to content
Snippets Groups Projects
Commit 13115594 authored by Olivier Sauter's avatar Olivier Sauter
Browse files

add spectrogram of sum(n=i) modes

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@5214 d63d8f72-b253-0410-a779-e742ad2e26cf
parent ec3c40f0
Branches
Tags
No related merge requests found
......@@ -80,16 +80,24 @@ if prod(isfield(gdat_data,{'data','t'})) && ~isempty(gdat_data.data) && ~isempty
delete(hylabel);
legend(gdat_data.label,2);
% add spectrogram per signal
mhd_sum_data = 0.;
nfft=1024;
tmhdm=mean(reshape(gdat_data.t(1:nfft*fix(length(gdat_data.t)/nfft)),nfft,fix(length(gdat_data.t)/nfft)));
for i=1:size(gdat_data.data,2)
nfft=1024;
tmhdm=mean(reshape(gdat_data.t(1:nfft*fix(length(gdat_data.t)/nfft)),nfft,fix(length(gdat_data.t)/nfft)));
[B,F,T]=specgram(gdat_data.data(:,i),nfft,1/mean(diff(gdat_data.t)),hanning(nfft),nfft/2);
figure
imagesc(T+tmhdm(1),F/1e3,20*log10(abs(B)));axis xy;colormap jet;
ylabel('freq')
xlabel(gdat_data.dimunits{1})
title(gdat_data.label{i})
mhd_sum_data = mhd_sum_data + gdat_data.data(:,i);
end
[B,F,T]=specgram(mhd_sum_data./size(gdat_data.data,2),nfft,1/mean(diff(gdat_data.t)),hanning(nfft),nfft/2);
figure
imagesc(T+tmhdm(1),F/1e3,20*log10(abs(B)));axis xy;colormap jet;
ylabel('freq')
xlabel(gdat_data.dimunits{1})
title(['sum of ' gdat_data.label{:}])
end
else
disp('cannot plot gdat_data, has empty data or t field')
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment