定点位宽对滤波器处理过程影响.pdf

%finite-word length effect on filter processing %y[n]=x[n]+0.875y[n-1]-0.75y[n-2] %% %%%%%%%%%%%%%%%%%%%%%%%%(i) %H(z)=Y(z)/X(z)=1/(1-0.875z^(-1)+0.75z^(-2)) %The we get parameter of b and a b=[1,0,0]; a=[1,-0.875,0.75]; %%%%(i-a) %plot zero-pole pattern figure(1) zplane(b,a);grid title((i)zero-pole pattern); %%%%(i-b) %plot frequency response figure(2) [h,w]=freqz(b,a,2001); plot(w/pi,20*log10(abs(h))),grid on xlabel(Normalized Frequency (\times\pi rad/sample)); ylabel(Magnitude (dB)); title((i)Frequency Response); %%%%(i-c) %x[n]=0.375*delta[n] %X(z)=0.375 for all z %Y(z)=X(z)H(z)=0.375/(1-0.875z^(-1)+0.75z^(-2)) b_1=0.375*b; a_1=a; disp(r,p,k of question(i-c):); [r,p,k]=residuez(b_1,a_1) %from Parameter r,p,k we know %Y(z)=r(1)/(1-p(1)z^(-1))+r(2)/(1-p(2)z^(-1)) %Then %y[n]=(r(1)(p(1))^n+r(2)(p(2))^n)u[n] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% %% %%%%%%%%%%%%%%%%%%%%%%%%(ii) %decimal to 2s complement %1sign+4fraction b_2=sfi(b,5,4); a_2=sfi(a,5,4); %%%%(ii-a) %plot zero-pole pattern figure(3) zplane(b_2.data,a_2.data);grid title((ii)zero-pole pattern); %%%%(ii-b) %plot frequency response figure(4) [h_2,w_2]=freqz(b_2.data,a_2.data,2001); plot(w_2/pi,20*log10(abs(h_2))),grid on xlabel(Normalized Frequency (\times\pi rad/sample)); ylabel(Magnitude (dB)); title((ii)Frequency Response); %%%%(ii-c) %using rouding method %b_2c=sfi(sfi(0.375,5,4)*b_2,5,4); b_2c_t=sfi(0.375,5,4)*b_2; b_2c=roundn(b_2c_t.data,-4); a_2c=a_2; [r_2c,p_2c,k_2c]=residuez(b_2c,a_2c.data); disp(r(ii-c):); r_2c=roundn(r_2c,-4) disp(p(ii-c):); p_2c=roundn(p_2c,-4) disp(k(ii-c):); k_2c=roundn(k_2c,-4) %from Parameter r(ii-c),p(ii-c),k(ii-c) we know %Y(z)=r(1)/(1-p(1)z^(-1))+r(2)/(1-p(2)z^(-1)) %Then %y[n]=(r(1)(p(1))^n+r(2)(p(2))^n)u[n] %%%%(ii-d) b_2d=fix(b_2c_t.data*10000)/10000; a_2d=a_2; [r_2d,p_2d,k_2d]=residuez(b_2d,a_2d.data); disp(r(ii-d):); r_2d=fix(r_2d*10000)/10000 disp(p(ii-d):); p_2d=fix(p_2d*