% 
% Written by:
% -- 
% John L. Weatherwax                2005-08-04
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
% Epage 
% 
%-----

close all; clc; clear; 

X = [ 1, 2, 1; 0, 1, 1; 2, 1, 2 ];

% Select a 2 x 2 matrix to transform
XHat = X(1:2,1:2);

% Generate some Hadamard matrices of various orders:
for n=1:5, 
  H_n = hadamard( 2^n )/(2^(n/2));
end

n=1; 
H_n = hadamard( 2^n )/(2^(n/2));

% Implement the 
Y = H_n * XHat * H_n;
Y