第三章习题答案1要点Word格式.docx
《第三章习题答案1要点Word格式.docx》由会员分享,可在线阅读,更多相关《第三章习题答案1要点Word格式.docx(29页珍藏版)》请在冰豆网上搜索。
figure
(2);
OMEGA=-pi:
pi/100:
pi;
H=freqz(b,a,OMEGA);
subplot(2,1,1),plot(OMEGA,abs(H));
subplot(2,1,2),plot(OMEGA,180/pi*unwrap(angle(H)));
题3-5(b)图
3-6给定,确定模拟滤波器的系统函数。
根据给定的平方幅度响应,得
与比较,得到。
取左半平面的三个极点,得
,极点;
因此由,得
对共轭极点,有
代入上式,得
3-7模拟低通滤波器的参数如下:
,,,,用巴特沃斯近似求。
已知,,,,确定巴特沃斯滤波器的阶数如下:
取。
本题由于正好是,故低通滤波器的截止频率为:
或者,由下式来求取。
将代入五阶巴特沃斯模拟低通传递函数
3-8已知,使用脉冲响应不变法和双线性方法分别设计数字低通滤波
器,使得截止频率为=。
(1)双线性变换法:
3截止频率为=,
于是
=0.2920
参数不参与设计
(2)脉冲响应不变法:
因为脉冲响应不变法是由下面的映射完成的:
所以
3-9用脉冲响应不变法将转换为,采样周期为
,其中为任意整数
上式递推可得:
3-10要求设计一个数字低通滤波器,在频率低于的范围内,低通幅度特性为常数,并且不低于0.75,在频率和之间,阻带衰减至少为20。
试求出满足这些指标的最低阶巴特沃斯滤波器的传递函数,采用双线性变换。
令为模拟滤波器的平方幅度函数,且由于采用双线性变换,
若,故我们要求
因此巴特沃斯滤波器的形式为:
因此:
指标放松一点,可以取,代入上式得
对于这个值,通带技术指标基本达到,阻带技术指标刚好满足,在s平面左半部由三个极点对,其坐标为。
极点对1:
;
极点对2:
极点对3:
以代入上式,最后可得
3-11试设计一巴特沃斯数字低通滤波器,设计指标为:
在0.3通带频率范围内,通带幅度波动小于1,在0.5阻带频率范围内,阻带衰减大于12。
由题意可以得出:
=0.3,=1
=0.5,=12
(1)频率预畸变
===1.019/
===2/
(2)确定滤波器阶数:
===0.1321
===1.9627
=-=-=3.002,取=3
(3)查表求归一化低通滤波器函数
=
(4)求模拟滤波器系统函数
()
==
(5)求系统函数
将代入得:
=
3-12用双线性变换法设计数字低通滤波器,等效模拟滤波器指标参数如下:
输入模拟信号的最高频率;
选用巴特沃斯滤波器,3截止频率,阻带截止频率,阻带最小衰减20。
(1)因为采用双线性变换法设计,数字频率与相应的模拟频率之间为非线性关系
但采样数字滤波器要求其中的数字滤波器与总等效模拟滤波器之间的频率映射关系为线性关系=。
所以,不能直接按等效模拟滤波器技术指标设计相应模拟滤波器,再将其双线性变换法映射成数字滤波器。
因此,我们必须按
将等效模拟滤波器指标参数转换成采样数字滤波器系统中数字滤波器指标参数,再用双线性变换法的一般步骤设计该数字滤波器。
通带边界频率:
===;
通带最大衰减:
=3
阻带截至频率:
阻带最小衰减:
=20
以下为双线性变换法设计数字滤波器的一般过程。
(2)预畸变校正,确定相应的模拟滤波器指标参数:
==800=800,=3
==800=1931.37,=20
(3)确定滤波器阶数:
===0.1003
===2.414
===2.609
取=3。
(4)查表求归一化低通滤波器函数
(5)求模拟滤波器系统函数
==
(6)求系统函数=
3-13试设计一个数字高通滤波器,要求通带下限频率rad。
阻带上限频率为,通带衰减不大于3,阻带衰减不小于20。
(1)已知数字高通滤波器指标:
=0.8rad,=3
=0.44rad,=20
(2)由于设计的是高通数字滤波器,所以采用双线性变换法,所以要进行预畸变校正,确定相应的模拟高通滤波器指标(为了计算方便,取T=2s):
===3.0777rad/s
===0.8273rad/s
(3)将高通滤波器指标转换成模拟低通指标。
高通归一边界频率为(本题=):
==1
低通指标为:
==3.7203
(4)设计归一化低通:
===0.1003
==3.7203
=-=1.75,取=2
查表可得
(5)频率变换,求模拟高通:
(6)用双线性变换法将转换成:
=|s==
3-14一个数字系统的抽样频率,试设计一个为此系统使用的带通数字滤波器
,希望采用巴特沃斯滤波器,通带范围为到,在带边频率处的衰减不大于3;
在以下和以上衰减不小18。
(1)确定带通滤波器技术指标:
因为,所以得出:
,
,
通带内最大衰减,阻带内最小衰减。
(2)确定相应模拟滤波器技术指标。
为计算简单,设。
通带中心频率
带宽
=-
将以上边界频率对归一化,得到相应归一化带通边界频率:
=,=
(3)由归一化带通指标确定相应模拟归一化低通技术指标。
归一化阻带截止频率为:
归一化通带截止频率为
(4)设计模拟归一化低通
==
--
取,因为1.9394很接近2,所以取基本满足要求,且系统简单。
查表可得归一化低通系统函数:
(5)频率变换,将转换成模拟带通
(6)用双线性变换公式将转换成
3-15一个数字系统的抽样频率为1000,已知该系统受到频率为100的噪声干扰,现设计一带阻滤波器去掉该噪声。
要求3的带边频率为95和105,阻带衰减不小于14,阻带的下边和上边频率分别为99和101。
(1)确定带阻滤波器技术指标:
阻带中心频率
带阻带宽
将以上边界频率对归一化,得到相应归一化带阻边界频率:
(3)由归一化带阻指标确定相应模拟归一化低通技术指标。
取,。
取,因为1.055很接近1,所以取,基本满足要求,且系统简单。
(5)频率变换,将转换成模拟带阻
=
3-16试用矩形窗口设计法设计一个FIR线形相位低通数字滤波器,已知,
画出和曲线,再计算正、负肩峰值的位置和过渡带宽度。
写出理想的频响为
求得理想冲激响应为
计算得
加矩形窗:
即
正、负肩峰值的位置如题3-16表所示。
题3-16表
轴上的位置
值
正肩峰(A点)
临界频率(B点)
负肩峰(C点)
过渡带宽度为
利用MATLAB演示其结果如题3-16图所示:
N=21;
a=(N-1)/2;
Wc=0.5*pi;
n=[0:
1:
(N-1)];
m=n-a+eps;
避免被零除
hd=sin(Wc*m)./(pi*m);
[H1,W]=freqz(hd,1);
figure
(1);
subplot(211);
stem(n,hd);
xlabel('
n'
);
ylabel('
h(n)'
subplot(212);
plot(W/pi,20*log10(abs(H1)/max(H1)));
频率'
幅频响应'
题3-16图加矩形窗时的脉冲响应及其频谱图
3-17试用窗函数法设计一个第一类线性相位FIR数字高通滤波器,已知
,,,对于矩形窗,过渡带宽度为。
(1)确定的长度;
(2)求的表达式;
(3)?
(1)偶数,取=65
(2)
(3)
3-18用矩形窗设计线性相位数字低通滤波器,理想滤波器传输函数为:
(1)求出相应的理想低通滤波器的单位脉冲响应。
(2)求出用矩形窗函数法设计的FIR滤波器的表达式。
(1)
(2)为满足线性相位条件,要求,为矩形窗函数的长度。
加矩形窗函数得。
3-19用矩形窗设计线性相位高通滤波器,逼近滤波器传输函数为
=
(1)求出相应于理想低通的单位脉冲响应;
(2)求出矩形窗设计法的表达式,确定与之间的关系;
(3)的取值有什么限制?
为什么?
(1)直接IDTFT[]计算:
表达式中第2项正好是截止频率为的理想低通滤波器的单位脉冲相应。
而对应于一个线性相位全通滤波器:
即高通滤波器可由全通滤波器减去低通滤波器实现。
(2)用N表示
为了满足线性相位条件:
要求满足=
(3)必须取奇数。
因为N为偶数时(情况2),=0,不能实现高通。
3-20使用频率取样设计法(第一种形式取样)设计一个FIR线性相位低通数字滤波器。
已知。
理想低通0到和到处幅度函数为1,其余为0。
采样频率间隔为
,的位置在,即和之间,其对称点位置是,即和之间。
对理想低通采样,可得
第一类FIRDF的相位特性为
综合幅度和相位,FIRDF的离散频域抽样值为
3-21用频率采样法设计第一类线性相位FIR低通滤波器,要求通带截止频率=,阻
带最大衰减25,过渡带宽度=,问滤波器长度至少为多少才可能满足要求.
因为阻带最大衰减为25所以需要一个过渡点即,所以得出滤波器长度为
3-22利用频率采样法设计线性相位FIR低通滤波器,设N=16,给定希望逼近滤波器的幅
度采样为
有希望逼近的滤波器幅度采样可构造出的:
=IDFT[]==
阻带最小衰减接近-40dB。
3-23一个IIR网络的差分方程为=,当输入序列