beamlet3dfig17.m

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

function Beamlet3Dfig17(OpenNewWindow,HaveTitle,LoadData)

%Auxillary functions: 3D beamlets

LoadData=0;
ty=8; tz=8; sy=8; sz=8;
n=16;
if nargin<2
    HaveTitle=1;
end
if nargin<1
    OpenNewWindow=1;
end
if OpenNewWindow
    figure
end

B=zeros(n,n,n,n);
B(ty,tz,sy,sz)=1;
vol=adj_newbt3aa(B);
BB=newbt3aaa(vol);
subplot(2,2,1);
imagesc(squeeze(-BB(8,8,:,:)))
axis square
subplot(2,2,2);
imagesc(squeeze(-BB(:,:,8,8)))
axis square
subplot(2,2,3);
axis square
sb=sort(BB(:));
sb=flipud(sb);
plot(log(sb(1:4000)))
set(gca,'Ylim',[-2 2]);
axis square
W=fwt4_po(BB,1,makeonfilter('Haar',2));
sw=sort(W(:));
sw=flipud(sw);
subplot(2,2,4);
plot(log(sw(1:4000)))
set(gca,'Ylim',[-2 2]);
axis square
colormap bone

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