%1.1
impz([1 0.74],[1],5,1000)
%1.2
stepz([1 -0.74],[1],5,1000)
%1.3
plot(abs(fft([1 0.74],8000)))
xlabel("Frequency [kHz]")
ylabel("Magnitude")
%1.4
plot(angle(fft([1 0.74],8000)))
xlabel("Frequency [kHz]")
ylabel("Phase [rad]")
%2.1
stem(0:20,filter(1,[1 -0.5],[1; zeros(20,1)]))
%2.2
stepz([1 0.77],[1],21,100)
%2.3
x = 0:8000;
plot(x,(1./(1-2*0.77*cos(2*pi*x*(1/8000))+(0.77).^2).^(1/2)))
%2.4
plot((0:(length(fft(filter([1],[1 -0.77],[1,zeros(1,7999)])))-1))*8000/length(fft(filter([1],[1 -0.77],[1,zeros(1,7999)]))),angle(fft(filter([1],[1 -0.77],[1,zeros(1,7999)]))))
impz([1 0.74],[1],5,1000)
%1.2
stepz([1 -0.74],[1],5,1000)
%1.3
plot(abs(fft([1 0.74],8000)))
xlabel("Frequency [kHz]")
ylabel("Magnitude")
%1.4
plot(angle(fft([1 0.74],8000)))
xlabel("Frequency [kHz]")
ylabel("Phase [rad]")
%2.1
stem(0:20,filter(1,[1 -0.5],[1; zeros(20,1)]))
%2.2
stepz([1 0.77],[1],21,100)
%2.3
x = 0:8000;
plot(x,(1./(1-2*0.77*cos(2*pi*x*(1/8000))+(0.77).^2).^(1/2)))
%2.4
plot((0:(length(fft(filter([1],[1 -0.77],[1,zeros(1,7999)])))-1))*8000/length(fft(filter([1],[1 -0.77],[1,zeros(1,7999)]))),angle(fft(filter([1],[1 -0.77],[1,zeros(1,7999)]))))