toon0821.m

来自「老外写的小波变换的工具箱」· M 代码 · 共 35 行

% toon0821 -- Some Cosine Packets in the Frequency Domain
%
	posarray = [ 1 0 64; 2 1 32; 3 2 16; 4 4 8 ; 5 8 4; 6 16 2; 7 32 1];
	sz = size(posarray);
	nr = sz(1);
	n  = 512;
	w = zeros(1,n);
	t = (.5:(n-.5)) ./n;
%
	clf; subplot(111);
	LockAxes([0 .5 0 (nr+1)]); 
	title('Some Cosine Packets in Frequency Domain')
	t = t - .5;
%
	for iter = 1:nr,
	    d = posarray(iter,1);
	    b = posarray(iter,2);
	    k = posarray(iter,3);
	    w = MakeCosinePacket(d,b,k,'Sine',7,n);
	    a = abs(fftshift(fft(w))); a= (.9 ./max(a)) .* a;
	    PlotSpikes(iter,t,a);
	    txt = sprintf('(%1.0f,%2.0f,%2.0f)',d,b,k);
	    text(.43,(iter)+.275,txt);
	end
    
    
 
 
%
%  Part of Wavelab Version 850
%  Built Tue Jan  3 13:20:43 EST 2006
%  This is Copyrighted Material
%  For Copying permissions see COPYING.m
%  Comments? e-mail [e-mail]