50t.m

包含50个数字信号处理方面的matlab程序

N=100;n=[0:1:N-1];Wc1=0.2*pi;Wc2=0.4*pi;
f1=50;f2=150;fs=1000;
hd1=ideal_lp(Wc1,N);
hd2=ideal_lp(Wc2,N);
hd0=ideal_lp(pi,N);
w_han=(blackman(N))';
h1=hd1.*w_han;
h2=hd2.*w_han;
h0=hd0.*w_han;
h=h0-h2+h1;
[db,mag,pha,grd,w]=freqz_m(h,1);
subplot(3,1,1);plot((w*1000)/(pi*2),db);axis([0,400,-100,5]);title('frequency domain');
n=0:240;t=0.001*n;
x=sin(2*f1*pi*t)+sin(2*f2*pi*t);
subplot(3,1,2);plot(n,x);axis([0,250,-2,2]);
y=filter2(h,x);
subplot(3,1,3);plot(n,y);axis([0,250,-1.2,1.2]);