wvd的matlab实现程序,tfrqview.m 源代码在线查看 - WVD变换的matlab工具包 资源下载 虫虫电子下载站...

'cancel'); if layout==1, issig=~issig; if issig==1, if isempty(sig), disp('Impossible action : the signal is unavailable'); issig=0; CallTfrView=0; else sigenveloppe=menu('SIGNAL REPRESENTATION','signal only','signal with enveloppe')-1; CallTfrView=1; end; else isgridsig=0; end; elseif layout==2, isspec=~isspec; if isspec==1, if isempty(sig), disp('Impossible action : the signal is unavailable'); isspec=0; CallTfrView=0; else linlogspec=menu('FREQUENCY REPRESENTATION','linear scale','log scale')-1; CallTfrView=1; end; else isgridspec=0; end; elseif layout==3, iscolorbar=~iscolorbar; CallTfrView=1; elseif layout==4, CallTfrView=0; end; elseif (choice==ChoiceColormap), % Change the color map colmap=menu('COLOR MAP :',... 'hsv','jet','cool','bone','gray','hot','prism',... 'brighten','darken','permute','spin','cancel'); if colmap==1, colormap(hsv); elseif colmap==2, colormap(jet); elseif colmap==3, colormap(cool); elseif colmap==4, colormap(bone); elseif colmap==5, colormap(gray); elseif colmap==6, colormap(hot); elseif colmap==7, colormap(prism); elseif colmap==8, brighten(0.25); elseif colmap==9, brighten(-0.25); elseif colmap==10, c = colormap; colormap(flipud(c)); elseif colmap==11, spinmap; end CallTfrView=0; elseif (choice==ChoiceSampling), % Change the sampling frequency fprintf(' Old sampling frequency: %f\n',fs); fsold=fs; fs=input(' New sampling frequency: '); if isempty(fs), fs=fsold; CallTfrView=0; else CallTfrView=1; end; elseif (choice==ChoiceThreshold), % Change the threshold fprintf(' Old threshold: %f\n', threshold); throld=threshold; threshold=input(' New threshold: '); if isempty(threshold), threshold=throld; CallTfrView=0; else CallTfrView=1; end elseif (choice==ChoiceLinlog), % Change the lin/log scale of tfr linlogtfr=1-linlogtfr; elseif (choice==ChoiceLevels), % Change the number of levels fprintf(' Old number of levels: %f\n',levelnumb); levelold=levelnumb; levelnumb=input(' New number of levels: '); if isempty(levelnumb), levelnumb=levelold; CallTfrView=0; else CallTfrView=1; end elseif (choice==ChoiceGrid), % Grids if ~issig & ~isspec, % No signal and no spectrum isgridtfr=1-isgridtfr; CallTfrView=1; elseif issig & ~isspec, % A signal, no spectrum Grid=1; if ~isgridsig, gridsigstr='add a grid on the signal'; else gridsigstr='remove the grid on the signal'; end if ~isgridtfr, gridtfrstr='add a grid on the TFR'; else gridtfrstr='remove the grid on the TFR'; end Grid=menu('GRID MENU :',gridsigstr,gridtfrstr,'cancel'); if Grid==1, isgridsig=1-isgridsig; CallTfrView=1; elseif Grid==2, isgridtfr=1-isgridtfr; CallTfrView=1; else CallTfrView=0; end elseif ~issig & isspec, % No signal, a spectrum Grid=1; if ~isgridspec, gridspestr='add a grid on the spectrum'; else gridspestr='remove the grid on the spectrum'; end if ~isgridtfr, gridtfrstr='add a grid on the TFR'; else gridtfrstr='remove the grid on the TFR'; end Grid=menu('GRID MENU :',gridspestr,gridtfrstr,'Close'); if Grid==1, isgridspec=1-isgridspec; CallTfrView=1; elseif Grid==2, isgridtfr=1-isgridtfr; CallTfrView=1; else CallTfrView=0; end else % A signal and a spectrum Grid=1; if ~isgridsig, gridsigstr='add a grid on the signal'; else gridsigstr='remove the grid on the signal'; end if ~isgridspec, gridspestr='add a grid on the spectrum'; else gridspestr='remove the grid on the spectrum'; end if ~isgridtfr, gridtfrstr='add a grid on the TFR'; else gridtfrstr='remove the grid on the TFR'; end Grid=menu('GRID MENU :',gridsigstr,gridspestr,gridtfrstr,'cancel'); if Grid==1, isgridsig=1-isgridsig; CallTfrView=1; elseif Grid==2, isgridspec=1-isgridspec; CallTfrView=1; elseif Grid==3, isgridtfr=1-isgridtfr; CallTfrView=1; else CallTfrView=0; end end elseif (choice==ChoiceRedraw), % redraw ? RefreshFigure=1-RefreshFigure; if RefreshFigure==1, CallTfrView=1; end; elseif (choice==ChoiceNewFigure), % new figure figure; CallTfrView=1; elseif (choice==ChoiceSaveResults), % Save the results f=freq*fs; Nmethod=length(method); if comp(1:2)=='PC', DefaultName=[method(4:Nmethod),num2str(Nsig),'.mat']; [name,PathWorkDir] = uiputfile(DefaultName, 'Save As'); else DefaultName=[method(4:Nmethod),num2str(Nsig)]; nameStr=[' Name of the MAT file [',DefaultName,'] : ']; name=input(nameStr,'s'); while (length(name)>8), disp(' The name must have less than 8 characters'); name=input(nameStr,'s'); end if isempty(name), name=DefaultName; end end linlog=linlogtfr+2*linlogspec+4*sigenveloppe; isgrid=isgridsig+2*isgridspec+4*isgridtfr; param = [display,linlog,threshold,levelnumb,Nf2,layout,fs,isgrid]; SavedColormap=colormap; if (nargin TfrQView=['colormap(SavedColormap); tfrqview(tfr,sig,t,method)']; TfrView =['clf;colormap(SavedColormap); tfrview(tfr,sig,t,method,param)']; eval(['save ',PathWorkDir,name,... ' tfr sig t f fs method param SavedColormap TfrView TfrQView']); elseif (nargin==5), TfrQView=['colormap(SavedColormap); tfrqview(tfr,sig,t,method,p1)']; TfrView =['clf; colormap(SavedColormap); tfrview(tfr,sig,t,method,param,p1)']; eval(['save ',PathWorkDir,name, ... ' tfr sig t f fs method param p1 SavedColormap TfrView TfrQView']); elseif (nargin==6), TfrQView=['colormap(SavedColormap); tfrqview(tfr,sig,t,method,p1,p2)']; TfrView =['clf; colormap(SavedColormap); tfrview(tfr,sig,t,method,param,p1,p2)']; eval(['save ',PathWorkDir,name,... ' tfr sig t f fs method param p1 p2 SavedColormap TfrView TfrQView']); elseif (nargin==7), TfrQView=['colormap(SavedColormap); tfrqview(tfr,sig,t,method,p1,p2,p3)']; TfrView =['clf; colormap(SavedColormap); tfrview(tfr,sig,t,method,param,p1,p2,p3)']; eval(['save ',PathWorkDir,name,... ' tfr sig t f fs method param p1 p2 p3 SavedColormap TfrView TfrQView']); elseif (nargin==8), TfrQView=['colormap(SavedColormap); tfrqview(tfr,sig,t,method,p1,p2,p3,p4)']; TfrView =['clf; colormap(SavedColormap); tfrview(tfr,sig,t,method,param,p1,p2,p3,p4)']; eval(['save ',PathWorkDir,name,... ' tfr sig t f fs method param p1 p2 p3 p4 SavedColormap TfrView TfrQView']); elseif (nargin==9), TfrQView=['colormap(SavedColormap); tfrqview(tfr,sig,t,method,p1,p2,p3,p4,p5)']; TfrView =['clf; colormap(SavedColormap); tfrview(tfr,sig,t,method,param,p1,p2,p3,p4,p5)']; eval(['save ',PathWorkDir,name,... ' tfr sig t f fs method param p1 p2 p3 p4 p5 SavedColormap TfrView TfrQView']); end; disp(' '); fprintf('The file is saved in the directory %s\n',PathWorkDir); fprintf('under the name %s\n',name); fprintf('If you want to find again the exact layout of this screen, do\n'); fprintf('load %s%s; eval(TfrView);\n\n', PathWorkDir,name); fprintf('If you want to restart the display session under tfrqview, do\n'); fprintf('load %s%s; eval(TfrQView);\n',PathWorkDir,name); CallTfrView=0; elseif (choice==ChoiceSaveOptions), % Save options save options fs threshold linlogtfr linlogspec levelnumb ... display layout colmap iscolorbar ... isgridsig isgridspec isgridtfr issig sigenveloppe isspec; fprintf('\n Options saved\n'); CallTfrView=0; elseif (choice==ChoicePrint), % Print the current figure Nmethod=length(method); if comp(1:2)=='PC', DefaultName=[method(4:Nmethod),num2str(Nsig),'.eps']; [name,PathWorkDir] = uiputfile(DefaultName, 'Save As'); else DefaultName=[method(4:Nmethod),num2str(Nsig)]; nameStr=[' Name of the EPS file [',DefaultName,'] : ']; name=input(nameStr,'s'); while (length(name)>8), disp('The name must have less than 8 characters'); name=input(nameStr,'s'); end if isempty(name), name=DefaultName; end end eval(['print -deps ',PathWorkDir,name]); fprintf(' The file is saved in the directory %s\n',PathWorkDir); fprintf('under the name %s\n',name); CallTfrView=0; end;end;% good bye. I hope that everything happened fine.