drdrawfig26.m

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

function dRDrawFig26(OpenNewWindow,HaveTitle,LoadData)

% if you don't want to process it, set LoadData=1

if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end

LoadData=0;
% DrawFig26.m: Inspect some Dual Ridgelets (Old Code)


parlist = [ 4  0 2  2 1 ; ...
			4  14 2  0 1 ];

%figure;
fignum = 11;
n=64;
npar = 2;
str = 'abcdfegijk';
for pp = 1:npar,
	par = parlist(pp,:);
	theta = zeros(n,2*n);
	theta(dyad2ix(par(1),par(2)),dyad2ix(par(3),par(4))+(par(5)-1)*n) =1;
	ridgelet = t_FastRidgeXform(theta);
	subplot(1,2,pp)
	imagesc(abs(ridgelet));
	axis('square');set (gca,'XTick',[]);
    set (gca,'YTick',[]);
    if HaveTitle
	title(sprintf('(%s): Dual Ridgelet s=%i;(j,k)=(%i,%i);(i,l)=(%i,%i)',...
		str(pp),par(5),par(1),par(2),par(3),par(4)));
end
	fignum = fignum+1;
end
%print -deps Fig_11_12.eps

% Copyright (c) 2002 David Donoho and 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]
%