-
Notifications
You must be signed in to change notification settings - Fork 1
/
my_dct2.m
39 lines (34 loc) · 918 Bytes
/
my_dct2.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
function b=my_dct2(arg1)
% beta(1)=1/sqrt(2);
% beta(2:128)=1;
% ImF(1:36,1:36)=0;
% N=36;
% for k=0:N-1
% for l=0:N-1
% for i=1:N
% for j=1:N
% ImF(k+1,l+1) = ImF(k+1,l+1) + (1/sqrt(2*N))*beta(k+1)*beta(l+1)*Im(i,j)*cos(pi*(2*j+1)*(k+1)/(2*N))*cos(pi*(2*i+1)*(l+1)/(2*N));
% end
% end
% end
[m, n] = size(arg1);
% Basic algorithm.
if (nargin == 1),
if (m > 1) && (n > 1),
b = dct(dct(arg1).').';
return;
else
mrows = m;
ncols = n;
end
end
% Padding for vector input.
a = arg1;
if nargin==2, ncols = mrows(2); mrows = mrows(1); end
mpad = mrows; npad = ncols;
if m == 1 && mpad > m, a(2, 1) = 0; m = 2; end
if n == 1 && npad > n, a(1, 2) = 0; n = 2; end
if m == 1, mpad = npad; npad = 1; end % For row vector.
% Transform.
b = dct(a, mpad);
if m > 1 && n > 1, b = dct(b.', npad).'; end