一、使用自编函数设计IIR滤波器

1、冲激响应法

(1)注给出的数字滤波器指标先化成模拟指标

(2)设计出模拟滤波器:

(3)使用冲激响应法转化成数字滤波器

(4)一个demo

clear all;
wp=0.2*pi; %数字指标
ws=0.3*pi;
Rp=1;
As=15;
T=1;Fs=1/T;
%冲激响应法
[cs,ds]=afd_butt(wp/T,ws/T,Rp,As);
[b,a]=imp_invr(cs,ds,T);
[C,B,A]=dir2par(b,a)
[db,mag,pha,grd,w]=freqz_m(b,a); figure(1)
subplot(2,2,1),plot(w/pi,mag);title('幅度');grid;axis([0,1,0,1.1]);
subplot(2,2,2),plot(w/pi,db);title('幅度 in db');grid;axis([0,1,-30,5]);
subplot(2,2,3),plot(w/pi,pha/pi);title('相位');grid;axis([0,1,-1,1]);
subplot(2,2,4),plot(w/pi,grd),title('群延时');grid;axis([0,1,0,10]);

 

2、双线性法设计IIR

(1)预畸变:

(2)设计低通滤波器

(3)转化成数字滤波器

(4)一个demo

%双线性法
OmegaP = (2/T)*tan(wp/2);
OmegaS = (2/T)*tan(ws/2); % Analog Butterworth Prototype Filter Calculation:
[cs,ds] = afd_butt(OmegaP,OmegaS,Rp,As);
% Bilinear transformation:
[b,a] = bilinear(cs,ds,Fs);
[C,B,A] = dir2cas(b,a)
[db,mag,pha,grd,w]=freqz_m(b,a);

 

二、MATLAB自带函数进行设计

1、模拟低通滤波器转数字滤波器设计

(1)根据数字指标设计出模拟滤波器的阶数n和截止频率Wn

(2)求出低通模拟原型滤波器

(3)由于所得的结果为零极点型,还必须转成b/a型,可用函数zp2tf

[bap,aap]=zp2tf(z,p,k)

(4)将模拟滤波器转换成数字滤波器

(低通,高通,带通,带阻)可用函数lp2lp, lp2hp,lp2bs, lp2bp。

(5)求出滤波器幅频,相频和冲激响应

(6)一个demo

Wp=0.2*pi;Rp=1;Ws=0.3*pi;As=15;
[n,Wn]=buttord(Wp,Ws,Rp,As,'s')
[z,p,k]=buttap(n);
[bap,aap]=zp2tf(z,p,k)
[b,a]=lp2lp(bap,aap,Wn);
[db,mag,pha,w]=freqs_m(b,a,0.5*pi);
[ha,x,t]=impulse(b,a);
figure(1)
subplot(2,2,1),plot(w/pi,mag);title('幅度');
subplot(2,2,2),plot(w/pi,db);title('幅度 in db');
subplot(2,2,3),plot(w/pi,pha/pi);title('相位');
subplot(2,2,4),plot(t,ha),title('冲激响应');
figure(2)
freqs(b,a);

 

2、直接进行数字滤波器设计

(1)根据给定指标得出,低通模拟原型滤波器的阶数和截止频率:

[n,wn]=buttord(wp,ws,Rp,As)
[n,wn]=cheb1ord(wp,ws,Rp,As)

(2)直接求出数字滤波器,利用函数butter,cheby1,cheby2。

[b,a]=butter(n,wn,’ftype’)
[b,a]=cheby1(n, Rp,wn,’ftype’) ;‘ftype’指的是:’low’,’bandpass’,’high’,’stop’

(3)求出滤波器的幅频,相频及群延时

[db,mag,pha,grd,w]=freqz_m(b,a);(P258)
直接画出幅频特性和相频特性:freqz(b,a) (matlab自带)

(4)一个demo

wp=0.2;Rp=1;ws=0.3;As=15;
[n,wn]=buttord(wp,ws,Rp,As)
[b,a]=butter(n,wn);
[db,mag,pha,grd,w]=freqz_m(b,a);
figure(3)
subplot(2,2,1),plot(w/pi,mag);title('幅度');grid; subplot(2,2,2),plot(w/pi,db);title('幅度 in db');grid; subplot(2,2,3),plot(w/pi,pha/pi);title('相位');grid;
subplot(2,2,4),plot(w/pi,grd),title('群延时');grid;
figure(4);freqz(b,a)

PS:

matlab中数字角频率进行了归一化,即0-1表示0-pi

频率以Hz给出,均表示模拟频率,因此利用matlab直接设计滤波器时,先要化成数字频率。

例:wp=[100 200]*2/fs

