Home > src > main > matlab > flaglet_demo1.m

flaglet_demo1

PURPOSE ^

flaglet_demo1

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

 flaglet_demo1
 Plot wavelets functions

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 % flaglet_demo1
0002 % Plot wavelets functions
0003 
0004 s = 0.3;
0005 
0006 R = 1.0 ;
0007 Rplot = 0.39 ;
0008 
0009 L = 92 ;
0010 N = 92 ;
0011 P = 128 ;
0012 
0013 scales_l = [4 5];
0014 scales_n = [5 6];
0015 
0016 J_l = numel(scales_l);
0017 J_n = numel(scales_n);
0018 
0019 
0020 B_l = 2 ;
0021 B_n = 2 ;
0022 J_min_l = 0 ;
0023 J_min_n = 0 ;
0024 
0025 
0026 % Translation
0027 K_n_s = zeros(N,L)  ;
0028 birs = [ s, R ];
0029 for n = 1:N
0030     fn = zeros(1,N) ;
0031     fn(n) = 1.0 ;
0032     [K_n_s_temp, ~] = slag_synthesis(fn, 'Nodes', birs) ;
0033     K_n_s(n,:) = K_n_s_temp(1) ;
0034 end
0035 
0036 [kappa kappa0] = flaglet_axisym_tiling(B_l, B_n, L, N, J_min_l, J_min_n);
0037 
0038 figure('Position',[100 100 750 800])
0039 for jl=scales_l
0040     for jn=scales_n
0041         ind = (find(jl == scales_l)-1)*(J_n) + find(jn == scales_n);
0042         subplot(J_l, J_n, ind)
0043         temp = kappa{jl,jn} ;
0044         temp = temp .* K_n_s ; 
0045         flaglet_plot_axisym_wavelet_kernel(temp, R, Rplot, P)
0046         %title(['jl=',int2str(jl),' jn=',int2str(jn)])
0047         caxis([-1 1])
0048         %xlabel(['x',int2str(ind)])
0049     end
0050 end
0051 colormap(hot(256))
0052 
0053

Generated on Mon 24-Sep-2012 12:26:33 by m2html © 2005