毕业论文论文范文课程设计实践报告法律论文英语论文教学论文医学论文农学论文艺术论文行政论文管理论文计算机安全
您现在的位置: 毕业论文 >> 报告 >> 正文

数字巴特沃思滤波器的设计

更新时间:2008-12-21:  来源:毕业论文

数字巴特沃思滤波器的设计
实验目的 通过本次实验掌握数字巴特沃思滤波器的设计。
实验要求 在掌握巴特沃斯滤波器的原理的基础上,通过MATLAB软件进行数字巴特沃斯滤波器的设计。

实验原理 1. 利用原型设计法,对满足指标要求的模拟低通原型滤波器,通过双线性变换法,得到满足指标要求的数字巴特沃思低通滤波器;
2. 双线性变换法原理参见《数字信号处理》教材。
实验步骤 1、编写程序。设计数字巴特沃思低通滤波器dbutter.m文件。
输入参数:  Fs
输出参数:系数
2、编写程序。通过调用dbutter.m函数,设计一个数字低通滤波器。
   绘出零极点分布图、单位取样响应(可调用库函数impz)、幅频响应、相位响应。
3、用help查看内部函数butter.m,了解调用格式,并用此函数重新设计2中要求的滤波器。
4、利用fdatool设计各种类型的数字滤波器。


实验内容 1、设计巴特沃思低通滤波器
2、利用MATLAB工具fdatool设计各种类型的滤波器。
实验数据 编写的设计数字巴特沃思低通
波器的dbutter.m文件如下:
冲激响应不变法:
function
[b,a]=dbutter1(wp,ws,ap,as,fs)
T=1/fs;
np=wp*pi*fs;
ns=ws*pi*fs;
N1=(log10((10^(ap/10)-1)/…
(10^(as/10)-1)))/…
(2*log10(np/ns));
N=ceil(N1);  
nc=np/((10^(ap/10)-1)^…
(1/(2*N)));
for k=1:N
s(k)=nc*exp(j*pi*((2*k+N-1)/…
(2*N)));
lwfree.cn   ,N);
for n=1:N
    A=1;
    for m=1:N
        if m~=n
            A=A*(s(n)-s(m));
            R(n)=(nc^N)/A;
       % else continue
        end
    end
end
[bb,aa]=residue(R,s,0);
[b,a]=impinvar(bb,aa,fs);
end
双线性变换法:
function
[b,a]=dbutter2(wp,ws,ap,as,fs)
T=1/fs;
np=(2*fs)*tan(wp*pi/2);
ns=(2*fs)*tan(ws*pi/2);
N1=(log10((10^(ap/10)-1)…
/(10^(as/…
10)-1)))/(2*log10(np/ns));
N=ceil(N1); 
nc=np/((10^(ap/10)-1)^…
(1/(2*N)));
for k=1:N
s(k)=nc*exp(j*pi*((2*k+N-1)/…
(2*N)));
end
R=zeros(1,N);
for n=1:N
    A=1;
    for m=1:N
        if m~=n
            A=A*(s(n)-s(m));
            R(n)=(nc^N)/A;
       % else continue
        end
    end
end
[bb,aa]=residue(R,s,0);
[b,a]=bilinear(bb,aa,fs);
end
>> [b1,a1]=dbutter2(0.2,0.3,1,15,100000);
lwfree.cn
>> [h1,w1]=freqz(b1,a1);
>> plot(w1,abs(h1),'LineWidth',3)
>> title('幅频响应'),grid on
>>  plot(w1,(180/pi)*angle(h1),'LineWidth',1.5)
>> grid on,title('相位响应');172

[1] [2] 下一页

数字巴特沃思滤波器的设计下载如图片无法显示或论文不完整,请联系qq752018766
设为首页 | 联系站长 | 友情链接 | 网站地图 |

copyright©lwfree.cn 六维论文网 严禁转载
如果本毕业论文网损害了您的利益或者侵犯了您的权利,请及时联系,我们一定会及时改正。