1.总结学过的滤波器设计方法,用matlab**例子分析不同设计方法的滤波器的性能及适应场合。
答:1.1模拟低通滤波器的设计方法。
1.1.1 butterworth滤波器设计步骤:
.确定阶次n
1 已知ωc、ωs和as求butterworth df阶数n
由: 求出n:
2 已知ωc、ωs和ω=ωp()的衰减ap 求butterworth df阶数n
3 已知ωp、ωs和ω=ωp的衰减ap 和as 求butterworth df阶数n
由 .用阶次n确定
根据公式:在左半平面的极点即为的极点,因而。
1.1.2 切比雪夫低通滤波器设计步骤:
.确定技术指标
归一化。.根据技术指标求出滤波器阶数n及:
其中 .求出归一化系统函数。
其中极点由下式求出:
或者由和s直接查表得。
.去归一化:
2.数字低通滤波器的设计步骤:
1) 确定数字低通滤波器的技术指标:通带截止频率、通带最大衰减系数、阻带截止频率、阻带最小衰减系数。
2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。
巴特沃斯:切比雪夫。
其中。3)把模拟滤波器变换成数字滤波器,即把模拟滤波器的系数映射成数字滤波器的系统函数。
实现系统传递函数s域至z域映射有脉冲响应不变法和双线性映射两种方法。
3.1)脉冲响应不变法。
按照技术要求设计一个模拟低通滤波器,得到模拟低通滤波器的传输函数转换成数字低通滤波器的系统函数h(z)。
设模拟滤波器的传输函数为,相应的单位冲激响应是, =lt,lt[.]代表拉氏变换,对进行等间隔采样,采样间隔为t,得到,将h(n)= 作为数字滤波器的单位取样响应,那么数字滤波器的系统函数h(z)便是h(n)的z变换。因此脉冲响应不变法是一种时域上的转换方法,它是h(n)在采样点上等于。
设模拟滤波器只有单阶极点,且分母多项式的阶次高于分子多项式的阶次,将用部分分式表示:
式中为的单阶极点。
将逆拉氏变换得到:
式中u(t)是单位阶跃函数。
对进行等间隔采样,采样间隔为t,得到。
对上式进行z变换,得到数字滤波器的系统函数h(z):
经过一系列变换得到:
3.2)双线性变换法。
这种变换方法,采用非线性频率压缩方法,将整个频率轴上的频率范围压缩到±π/t之间,再用转换到z平面上。设ha(s),s=jω,经过非线性频率压缩后用ha(s1),=jω1表示,这里用正切变换实现频率压缩:
式中t仍是采样间隔,当ω1从-π/t经过0变化到π/t时,ω则由-∞经过0变化到+∞,实现了s平面上整个虚轴完全压缩到平面上虚轴的±π/t之间的转换。这样便有。
再通过转换到z平面上,得到:
令,有。两种方法比较:
脉冲响应不变法的优点:
1)模拟频率到数字频率的转换时线性的;
2)数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应,因此时域特性逼近好。
缺点:会产生频谱混叠现象,只适合带限滤波器。
双线性变换法优点:
克服多值映射得关系,可以消除频率的混叠。
缺点:时域到频域的变换是非线性的,在高频处有较大的失真。
3.数字高通滤波器的设计步骤:
数字高通滤波器的技术指标为:通带截止频率ωp 阻带截止频率ωs
通带衰减频率阻带衰减频率
预畸变处理,将数字高通指标转换为模拟低通指标。
4 确定阶数n由 (可由模拟低通滤波器设计方法可得)
归一化及去归一化查表令s=s/ω 归一化模拟低通圆型系统函数。
低通向高通转化令s1=1/s 由频率变换公式。
即可得。 滤波器数字化令利用双线性变换化。
带入数据可得数字高通h(z)
数字高通不能采用脉冲响应不变法原因是:脉冲响应不变法有频谱周期延拓效应,因此只能用于带限的频响特性,如衰减特性很好的低通或带通。而高频衰减越大,频响的混淆效应越小,至于高通滤波器,由于它们在高频部分不衰减,因此将完全混淆在低频响应中。
4.数字带通滤波器的设计:
步骤:1) 确定性能指标: 在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:通带截止频率wc1,wc2、阻带截止频率wr1,wr2、阻带最小衰减αs通带最大衰减αp
2) 对带通数字滤波器h(z)的数字边界频率预畸变。
得到带通模拟滤波器h(s)的边界频率。
主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1,ωs2的转换。
对双线性变换法一般t=2s
通带截止频率wc1=(2/t)*tan(wp1/2) 、wc2=(2/t)*tan(wp2/2)
阻带截止频率wr1=(2/t)*tan(ws1/2)、wr2=(2/t)*tan(ws2/2)
阻带最小衰减αs
通带最大衰减αp
3)低通到带通频率变换。
将模拟带通滤波器指标转换为模拟低通滤波器指标。
b=wc2-wc1
normwr1=((wr1^2)-(w0^2))/b*wr1))
normwr2=((wr2^2)-(w0^2))/b*wr2))
normwc1=((wc1^2)-(w0^2))/b*wc1))
normwc2=((wc2^2)-(w0^2))/b*wc2))
模拟低通滤波器指标:
normwc,normwr,αp,αs
4)设计模拟低通原型滤波器。查表得到归一化低通传输函数g(p):
用模拟低通滤波器设计方法(由巴特沃斯设计步骤或切比雪夫设计步骤)得到模拟低通滤波器的传输函数ha(s)
5)模拟低通滤波器转化为模拟带通滤波器。
6)利用双线性变换法将模拟带通滤波器ha(s)转换成数字带通滤波器h(z)
由得到 也可以用脉冲响应不变法:
两种方法比较比较:
脉冲响应不变法数字滤波器单位脉冲响应的数字表示近似原型的模拟滤波器单位脉冲响应,因此时域特性逼近好。但会产生频谱混叠现象,只适合带限滤波器。
双线性变换法可以克服多值映射得关系,可以消除频率的混叠但时域到频域的变换是非线性的,在高频处有较大的失真。
5.数字带阻滤波器的设计:
步骤:1) 确定性能指标:
通带截止频率wc1,wc2、阻带截止频率wr1,wr2、阻带最小衰减αs通带最大衰减αp
2) 对带通数字滤波器h(z)的数字边界频率预畸变。
主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1,ωs2的转换。
对双线性变换法一般t=2s
通带截止频率wc1=(2/t)*tan(wp1/2) 、wc2=(2/t)*tan(wp2/2)
阻带截止频率wr1=(2/t)*tan(ws1/2)、wr2=(2/t)*tan(ws2/2)
阻带最小衰减αs 通带最大衰减αp
由模拟低通滤波器设计方法可得)
模拟低通滤波器确定模拟带阻滤波器。
由模拟低通到模拟带阻的变换这一模拟低通到带阻的变换关系为
式中s为模拟低通原型拉普拉斯变量(s=σ+jω),为模拟带阻的拉普拉斯变量(=σjω )是模拟带阻滤波器的几何中心频率。令=jω可得
故平面的虚轴与s平面的虚轴相对应,代入s=jω,消去j,可得。
低通的阻带映像到带阻的阻带。
化简得到:在h()中的变换关系,可得到带阻滤波器系统函数。
由模拟带阻到数字带阻的变换仍利用双线性变换。
模拟低通原型滤波器的平面变换成数字带阻滤波器的平面的表达式
可得。从模拟低通系统函数h(),转换数字带阻系统函数
数字带阻滤波器不能用脉冲响应不变法:原因是脉冲响应不变法有频谱周期延拓效应,因此只能用于带限的频响特性,如衰减特性很好的低通或带通。而高频衰减越大,频响的混淆效应越小,至于带阻滤波器,由于它们在高频部分不衰减,因此将完全混淆在低频响应中。
matlab中程序运行:
> fp=2100; fs=8000;
>fs=20000;
>rp=0.5; rs=30;
>t=1/fs; %设计指标
>w1p=fp/fs*2;w1s=fs/fs*2;
> [n,wn] =buttord (w1p,w1s,rp,rs,'s');
> [z,p,k]=buttap(n);
> [bp,ap]=zp2tf(z,p,k);
> [bs,as]=lp2lp(bp,ap,wn*pi*fs);
> [bz,az]=impinvar(bs,as,fs); 用脉冲响应不变法进行模数变换
>sys=tf(bz,az,t
> [h,w]=freqz(bz,az,512,fs);
> plot(w,abs(h));grid on; xlabel('频率/hz');ylabel('振幅/h');title('巴特沃斯滤波器')
现代信号处理作业
中国矿业大学。20 级硕士研究生课程考试试卷。考试科目现代信号处理基础 考试时间。学生姓名。学生学号。所在院系信电学院 任课教师。中国矿业大学研究生院培养管理处印制。一 自适应滤波。1 自适应df的工作原理 自适应滤波器是以最小均方误差为准则的最佳滤波器,它能自动调节其本身的单位脉冲响应h n 特性...
现代信号处理作业
中北大学。研究生 现代信号处理 作业。所在院 系 信息与通信工程学院 专业 测试计量技术及仪器 姓名王栋。学号 s20100329 班级 y100503 156 小波变换与matlab相关程序的学习 小波变换是在傅里叶变换的基础上发展起来的,它优于傅里叶分析的地方是它在时域与空域都是局部化的,其局部...
现代信号处理作业
2014年6月。cramer rao不等式的证明。定义 设为参数的估计子,参数的估计子的偏差定义为该估计子误差的期望值,即。则估计子称为无偏估计子。若变差等于零或者,即估计子的期望值等于真实参数,则成估计子的渐进无偏估计子。无偏性反映了参数估计量的取值在真值周围摆动程度。一个参数往往具有不止一个无偏...