% SOE313 Digital signalbehandling  Terje Natås
% Høsten 2004  
% Sinc-funksjon Fig. 3.11 ---------------
N=40;		   	        % Antall stolper 
d=0.2;                  % Nyttefaktor
p=-N/2:1:N/2-1;
y=sinc(d*p);
bar(p,y,0.1)            %Plotting. Stolpefunksjon
axis([-20 20 -0.25 1])
xlabel('p')
grid; pause
% Likerettet Sinc-funksjon  Fig 3.12---------
bar(p,abs(y),0.1)
axis([-20 20 -0.25 1])
xlabel('p')
grid; pause
%---------------------------------
% Sinc-funksjon. Diskret funksjon av f.  Fig. 3.19
A=5; tau=1; T=5;
d=tau/T;                % Nyttefaktor    
p=-N/2:1:N/2-1;
y=A*d*sinc(p*d);
f=p/T;                  % Punkter på frekvensaksen
bar(f,y,0.1)            % Plotting. Stolpefunksjon
axis([-4 4 -0.25 1])
xlabel('Frekvens f [Hz]')
grid; pause
% Sinc. Diskret funksjon av f.  Fig. 3.20 
T=10;                   % Dobler perioden
N=80;                   % Dobler antall plottinger
d=tau/T;
p=-N/2:1:N/2-1;
y=A*d*sinc(p*d);
f=p/T;
bar(f,y,0.1)
axis([-4 4 -0.25 1])
xlabel('Frekvens f [Hz]')
grid; pause
% Sinc. Diskret funksjon av f Fig 3.21
T=20                    % Dobler perioden
N=160;                  % Dobler antall plottinger
d=tau/T;
p=-N/2:1:N/2-1;
y=A*d*sinc(p*d);
f=p/T;
bar(f,y,0.1)
axis([-4 4 -0.25 1])
xlabel('Frekvens f [Hz]')
grid
pause
%----------------------
% Sinc. Kontinuerlig funksjon av f. Fig 3.25
N=160;tau=1;A=1;
T=20; d=tau/T;
p=-N/2:1:N/2-1;
y=A*tau*sinc(p*d);
f=p/T;
plot(f,y)
axis([-4 4 -0.25 1])
xlabel('Frekvens f [Hz]')
grid
pause
