
% provede konvoluci dvou obdelniku a konvoluci obdelnikove vlny s Gauss. fci

function sfa2_1


L=10, N=128
del=L/N;
delf=1/N/del;
x=(-N/2:N/2-1)*del; 
f=(-N/2:N/2-1)/N/del;

% konvoluce stejnych obdelniku
g=(abs(x)<1.5);
h=g;

figure(1)
plot(x,g,"*")
xlim([-6,6])
ylim([-0.1,1.1])


%spektra
G=ft(g,del);
H=ft(h,del);

%konvoluce
S=G.*H;
s=ift(S,delf);

figure(2)
plot(x,real(s),"*")
xlim([-6,6])
ylim([-0.3,3.3])






% ted konvoluce s Gauss. fci


% obdelnikova vlna
g=(-1).^floor(x);
figure(1)
plot(x,g,".-")
xlim([-6,6])
ylim([-1.5,1.5])

%Gauss
h=exp(-15*x.^2)/sqrt(pi/15);
figure(2)
plot(x,h,".")
xlim([-3,3])
ylim([-0.3,2.5])

%konvoluce
G=ft(g,del);
H=ft(h,del);
S=G.*H;
s=ift(S,delf);

figure(3)
plot(x,real(s),".-")
xlim([-6,6])
ylim([-1.5,1.5])


endfunction



function G=ft(g,del)
 G=fftshift(fft(ifftshift(g)))*del;
end


function g=ift(G,delf)
 g=fftshift(ifft(ifftshift(G)))*delf*length(G);
end
