程序如下：yFsbitswavreadd声音片段01wav读出信号，采样率和采样位数Sou
dyYffty4096figure1subplot221plotytitle原时域波形ylabelamplitudexlabel
subplot222plotabsY对频域取模axis0409603横纵坐标确定title原频谱特性标题ylabelamplitudeY轴显示xlabelfreque
cyhzX轴显示程序运行的结果如下图所示：
famplitudeamplitude
105
005
10
原时域波形
原频谱特性3
2
1
50001000015000

001000200030004000freque
cyhz
（3）设计数字滤波器和画出频率响应根据语音信号的特点给出有关滤波器的性能指标：
1低通滤波器性能指标，fp1000Hzfc2000HzAs100db
Ap1db
用巴特沃斯滤设计的低通滤波器的程序如下：
wp212008000
归一化通带数字频率wp
ws220008000
归一化阻带数字截止频率ws
Rp05
通带波纹系数Rp
Rs60
最小阻带衰减Rs
NW
buttordwpwsRpRs求阶数N和滤波器截止频率W

umde
butterNRp
传输分子和分母的系数
hwfreqz
umde

求频率响应
fsubplot211
窗口分成2行1列绘图区第
1个区为当前活动区
plotwpiabsh
二维连续图形
grid
添加网格
title巴特沃思型低通滤波器的幅频响应加图形说明
xlabelomegapi
加X轴说明（后表示转义字符）
ylabel振幅幅值
加Y轴说明
subplot212
窗口分成2行1列绘图区第2
个区为当前活动区
plotwpi20log10absh
二维连续图形
grid
添加网格
title巴特沃思型低通滤波器的幅频响应加图形说明
xlabelomegapi
加X轴说明（后表示转
义字符）
ylabel振幅分贝
加Y轴说明
程序运行
f巴特沃思型低通滤波器的幅频响应15
振幅幅值
1
05
000102030405060708091
巴特沃思型低通滤波器的幅频响应200
振幅分贝
0
200
40000102030405060708091
（4）用滤波器对信号进行滤波
比较两种滤波器的性能，然后用性能较好的滤波器对采集的
信号进行滤波。在MATLAB中，FIR滤波器利用函数fftfilt对信号进行滤波，IIR滤波器利用函数filter对信号进行滤波。应用以上滤波器对信号进行滤波：
fs8000ywavreaddji
gjia
gxueyua
wavffilter
umde
y
读取原始语音信号滤波
f1fftf7000
subplot211plotabsf1图
画出滤波后的时域
ftitle滤波后的幅度谱
grid
subplot212
plota
glef1
画出滤波后的时域
图
title滤波后的相位谱
grid
wavwriteffsdji
gjia
gxueyua
1wav
其图形如图：
250200150100
5000
42024
0
滤波后的幅度谱
1000
2000
3000
4000
滤波后的相位谱
5000
6000
7000
1000
2000
3000
r