语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

subplot(2,1,1) plot(yn1);

title('?¤?°ìa1-1??D??í?y') xn1=[2,0,0,8,8,4,2,5,0,1,1,0];

hn=[1, 2.43, 6.17,12.93,22.17,32.25,40.88, 45.87, 45.87, 40.88, 32.25, 22.17, 12.93, 6.17, 2.43,1.0000] yn2=fftfilt(xn1,hn); subplot(2,1,2); plot(yn2);

title(' ?¤?°ìa1-2?2?ü?í?y')

实现结果绘图如图7.1所示。

图7.1 预习题一结果图

(2)编写程序演示采样定理（时域采样、频谱周期延拓），同时演示采样频率小于2fc时，产生的混叠效应： ① 对下面连续信号进行采样：

a为衰xa(t)?Ae?atsin(?0t)u(t),A?n，a??0?n2?,n?学号， A为幅度因子，减因子，?0为模拟角频率，其中n为学号(例如，王墨同学n=23)

第 41 页 共 56 页

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

② 要求输入采样频率fs（根据程序处理需要指定范围）后，在时域演示信号波形、采样脉冲及采样后信号；在频域演示不同采样频率下对应信号的频谱。

MATLAB程序如下： clf;

A=27;a=27*sqrt(2)*pi;w0=27*sqrt(2)*pi; %xt=A*exp(-a*t)*sin(w0*t); n=0:1:100;T1=w0/100000;

xn=A*exp(-a*n*T1).*sin(w0*n*T1); subplot(2,1,1); stem(n,xn);

title('预习题2-1时域采样') T2=w0/100000;f=1/T2; t=0:T2:1;%采样周期为T xt=A*exp(-a*t).*sin(w0*t); y=fft(xt,512);

f1=f*(0:256)/512;%设置频率轴 subplot(2,1,2);

stem(f1,y(1:257));%画出频域内的信号 title('预习题2-2频域采样') 实现结果绘图如图7.2所示。

第 42 页 共 56 页

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

图7.2 预习题二结果图

参考文献

[1]《数字信号处理（第二版）》，丁玉美等，西安电子科技大学出版社

[2]《数字信号处理试验指导书》王创新、文卉编 长沙理工大学印刷（内部使用） [3]《数字信号处理及其MATLAB实现》，陈怀琛等译，电子工业出版社 [4]《MATLAB及在电子信息课程中的应用》，陈怀琛等，电子工业出版社 [5]《数字信号处理》A.V.奥本海姆，R.W.谢弗著，北京：科学出版社 [6]《数字信号处理——理论、算法与实现（第二版）》胡广书编著，北京：电子

第 43 页 共 56 页

语音信号的数字滤波处理——巴特沃思、bartlett窗滤波器设计

工业出版社

附录

附录A巴特沃思低通滤波器仿真程序

%语音信号 Fs=22050;

[x,FS,bits]=wavread('D:\\ding.wav'); x=x(:,1); figure(1);

第 44 页 共 56 页