%-- 20.10.2010 14:04 --%
imshow(kvetsctsc.m)
imshow(kvet,[0,255],'InitialMagnification','fit');
imshow(kvetsc,[0,255],'InitialMagnification','fit');
imshow(flower,[0,255],'InitialMagnification','fit');
imshow(Cmy343,[0,255],'InitialMagnification','fit');
imshow(Cmy343.bmp,[0,255],'InitialMagnification','fit');
imshow(kvet,[0 255],'InitialMagnification','fit');
imshow(kvetsc,[0 255],'InitialMagnification','fit');
imshow(flower,[0 255],'InitialMagnification','fit');
kvetsc
imshow(kvet,[0 255],'InitialMagnification','fit');
kvetn = kvet/256;
imwrite(kvetn, 'kvet.bmp');
I = imread('kvet.bmp');
I = im2double(I);
[U,S,V] = svd(I)
eol = S(1,1)*U(:,1)*V(:,1)';
imshow(eol,'InitialMagnification','fit');
eo1 = S(1,1)*U(:,1)*V(:,1)';
imshow(eo1,'InitialMagnification','fit');
eo2=S(2,2)*U(:,2)*V(:,2)';
eo3=S(3,3)*U(:,3)*V(:,3)';
imshow(eo2,'InitialMagnification','fit');
imshow(eo3,'InitialMagnification','fit');
eo4=S(4,4)*U(:,4)*V(:,4)';
eo5=S(5,5)*U(:,5)*V(:,5)';
imshow(eo5,'InitialMagnification','fit');
kveteo = eo1+eo2+eo3+eo4+eo5;
imshow(kveteo,'InitialMagnification','fit');
clear all
X = [1;1;1;1;0;0;0;0];
N = 8; % možná 0
Fx = 1 / sqrt(N) * fft(X);
Fx
plot(abs(Fx));
x=[1;1;1;1;1;1;1;1]
Fx = 1 / sqrt(N) * fft(x)
f = ones(30,30);
f(:,13:19) = 0;
imshow(f);
F = fft(f);
imshow(abs(Fs),[0 10],'InitialMagnification', 'fit');
Fs = fftshift(F);
imshow(abs(Fs),[0 10],'InitialMagnification', 'fit');
F = 1/sqrt(30) * fft(f);
Fs = fftshift(F);
imshow(abs(Fs),[0 10],'InitialMagnification', 'fit');
imshow(abs(Fs),[0 50],'InitialMagnification', 'fit');
F = 1/sqrt(30) * fft2(f);
Fs = fftshift(F);
imshow(abs(Fs),[0 10],'InitialMagnification', 'fit');
U = (0.5, 0.5, 0.5, 0.5; 0.5, -0.5i, -0.5, 0.5i; 0.5, -0.5, 0.5, -0.5; 0.5, 0.5i, -0.5, -0.5i);
U = [0.5, 0.5, 0.5, 0.5; 0.5, -0.5i, -0.5, 0.5i; 0.5, -0.5, 0.5, -0.5; 0.5, 0.5i, -0.5, -0.5i]
f=[0,0,0,0;0,1,1,0;0,1,1,0;0,0,0,0]
iFf = conj(U)*Ff*conj(U)
Ff = U*f*U
iFf = conj(U)*Ff*conj(U)
