% Lab5c.m
Tellerkoeff=[1 1]	;		% Koeffisienter i teller [a0 a1]
Nevnerkoeff=[1 -0.9];	    % Koeffisienter i nevner [ b0 b1]
Fs=8000;				      % Samplingsfrekvens		
poler=roots(Nevnerkoeff)			
nullpunkt=roots(Tellerkoeff)
zplane(Tellerkoeff,Nevnerkoeff)	% Poler og nullpunkt
pause
% Oversiktplot av forsterkning/fase
N=128;		                           % Antall punkter	
freqz(Tellerkoeff,Nevnerkoeff,N,Fs)
pause
% Plotting av forsterkningskurve mellom 0 og 500 Hz
f=[0:500];
[H,f]=freqz(Tellerkoeff,Nevnerkoeff,f,Fs);
plot(f,20*log10(abs(H)))
xlabel('Frekvens [Hz]'),ylabel('Forsterkning [dB]')
pause
%Utskrift av enkeltverdier av forsterkning og fase
f=[0 500 1000 2000 3000 4000]                             % Frekvensene 0  og 1000 Hz 
H=freqz(Tellerkoeff,Nevnerkoeff,f,Fs); % Regner ut forst/fase ved frekvensene i f
Forsterkning=20*log10(abs(H))          % Regner om til dB.  NB:skrivemåten log10!
Fase=angle(H)
pause
