drdrawfig33.m

来自「beamlet变化的工具箱」· M 代码 · 共 104 行

function dRDrawFig33(OpenNewWindow,HaveTitle,LoadData)

%In the second row, four subpanels: first 3 are easy to
%make: they show: a ridgelet, a ridgelet on n*2n grid,
%a horizontal ridge function.

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end
LoadData=0;

subplot(1,4,1);
subplot(1,4,2);
subplot(1,4,3);
subplot(1,4,4);


%In the second row, four subpanels: first 3 are easy to
%make: they show: a ridgelet, a ridgelet on n*2n grid,
%a horizontal ridge function.
n1 = 64;
n2 = 128;

%We start with a Meyer wavelet on 2n points. 

wave = MakeWavelet(4,8,'Meyer',2,'Mother',128);
BigWave=wave';
%We then create an 2n*n array where each column contains a copy of
%that wavelet. That gives us panel (c). 
for k=1:63;
    
    BigWave=[BigWave wave'];
end
subplot(1,4,3)
imagesc(abs(BigWave));axis image;set (gca,'XTick',[]);
    set (gca,'YTick',[]);


     
%We then create an n*2n array in which Panel (c) has been sheared 23 degrees.
%That gives us Panel(b).  

theta2=23*pi/180;
%

%
% Build matrices of Difference
for i=1:n1,
	vshift2 = Shift1dSignal(BigWave(:,i),(i-n1/2-1).*sin(theta2));
	ExactShift2(:,i) = real(vshift2(:));
end

subplot(1,4,2)
imagesc(abs(ExactShift2));axis image;set (gca,'XTick',[]);
    set (gca,'YTick',[]);



%We then create an n*n array by mutilation,
%i.e. dropping the outer n rows -- n/2 on top and n/2 on bottom.

ExactShift2_crop=ExactShift2(n2/4+1:3*n2/4,:);

subplot(1,4,1)
imagesc(abs(ExactShift2_crop));axis image;set (gca,'XTick',[]);
    set (gca,'YTick',[]);


S2=sum(BigWave,2);

x2=1:length(S2);
subplot(1,4,4)
plot(S2,x2, 'k','LineWidth',2);ylim([1,128])
axis off
set(gca,'PlotBoxAspectRatio',[1 2 1]);set(gca,'YDir','reverse')


% Copyright (c) 2002 Ana Georgina Flesia
%
% Part of BeamLab Version:200
% Built:Friday,23-Aug-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%
%
% Part of BeamLab Version:200
% Built:Saturday,14-Sep-2002 00:00:00
% This is Copyrighted Material
% For Copying permissions see COPYING.m
% Comments? e-mail [e-mail]
%