FINAL EXAM.docVIP

FINAL EXAM PROBLEM 1a: Using the following code in MatLab: Num = 105; Den = conv([1, 5.792, 9.14], [1, 4.2, 11.847]) % Part A: Plot the MAG frequency response of Ha(s) % At what frequency will the magnitude drop by 1dB from its dc value W =[0: 0.001: 4];% Omega vector Ha=freqs(Num, Den, W);% Finds the frequency response magHa=abs(Ha); magHadB=20*log10(magHa); figure(1) plot(W, magHadB); % Plots the Magnitude response in dB xlabel(Frequency in rad/sec); ylabel(Amplitude in dB); title(Final - Part 1a: Magnitude FREQ Response of the Filter in dB) omega1 = find(magHadB -1.001 magHadB -1) omega1value = W(omega1) I obtained figure : I also obtained the following Outputs: Den = 1.0000 9.9920 45.3134 107.0058 108.2816 omega1 = 1093 omega1value = 1.0920 PROBLEM 1b: Using the following code in MatLab: figure(2) plot(W,angle(Ha)); xlabel(Freq in rad/sec); ylabel(AMP); title(Phase response); I obtained the following figure. PROBLEM 1c: Using the following code in MatLab: %problem 1c % first find T % Find T by Hand Calculations % W1 = pi/4; % Given value % 1.092 = (2/T)*tan(.25*pi/2) % T = (2/1.092)*tan(.25*pi/2) = .76 T = 0.76; % Use bilinear transformation to find x(n) [b, a] = bilinear(Num, Den, 1/T); w = [0: pi/1024: pi]; tranfun=tf(b,a) I was able to deterimine T and obtain H(z) below; Transfer function: 0.1125 s^4 + 0.4498 s^3 + 0.6747 s^2 + 0.4498 s + 0.1125 -------------------------------------------------------- s^4 + 0.4715 s^3 + 0.3318 s^2 + 0.0453 s + 0.006805 PROBLEM 1d Using the following code in MatLab: %%problem 1d Htran = freqz(b, a, w); magH = abs(Htran); magHdB = 20*log10(magH); figure(3) plot(w, magHdB); % Plots the Magnitude response in dB xlabel(Frequency in rad/samp); ylabel(Amplitude in db); title(Final - Part 1d: Magnitude Response of the Digital Filter in dB) figure(4) plot(w,angle(Htran)); xlabel(F