dsp-1.rar_Sound quality

% Program P2_6 % Cascade Realization clf; n = 0:40; x = [1 zeros(1,40)];% Generate the input % Coefficients of 4th-order system den = [1 1.6 2.28 1.325 0.68]; num = [0.06 -0.19 0.27 -0.26 0.12]; % Compute the output of 4th-order system ic=[-1 1 0 1]; y = filter(num,den,x,ic); % Coefficients of the two 2nd-order systems num1 = [0.3 -0.2 0.4];den1 = [1 0.9 0.8]; num2 = [0.2 -0.5 0.3];den2 = [1 0.7 0.85]; % Output y1[n] of the first stage in the cascade ic1=[-1 1]; ic2=[0 1]; y1 = filter(num2,den2,x,ic1); % Output y2[n] of the second stage in the cascade y2 = filter(num1,den1,y1,ic2); % Difference between y[n] and y2[n] d = y - y2; % Plot output and difference signals subplot(3,1,1); stem(n,y); ylabel('Amplitude'); title('y[n]:Output of 4th-order Realization');grid; subplot(3,1,2); stem(n,y2) ylabel('Amplitude'); title('y2[n]:Output of Cascade Realization');grid; subplot(3,1,3); stem(n,d) xlabel('Time index n');ylabel('Amplitude'); title('d[n]:Difference Signal');grid;