s2let_plot_wavs

 s2let_plot_wavs

0001 % s2let_plot_wavs
0002 
0003 load('kappas_spline');
0004 load('kappas_s2dw');
0005 load('kappas_need');
0006 
0007 B = 3;
0008 J_min = 2;
0009 L = 128;
0010 J = s2let_jmax(L, B);
0011 Jmax = 3;
0012 
0013 ns = ceil(sqrt(2+J-J_min+1)) ;
0014 nx = 1;
0015 ny = 3;
0016 
0017 
0018 
0019 J = s2let_jmax(L, B);
0020 xi = 0:0.01:L-1;
0021 x = 0:L-1;
0022 
0023 
0024 figure('Position',[100 100 900 450])
0025 yi = interp1(x,kappa0_spline,xi,'pchip');
0026 semilogx(xi, yi, '-.r', 'LineWidth', 2);
0027 hold on;
0028 yi = interp1(x,kappa0_s2dw,xi,'pchip');
0029 plot(xi, yi, '-k', 'LineWidth', 2);
0030 yi = interp1(x,kappa0_need,xi,'pchip');
0031 plot(xi, yi, '--b', 'LineWidth', 2);
0032 for j = J_min:J  
0033   colour = rand(1,3)*0.9;
0034   yi = interp1(x, kappa_spline(j+1,:), xi,'pchip');
0035   plot(xi, yi, '-.r', 'LineWidth', 2)%, 'Color', colour);
0036   yi = interp1(x, kappa_s2dw(j+1,:), xi,'pchip');
0037   plot(xi, yi, '-k', 'LineWidth', 2)%, 'Color', colour);
0038   yi = interp1(x, kappa_need(j+1,:), xi,'pchip');
0039   plot(xi, yi, '--b', 'LineWidth', 2)%, 'Color', colour);
0040 end
0041 axis([1 L -0.05 1.15]);
0042 set(gca,'XTick',2.^[0:(J+2)]);
0043 hleg1 = legend('B-Spline', 'SD', 'Needlet');
0044 set(hleg1, 'Position', [.15,.25,.1,.2]);
0045 
0046 
0047 [thetas, phis, n, ntheta, nphi] = ssht_sampling(L);
0048 figure('Position',[100 100 900 200]) 
0049 
0050 h = subplot(nx, ny, 1);
0051 hold on
0052 flm = zeros(L^2,1);
0053 for l = 0:L-1
0054     flm(l^2+l+1,1) = kappa0_spline(l+1);
0055 end     
0056 f = ssht_inverse(flm, L, 'Reality', true);
0057 plot(thetas, f(:,1), '-.r', 'LineWidth', 2)
0058 mx = 1.1*max(f(:,1));
0059 axis([0 2. -mx/8 mx ])
0060 flm = zeros(L^2,1);
0061 for l = 0:L-1
0062     flm(l^2+l+1,1) = kappa0_s2dw(l+1);
0063 end     
0064 f = ssht_inverse(flm, L, 'Reality', true);
0065 plot(thetas, f(:,1), '-k', 'LineWidth', 2)
0066    
0067 flm = zeros(L^2,1);
0068 for l = 0:L-1
0069     flm(l^2+l+1,1) = kappa0_need(l+1);
0070 end     
0071 f = ssht_inverse(flm, L, 'Reality', true);
0072 plot(thetas, f(:,1), '--b', 'LineWidth', 2)
0073    
0074    
0075 for j = J_min:Jmax
0076    h = subplot(nx, ny, j-J_min+2);
0077    hold on
0078    flm = zeros(L^2,1);
0079     for l = 0:L-1
0080         flm(l^2+l+1,1) = kappa_spline(j+1,l+1);
0081     end  
0082    f = ssht_inverse(flm, L, 'Reality', true);
0083    plot(thetas, f(:,1), '-.r', 'LineWidth', 2) 
0084     mx = 1.1*max(f(:,1));
0085     axis([0 2. -mx/7 mx ])
0086    flm = zeros(L^2,1);
0087     for l = 0:L-1
0088         flm(l^2+l+1,1) = kappa_s2dw(j+1,l+1);
0089     end  
0090    f = ssht_inverse(flm, L, 'Reality', true);
0091    plot(thetas, f(:,1), '-k', 'LineWidth', 2) 
0092    flm = zeros(L^2,1);
0093     for l = 0:L-1
0094         flm(l^2+l+1,1) = kappa_need(j+1,l+1);
0095     end  
0096    f = ssht_inverse(flm, L, 'Reality', true);
0097    plot(thetas, f(:,1), '--b', 'LineWidth', 2) 
0098 end