器的参数显示滤波器的频率响应选择信号等。5滤波器设计方法综述滤波器的设计步骤1）给出系统的性能指标2）用一个离散的时间系统逼近这些性能指标3）实现该系统一般我们利用数字计算的方法实现系统所以将该离散时间滤波器称为数字滤波器滤波器的指标往往是以频域的形式给出的尤其是低通带通高通和带阻这些选频滤波器如图所示的一个线性时不变离散系统如果输入是带限的且采样率满足奈奎斯特采样率这系统是一个线性时不变的连续系统
HejΩTΩπTωHeffjΩHejωHeffjωπT
6实验程序
fs22050
语音信号采样频率为22050
x1wavreadWi
dowsCriticalStopwav读取语音信号的数据，赋给变量x1
sou
dx122050
播放语音信号
y1fftx11024
对信号做1024点FFT变换
ffs05111024
figure1plotx1title原始语音信号
做原始语音信号的时域图形
xlabeltime

ylabelfuzhi

figure2freqzx1title频率响应图
绘制原始语音信号的频率响应图
figure3
fsubplot211
plotabsy11512
做原始语音信号的FFT频谱图
title原始语音信号FFT频谱
subplot212
plotfabsy11512title原始语音信号频谱
xlabelHz
ylabelfuzhi
程序2：
fs22050
语音信号采样频率为22050
x1wavreadWi
dowsCriticalStopwav读取语音信号的数据，赋给变量x1
t0122050sizex1122050
y1fftx11024
对信号做1024点FFT变换
ffs05111024x2ra
d
1le
gthx1产生一与x长度一致的随机信号
sou
dx222050
figure1plotx2title高斯随机噪声
做原始语音信号的时域图形
xlabeltime

ylabelfuzhi

ra
d
state0mra
d
sizex1x201mx1
sou
dx222050播放加噪声后的语音信号y2fftx21024figure2plottx2title加噪后的语音信号xlabeltime
ylabelfuzhi
figure3subplot211plotfabsy21512title原始语音信号频谱xlabelHzylabelfuzhisubplot212
fplotfabsy21512title加噪后的语音信号频谱xlabelHzylabelfuzhi
根据以上代码，你可以修改下面有错误的代码程序3双线性变换法设计Butterworth滤波器
fs22050x1wavreadh课程设计2shuziwav
t0122050sizex1122050
Au003
dAucos2pi5000t
x2x1d
wp025pi
ws03pi
Rp1
Rs15
Fs22050
Ts1Fswp12Tsta
wp2
将模拟指标转换成数字指标
ws12Tsta
ws2
NW
buttordwp1ws1RpRss选择滤波器的最小阶数
ZPKbuttapN
创建butterworth模拟滤波器
BapAapzp2tfZPK
balp2lpBapAapW
bzazbili
earbaFs滤波器的转换
HWfreqzbzaz
用双线性变换法实现r