Homec4science
Diffusion Research Scripts (Thomas Bolton) (master) Browse

Scripts_Sharing/ISFC/Utils_Graphs/gspbox/plotting/gsp_plot_signal_spectral.m
8897434f9116master

Last Change
Raw File

gsp_plot_signal_spectral.m
View Options

function gsp_plot_signal_spectral(G,f_hat,param)
%GSP_PLOT_SIGNAL_SPECTRAL Plot a graph signal in the graph spectral domain
% Usage: gsp_plot_signal_spectral(G,signal);
% gsp_plot_signal_spectral(G,signal,param);
%
% Input parameters:
% G : Graph or laplacian eigenvalues
% f_hat : Graph signal in the spectral domain.
% param : Optional variable containing additional parameters.
% Output parameters:
% none
%
% 'gsp_plot_signal_spectral(G,signal)' plots a graph signal in the graph
% Fourier domain.
%
% Warning: to use this function, the Fourier basis of the graph should be
% first computed. You can do it with:
%
% G = gsp_full_eigen(G);
%
% Example:
%
% N = 32;
% G = gsp_path(N);
% G = gsp_compute_fourier_basis(G);
% f = sin((1:N)'*2*pi*4/N);
% fhat = gsp_gft(G,f);
% gsp_plot_signal_spectral(G,fhat);
%
% Additional parameters
% ---------------------
%
% param.highlight : Highlight one frequency component.
% param.plot_abs : Option to plot the absolute value of f_hat
%
% See also: gsp_plot_graph gsp_plot_signal
%
%
% Url: http://lts2research.epfl.ch/gsp/doc/plotting/gsp_plot_signal_spectral.php
% Copyright (C) 2013-2016 Nathanael Perraudin, Johan Paratte, David I Shuman.
% This file is part of GSPbox version 0.5.2
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% This program is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with this program. If not, see <http://www.gnu.org/licenses/>.
% If you use this toolbox please kindly cite
% N. Perraudin, J. Paratte, D. Shuman, V. Kalofolias, P. Vandergheynst,
% and D. K. Hammond. GSPBOX: A toolbox for signal processing on graphs.
% ArXiv e-prints, Aug. 2014.
% http://arxiv.org/abs/1408.5781
% Author : David I Shuman, Nathanael Perraudin
% Testing: test_plotting
% Read input parameters
if nargin < 3
param = struct;
end
if ~isfield(param, 'highlight'), param.highlight = 0; end
if ~isfield(param, 'plot_abs'), param.plot_abs = 1; end
if param.plot_abs
f_hat=abs(f_hat);
end
if isstruct(G)
E = G.e;
end
cla;
hold on;
stem(E,f_hat,'LineWidth',1)
if param.highlight>0
stem(E(highlight),f_hat(highlight),'filled','r','LineWidth',1);
end
tic_int=ceil(max(E)/10);
set(gca,'XTick',0:tic_int:ceil(max(E)));
set(gca,'color','none')
if sign(min(f_hat))>0
aymin = 0.9*min(f_hat);
else
aymin = 1.1*min(f_hat);
end
aymax = max(1.1*max(f_hat),1.1*min(f_hat)+eps);
axis([-.1 ceil(max(E)) aymin aymax]);
% set(gca,'fontsize',18)
box on;
hold off;
end
c4science · Help