资源描述:
《matlab实验窗函数》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验程序及运行结果:1.矩形窗N=32;w=boxcar(N);nn=0:(N-1);stem(nn,w)Hanning窗N=32;w=hanning(N);nn=0:(N-1);stem(nn,w)Hamming窗N=32;w=hamming(N);nn=0:(N-1);stem(nn,w)14Bartelett窗N=32;w=bartelett(N);nn=0:(N-1);stem(nn,w)Blackman窗N=32;w=blackman(N);nn=0:(N-1);stem(nn,w)Triang窗N=32;w=blackman(N);nn=0:(N-1);s
2、tem(nn,w)14Kaiser窗N=32;b=8;w=kaiser(N,b);nn=0:(N-1);stem(nn,w)Chebwin窗N=32;r1=30;r2=50;w1=chebwin(N,r1);w2=chebwin(N,r2);nn=0:(N-1);subplot(211);stem(nn,w1),title('r=30');subplot(212);stem(nn,w2),title('r=50');2.矩形窗N=16;w=boxcar(N);[H,W]=dtft(w,2048);plot(W/pi,abs(H))14Hanning窗N=16;w=ha
3、nning(N);[H,W]=dtft(w,2048);plot(W/pi,abs(H))Hamming窗N=16;w=hamming(N);[H,W]=dtft(w,2048);plot(W/pi,abs(H))Bartlett窗N=16;w=bartlett(N);[H,W]=dtft(w,2048);plot(W/pi,abs(H))14Blackman窗N=16;w=blackman(N);[H,W]=dtft(w,2048);plot(W/pi,abs(H))Triang窗N=16;w=triang(N);[H,W]=dtft(w,2048);plot(W/
4、pi,abs(H))Kaiser窗N=16;b=8;w=kaiser(N,b);[H,W]=dtft(w,2048);plot(W/pi,abs(H))14Chebwin窗N=16;b=50;w=kaiser(N,b);[H,W]=dtft(w,2048);plot(W/pi,abs(H))3.N=10;w=boxcar(N);[H,W]=dtft(w,1024);subplot(211);plot(W/pi,abs(H)),grid;xlabel('frequency'),ylabel('magnitude');subplot(212);plot(W/pi,angl
5、e(H)),grid;xlabel('frequency'),ylabel('phase');14N=20;w=boxcar(N);[H,W]=dtft(w,1024);subplot(211);plot(W/pi,abs(H)),grid;xlabel('frequency'),ylabel('magnitude');subplot(212);plot(W/pi,angle(H)),grid;xlabel('frequency'),ylabel('phase');N=50;w=boxcar(N);[H,W]=dtft(w,1024);subplot(211);plo
6、t(W/pi,abs(H)),grid;xlabel('frequency'),ylabel('magnitude');subplot(212);plot(W/pi,angle(H)),grid;xlabel('frequency'),ylabel('phase');14N=100;w=boxcar(N);[H,W]=dtft(w,1024);subplot(211);plot(W/pi,abs(H)),grid;xlabel('frequency'),ylabel('magnitude');subplot(212);plot(W/pi,angle(H)),grid;
7、xlabel('frequency'),ylabel('phase');144.矩形窗fs=10;Tp=4;f=25/16;N=0.9*Tp*fs;n=[0:N-1];w=boxcar(N);x=0.75+3.4*cos(2*pi*f*n/fs)+2.7*cos(4*pi*f*n/fs)+1.5*sin(3.5*pi*f*n/fs)+2.5*sin(7*pi*f*n/fs);y=w.*x';[H,W]=dtft(y,1024);subplot(211);plot(W/2/pi,abs(H)),grid;xlabel('frequency'),ylab