Použité podprogramy pro MATLAB:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=celek(x)
%
% Demonstrace zlepseni pomeru signal/sum
% upravou jeho spektra
%
% generovani zasumeneho signalu pomoci generatoru
% nahodnych cisel:
sig=gener(0);
figure(2);
% vytvoreni spektra signalu:
spektr=spektrum(sig);
h=[1:1:2];
a(1)=3;
figure(3);
% Uprava spektra odstranenim spodni casti
orezanej=orez(spektr,20);
figure(4);
% zpetna transformace spektra
y=fzpet(orezanej,10,1,512,100);
t=[1:1:50];
plot(t,y(1:50));
xlabel('time (milliseconds)');
title('Output signal');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=gener(x)
%
%Generuje a nakresli zasumely signal
% pomoci generatoru nahodnych cisel
%
t = 0:0.001:0.6;
x = sin(2*pi*50*t);
%+sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(1000*t(1:50),y(1:50));
title('Signal Corrupted with Zero-Mean Random Noise');
xlabel('time (milliseconds)');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=spektrum(x)
%
% vygeneruje a zobrazi spektrum signalu
%
Y = fft(x,512);
Pyy = Y.* conj(Y) / 512;
f = 1000*(0:256)/512;
plot(f,Pyy(1:257));
title('Frequency content of y');
xlabel('frequency (Hz)');
y=Pyy; 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=orez(x,h)
%
% orizne spodek koeficientu podle x
%
y = [1:1:512];
for i=1:512
    if x(i)<h y(i)=0; 
    else y(i)=x(i)-h;
    end;
end;

t = [1:1:256];
plot(t,y(1:256));
title('Frequency content of y');
xlabel('frequency (Hz)');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function y=fzpet(b,omega,min,max,kmax)
%sestavi puvodni funkci sectenim cosinove F. rady
% b je matice koeficientu
% omega 
% min je pocetatek casu - setiny vteriny (casu)
% max je konec casove osy - setiny vteriny (casu)
% celkem je m=max-min casovych bodu

m=max-min;
for i=1:m
% y(i)=b(1)/2;
y(i)=0;
for j=1:2*kmax+1
    n=j-kmax;
% y(i)=y(i)+b(j+1)*cos(j*(min+i)*omega/100);
y(i)=y(i)+b(j)*cos(n*(min+i)*omega/100);

end
end