drpdrawfig10.m

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

function dRPDrawFig10(OpenNewWindow,HaveTitle,LoadData)

%DrawFig 9:   basis functions derived from tiles of different heights


if nargin<3
    LoadData=0;
end

if nargin<2
    HaveTitle=1;
end

if nargin<1
    OpenNewWindow=1;
end

if OpenNewWindow
    figure
end


LoadData=0;
n=64;
wplist = [2 3 2; 2 2 1 ; 3 2 1];% 4 3 2 ];
sz = size(wplist); nr = sz(1); 
%figure;
s = [19 -19  -12];
for i=1:nr,
	d = wplist(i,1);
    b = wplist(i,2);
    k = wplist(i,3);	
    wp = MakeWaveletPacket(d,b,k,'Coiflet',3,n);
    wp = circshift(wp,[1 s(1,i)]);
    t = (0:(n-1)) ./n;
    wFather = zeros(1,64);
    wFather([1 9]) =[1 0];
    QMFSymm8 = MakeONFilter('Symmlet',8);
    x = IWT_PO(wFather,4,QMFSymm8);
    B = repmat(wp',1,64);
    A = repmat(x,64,1);
    C=A.*B;
    S=zeros(128,128);
    S(32:95,32:95) =C; 
    X =  Inv_FastSlantStack(S);
    subplot(2,3,i)
	plot(t,wp,'LineWidth',2);axis square;set(gca,'YDir','reverse');axis off
    subplot(2,3,i+3);
    imagesc(abs(X));axis image;set (gca,'XTick',[]);
   set (gca,'YTick',[]);

end
%
% 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]
%