Matlab fir iir滤波器的比较设计程序及图片原创文档格式.docx
- 文档编号:1277662
- 上传时间:2023-04-30
- 格式:DOCX
- 页数:19
- 大小:92.58KB
Matlab fir iir滤波器的比较设计程序及图片原创文档格式.docx
《Matlab fir iir滤波器的比较设计程序及图片原创文档格式.docx》由会员分享,可在线阅读,更多相关《Matlab fir iir滤波器的比较设计程序及图片原创文档格式.docx(19页珍藏版)》请在冰点文库上搜索。
IIR
fpl=10000;
wpl=2*pi*fpl/Fs;
wsl=2*pi*fsl/Fs;
wsu=2*pi*fsu/Fs;
wpu=2*pi*fpu/Fs;
wp=[wpl/pi,wpu/pi];
ws=[wsl/pi,wsu/pi];
[N,Wpo]=ellipord(wp,ws,Rp,Rs);
%求模拟带阻滤波器的阶数和通带截止频率
[b,a]=ellip(N,Rp,Rs,Wpo,'
stop'
);
w=0:
0.005*pi:
pi;
[h,w]=freqz(b,a,w);
阶数N=6
clc;
clear;
hw1=20*log10(abs(h));
plot((0:
199)/200,hw1(1:
200));
gridon;
xlabel('
\omega/\pi'
ylabel('
幅度/dB'
%也可用plot(20*log10(abs(h)));
thet=unwrap(angle(h));
figure;
199)/200,thet(1:
相频响应'
x=sin(2*pi*(1:
256)*7000/Fs);
%输入信号:
0<
=n<
=255,f0=7kHz;
yi=filter(b,a,x);
plot(yi(1:
250));
n'
y(n)'
幅频响应
相频响应
系统输出
f0=12kHz(在阻带)
f0=7kHz(在通带)
FIR
B1=wpu-wsu;
B2=wsl-wpl;
B=min(B1,B2);
NF=ceil(6.6*pi/B);
wp1=[(fpl+fsl)/Fs,(fpu+fsu)/Fs];
h1=fir1(NF-1,wp1,'
[hf,w]=freqz(h1,1,w);
freqz(h1,1,w);
hw2=20*log10(abs(hf));
199)/200,hw2(1:
thet2=unwrap(angle(hf));
199)/200,thet2(1:
256)*12000/Fs);
yi=filter(h1,1,x);
阶数NF=166
(1)f0=12kHz
(2)f0=7kHz(在通带中)
Fs与滤波器阶数的关系
(1)Fs=50kHz
N=6
(2)Fs=100kHz
N=7
NF=331
(3)Fs=150kHz
NF=496
(4)Fs=2MHz
NF=6601
验证了书上
(1)P251中间几行“无论是设计IIR数字滤波器还是FIR。
。
”
(2)P251第三行IIR数字滤波器随着Fs的提高,阶数提高不明显。
IIR滤波器频率选择性越好,其相位的非线性越严重。
fsl=10250;
fsu=20750;
fsl=12000;
fsu=19000;
fsl=13000;
fsu=18000;
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab fir iir滤波器的比较设计程序及图片原创 iir 滤波器 比较 设计 程序 图片 原创