frt_zp2.m

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

function WC = FRT_Zp2(x)
% FRT_Zp2: Fast Ridgelet Transform for Zp2
% Usage
%	WC = FRT_Zp2(x)
% Inputs
%	x	p by p array, p prime
% Outputs
%   WC  p by (p+1) array
%
p = size(x,1);
Ave = mean(x(:));
Rad = Radon_Zp2_Fast(x-Ave);
Rad(:,1) = Rad(:,1) + Ave;
WC = zeros(p+1,p+1);
S6 = MakeONFilter('Symmlet',6);
Temp = zeros(p+1,1);
for i=1:(p+1),
  Temp(1:p,1) = Rad(:,i);
  WC(:,i) = FWT_PO(Temp,3,S6); % Emmanuel's Favorite
end

% Copyright (c) 2000 David L. Donoho
  

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