实验2数字信号处理Word文档格式.docx
- 文档编号:4400639
- 上传时间:2023-05-03
- 格式:DOCX
- 页数:15
- 大小:394.88KB
实验2数字信号处理Word文档格式.docx
《实验2数字信号处理Word文档格式.docx》由会员分享,可在线阅读,更多相关《实验2数字信号处理Word文档格式.docx(15页珍藏版)》请在冰点文库上搜索。
ylabel('
幅度响应'
wp=2*1000*tan(2*pi*200/(2*1000));
ws=2*1000*tan(2*pi*300/(2*1000));
Ap=1;
As=25;
[n,Wn]=buttord(Wp,Ws,Ap,As,'
[b,a]=butter(n,Wn,'
[bn1,an1]=bilinear(b,a,1000);
[H1,W]=freqz(bn1,an1);
[bn2,an2]=impinvar(b,a,1000);
[H2,W]=freqz(bn2,an2);
f=W/(2*pi)*1000;
plot(f,abs(H1),'
-.'
f,abs(H2),'
-'
gridon
legend('
双线性变法'
'
脉冲响应不变法'
wp=2*8000*tan(2*pi*1200/(2*8000));
ws=2*8000*tan(2*pi*2000/(2*8000));
[N1,wn1]=buttord(wp,ws,0.5,40,'
[B1,A1]=butter(N1,wn1,'
[bn1,an1]=bilinear(B1,A1,8000);
[h1,w1]=freqz(bn1,an1);
f1=w1/(2*pi)*8000;
subplot(311);
plot(f1,20*log10(abs(h1)));
axis([0,4000,-80,10]);
频率1'
幅度响应1'
title('
butterworth滤波器'
[N2,wn2]=cheb1ord(wp,ws,0.5,40,'
[B2,A2]=cheby1(N2,0.5,wn2,'
[bn2,an2]=bilinear(B2,A2,8000);
[h2,w2]=freqz(bn2,an2);
f2=w2/(2*pi)*8000;
subplot(312);
plot(f2,20*log10(abs(h2)));
频率2'
幅度响应2'
chebyshev滤波器'
[N3,wn3]=ellipord(wp,ws,0.5,40,'
[B3,A3]=ellip(N3,0.5,40,wn3,'
[bn3,an3]=bilinear(B3,A3,8000);
[h3,w3]=freqz(bn3,an3);
f3=w3/(2*pi)*8000;
subplot(313);
plot(f3,20*log10(abs(h3)));
频率3'
幅度响应3'
椭圆滤波器'
w1=2*30000*tan(2*pi*2000/(2*30000));
w2=2*30000*tan(2*pi*3000/(2*30000));
wr1=2*30000*tan(2*pi*1500/(2*30000));
wr2=2*30000*tan(2*pi*6000/(2*30000));
[N,wn]=buttord([w1w2],[wr1wr2],3,20,'
[B,A]=butter(N,wn,'
[bn1,an1]=bilinear(B,A,30000);
f=w/(2*pi)*30000;
axis([0,15000,-80,10]);
N=45,计算并画出矩形窗、汉明窗、布莱克曼窗的归一化的幅度谱,并比较各自的主要
特点。
N=45;
wn1=kaiser(N,0);
wn2=hamming(N);
wn3=blackman(N);
[h1,w1]=freqz(wn1,N);
[h2,w2]=freqz(wn2,N);
[h3,w3]=freqz(wn3,N);
plot(w1/pi,20*log10(abs(h1)),'
w2/pi,20*log10(abs(h2)),'
--'
w3/pi,20*log10(abs(h3)),'
:
'
axis([0,1,-120,10]);
grid;
归一化频率/\pi'
幅度/dB'
三种窗口函数'
矩形窗'
汉明窗'
布莱克曼窗'
3);
N=15;
h=fir1(N-1,[0.30.5],'
bandpass'
hanning(N));
[h1,w1]=freqz(h,1);
subplot(2,1,1);
plot(w1/pi,20*log10(abs(h1)));
axis([0,1,-80,10]);
频率/\pi'
N=15,矩形窗'
subplot(2,1,2);
N=45,矩形'
figure
(1)
freqz(h,1)
N=15,汉宁窗'
figure
(2)
N=45,汉宁窗'
将上式的hanning改为下面的形式
%矩形窗
boxcar(N));
blackman(N));
N=15,布莱克曼窗'
N=45,布莱克曼窗'
如:
N=40;
f=[00.20.20.40.40.60.60.80.81];
a=[0011001100];
beta=4;
h=fir2(N-1,f,a,kaiser(N,beta));
subplot(3,1,1);
beta=4时凯泽窗专用线性相位滤波器'
beta=6;
subplot(3,1,2);
beta=6时凯泽窗专用线性相位滤波器'
beta=10;
subplot(3,1,3);
beta=10时凯泽窗专用线性相位滤波器'
第四题
f=[00.50.51];
m=[1100];
b=fir2(63,f,m);
[h,w]=freqz(b,1,128);
plot(f,m,w/pi,abs(h))
Ideal'
fir2Designed'
)
频率采样的低通滤波器'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 数字信号 处理