Курсовая работа: Проектування блоку обробки сигналів
www=250/f_discr;
//zz_t=iir(4,'lp','butt',[www www],[0 0]);
zz_t=iir(4,'lp','cheb1',[www www],[0.1 0]);
//zz_t=iir(4,'lp','cheb2',[www www],[0 0.1]);
//zz_t=iir(4,'lp','ellip',[www www],[0.1 0.1]);
sz_t=syslin('d',zz_t);
Hz_t=freq(sz_t.num,sz_t.den,exp(%i*2*%pi*ft));
ampz_t=abs(Hz_t);
phiz_t=atan(imag(Hz_t),real(Hz_t));
scf(10)
subplot(211)
plot2d(frq_s,ampz_t,logflag="ln")
xgrid()
subplot(212)
plot2d(frq_s,phiz_t,logflag="ln")
xgrid()
yf_out=flts(s_in,sz_t);
scf(11)
subplot(211), plot(yf_out), xgrid()
subplot(212), plot(s_in), xgrid()
scf(12)
ach=abs(fft(yf_out));
plot2d(frq_s,ach(1:641)), xgrid()
err1=0;
for n=3:(N-2),
err1=err1+(ys(n-2)-yf_out(n+2))^2;
end
err1=sqrt(err1/n);
//Режекторный фильтр;