实验七华北电力大学数字信号处理实验.docx
- 文档编号:1505609
- 上传时间:2022-10-22
- 格式:DOCX
- 页数:22
- 大小:538.91KB
实验七华北电力大学数字信号处理实验.docx
《实验七华北电力大学数字信号处理实验.docx》由会员分享,可在线阅读,更多相关《实验七华北电力大学数字信号处理实验.docx(22页珍藏版)》请在冰豆网上搜索。
实验七华北电力大学数字信号处理实验
实验报告
实验名称________________
课程名称________________
院系部:
专业班级:
学生姓名:
学号:
同组人:
实验台号:
指导教师:
成绩:
实验日期:
华北电力大学
1.实验目的
加深理解IIR数字滤波器的特性,掌握IIR数字滤波器的设计原理与设计方法,以及IIR数字滤波器的应用。
2.实验原理
3.实验内容及步骤
4.数据处理与总结
1.
(1)
clear;
fsam=100;
t0=1/fsam;t=6;k=0:
t0:
t;
N=512;
x=1+cos(pi/4.*k/t0)+cos(2*pi/3.*k/t0);
f=(-N/2:
(N/2-1))/N*2*pi;
[N1,Wc]=buttord(1/4,1/2,3,60);
[b,a]=butter(N1,Wc,'low');
freqz(b,a);
figure;
axis([0,1,-120,0]);
y=filter(b,a,x);
Y=fftshift(fft(y,N));
stem(f,abs(Y));
右图一为数字低通滤波器的幅频响应
和相频响应。
右图二为信号经过滤波器后输出信号的频谱。
(2)
clear;
fsam=100;
t0=1/fsam;t=6;k=0:
t0:
t;
N=512;
x=1+cos(pi/4.*k/t0)+cos(2*pi/3.*k/t0);
f=(-N/2:
(N/2-1))/N*2*pi;
[N1,Wc]=buttord(2/3,10/24,3,100);
[b,a]=butter(N1,Wc,'high');
freqz(b,a);
figure;
axis([0,1,-120,0]);
y=filter(b,a,x);
Y=fftshift(fft(y,N));
stem(f,abs(Y));
右图一为数字高通滤波器的幅频响应
和相频响应。
右图二为信号经过滤波器后输出信号的频谱。
(3)
clear;
fsam=100;
t0=1/fsam;t=6;k=0:
t0:
t;
N=512;
x=1+cos(pi/4.*k/t0)+cos(2*pi/3.*k/t0);
f=(-N/2:
(N/2-1))/N*2*pi;
[N1,Wc]=buttord([1/88/24],[1/322/3],3,60);
[b,a]=butter(N1,Wc,'bandpass');
freqz(b,a);
figure;
axis([0,1,-120,0]);
y=filter(b,a,x);
Y=fftshift(fft(y,N));
stem(f,abs(Y));
右图一为数字带通滤波器的幅频响应
和相频响应。
右图二为信号经过滤波器后输出信号的频谱。
(4)
clear;
fsam=100;
t0=1/fsam;t=6;k=0:
t0:
t;
N=512;
x=1+cos(pi/4.*k/t0)+cos(2*pi/3.*k/t0);
f=(-N/2:
(N/2-1))/N*2*pi;
[N1,Wc]=buttord([1/322/3],[0.210/27],3,60);
[b,a]=butter(N1,Wc,'stop');
freqz(b,a);
figure;
axis([0,1,-120,0]);
y=filter(b,a,x);
Y=fftshift(fft(y,N));
stem(f,abs(Y));
右图一为数字带阻滤波器的幅频响应
和相频响应。
右图二为信号经过滤波器后输出信号的频谱。
2.
(1).脉冲响应不变法
clear;
omegas=[0.20.72]*pi;
omegap=[0.30.6]*pi;
Ap=1;As=42;
Fs=1;
ws=omegas*Fs;
wp=omegap*Fs;
[N,Wc]=buttord(wp,ws,Ap,As,'s');
[b,a]=butter(N,Wc,'bandpass','s');
%[N,Wc]=cheb1ord(wp,ws,Ap,As,'s');
%[b,a]=cheby1(N,Ap,Wc,'bandpass','s');
%[N,Wc]=cheb2ord(wp,ws,Ap,As,'s');
%[b,a]=cheby2(N,As,Wc,'bandpass','s');
%
%[N,Wc]=ellipord(wp,ws,Ap,As,'s');
%[b,a]=ellip(N,Ap,As,Wc,'bandpass','s');
[numd,dend]=impinvar(b,a,Fs);
w=linspace(0,pi,512);
h=freqz(numd,dend,w);
norm=max(abs(h));
numd=numd/norm;
plot(w/pi,20*log10(abs(h)/norm));
gridon
右图从上到下分别是buttord
Cheby1cheby2ellipord型滤波器的
频率特性(脉冲响应不变法)
(2)双线性变换法
clear;
omegas=[0.20.72]*pi;
omegap=[0.30.6]*pi;
Ap=1;As=42;
Fs=1;T=1/Fs;
wp=2*tan(omegap/2)/T;
ws=2*tan(omegas/2)/T;
%[N,Wc]=buttord(wp,ws,Ap,As,'s');
%[b,a]=butter(N,Wc,'bandpass','s');
%[N,Wc]=cheb1ord(wp,ws,Ap,As,'s');
%[b,a]=cheby1(N,Ap,Wc,'bandpass','s');
%
%[N,Wc]=cheb2ord(wp,ws,Ap,As,'s');
%[b,a]=cheby2(N,As,Wc,'bandpass','s');
[N,Wc]=ellipord(wp,ws,Ap,As,'s');
[b,a]=ellip(N,Ap,As,Wc,'bandpass','s');
[numd,dend]=bilinear(b,a,Fs);
w=linspace(0,pi,512);
h=freqz(numd,dend,w);
norm=max(abs(h));
numd=numd/norm;
plot(w/pi,20*log10(abs(h)/norm));
gridon
title('ellipord')
右图从上到下分别是buttord
Cheby1cheby2ellipord型滤波器的
频率特性(双线性变换法)
3.
(1)
clear;
k=0:
204;N=205;
x=sin(2*pi*770/8000.*k)+sin(2*pi*1209/8000.*k);
stem(k,x);
title('时域波形')
figure;
X=fftshift(fft(x));
omega=((-(N-1)/2):
((N-1)/2))/N*2*pi;
stem(omega,abs(X));
title('频谱图')
(2)
clear;
k=0:
204;N=205;
x=sin(2*pi*770/8000.*k)+sin(2*pi*1209/8000.*k);
omega=((-(N-1)/2):
((N-1)/2))/N*2*pi;
%低通
wp=2*941/8000;ws=2*1200/8000;Ap=3;As=60;
[N,wc]=buttord(wp,ws,Ap,As);
[b,a]=butter(N,wc);
freqz(b,a);
figure
%滤波输出
y=filter(b,a,x);
Y=fftshift(fft(y));
stem(omega,abs(Y))
clear;
k=0:
204;N=205;
x=sin(2*pi*770/8000.*k)+sin(2*pi*1209/8000.*k);
omega=((-(N-1)/2):
((N-1)/2))/N*2*pi;
%高通
wp=2*1200/8000;ws=2*1000/8000;Ap=3;As=60;
[N,wc]=buttord(wp,ws,Ap,As);
[b,a]=butter(N,wc,'high');
freqz(b,a);
figure
%滤波输出
y=filter(b,a,x);
Y=fftshift(fft(y));
stem(omega,abs(Y))
(3)
clear;
k=0:
204;N=205;
x=sin(2*pi*770/8000.*k)+sin(2*pi*1209/8000.*k);
omega=((-(N-1)/2):
((N-1)/2))/N*2*pi;
wp1=[692702]./4000;ws1=[687707]./4000;
wp2=[765775]./4000;ws2=[760780]./4000;
wp3=[847857]./4000;ws3=[842862]./4000;
wp4=[936946]./4000;ws4=[931951]./4000;
Ap=3;As=30;
[N1,wc1]=buttord(wp1,ws1,Ap,As);
[b1,a1]=butter(N1,wc1,'bandpass');
[N2,wc2]=buttord(wp2,ws2,Ap,As);
[b2,a2]=butter(N2,wc2,'bandpass');
[N3,wc3]=buttord(wp3,ws3,Ap,As);
[b3,a3]=butter(N3,wc3,'bandpass');
[N4,wc4]=buttord(wp4,ws4,Ap,As);
[b4,a4]=butter(N4,wc4,'bandpass');
figure
(1)
freqz(b1,a1);
figure
(2)
freqz(b2,a2);
figure(3)
freqz(b3,a3);
figure(4)
freqz(b4,a4);
以上分别是4个带通滤波器的幅频响应.
clear;
k=0:
204;N=205;
x=sin(2*pi*770/8000.*k)+sin(2*pi*1209/8000.*k);
omega=((-(N-1)/2):
((N-1)/2))/N*2*pi;
wp1=[692702]./4000;ws1=[687707]./4000;
wp2=[765775]./4000;ws2=[760780]./4000;
wp3=[847857]./4000;ws3=[842862]./4000;
wp4=[936946]./4000;ws4=[931951]./4000;
Ap=3;As=20;
[N1,wc1]=buttord(wp1,ws1,Ap,As);
[b1,a1]=butter(N1,wc1,'bandpass');
[N2,wc2]=buttord(wp2,ws2,Ap,As);
[b2,a2]=butter(N2,wc2,'bandpass');
[N3,wc3]=buttord
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 华北电力 大学 数字信号 处理