function [margt,margf,E]=margtfr(tfr,t,f) %MARGTFR Marginals and energy of a time-frequency representation.%	[MARGT,MARGF,E]=MARGTFR(TFR,T,F) calculates the time and%	frequency marginals and the energy of a time-frequency%	representation. %%	TFR : time-frequency representation (M,N)%	T   : vector containing the time samples in sec. %	       (default : (1:N))%	F   : vector containing the frequency samples in Hz, not%	       necessary uniformly sampled. (default : (1:M))%	MARGT : time marginal%	MARGF : frequency marginal%	E     : energy of TFR%%	Example :    %	 S=altes(128,0.05,0.45); TFR=tfrscalo(S,1:128,8,'auto');%	 [MARGT,MARGF,E] = margtfr(TFR); %	 subplot(211); plot(T,MARGT); subplot(212); plot(F,MARGF);%%	See also MOMTTFR.%	P. Goncalves, October 95%	Copyright (c) 1995 Rice University%%	------------------- CONFIDENTIAL PROGRAM -------------------- %	This program can not be used without the authorization of its%	author(s). For any comment or bug report, please send e-mail to %	f.auger@ieee.org [M,N] = size(tfr) ;if nargin == 1, t = (1:N); f = (1:M)';elseif nargin == 2, f = (1:M)';end  E     = real(integ2d(tfr,t,f)/M) ;margt = real(integ(tfr.',f')/M) ;margf = real(integ(tfr,t)/N) ;