matlab实验窗函数
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验程序及运行结果: 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)
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); stem(nn,w)
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))
05101520253035
05101520253035
Hanning窗
N=16;
w=hanning(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))
-1-0.8-0.6-0.4-0.200.20.40.60.81
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/pi,abs(H))
Kaiser窗
N=16;
b=8;
w=kaiser(N,b); [H,W]=dtft(w,2048); plot(W/pi,abs(H))
-1
-0.8-0.6-0.4-0.200.20.40.60.81
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,angle(H)),grid;
xlabel('frequency'),ylabel('phase');
-1
-0.8-0.6-0.4-0.2
00.20.40.60.81
05
10
frequency
m a g n i t u d
e
-1
-0.8-0.6-0.4-0.2
00.20.40.60.81
-4-202
4frequency
p h a s e
N=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');
-1
-0.8-0.6-0.4-0.2
00.20.40.60.81
051015
20frequency
m a g n i t u d e
-1
-0.8-0.6-0.4-0.2
00.20.40.60.81
-4-202
4frequency
p h a s e
N=50;
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');