rand_wpe.m

function f=rand_wpe(varargin)
%RAND_HPE  Random WPE even function.
%   Usage:  f=rand_wpe(s);
%
%   RAND_WPE(s) generates an array of size s, which is WPE along the
%   first dimension.

if isnumeric(varargin);
  s=varargin;
else
  s=cell2mat(varargin);
end;


if length(s)==1
  error('To avoid confusion, the size must be at least two-dimensional.');
end;

shalf=s;

shalf(1)=floor((shalf(1)-1)/2);

f=(randn(shalf)-.5)+(randn(shalf)-.5)*i;

if rem(s(1),2)==0
  f=[randn([1 s(2:end)])-.5; ...
     f; ...
     randn([1 s(2:end)])-.5; ...
     flipud(conj(f))];
else
  f=[randn([1 s(2:end)])-.5; ...
     f; ...
     flipud(conj(f))];
end;