Курсовая работа: Проектування блоку обробки сигналів

www2sr=48/f_discr;

//zz_tsr=iir(4,'sb','butt',[www1sr www2sr],[0.5 0]);

//zz_tsr=iir(4,'sb','cheb1',[www1sr www2sr],[0.5 0]);

zz_tsr=iir(4,'sb','cheb2',[www1sr www2sr],[0 0.5]);

//zz_tsr=iir(4,'sb','ellip',[www1sr www2sr],[0.5 0.3]);

sz_tsr=syslin('d',zz_tsr);

Hz_tsr=freq(sz_tsr.num,sz_tsr.den,exp(%i*2*%pi*ft));

ampz_tsr=abs(Hz_tsr);

phiz_tsr=atan(imag(Hz_tsr),real(Hz_tsr));

scf(13)

subplot(211)

plot2d(frq_s,ampz_tsr,logflag="ln")

xgrid()

subplot(212)

plot2d(frq_s,phiz_tsr,logflag="ln")

xgrid()

yf_out2=yf_out;

yf_out2=flts(yf_out2,sz_tsr);

scf(14)

subplot(211), plot(yf_out2), xgrid()

subplot(212), plot(s_in), xgrid()

err2=0;

for n=4:(N-3),

err2=err2+(ys(n-2)-yf_out2(n+2))^2;

end

err2=sqrt(err2/n);

scf(16)

ach2=abs(fft(yf_out2));

plot2d(frq_s,ach2(1:641)), xgrid()