Mercurial > hg > ishara

view dsp/rot45.m @ 34:c75bb62b90a9

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Imported audio synthesis tools.
author samer
date Sun, 20 Jan 2013 19:05:05 +0000
parents db7f4afd27c5
children
line wrap: on
line source
function B = rot45(A)
% rot45 - rotate square array by 45 degrees
%
% rot45 :: [[N,N]] -> [[N,N]].
%
% Array is SUBSAMPLED, so, you should pre-filter to avoid aliasing.
% Result is padded with zeros in the corners.

n = length(A);

B = diag(A,0)';
pad = [0];

for k=2:2:n-1
	B = [ pad diag(A,k)' pad; B; pad diag(A,-k)' pad];
	pad = [pad 0];
end