scinit.m

老外写的小波变换的工具箱

function SCInit
% SCInit -- Setup global datastructures for ShortCourse
%
global Blocks Bumps HeaviSine Doppler
global xblocks xbumps xheavi xDoppler
global yblocks ybumps yheavi yDoppler
global t N rho alf
%
	if isempty(N),
	     N = 2048;
	end
	if isempty(rho),
	     rho = 7;
	end
	if isempty(alf),
	      alf = .9;
	end
	t = (0:(N-1))/N;
%
	disp('Initializing Parameters:')
	disp(sprintf('N = %g',N))
	disp(sprintf('rho = %g',rho))
	disp(sprintf('alf = %g',alf))
%
	HeaviSine = MakeSignal('HeaviSine',N);
	Bumps     = MakeSignal('Bumps',N);
	Blocks    = MakeSignal('Blocks',N);
	Doppler   = MakeSignal('Doppler',N);
%
	[xblocks,yblocks]   = NoiseMaker(Blocks,rho);
	[xbumps,ybumps]     = NoiseMaker(Bumps,rho);
	[xheavi,yheavi]     = NoiseMaker(HeaviSine,rho);
	[xDoppler,yDoppler] = NoiseMaker(Doppler,rho);
%
global yQuad yMish vBumps vDoppler tt
global xQuad xMish uBumps uDoppler
global rhoTFA
global nTFA
%
	if isempty(nTFA),
	     nTFA = 1024;
	end
	if isempty(rhoTFA),
	     rhoTFA = 7;
	end
	tt = (0:(nTFA-1))/nTFA;
%
	disp(sprintf('nTFA = %g',nTFA))
	disp(sprintf('rhoTFA = %g',rhoTFA))
%
	QuadChirp = MakeSignal('QuadChirp',nTFA);
	MishMash  = MakeSignal('MishMash' ,nTFA);
	wBumps    = MakeSignal('Bumps'    ,nTFA);
	wDoppler  = MakeSignal('Doppler'  ,nTFA);
%
	[xQuad,yQuad]       = NoiseMaker(QuadChirp,rhoTFA);
	[xMish,yMish]       = NoiseMaker(MishMash,rhoTFA);
	[uBumps,vBumps]     = NoiseMaker(wBumps,rhoTFA);
	[uDoppler,vDoppler] = NoiseMaker(wDoppler,rhoTFA);

    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:42 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]