哈工大数字信号处理程序清单
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验程序及结果：
（1）矩形窗
clear all;
n=0:1:14;
wR=ones(1,15);
hd=sin(0.25*pi*(n-7+eps))./(pi*(n-7+eps));
h1=hd.*wR; N=64;
H1=fft(h1,N);
n=0:N-1;
w=2*pi/64*n;subplot(2,2,1);plot(w,fftshift(20*log10((abs(H1))))); grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB')
title('幅度曲线和相频曲线（n=15）');
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,3);plot(w,unwrap(phase(H1)));
xlabel('w/rad')
clear all;
n=0:1:32;
wR=ones(1,33);
hd=sin(0.25*pi*(n-16+eps))./(pi*(n-16+eps));
h1=hd.*wR;
N=64;
H1=fft(h1,N);
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,2);plot(w,fftshift(20*log10((abs(H1)))));
grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB')
title('幅度曲线和相频曲线（n=33）');
n=0:N-1;w=2*pi/64*n;subplot(2,2,4);plot(w,unwrap(phase(H1))); xlabel('w/rad')
（2）汉宁窗
clear all;
n=0:1:14;
wH=0.5*(1-cos(2*pi/14*n));
hd=sin(0.25*pi*(n-7+eps))./(pi*(n-7+eps)); h1=hd.*wH;
N=64;
H1=fft(h1,N);
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,1);
plot(w,fftshift(20*log10((abs(H1)))));
grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB');
title('幅度曲线和相频曲线（n=15）');
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,3);
plot(w,unwrap(phase(H1)));
grid
xlabel('w/rad')
n=0:1:32;
wH=0.5*(1-cos(2*pi/32*n));
hd=sin(0.25*pi*(n-16+eps))./(pi*(n-16+eps)); h1=hd.*wH;
N=64;
H1=fft(h1,N);
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,2);
plot(w,fftshift(20*log10((abs(H1)))));
grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB')
title('幅度曲线和相频曲线（n=33）');
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,4);plot(w,unwrap(phase(H1))); grid
xlabel('w/rad')
3.海明窗
n=0:1:14;
wH=0.54-0.46*cos(2*pi*n/(14+eps));
hd=sin(0.25*pi*(n-7+eps))./(pi*(n-7+eps));
h1=hd.*wH;
N=64;
H1=fft(h1,N);
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,1);
plot(w,fftshift(20*log10((abs(H1)))));
grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB')
title('幅度曲线和相频曲线（n=15）');
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,3);
plot(w,unwrap(phase(H1)));
grid
xlabel('w/rad')
n=0:1:32;
wH=0.54-0.46*cos(2*pi*n/(32+eps));
hd=sin(0.25*pi*(n-16+eps))./(pi*(n-16+eps)); h1=hd.*wH; N=64;
H1=fft(h1,N);
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,2);plot(w,fftshift(20*log10((abs(H1))))); grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB')
title('幅度曲线和相频曲线（n=33）');
n=0:N-1;
w=2*pi/64*n;
subplot(2,2,4);
plot(w,unwrap(phase(H1)));
grid
xlabel('w/rad')
n=0:1:14;
wH=0.54-0.46*cos(2*pi*n/(14+eps));
hd=sin(0.25*pi*(n-7+eps))./(pi*(n-7+eps));
h1=hd.*wH;
N=64;
H1=fft(h1,N);
n=0:N-1;w=2*pi/64*n;
subplot(2,2,1);
plot(w,fftshift(20*log10((abs(H1)))));
grid
xlabel('w/rad')
ylabel('20lg|H(jw)|/dB')
title('幅度曲线和相频曲线（n=15）');