数字信号处理实验报告电子版Word下载.docx
《数字信号处理实验报告电子版Word下载.docx》由会员分享，可在线阅读，更多相关《数字信号处理实验报告电子版Word下载.docx（16页珍藏版）》请在冰豆网上搜索。
2.MATLAB7.0软件
四.实验原理
1.学习估计巴特沃兹滤波器的阶数的MATLAB命令。
2.因果无限冲激响应和有限冲激响应传输函数的结构表示法，的基础上学会设计数字滤波器。
3.运用理论课知识求解问题。
五、实验步骤
1.根据实验题目要求进行分析；
2.运用所学知识用MATLAB编程实现题目要求；
3.对结果进行分析总结。
六、实验记录（数据、图表、波形、程序等）
7.1：
[N1,Wn1]=buttord（0.1,0.2,0.5,40）
[N2,Wn2]=cheb1ord（0.1,0.2,0.5,40）
[N3,Wn3]=cheb2ord（0.1,0.2,0.5,40）
[N4,Wn4]=ellipord（0.1,0.2,0.5,40）
N1=8Wn1=0.1151
N2=5Wn2=0.1000
N3=5Wn3=0.2000
N4=4Wn4=0.1000
（1）
使用这些值得到巴特沃斯低通滤波器最低阶数N=8，相应的
标准通带边缘频率Wn是0.1151.
（2）
使用这些值得到切比雪夫1型低通滤波器最低阶数N=5，相
应的标准通带边缘频率Wn是0.1000.
（3）
使用这些值得到切比雪夫2型低通滤波器最低阶数N=5，相
应的标准通带边缘频率Wn是0.2000.
（4）
使用这些值得到椭圆低通滤波器最低阶数N=4，相应的标准
通带边缘频率Wn是0.1000.
从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求
7.2：
[N1,Wn1]=buttord（0.3,0.17,1,50）
[N2,Wn2]=cheb1ord（0.3,0.17,1,50）
[N3,Wn3]=cheb2ord（0.3,0.17,1,50）
[N4,Wn4]=ellipord（0.3,0.17,1,50）
N1=11Wn1=0.2754
N2=6Wn2=0.3000
N3=6Wn3=0.1700
N4=5Wn4=0.3000
使用这些值得到巴特沃斯低通滤波器最低阶数N=11，相应的
标准通带边缘频率Wn是0.2754.。
使用这些值得到切比雪夫1型低通滤波器最低阶数N=6，相
应的标准通带边缘频率Wn是0.3000。
使用这些值得到切比雪夫2型低通滤波器最低阶数N=6，相
应的标准通带边缘频率Wn是0.1700.。
使用这些值得到椭圆低通滤波器最低阶数N=5，相应的标准
通带边缘频率Wn是0.3000.。
7.5：
Ws=[0.40.6];
Wp=[0.20.8];
Rp=0.4;
Rs=50;
[N1,Wn1]=buttord（Wp,Ws,Rp,Rs）;
[num,den]=butter（N1,Wn1,'
stop'
）;
disp（'
分子系数是'
disp（num）;
分母系数是'
disp（den）;
[g,w]=gain（num,den）;
h=freqz（num,den,w）;
[gd,w1]=grpdelay（num,den）;
plot（w/pi,g）
grid
axis（[01-605]）;
xlabel（'
\omega?
/\pi'
ylabel（'
增益dB'
title（'
Gain?
Responseofa?
ButterworthBandstopFilter'
pause
plot（w/pi,unwrap（angle（h）））;
Phase?
Spectrum?
arg[H（e^{j\omega}）]'
）
in?
radians'
plot（w1/pi,unwrap（gd））;
Delay?
Characteristic'
Group?
Delay'
分子系数是
Columns1through10
0.04930.00000.24650.00000.49300.00000.49300.00000.24650.0000
Column11
0.0493
分母系数是
1.00000.0000-0.08500.00000.63600.0000-0.02880.00000.05610.0000
-0.0008
7.6：
[num,den]=cheby1（N1,Rp,Wn1,'
dB'
Response?
of?
a?
Butterworth?
Bandstop?
Filter'
0.02760.00000.13790.00000.27570.00000.27570.00000.13790.0000
0.0276
1.00000.0000-0.9675-0.00001.44010.0000-0.95370.00000.5123-0.0000
7.20：
Ft=10000;
Fp=2000;
Fs=2500;
dev=[0.0050.005];
fedge=[FpFs];
aval=[10];
[N,Wn,beta,ftype]=kaiserord（fedge,aval,dev,Ft）;
b=fir1（N,Wn）;
[g,w]=gain（b,1）;
h=freqz（b,1,w）;
plot（w/pi,g）;
Gain'
移除跳变后的相位。
7.23：
w1=kaiser（184,4.6）;
b=fir1（183,0.31,w1）;
[g,w]=Gain（b,1）;
由图可知，该滤波器的设计满足指标。
7.24：
Wp=0.31;
Ws=0.41;
As=50;
Ds=10^（-As/20）;
F=[WpWs];
A=[10];
DEV=[DsDs];
[N,Wn,BTA,Ftype]=kaiserord（F,A,DEV）;
Win=kaiser（N+1,BTA）;
h=fir1（N,Wn,Ftype,Win）;
Numerator?
Coefficients?
are?
'
disp（h）;
[g,w]=gain（h,[1]）;
figure
（1）;
grid;
axis（[01-805]）;
Response'
w2=0:
pi/511:
pi;
Hz=freqz（h,[1],w2）;
figure
（2）;
Phase=angle（Hz）;
plot（w2/pi,Phase）;
（rad）'
figure（3）;
UPhase=unwrap（Phase）;
plot（w2/pi,UPhase）;
Unwrapped?
从图中可以看出设计的滤波器满足要求。
7.25
fpts=[00.250.30.450.51];
mval=[0.40.4110.80.8];
b=fir2（95,fpts,mval）;
plot（w/pi,abs（h））;
Amplitude?
response'
Amplitude'
从幅度响应中可以看出，此滤波器满足指标。
7.27
fedge=[1500180030004200];
aval=[010];
dev=[0.10.020.1];
FT=12000;
[NWnBetaftype]=kaiserord（fedge,aval,dev,FT）;
fpts=fedge/FT;
fpts=[0fpts1];
mval=[010010];
b=remez（N,fpts,mval）;
[g,w]=gain（b,1）;
观察上图可知，其阻带波纹过大，不符合题意。
七、实验思考题及解答
7.1
答：
从以上结果中观察到椭圆滤波器的阶数最低，并且符合要求。
7.2
7.20
是移除跳变后的相位。
7.23
滤波器的设计满足指标。
7.24
八、实验结果分析与总结
通过实验加深了对理论知识的理解。
对线性相位传输函数的判断方法有了进一步的掌握，对稳定传输函数的判断也有了更深的认识。
若传输函数的分子系数是对称的，则该传输函数是一个线性相位传输函数，若不对称，则不是线性的。
掌握了有限冲激响应传输函数以及无限冲激响应传输函数的实现。
同时对级联框图的画法更加熟悉了。
此外，对MATLAB更加熟悉了。
升级会员 