50.m

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

f1=50;f2=150;
wc1=0.1*pi;wc2=0.2*pi;
N=100;n=[0:1:N-1];
hd=ideal_lp(wc1,N)+ideal_lp(pi-wc2,N).*exp(j*pi*n);
w_han=(hanning(N))';
h=hd.*w_han;
[db,mag,pha,grd,w]=freqz_m(h,1);
subplot(3,1,1);plot(w/pi,db);xlabel('w/pi');ylabel('pha');
axis([0,1,-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=filter(h,1,x);
subplot(3,1,3);plot(n,y);axis([0,250,-1.2,1.2]);