VISI Lab 实验报告-----澳洲Swinbourne大学.docx

PAGE \* MERGEFORMAT29 Swinburne University of Technology Faculty of Engineering and Industrial Sciences EEE80003 Lab 1 \ Cascade Realization The factored form of a casual FIR transfer function H(z) of order M-1, as given below: (1) can be determined from its polynomial form representation given by: (2) This factored form given in (1) can then be utilized to realize H(z) in a cascade form. The Matlab function below achieves this conversion. For this first part, set den to 1 in accordance to equation (2). function [sos] = rat2sos(num,den) % [sos] = rat2sos(num,den) % Conversion of a rational transfer function % to its factored form % Adapted from Mitra 2000 % Output % ===== % sos : second order systems % Input %===== % num: numerator of the rational function % den: denominator of the rational function % Find the zeros and poles [z,p,k] = tf2zp(num,den); % then find the SOS matrix, you may also use the tf2sos function directly, for more recent versions of MATLAB. sos = zp2sos(z,p,k); The output is arranged in the follow fashion: Q1: Using the function rat2sos develop a cascade realization of the following FIR transfer function: Sketch the block diagram of the cascade realization. Is H1(z) a linear-phase transfer function? Firstly, using the function rat2sos find the SOS matrix of : Q1.m: num=[2 10 23 34 31 16 4]; % num : numerator of the rational function den=[1 0 0 0 0 0 0]; % den : denominator of the rational function sos=rat2sos(num,den) % Conversion of a rational transfer function % to its factored form % sos : second order systems h=impz(num,den); figure(1);stem(h) % make the Impulse Response plot of h(n) in figure (1) xlabel(n);ylabel(h(n)) % the xlabel of figure(1) is n, and ylable is h(n) title(Impulse Respone) % the title of figure(1) is Impulse Response. Run Q1.m, then we can get the SOS matrix and impulse response plot as following: Q1 sos = 2.0000 6.0000 4.0000 1.0000 0 0 1.0000 1.0000