filter_parameter.m

这里面有很多dsp应用程序源码

clear all
%  clf
N=256;
fs=600;
dt=1/fs;
for k=1:N
    f1=100;f2=250;f3=270;
    pi_k_dt=2*pi*k/fs;
    y(k)=sin(pi_k_dt*f1)+sin(pi_k_dt*f2)+sin(pi_k_dt*f3);
end
y=round((y/3)*2^15);
%  inputdata=(buffer(int16(y),8))'

lp=150;     
wnl=2*lp/fs; 
[B,A]=CHEBY1(3,0.5,wnl)
[Z,P,K]=CHEBY1(3,0.5,wnl)
YY=filter(B,A,y);
B=[0.1589 0.4768 0.4768 0.1589];
A=[1 0.1268 -0.5239 0.1257];
YY2=filter(B,A,y);
subplot(3,1,1)
plot([1:100],y(1:100))
subplot(3,1,2)
plot([1:100],YY(1:100));
subplot(3,1,3)
plot([1:100],YY2(1:100));

% repmat('.WORD',32,1)we