calcthermopartition.sci

来自「小波分解源代码」· SCI 代码 · 共 39 行

function z = CalcThermoPartition(cw,mm,q)
// CalcThermoPartition -- Build Thermodynamic Partition Function
//  Usage
//    z = CalcThermoPartition(cw,mm,q)
//  Inputs
//    cw        matrix output by CWT
//    mm        maxmap output by MMWT
//    q         optional list of exponents default linspace(-2,5,61)
//  Outputs
//    z         matrix nexp by nscale of z(q,a)
//
//  Description
//    z(q,a) = sum_i( |CWT(a,b(i))|^q ),  where b = (b(i)) is a list
//    of wavelet transform maxima  
//
//  See Also
//    CWT, WTMM
//
//  Copyright Aldo I Maalouf

  [lhs,rhs]=argn();
	if rhs < 3,
	   q = linspace(-2,5,61);
	end
	
	[n,nscale] = size(cw);
	for k=1:nscale,
	    j = find(mm(:,k));
		if ~isempty(j),
			C = abs(cw(j,k));
			for i=1:length2(q),
				z(i,k) = sum(C .^ q(i));
			end
		else
			z(:,k) = eps^q;
		end
	end
endfunction