数字信号处理实验(五)——IIR滤波器的设计的更多相关文章

  1. 2018-2019-1-20165221&20165225 《信息安全系统设计》实验五:通讯协议设计

    2018-2019-1-20165221&20165225 <信息安全系统设计>-实验五:通讯协议设计 OpenSSL学习: 简介: OpenSSL是为网络通信提供安全及数据完整性 ...

  2. 数字信号处理实验(四)——数字滤波器结构

    一.滤波器结构 1.IIR滤波器 (1)系统函数   (2)差分方程   (3)级联形式:   (4)并联形式   2.FIR滤波器 (1)系统函数   (2)差分方程   (3)级联形式:   (4 ...

  3. 数字信号处理实验(六)——FIR滤波器的设计

    一.四种线性相位FIR滤波器的振幅响应 1.自编函数 [Hr,w,a,L]=-n) [Hr,w,a,L]=-n) [Hr,w,a,L]=-n) [Hr,w,a,L]=-n) 2.一个demo clea ...

  4. 数字信号处理实验(一)——DTFT

    1.MATLAB自编绘图函数 function [] = signal_write(X,w,flag) % X:数据 % w:频率向量 magX=abs(X);angX=angle(X); realX ...

  5. 数字信号处理实验(零)—— 一维声音信号处理和二维图像处理

    一.在matlab下声音信号的I/O 1.读wav文件函数 •y = wavread('filename') •[y,Fs,bits] = wavread('filename') •[...] = w ...

  6. IIR滤波器和FIR滤波器的区别与联系zz

      -------------------------------------------------------------------------------------------------- ...

  7. 《数字信号处理》课程实验2 – FIR数字滤波器设计

    一.FIR数字滤波器设计原理  本实验采用窗函数法设计FIR数字低通滤波器.我们希望设计的滤波器系统函数如下: \(H_{d}\left( e^{jw} \right) = \left\{ \begi ...

  8. 手把手教系列之IIR滤波器设计

    [导读]:在嵌入式系统中经常需要采集模拟信号,采集模拟信号的信号链中难免引入干扰,那么如何滤除干扰呢?今天就来个一步一步描述如何设计部署一个IIR滤波器到你的系统. 何为IIR滤波器? 无限冲激响应( ...

  9. IIR滤波器设计(调用MATLAB IIR函数来实现)

    转载请注明文章来源 – http://blog.csdn.net/v_hyx ,请勿用于任何商业用途         对于滤波器设计,以前虽然学过相关的理论(现代数字信号处理和DSP设计),但一直不求 ...

随机推荐

  1. Spring MVC配置静态资源的正常访问

    SpringMVC如果过滤器过滤范围配置了/或者/*,那么框架会过滤所有请求,包括自己写的请求和静态资源请求,这样静态资源就不能正常加载,包括js文件.css文件.图片资源访问的时候都会出现404页面 ...

  2. ACM/ICPC 之 昂贵的聘礼-最短路解法(POJ1062)

    //转移为最短路问题,枚举必经每一个不小于酋长等级的人的最短路 //Time:16Ms Memory:208K #include<iostream> #include<cstring ...

  3. 【leetcode】Happy Number(easy)

    Write an algorithm to determine if a number is "happy". A happy number is a number defined ...

  4. 【编程题目】在从 1 到 n 的正数中 1 出现的次数

    30.在从 1 到 n 的正数中 1 出现的次数(数组)题目:输入一个整数 n,求从 1 到 n 这 n 个整数的十进制表示中 1 出现的次数.例如输入 12,从 1 到 12 这些整数中包含 1 的 ...

  5. ajax,下拉框级联

    js代码: $(document).ready(function() { $("#type1").change(function(){ var type1Code=$(" ...

  6. yum Error: Cannot retrieve metalink for repository: epel. Please verify its path and try again

    今天使用yum时出现如上问题,解决办法: 1.编辑:/etc/yum.repos.d/epel.repo 将baseurl的注释取消, mirrorlist注释掉,如图:

  7. [Android Pro] Gradle Tips#1-tasks

    reference to : http://trickyandroid.com/gradle-tip-1-tasks/ http://blog.csdn.net/lzyzsd/article/deta ...

  8. ajax实例1

    前台: function getDetail(index){ $.post("<%=request.getContextPath() %>/member/dbcenter!get ...

  9. Jmeter 提取http请求返回值里json数据参数化方法

    第三方插件下载地址:http://jmeter-plugins.org/downloads/all/ 插件下载后解压:找到JMeterPlugins-Extras.jar,把JMeterPlugins ...

  10. 【Java EE 学习 19】【使用过滤器实现全站压缩】【使用ThreadLocal模式解决跨DAO事务回滚问题】

    一.使用过滤器实现全站压缩 1.目标:对网站的所有JSP页面进行页面压缩,减少用户流量的使用.但是对图片和视频不进行压缩,因为图片和视频的压缩率很小,而且处理所需要的服务器资源很大. 2.实现原理: ...