数字信号处理实验三
数字信号处理实验报告 3
数字信号处理实验报告姓名：班级：通信学号：实验名称：频域抽样定理验证实验类型：验证试验指导教师:实习日期：2013.频域采样定理验证实验一． 实验目的：1. 加深对离散序列频域抽样定理的理解2.了解由频谱通过IFFT 计算连续时间信号的方法3.掌握用MATLAB 语言进行频域抽样与恢复时程序的编写方法 4、用MATLAB 语言将X(k)恢复为X(z)及X(e jw )。
二． 实验原理：1、1、频域采样定理： 如果序列x(n)的长度为M ，频域抽样点数为N ，则只有当频域采样点数N ≥M 时，才有x N (n)=IDFT[X(k)]=x(n),即可由频域采样X(k)无失真的恢复原序列 x(n)。
2、用X(k)表示X(z)的内插公式：∑-=-----=10111)(1)(N k kNNzWz k X Nz X内插函数： zWzkNNN z 1k111)(-----=ϕ频域内插公式：∑-=-=10)2()()(N K j k Nk X e X πωϕω频域内插函数：e N j N N )21()2sin()2sin(1)(--=ωωωωϕ三． 实验任务与步骤：实验一：长度为26的三角形序列x(n)如图(b)所示，编写MATLAB 程序验证频域抽样定理。
实验二：已知一个时间序列的频谱为X(e jw )=2+4e -jw +6e -j2w +4e -j3w +2e -j4w分别取频域抽样点数N为3、5和10，用IPPT计算并求出其时间序列x(n)，用图形显示各时间序列。
由此讨论原时域信号不失真地由频域抽样恢复的条件。
实验三：由X32(k)恢复X(z)和X(e jw)。
四．实验结论与分析：实验一：源程序：M=26;N=32;n=0:M; %产生M长三角波序列x(n)xa=0:floor(M/2);xb= ceil(M/2)-1:-1:0; xn=[xa,xb];Xk=fft(xn,512); %1024点FFT[x(n)], 用于近似序列x(n)的TFX32k=fft(xn,32); %32点FFT[x(n)]x32n=ifft(X32k); %32点IFFT[X32(k)]得到x32(n)X16k=X32k(1:2:N); %隔点抽取X32k得到X16(K)x16n=ifft(X16k,N/2); %16点IFFT[X16(k)]得到x16(n)subplot(3,2,2);stem(n,xn,'.');box ontitle('(b) 三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,20])k=0:511;wk=2*k/512;subplot(3,2,1);plot(wk,abs(Xk));title('(a)FT[x(n)]');xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|');axis([0,1,0,200])k=0:N/2-1;subplot(3,2,3);stem(k,abs(X16k),'.');box ontitle('(c) 16点频域');xlabel('k');ylabel('|X_1_6(k)|');axis([0,8,0,200])n1=0:N/2-1;subplot(3,2,4);stem(n1,x16n,'.');box ontitle('(d) 16点IDFT[X_1_6(k)]');xlabel('n');ylabel('x_1_6(n)');axis([0,32,0,20])k=0:N-1;subplot(3,2,5);stem(k,abs(X32k),'.');box ontitle('(e) 32点频域采样');xlabel('k');ylabel('|X_3_2(k)|');axis([0,16,0,200])n1=0:N-1;subplot(3,2,6);stem(n1,x32n,'.');box ontitle('(f) 32点IDFT[X_3_2(k)]');xlabel('n');ylabel('x_3_2(n)');axis([0,32,0,20])结果如下所示：实验一分析：序列x(n)的长度M=26，由图中可以看出，当采样点数N=16<M时，x16(n)确实等于原三角序列x(n)以16为周期的周期延拓序列的主值序列。
202数字信号处理实验讲义 - 实验三
实验三 信号检测与分析一、实验目的1．掌握FFT 函数的用法。
2. 利用FFT 进行信号检测及谱分析。
3．了解信号截取长度对谱分析的影响。
二、实验内容1．利用FFT 计算信号功率谱。
2. 进行信号检测。
分析信号频谱所对应频率轴的数字频率和频率之间的关系。
3. 对例2，进一步增加截取长度和FFT 点数，如N 加大到256，观察信号频谱的变化，分析产生这一变化的原因。
在截取长度不变的条件下改变采样频率，观察信号频谱的变化，分析产生这一变化的原因。
4. 对例3，加大噪声到2*randn(1,N)和8*randn(1,N)，画出并比较不同噪声下时域波形和频谱。
5. 用一个N 点FFT 计算两个长度为N 的实序列N 点离散傅里叶变换，并将结果和直接使用两个N 点FFT 得到的结果进行比较。
6．比较DFT 和FFT 的运算时间。
(计时函数 tic ， toc)7. 比较利用FFT 计算线性卷积的和直接计算线性卷积的时间及结果。
(选作)8．对给定语音信号进行谱分析，写出采样频率，画出语音信号的波形及频谱，并分析语音信号的频率分布特点。
三、 实验报告要求1．记录例题程序的实验结果、图形。
2．写出自己编写的程序并记录结果、图形。
3．对实验结果进行分析。
% fft 一维快速傅立叶变换函数% 格式： y=fft(x) FFT 算法计算矢量x 的离散傅立叶变换，当x 为矩阵时，y 为矩阵x 的每一列 % 的FFT 。
%当x 的长度为2的幂次方时，则fft 采用基2的FFT 算法，否则采用稍慢的混合基算法。
% y=fft(x,n) 实现n 点FFT 。
当x 的长度小于n 时，fft 函数在x 的尾部补零，以构成n 点数据；% 当x 的长度大于n 时，fft 函数会截断序列x 。
当x 为矩阵时，fft 函数按类似的方法处理列长度。
% ifft 一维快速傅立叶反变换% 格式：y=ifft(x) 用于计算矢量x 的IFFT 。
数字信号处理实验3
实验3 离散时间系统的频域分析一、实验目的（1）了解DFS 、DFT 与DTFT 的联系；加深对FFT 基本理论的理解；掌握用MATLB 语言进行傅里叶变换时常用的子函数；（2）了解离散系统的零极点与系统因果性和稳定性的关系；加深对离散系统的频率响应特性基本概念的理解；熟悉MATLAB 中进行离散系统零极点分析的常用子函数；掌握离散系统幅频响应和相频响应的求解方法。
二、实验内容1. 已知离散时间系统函数为 用matlab 中的函数()432143213.07.05.11.112.01.03.01.02.0--------+-+-++++=zz z z z z z z z H 求该系统的零极点及零极点分布图，并判断系统的因果稳定性。
方法一：利用tf3zp 函数b=[0.2 0.1 0.3 0.1 0.2]; a=[1 -1.1 1.5 -0.7 0.3]; [z,p,k]=tf2zp(b,a); c1=abs(z);c2=angle(z); c3=abs(p);c4=angle(p); polar(c4,c3,'bx') hold onpolar(c2,c1,'ro') disp(z) disp(p)disp(abs(z)) disp(abs(p))90270方法二：利用zplaneb=[0.2 0.1 0.3 0.1 0.2];a=[1 -1.1 1.5 -0.7 0.3];z=roots(b);p=roots(a);zplane(b,a)disp(z)disp(p)disp(abs(z))disp(abs(p))-1-0.500.51-1-0.8-0.6-0.4-0.200.20.40.60.81Real PartI m a g i n a r y P a r t由于极点都在单位圆内，故该系统稳定。
若其收敛域为圆外区域，则系统是因果系统。
2. 已知离散时间系统的系统函数为()432143213.07.05.11.112.01.03.01.02.0--------+-+-++++=z z z z z z z z z H求该系统在π~0频率范围内的绝对幅频响应、相频响应。
(完整版)数字信号处理实验三
3.41；3.42 由教材可知： ，即序列的偶部分的傅立叶变换是序列的傅立叶变换的实部。
5、实验步骤
1、进行本实验，首先必须熟悉matlab的运用，所以第一步是学会使用matlab。
2、学习相关基础知识，根据《数字信号处理》课程的学习理解实验内容和目的。
plot(w/pi,angle(h1));grid
xlabel('\omega/\pi');ylabel('以弧度为单位的相位');
title('原序列的相位谱')
subplot(2,2,4)
plot(w/pi,angle(h2));grid
xlabel('\omega/\pi');ylabel('以弧度为单位的相位');
grid;
title('相位谱arg[H(e^{j\omega})]');
xlabel('\omega/\pi');
ylabel('以弧度为单位的相位');
3.4
clf;
w=-4*pi:8*pi/511:4*pi;
num1=[1 3 5 7 9 11 13 15 17];
h=freqz(num,1,w);
Q3.32 通过加入合适的注释语句和程序语句，修改程序P3.8，对程序生成的图形中的两个轴加标记。时移量是多少？
Q3.33 运行修改后的程序并验证离散傅里叶变换的圆周时移性质。
Q3.36 运行程序P3.9并验证离散傅里叶变换的圆周卷积性质。
Q3.38 运行程序P3.10并验证线性卷积可通过圆周卷积得到。
数字信号处理实验3
实验报告课程名称数字信号处理实验名称利用频率抽样法设计线性相位数字FIR 带通滤波器指导教师罗倩实验日期2013年12月15日学院信息与通信工程专业电子信息工程学生姓名于卓立班级/学号电信1103/ 2011010533 成绩指导老师签字实验三、利用频率抽样法设计线性相位数字 FIR 带通滤波器一、实验目的：1、掌握利用频率抽样法设计FIR数字滤波器的原理和方法；2、了解频率抽样对滤波特性的影响。
二、题目：设计数字带通滤波器，指标为低阻带：ω1s=0.1π， As=60dB低通带：ω1p=0.3π， Rp=0.1dB高通带：ω2p=0.4π， Rp=0.1dB高阻带：ω2s=0.6π， As=60dB画出滤波器的频率响应曲线和单位冲激响应曲线，给出实际滤波器长度、通带衰减和阻带衰减。
（选作：）仿真0-9 DTMF 信号；播放拨号音；在所设计的带通滤波器基础上（或用Goertzel算法），改变截止频率，构成一组带通滤波器组，识别拨号数字，并显示出来；考虑到线路质量，在双音多频信号上加入噪声，再次识别。
设计步骤参见p268图7-14。
三、实验流图四、实验结果五、参考文献[1] 焦瑞莉,罗倩,汪毓铎,顾奕，信号与系统引论，北京：高等教育出版社，2009年3月.。
[2] 胡广书，数字信号处理——理论、算法与实现，北京：清华大学出版社，2003年8月。
六、设计结果分析实验代码在修改后较符合要求，在0.1之前没有超过-20DB的信号，通带、阻带、过渡带都与设计要求相差不大。
七、实验心得做完这次实验，实验者进一步的掌握了matlab设计程序的一些方法，虽然由于审题疏忽，实验者吧频率抽样法一开始做成了窗函数法，但是在老师的现场讲解下下，还是及时改成了错误，修改代码，完成实验。
附录clear allclcclfN=1:40;H(1:3)=0;H(4)=0.1;H(5)=0.55;H(6:9)=1;H(10)=0.5;H(11)=0.1;H(12:31)=0;H(32)=0.1;H(33)=0.5;H(34:37)=1;H(38)=0.5;H(39)=0.1;H(40)=0;figure(1)stem(N*2/40,H);phi=-1*(40-1)*pi*N/40;c=1*1j.*phi;Y=H.*exp(c);y=ifft(Y,40);%M=1:400;w=0.01:0.001:pi;HW=freqz(y,[1],w);Y1=abs(HW);Y2=20*log10(Y1);figure(2);plot(w/pi,Y2);title('频率响应曲线');xlabel('w/pi' )ylabel('db') grid onfigure(3);t=1:40;stem(t,y);title('冲击响应h(n)曲线');xlabel('t' )ylabel('h(n)') grid on。
数字信号处理实验报告 3
数字信号处理实验报告姓名：班级：通信学号：实验名称：频域抽样定理验证实验类型：验证试验指导教师:实习日期：2013.频域采样定理验证实验一． 实验目的：1. 加深对离散序列频域抽样定理的理解2.了解由频谱通过IFFT 计算连续时间信号的方法3.掌握用MATLAB 语言进行频域抽样与恢复时程序的编写方法 4、用MATLAB 语言将X(k)恢复为X(z)及X(e jw )。
二． 实验原理：1、1、频域采样定理： 如果序列x(n)的长度为M ，频域抽样点数为N ，则只有当频域采样点数N ≥M 时，才有x N (n)=IDFT[X(k)]=x(n),即可由频域采样X(k)无失真的恢复原序列 x(n)。
2、用X(k)表示X(z)的内插公式：∑-=-----=10111)(1)(N k kNN zWz k X Nz X内插函数： zWzkNNN z 1k111)(-----=ϕ频域内插公式：∑-=-=10)2()()(N Kj k Nk X e X πωϕω频域内插函数：e N j N N )21()2sin()2sin(1)(--=ωωωωϕ三． 实验任务与步骤：实验一：长度为26的三角形序列x(n)如图(b)所示，编写MATLAB 程序验证频域抽样定理。
实验二：已知一个时间序列的频谱为X(e jw )=2+4e -jw +6e -j2w +4e -j3w +2e -j4w分别取频域抽样点数N为3、5和10，用IPPT计算并求出其时间序列x(n)，用图形显示各时间序列。
由此讨论原时域信号不失真地由频域抽样恢复的条件。
实验三：由X32(k)恢复X(z)和X(e jw)。
四．实验结论与分析：实验一：源程序：M=26;N=32;n=0:M; %产生M长三角波序列x(n)xa=0:floor(M/2);xb= ceil(M/2)-1:-1:0; xn=[xa,xb];Xk=fft(xn,512); %1024点FFT[x(n)], 用于近似序列x(n)的TF X32k=fft(xn,32); %32点FFT[x(n)]x32n=ifft(X32k); %32点IFFT[X32(k)]得到x32(n)X16k=X32k(1:2:N); %隔点抽取X32k得到X16(K)x16n=ifft(X16k,N/2); %16点IFFT[X16(k)]得到x16(n)subplot(3,2,2);stem(n,xn,'.');box ontitle('(b) 三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,20])k=0:511;wk=2*k/512;subplot(3,2,1);plot(wk,abs(Xk));title('(a)FT[x(n)]');xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|');axis([0,1,0,200])k=0:N/2-1;subplot(3,2,3);stem(k,abs(X16k),'.');box ontitle('(c) 16点频域');xlabel('k');ylabel('|X_1_6(k)|');axis([0,8,0,200])n1=0:N/2-1;subplot(3,2,4);stem(n1,x16n,'.');box ontitle('(d) 16点IDFT[X_1_6(k)]');xlabel('n');ylabel('x_1_6(n)');axis([0,32,0,20]) k=0:N-1;subplot(3,2,5);stem(k,abs(X32k),'.');box on title('(e) 32点频域采样');xlabel('k'); ylabel('|X_3_2(k)|');axis([0,16,0,200]) n1=0:N-1;subplot(3,2,6);stem(n1,x32n,'.');box on title('(f) 32点IDFT[X_3_2(k)]');xlabel('n'); ylabel('x_3_2(n)');axis([0,32,0,20])结果如下所示：实验一分析：序列x(n)的长度M=26，由图中可以看出，当采样点数N=16<M 时，x 16(n)确实等于原三角序列x(n)以16为周期的周期延拓序列的主值序列。
数字信号处理实验报告三
实验三 抽样一、实验原理抽样过程基于两个基本原理：混叠合重建，涉及正弦波合现行调频信号的混叠。
二、 实验内容.3.3.1抽样引起的混叠 由于在MATLAB 中不能产生模拟信号，实验需要做实时t 轴的仿真。
因此，把仿真时的△t 与所研究的抽样周期Ts 明确地区分开始很重要的。
1. 正弦信号混叠对连续时间正弦信号考虑下面表达式：()()φπ+=t f t 02sin x可以按抽样频率ST s f 1=对()t x 抽样来获得离散时间信号[]()()⎪⎪⎭⎫ ⎝⎛+=====φπn f f t x t x s f n t nT t ss 0/2sin n x 2. 实验内容以不同组合的0f 和s f 绘出[]n x ，可以说明混叠问题。
以下，取抽样频率s f =8kHz.a. 首先，绘出一个被抽样的正弦波的单图。
令正弦波的频率为300Hz ，然后在10ms 长间隔上抽样。
相位φ可以任意指定。
使用stem 绘出产生的离散时间信号。
因为是用眼睛实现重建可视化信号包络，应该很容易看到正弦信号的轮廓。
b. 如果必要，使用plot 绘图。
在这种情况下，点用直线段连接 起来， 这样正弦信号的特点应该是明显的。
用直线段连接信号 样本是自离散时间样本产生连续时间信号的“信号重建”的 一 种方式。
它 不是抽样定理所说的理想重建，但对大多数的情形它已经是足够好，很有用。
c. 把正弦的频率从100Hz 变至475Hz,每次增加125Hz ，如在()a 部 分中那样，绘出一系列相应的图。
注意，正如所预期的那样， 显 现的正弦信号的频率在逐渐增加。
最好用subplot 指令把四 个图放在同一屏上。
d. 把正弦的频率从7525Hz 变至7900Hz ，每次增加125Hz ，正 如在()c 中那样，另外绘出一系列相应的图。
注意现在显现的正弦信号的频率在逐渐减少。
解释这一现象。
e.把正弦的频率从32100Hz变至32475Hz，每次增加125Hz，再次绘出一些列类似的图。
数字信号处理--实验三
一、实验目的1．了解工程上两种最常用的变换方法：脉冲响应不变法和双线性变换法。
2．掌握双线性变换法设计IIR 滤波器的原理及具体设计方法，熟悉用双线性设计法设计低通、带通和高通IIR 数字滤波器的计算机程序.3．观察用双线性变换法设计的滤波器的频域特性，并与脉冲响应不变法相比较,了解双线性变换法的特点。
4．熟悉用双线性变换法设计数字Butterworth 和Chebyshev 滤波器的全过程。
5．了解多项式乘积和多项式乘方运算的计算机编程方法。
二、实验原理与方法从模拟滤波器设计IIR 数字滤波器具有四种方法：微分－差分变换法、脉冲响应不变法、双线性变换法、z 平面变换法。
工程上常用的是其中的两种:脉冲响应不变法、双线性变换法。
脉冲响应不变法需要经历如下基本步骤：由已知系统传输函数H(S)计算系统冲激响应h(t)；对h(t)等间隔采样得到h （n ）=h （n T)；由h （n ）获得数字滤波器的系统响应H （Z)。
这种方法非常直观，其算法宗旨是保证所设计的IIR 滤波器的脉冲响应和模拟滤波器的脉冲响应在采样点上完全一致。
而双线性变换法的设计准则是使数字滤波器的频率响应与参考模拟滤波器的频率响应相似。
脉冲响应不变法一个重要的特点是频率坐标的变换是线性的（），其确定是有频谱的周期延拓效应，存在频谱混叠的现象。
为了克服脉冲响应不变法可能产生的频谱混叠，提出了双线性变换法，它依靠双线性变换式：， ， 其中 ,建立其S 平面和Z 平面的单值映射关系，数字域频率和模拟域频率的关系是： , （3－1） 由上面的关系式可知,当时，终止在折叠频率处,整个轴单值的对应于单位圆的一周。
因此双线性变换法不同于脉冲响应不变法，不存在频谱混叠的问题。
从式（3－1）还可以看出，两者的频率不是线性关系。
这种非线性关系使得通带截至频率、过渡带的边缘频率的相对位置都发生了非线性畸变。
这种频率的畸变可以通过预畸变来校正。
用双线性变换法设计数字滤波器时，一般总是先将数字滤波器的个临界频率经过式（3－1）的频率预畸变,求得相应参考模拟滤波器的个临界频率，然后设计参考模拟滤波器的传递函数，最后通过双T Ω=ω1111--+-=z z s s s z -+=11Ω+=j s σωj re z =)2/(ωtg =Ω)(2Ω=arctg ω∞→Ωωπω=Ωj线性变换式求得数字滤波器的传递函数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%设计butterworth模拟低通滤波器
③设计chebyshevI型模拟低通滤波器
[z,p,k]=cheb1ap(n,Rp) %设计chebyshevI型模拟低通滤波器
%Rp：通带内的波纹系数，单位分贝
④设计chebyshevII型模拟低通滤波器
[z,p,k]=cheb2ap(n,Rs) %设计chebyshevII型模拟低通滤波器
数字域：[n,Wn]=ellipord(Wp,Ws,Rp,Rs) 模拟域：[n,Wn]=ellipord(Wp,Ws,Rp,Rs,’s’) 注意：
n：返回符合要求性能指标的数字滤波器或模拟滤波器的最小阶数 Wn：滤波器的截至频率（即3db频率） Wp：通带的截至频率，Ws：阻带的截至频率，单位rad/s。且均为归一化频率，即。1对
实验三 用双线性变换法设计 IIR数字滤波器
实验目的
1. 掌握用双线性变换法设计IIR DF的原理 及具体设计方法，熟悉用双线性变换法设 计IIR DF的计算机编程。
2. 观察用双线性变换法设计的DF的频响特 性，了解双线性变换法的特点。
3. 熟悉用双线性变换法设计BW和CB型DF的 全过程。
[z,p,k]=butter(n,Wn)
[zp,k]=buter(n,Wn,’ftype’)
[A,B,C,D]=butter(n,Wn)
[A,B,C,D]=butter(n,Wn,’ftype’) 模拟域：[b,a]=butter(n,Wn,’s’)可设计出截止频率为Wn的n阶模拟butterworth
通、带通、带阻滤波器 运用脉冲响应不变法或双线性变换法把模拟滤波器转
换成数字滤波器
归一化频率
1)模拟滤波器的设计函数
①设计bessel模拟低通滤波器
[z,p,k]=besselap(n) %设计bessel模拟低通滤波器
②设计butterworth模拟低通滤波器
[z,p,k]=buttap(n)
%Rs：阻带内的波纹系数低于通带Rs分贝
⑤设计椭圆模拟滤波器
[z,p,k]=ellipap(n,Rp.Rs) %设计椭圆模拟滤波器
2）滤波器阶数的选择 下列函数除了能选择模拟滤波器的阶数外，同时也能选择数字滤波器的阶数。
①选择butterworth滤波器阶数 数字域：[n,Wn]=buttord(Wp,Ws,Rp,Rs) 模拟域：[n,Wn]=buttord(Wp,Ws,Rp,Rs,’s’) ②选择chebyshevI型滤波器阶数 数字域：[n,Wn]=cheb1ord(Wp,Ws,Rp,Rs) 模拟域：[n,Wn]=cheb1ord(Wp,Ws,Rp,Rs,’s’) ③选择chebyshevII型滤波器阶数 数字域：[n,Wn]=cheb2ord(Wp,Ws,Rp,Rs) 模拟域：[n,Wn]=cheb2ord(Wp,Ws,Rp,Rs,’s’) ④选择椭圆滤波器阶数
（4
H z Ha s s 2 1z1 T 1z1
IIR数字滤波器设计流程
设计IIR数字滤波器的一般步骤：
把给出的数字滤波器的性能指标转换为模拟滤波器的 性能指标
根据转换后的性能指标，通过滤波器阶数选择函数， 来确定滤波器的最小阶数N和固有频率Wn
由最小阶数N得到低通滤波器原型 由固有频率Wn把模拟低通滤波器原型转换为低通、高
Wn的高通滤波器；当ftype＝stop时，可设计出带阻滤波器
[z,p,k]=cheby1(n,Rp,Wn)
[zp,k]= cheby1 (n,Rp,Wn,’ftype’)
[A,B,C,D]= cheby1 (n,Rp,Wn)
[A,B,C,D]= cheby1 (n,Rp,Wn,’ftype’) 模拟域：[b,a]= cheby1 (n,Rp,Wn,’s’)可设计出截止频率为Wn的n阶
滤波器，
其余形式类似于数字域的。
2）chebyshevI型滤波器（通带等波纹）设计 数字域：[b,a]=cheby1(n,Rp,Wn)可设计出n阶chebyshevI滤波器,其截止频率由
Wn确定，通带内的波纹由Rp确定 [b,a]=cheby1(n,Rp,Wn,’ftype’)当ftype=high时，可设计出截止频率为
应π弧度。
5．直接设计IIR数字滤波器 1）Butterworth模拟和数字滤波器设计 数字域：[b,a]=butter(n,Wn)可设计出截止频率为Wn的n阶butterworth滤波器
[b,a]=butter(n,Wn,’ftype’)当ftype=high时，可设计出截止频率为Wn的 高通滤波器；当ftype＝stop时，可设计出带阻滤波器
实验原理及方法
为了克服冲激响应不变法产生的频率混叠现象，这是从S 平面到Z平面的标准变换z＝esT的多值对应关系导致的，为了 克服这一缺点，产生了双线性变换法。
j
j Im[z]
-1
1 Re[z]
0
0
s平面
z平面
实验原理及方法
双线性变换法的映射函数：
s
2 T
1 1
z 1 z字域的。
3）chebyshevII型滤波器（阻带等波纹）设计 数 字 域 ： [b,a]=cheby2(n,Rs,Wn) 可 设 计 出 n 阶 chebyshevI滤波器,其截止频率由Wn确定，阻带内的 波纹由Rs确定
[b,a]=cheby2(n,Rs,Wn,’ftype’) 当ftype=high时，可设计出截止频率为Wn的高通滤波 器；当ftype＝stop时，可设计出带阻滤波器
z
1
s
2 T
2
分
式
展
开
T sz T s z 1 0 22
固定其中一个变量，则另一个是线性的。或者说，这种变 换对于s和z是双线性的。
实验原理及方法
2 tan1
T 2
2 T
tan
2
这说明模拟角频率Ω和数字角频率ω的关系是非线性的，但是没出 现混叠。在把Ω变换为ω时产生了非线性畸变。为了克服它带来的问题 ，通常要使Ω按上式预修正，以抵消畸变的影响。