基于matlab的带通带阻滤波器设计实例
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
基于的带通、带阻滤波器设计实例
以下两个滤波器都是切比雪夫I型数字滤波器,不是巴特沃尔滤波器,请使用者注意!
1.带通滤波器
(13)
%带通滤波
%使用注意事项:通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半%即,f131,的值小于 2
:需要带通滤波的序列
% f 1:通带左边界
% f 3:通带右边界
% 1:衰减截止左边界
% :衰变截止右边界
:边带区衰减数设置
:截止区衰减数设置
:序列x的采样频率
% f1=3003=500通带截止频率上下限
% 200600阻带截止频率上下限
% 0.130通带边衰减值和阻带边衰减值
% 2000采样率
%
1=2**f1;
3=2**f3;
2**;
2**;
[1 3];
[ ];
%
% 设计切比雪夫滤波器;
[]1();
[11]1();
%查看设计滤波器的曲线
[](11,256);
20*10((h));
()('所设计滤波器的通带曲线') ;
(11);
带通滤波器使用例子
%带通滤波器测试程序
2000;
(1);
1=100;
2=400;
3=700;
(2**1*t)(2**2*t)(2**3*t);
;
(211)();
(212)(,1);
% (11);
(x,300,500,200,600,0.1,30);
;
(211)();
(212)(,1);
%调用到的()函数代码如下
()
%当1,画幅值谱;当2,画功率谱;当其他的,那么花幅值谱和功率谱
%当1时,还可以多输入2个可选参数
%可选输入参数是用来控制需要查看的频率段的
%第一个是需要查看的频率段起点
%第二个是需要查看的频率段的终点
%其他不具备可选输入参数,如果输入发生位置错误
2^2((y))找出大于y的个数的最大的2的指数值(自动进算最佳步长)1024人为设置的步长
(y)去除直流分量
()对y信号进行,得到频率的幅值分布
.*()()函数是求y函数的共轭复数,实数的共轭复数是他本身。
*(02-1);�T变换后对应的频率的序列
% .*()()函数是求y函数的共轭复数,实数的共轭复数是他本身。
1
3
(,2*((12))(y))的帮助里画的方法
('幅值')('频率')('信号幅值谱');
(((12)))论坛上画的方法
f1{1};
{2};
(f1 * 1);
( * 1);
(()(()*2));
2
((12));
('功率谱密度')('频率')('信号功率谱');
(211)(,2*((12))(y));
('幅值')('频率')('信号幅值谱');
(212)((12));
('功率谱密度')('频率')('信号功率谱');
运行结果如下图,第一幅是滤波前测试信号的时频图,第二幅是滤波器的滤波曲线图,第三幅是经滤波后的测试信号时频图。
2.带阻滤波器
(13)
%带阻滤波
%使用注意事项:通带或阻带的截止频率与采样率的选取范围是不能超过采样率的一半%即,f131,的值小于 2
:需要带通滤波的序列
% f 1:通带左边界
% f 3:通带右边界
% 1:衰减截止左边界
% :衰变截止右边界
:边带区衰减数设置
:截止区衰减数设置
:序列x的采样频率
% f1=3003=500通带截止频率上下限
% 200600阻带截止频率上下限
% 0.130通带边衰减值和阻带边衰减值
% 2000采样率
%
1=2**f1;
3=2**f3;
2**;
2**;
[1 3];
[ ];
%
% 设计切比雪夫滤波器;
[]1();
[11]1(,'');
%查看设计滤波器的曲线
[](11,256);
20*10((h));
()('所设计滤波器的通带曲线') ;
(11);
使用例子
%带阻滤波器测试
1000;
(1);
(2**100*t)(2**150*t)(2**200*t);
(,1);
(y,110,190,140,160,0.1,30);
(,1);
运行结果如下图,第一幅是滤波前测试信号的频谱图,第二幅是滤波器的滤波曲线图,第三幅是经滤波后的测试信号频谱图。