MATLAB信号处理

-实验一 MATLAB入门（1）1.实验目的： （1）了解MATLAB的体系结构与特点，熟悉其集成开发环境。 （2）熟悉MATLAB界面窗口的功能和使用方法。 （3）熟悉MATLAB的帮助系统及使用方法。 （4）了解MATLAB的的数据类型、基本形式和数组的产生方法。 （5）掌握MATLAB基本的数学运算操作。2.实验原理 （1）MATLAB简介MATLAB是美国MathWorks公司开发的高性能的科学与工程计算软件。它在数值计算、自动控制、信号处理、神经网络、优化计算、小波分析、图像处理等领域有着广泛的用途。近年来，MATLAB在国内高等院校、科研院所的应用逐渐普及，成为广大科研、工程技术人员必备的工具之一。MATLAB具有矩阵和数组运算方便、编程效率极高、易学易用、可扩充性强和移植性好等优点，俗称为“草稿纸式的科学计算语言”。它把工程技术人员从繁琐的程序代码编写工作中解放出来，可以快速地验证自己的模型和算法。 经过几十年的扩充和完善,MATLAB已经发展成为集科学计算、可视化和编程于一体的高性能的科学计算语言和软件开发环境，整套软件由MATLAB开发环境、MATLAB语言、MATLAB数学函数库、MATLAB图形处理系统和MATLAB应用程序接口（API）等五大部分组成。MATLAB的主要特点包括强大的计算能力（尤其是矩阵计算能力）、方便的绘图功能及仿真能力、极高的编程效率。另外，MATLAB还附带了大量的专用工具箱，用于解决各种特定领域的问题。通过学习软件的基本操作及其编程方法，体会和逐步掌握它在矩阵运算、信号处理等方面的功能及其具体应用。通过本课程实验的学习，要求学生初步掌握MATLAB的使用方法，初步掌握M文件的编写和运行方法，初步将MATLAB运用于数字信号处理中。循序渐进地培养学生运用所学知识分析和解决问题的能力。 （2）MATLAB的工作界面（Desktop）与操作MATLAB安装成功后，第一次启动时，主界面如下图（不同版本可能有差异）所示： 其中 是命令窗口（Command Window），是MATLAB的主窗口，默认位于MATLAB界面的右侧，用于输入命令、运行命令并显示运行结果。 是历史命令窗（Command History），位于MATLAB界面的左下侧，默认为前台显示。历史命令窗用于保存用户输入过的所有的命令，为用户下一次使用同一个命令提供方便。 是当前目录浏览器（Current Directory），位于MATLAB 界面的左上侧，默认为前台显示。该窗口用于显示当前目录和目录中的所有文件及文件的相关信息，如类型，大小等。 是工作空间窗（Workspace），默认位于MATLAB界面左上侧的 后台。可以通过点击左上方的“Workspace”按钮使之在前台展现。MATLAB提供了方便实用的功能键用于编辑、修改命令窗口中当前和以前输入的命令行。这些功能键如下表所示：命令窗口中常用的功能键功能键功能功能键功能重新调入上一命令行End光标移到行尾重新调入下一命令行Ctrl+Home光标移到命令窗顶部光标左移一个字符Ctrl+End光标移到命令窗底部光标右移一个字符Esc清除命令行Ctrl+光标左移一个字Delete删除光标处字符Ctrl+光标右移一个字Backspace删除光标处左边字符Home光标移到行首clc清除命令窗显示内容clear清除工作空间中的变量和函数（3）MATLAB的帮助系统MATLAB提供了强大而完善的帮助系统，包括命令行帮助、联机帮助和演示帮助。要学会使用MATLAB，必须充分利用其帮助系统，尤其是命令行帮助功能。命令行帮助可以通过help命令获得，其格式为：help 或help 目录名目录名函数名主题名数学符号（全英文）第一种格式在命令窗直接输入help，不带任何参数，此时，将显示MATLAB的分类目录和对目录的简要说明；第二种格式可以显示具体目录所包含的命令和函数，或者具体的命令、函数、符号或某个主题的详细信息。例如，在命令窗口键入：help sin，将会显示关于正弦函数sin的详细的帮助信息，通常包括函数的输入输出参数的含义、函数的调用格式以及函数使用实例等。（4）MATLAB的数据类型、基本形式及产生方法 i) 数据类型：MATLAB中的数据类型主要有：数值型（numerical），字符型（char），元胞型（cell），结构体型（structure）等，实验中涉及的主要是数值型和字符型。与C语言不同的是，在MATLAB中使用变量时，不需要事先定义其数据类型，MATLAB会自动判定。 ii) 数据基本形式：数据、变量在MATLAB中的基本形式是数组（array），其中一维数组又称为向量（或矢量，vector），有行向量和列向量之分；二维数组又称为矩阵（matrix），而标量则视为11矩阵。 iii）数组的产生和引用：产生数组的方法很多，下面介绍两种基本的方法。一维数组的产生：st:step:en,其中st和en分别为数组第一个和最后一个元素，step为相邻两个元素之间的差值。例如 x=1:2:5 产生数组x=1 3 5。矩阵的产生：矩阵元素放在一对方括号中，同一行元素用逗号或空格隔开，不同行用分号隔开。例如 A=1 2;3,4。iiii）数组的引用：一维数组用单下标引用，比如x(2);矩阵用双下标或单下标引用均可，比如 A(1,2)与A（3）指的是同一个元素。除此以外，还有其他引用方式。3.实验内容及其步骤本次实验需要完成以下内容：（1） 常用窗口的基本操作注意：1）所有输入在英文输入状态下进行；2）表达式后加上分号“;”，则计算结果不显示，这种方法可以加速程序运行速度；表达式后无分号，计算结果马上显示在命令窗口中。在命令窗输入一条算式并计算，将计算结果赋给变量x，执行以下操作： 查看变量x的值； 清除命令窗口的显示（清屏）； 调出刚才输入的算式，修改并重新计算； 清除变量x参考：x =2*4+3/4-3.2（2） 向量产生产生一个行向量，元素全为1；产生一个列向量，元素全为0；产生一个元素为偶数的向量；参考：a=ones(1,5); b=zeros(6,1); c=2:2:12;（3） 输入简单的矩阵并计算：要点：英文标点；矩阵元素放在“”之中，每一行的元素间用逗号或空格隔开，不同的行以分号“;”或回车来分割；，并计算A*B，A.*B，A/B,A./B，比较结果有何不同。参考：A = 1,2;3,4; B = 2,4;1,3;A*B; A.*B; A/B;A./B; (4) 复数产生和运算：已知，求的模和幅角，计算。参考：z1 = 3+4i; z2 = 4+4i; mg =abs(z1);ang=angle(z2); z=z1/z2;4. 实验用MATLAB函数介绍在实验过程中用到的MATLAB指令（函数）有： 函数类：ones,zeros,abs, angle等，具体调用格式参看“help”或者查阅相关书籍。 命令或运算符类：clc，clear，help，lookfor，*，.*，.，/，./ , : 。5.思考练习题 （1）MATLAB在处理数据运算过程中与其他高级语言有何不同。 （2）熟悉MATLAB的常用数学运算符号，比较几种乘法的运算符号有何不同。（3）用help 命令查看以下函数的功能和用法：plot，stem，figure，subplot，axis，xlabel，ylabel，title，text，legend，sum。（4）已知，计算C=A*B，求C的实部，虚部，模和相角。（real , imag , abs , angle）实验二 MATLAB入门（2）1.实验目的：（1）掌握基本的绘图方法。（2）掌握M文件的编写和运行方法。（3）熟悉MATLAB的搜索路径及添加搜索路径的方法。2.实验原理 （1）二维图形的绘制与标注MATLAB提供了强大的绘图功能，便于数据的可视化操作。其中，二维图形的绘制与标注是最基本也是最常用的绘图操作。下面介绍所涉及的函数及其用法。1） plot 函数plot用于绘制连续函数的波形，其功能非常强大，可以同时绘制多条曲线，可以指定曲线的颜色、线型和数据点的形状。调用格式如下：plot(x,y) 绘制x为横坐标，y为纵坐标的曲线，x、y可以是向量，也可以是矩阵；例1：以下两行命令运行后，将绘制sin函数在0之间的波形t=0:0.01:pi;plot(t,sin(t)plot(x,y,s) 绘制曲线，并指定线型、颜色、数据点形状,s 为字符串例2：以下两行命令运行后，将绘制sin函数在0之间的波形,颜色为绿色，线型为虚线。t=0:0.01:pi;plot(t,sin(t),g-);plot(x1,y1,x2,y2) 于同一坐标系内分别x1,y1和x2，y2画波形；例3：以下命令在同一坐标系内画sin和cos函数的波形并指定颜色和线型。t=0:0.01:pi;plot(t,sin(t),g-,t,cos(t), r-);plot(x1,y1,s1,x2,y2,s2) 于同一坐标系内分别x1,y1和x2，y2画波形，并指定线型、颜色等例4：以下命令在同一坐标系内画sin和cos函数的波形，指定颜色和线型。t=0:0.01:pi;plot(t,sin(t),t,cos(t);2） stem 函数stem 用于绘制离散序列的波形。一般调用格式为：stem(n,x) 画出以n、x为坐标的图形。例如 n=0:10;x = cos(pi*n/3);stem(n,x);3） figure 和subplot函数figure 用来产生新的图形窗口。函数subplot用来将图形窗口切分为若干子窗口并激活其中某个子窗口用于当前绘图操作，一般和plot、stem等画图函数成对使用。调用格式为figure 生成一个新的图形窗口。subplot(mnk)或subplot(m,n,k) 将图形窗口切分为mn的子图，并激活第k个用于绘图。例5 以下命令在同一图形窗口的三个子图中分别画图：t=0:0.01:pi;n=0:10;subplot(311);plot(t,sin(t);subplot(312);plot(t,cos(t);subplot(313);stem(n,cos(pi*n/3);4） 用于图形标注的函数 xlabel，ylabel，title 坐标轴与图形名称标注 调用格式为：xlabel(s),ylabel(s),title(s),s为字符变量或常数。例如：xlabel(t),ylabel(y),title(waveform of cos) axis 坐标范围控制，调用格式为： axis（xmin xmax ymin ymax）例如：axis(-2 3 -3 3) text 添加文本注释，格式为text(x,y,s) 在坐标(x,y)处添加文本注释s，s是字符串。例如：text(0.5,1,maxpoint); legend 给多条曲线标注，同一坐标系中有多条曲线时，用legend区分它们。 调用格式：legend(str1,str2,),str1,str2,均为字符串例如：legend(sin,cos);（2）M文件的编写与运行MATLAB有两种运行方式：命令行运行方式和M文件运行方式。前者是在命令窗的提示符“”后直接输入命令或计算表达式，回车后MATLAB即执行运算并可以显示运行结果，这种方式适合实现一些简单的功能，如简单计算和绘图等；后一种方式中，M文件是用MATLAB语言编写的程序文件，其扩展名为.m，在MATLAB的M文件编辑器中输入、编辑和调试，生成后在命令窗口中运行该程序，与运行MATLAB其他命令相似。 从形式上来讲，M文件可以分成M脚本（M script）文件和M函数（M function）两类。脚本文件就是一系列MATLAB语句、命令的罗列、组合，运行脚本文件时，MATLAB依次执行文件中的每一行语句，与操作系统中批处理文件的运行方式相似。M函数文件和脚本文件的主要区别在于格式的不同：M函数文件第一行为函数定义行，以关键字“function”（通常为蓝色字体）开头，而脚本文件无此定义行；函数文件定义行的一般格式如下：function y1,y2,y3, = myfun(x1,x2,x3,)即定义行由关键字function，输出变量y1，y2，y3，,函数名myfun和输入变量x1，x2，x3,四部分组成。其中关键字和函数名不能省略，输入变量和输出变量可以有多个，也可以一个都没有。没有输入输出变量时，中括号，等号和小括号均可以省略。以下是M脚本文件和M函数文件的示例： %This is an M script sample。clear;x = -pi:0.1:pi;y1 = sin(x);y2 = cos(x);plot(x,y1,x,y2);title(cosine and sine functions);xlabel(time);ylabel(Amplitude);legend(y = cos(x),y = sin(x);grid on;%This is an M function sample.function y,pos = findmax(a)% findMax find the max value of matrix a togeter with its position% y = findmax(a):find the max value of matrix a% y,pos = findmax(a):find the max value of matrix a,get its position % as well in posy,p = max(a(:);r,c =ind2sub(size(a),p);pos = r,c;1）M文件的编写与保存将上述两段代码在MATLAB自带的文本编辑器里输入，并保存，就分别得到脚本文件和函数文件。保存时要注意几点，一是文件名必须以字母开头，二是函数文件名应与函数名字相同，三是名字不能与MATLAB关键字相同，四是扩展名为“.m”。2）M文件的运行脚本文件的运行很简单，只要在命令窗口中输入文件名（不含扩展名）并回车即可；比如，某脚本文件名为myfile.m，那么命令窗中输入myfile，并回车即可。对函数文件，在命令窗中调用形式与该函数定义行的格式相同。比如，对上述函数文件findmax，可以在命令窗口中输入: a=1,2;3,4;% 调用前需要给输入变量赋具体值mx,ps=findmax(a)完成对findmax的调用，两个输出值分别赋给mx和ps两个变量。（3）MATLAB的搜索路径及设置M文件运行时，MATLAB按照事先设置好的路径搜索顺序来寻找M文件，如果要执行的文件不在MATLAB的搜索路径中，就无法执行。利用MATLAB主界面File菜单中的“set path”可以将需要的目录、文件夹添加到MATLAB的搜索路径中，也可以设置新的搜索顺序。或者利用函数addpath，path，和pathtool设置搜索路径。3.实验内容及其步骤（4） 基本绘图操作将一个图形窗口分为三个子窗口，在三个字窗口中分别画y=sin(t),y=cos(t)和y=cos(n/3) 的波形。前两个为连续波形，第三个为离散波形。对第二个图的横坐标、纵坐标进行标注。（5） M文件的编写、调试及运行M 文件编写及运行：分别用编写脚本文件和函数文件实现 1+2+3+50。提示：150放在数组x中，x的生成语句为：n = 50; x = 1:n;求和用sum函数，函数用法请利用help命令获得。4. 实验用MATLAB函数介绍在实验过程中用到的MATLAB命令（函数）有： plot, stem, subplot, xlabel, ylabel, title, text，legend，figure，axis，addpath，path，pathtool，sum，sin，cos，exp等，具体调用格式参看“help”或者查阅相关书籍。另外，在具体的实验过程中也可以根据实际需要自己定义函数。 5.思考题 （1）用help命令查看plot的用法，回答以下问题： 1）plot函数中两个输入X和Y代表什么？他们是标量、向量还是矩阵？2）用plot绘图时，可以指定哪些线型和颜色？（2）绘制振荡曲线及其包络线的图形，其中t的取值范围是，并进行相应的标注。（3）编写一个函数文件，实现以下功能：输入一个数组，对其元素进行从大到小的排列，将排序后的数组输出；求数组元素的最大值并求出它在原数组中的位置。实验三 连续信号的频域分析1.实验目的： （1）掌握周期信号分解和合成的方法。 （2）掌握非周期信号频谱分析方法。 2.实验原理 （1）周期信号的分析与合成 周期为T的周期信号的傅里叶级数有三角型与指数型两种形式，分别表示为： 和 式子中，各系数的计算公式和相互关系如下：，傅里叶级数表明周期信号可以分解为正弦信号或虚指数的线性组合。由三角型的系数可画出周期信号的单边幅度谱和相位谱，由指数型的系数可画出周期信号的双边谱，它们都是离散谱。上述系数可以用MATLAB函数quad或quadl计算。它们的调用格式为：y = quad(FUN,A,B), 和 y = quadl(FUN,A,B)其中，FUN是被积函数名或函数句柄；A和B分别是积分区间的下限和上限。比如，系数可以如下计算：，(T为信号周期)PS: quad:采用递推自适应的Simpson法来计算,在低精度的非光滑曲线计算中是最有效的；quad1:采用递推自适应的Lobatto法来计算,在高精度的光滑曲线计算中更为高效;例1 求周期为4，幅度为1、脉冲宽度为2的对称矩形脉冲信号的三角型傅里叶级数。参考：T = 4;A = 1;tao = 2;w0 = 2*pi/T;N = 6;f = (n) (2*quad (t)(A*rectpuls(t,tao).*cos(n*w0*t),-T/2,T/2)/T);%(n)：匿名函数，自定义matlab中的函数，表示随后跟随的是属于变量n的函数。Matlab7以上版本才支持。for k = 1:N an(k) = f(k);endfun = (t)(A*rectpuls(t,tao);a0 = quad (fun,-T/2,T/2)/T;an = a0,an;nw0=w0*(0:N);stem(nw0,abs(an);例2 画出例1中前6次谐波合成的信号波形，要求画2个周期。t = -T:0.01:T;N=6;n=(0:N).;ft = an*cos(w0*n*t);plot(t,ft);（2）非周期信号的频谱分析非周期的频谱是连续谱，MATLAB通常用数值方法求其频谱，一种方法是利用DFT变换，另一种方法是利用函数quad和quadl。调用格式参见上文。其中，傅里叶变换为，数值计算时，频率处的频谱值为。例3求门宽为2的门函数的傅里叶变换，画出其幅度谱。参考：w1 = -2*pi;w2 = -w1;t1 = -tao/2;t2 = -t1;N = 500;wk = linspace(w1,w2,N);%线性间隔向量生成函数F = zeros(1,N);Fw = (w)(quad(t)(rectpuls(t,tao).*exp(-j*w*t),t1,t2);for k = 1:N F(k) = Fw(wk(k);endplot(wk,abs(F);3.实验内容及其步骤本次实验需要完成以下内容：（6） 复习有关傅里叶级数和傅里叶变换的有关内容。（7） 求周期为5，幅度为2、脉冲宽度为2的对称矩形脉冲信号的三角型傅里叶级数，i. 画出其幅度谱（取前7次谐波）。2）将前7次谐波叠加，画出合成的波形（画至少2个周期）。 提示：参考例1，例2.（8） 求单边指数信号的傅里叶变换，画出其在区间-,之间的幅度谱。提示：阶跃信号可使用heaviside(t)函数4. 实验用MATLAB函数介绍在实验过程中用到的MATLAB指令（函数）有： rectpuls, cos,quad,quadl, plot,linspace,zeros,abs,exp.5.思考练习题 （1）分析、对比实验内容1中合成波形和原周期信号波形之间的不同，简要说明原因，增多谐波个数，合成波形有什么变化？（2）对实验过程中所涉及的问题进行分析，试编写和修改相应的程序，得出最终正确的结果和波形图，并对实验报告进行整理分析。6.实验报告要求：（1）明确实验目的以及实验的原理。 （2）通过实验总结利用MATLAB分析信号频谱的一般方法。 （3）完成思考题的内容，对实验结果及其波形图进行分析对比，总结主要结论。实验四 离散信号的时域分析1. 实验目的（1） 学习MATLAB在信号分析和系统分析中的应用；（2） 掌握MATLAB表示离散序列的方法、序列运算的MATLAB实现；2. 实验原理（1） 离散信号的表示：在MATLAB中，可以用向量（数组）来表示有限长序列，不过这样的向量并没有包含样本位置的信息，因此，需要用两个向量（数组）来表示一个有限长序列，其中一个向量表示序列的值，称为样本向量；另一个向量表示样本的位置（即序列的序号），称为位置向量。比如序列在MATLAB中需要用以下两个向量来表示：n = -2,-1,0,1,2,3; x = 2,1,-1,5,2; 其中 向量n表示序号，向量x表示序列的值。当不需要样本位置信息或者序列从n = 0 开始时，可以只用样本向量表示。另外由于内存的限制，MATLAB无法表示无限长序列。（2） 离散信号的基本运算 信号相加：序列相加是对应样本相加，如果两序列长度不等或者位置向量不同，则不能用算数运算符“+”直接实现相加，必须对位置向量和长度统一处理后再相加。以下M函数可以实现任意两序列的相加运算： function y, n = sigadd(f1,n1,f2,n2)% y n = sigadd(f1,n1,f2,n2),Add two sequences.% Inputs:% f1 - the first sequence% n1 - index vector of f1% f2 - the second sequence% n2 - index vector of f2% Outputs:% y - the output sequence% n - index vector of y n = min(n1(1),n2(1):max(n1(end),n2(end); % index vector of y(n)y1 = zeros(1,length(n); y2 = y1; % initializationy1(n = n1(1) & n = n2(1) & n = n2(end) = f2;y = y1 + y2; 序列的移位：序列移位后，样本向量没有变化，只是位置向量变了。任意序列的移位可以用以下的M函数实现： function y,n = sigshift(x,m,n0)% y,n = sigshift(x,m,n0),result of y = x(n-n0);% Inputs:% x - sequence to be shifted% m - the index vector of x% n0 - shift amount % Outputs:% y - the output sequence% n - the index vector of yy = x;n = m + n0; 序列的反折：序列反折后，样本向量倒置，位置向量倒置且变号。以下函数可以实现序列的反折运算：function y,n = sigfold(x,m)% y,n = sigfold(x,n),result of y(n) = x(-n);% Inputs:% x - sequence to be folded% m - the index vector of x% Outputs:% y - the output sequence% n - the index vector of y y = fliplr(x);n = -fliplr(m);3. 实验内容和步骤（1） 序列产生及表示a 在MATLAB中产生有限区间上的（n）或（n-n0），画其波形。b 产生两个随机序列：x1序号始于-2，结束于7，样本值为整数；x2：序号始于0，结束于15，样本值介于 -2与2之间；对产生的序列画其波形图。参考：a. %Generate the delta sequence function x,n = deltaN(n0,n1,n2)% x,n deltaN(n0,n1,n2),generate delta(n-n0),n1 = n = n2;draw % its stem graph if needed% Inputs:% n0 - sample position% n1 - beginning index of x% n2 - end index of x% Outputs:% x - the delta sequence% n - index vector of xn = n1:n2;x = n = n0;if nargout =，=，=L y = y1,zeros(1,N-L);else m=L-N; y2 = circshift(y1,m); y3 = y2(1:m),zeros(1,N-m); y = y1(1:N)+y3;end函数中调用的circshift代码如下：function y=circshift(x,m)% y=x(n-m)N = length(x)if abs(m)N m= rem(m,N);end n=0:N-1;y=x(mod(n-m,N)+1);提示：分别计算x1、x2的DFT，相乘，求逆(用ifft)，将结果与圆周卷积的序列比较。4.实验用MATLAB函数介绍 用到的函数有：freqz，linspce,fft,ifft,subplot,plot,abs,angle,xlabel,title,hold on,conv,length,rem,mod等。5. 思考题 （1）如果，如何比较精确的求出其DTFT？ （2） 利用圆周卷积的矩阵形式，重新写一个计算圆周卷积的函数。实验八 离散时间系统的频域分析1.实验目的 （1）掌握系统频率响应的计算方法。（2）验证线性相位FIR滤波器的特性。2.实验原理（1）系统频响的计算 用freqz计算，格式是：X,w = freqz(b,a,N);可以计算出0区间N个频率点上的频响，结果存在X中，而w是这N个点的频率值构成的向量，b和a分别是系统传输函数的分子系数和分母系数。（2）线性相位FIR滤波器的特性 线性相位FIR滤波器共有4种类型，不同类型具有不同的零点特性。分别是：1 型：z=1 和z=-1处有偶数个零点或者没有零点；2型：z=1处有偶数个零点或无零点，z=-1处有奇数个零点；3型：z=1和z=-1处有奇数个零点；4型：z=1处有奇数个零点，z=-1处有偶数个零点或无零点；3.实验内容及其步骤（1）计算一个滑动平均滤波器的频响，取M=5，即,画出幅度和相位响应：提示：用freqz，abs和angle函数部分参考代码：M=5;b =ones(1,5)/M;a=1;H,w=freqz(b,a,500);（2）利用向量b=1 -8.5 12.4 -21,构造4种类型的线性相位滤波器，分别画出各自的单位样值响应，幅度响应、相位响应和零极点图。提示：用fliplr,freqz,abs,angle,zplane 函数部分参考代码：b1 =b,37,fliplr(b); %Design type1 FIR Filtern1=length(b1);subplot(221);stem(n1,b1);H,w=freqz(b1,1,500);%求H(e.jw)zplane(b1,1);%绘制零极点图4.实验用MATLAB函数介绍用到的函数有：freqz,stem,abs,angle,subplot,xlabel,title,zplane5.思考题（1）设系统函数 用MATLAB验证a0时是什么滤波器。（2）理想低通滤波器的单位样值响应可以用下式近似：，编写程序，画出的冲激响应，和幅度响应。选取不同的N值，比较所得结果。实验9 IIR数字滤波器设计1.实验目的 （1）掌握模拟滤波器设计方法（2）掌握双线性变换法和脉冲响应不变法设计滤波器；（3）掌握低通、高通、带通和带阻滤波器的设计。(4) 掌握信号滤波过程的具体实现2实验原理（1）IIR数字滤波器的设计方法目前，设计IIR数字滤波器的通用方法是先设计相应的模拟滤波器，然后再通过双线性变换法或脉冲响应不变法变换到相应的数字滤波器。模拟滤波器从功能上分有低通、高通、带通及带阻四种，从逼近方法上分有巴特沃兹（Butterworth）滤波器、切比雪夫（Chebyshev）I型滤波器、切比雪夫II型滤波器、椭圆（Elliptic）滤波器等。(2) 模拟滤波器设计第一步：求阶数及3dB截止频率巴特沃兹型：N,Wn=buttord(Wp,Ws,Ap,As,s)切比雪夫I型：N,Wn=cheb1ord(Wp,Ws,Ap,As,s)切比雪夫II型：N,Wn=cheb2ord(Wp,Ws,Ap,As,s)椭圆型：N,Wn=ellipord(Wp,Ws,Ap,As,s)第二步：求系统函数Ha(s)巴特沃兹型：b,a=butter(N,Wn,ftype,s)切比雪夫I型：b,a=cheby1(N,Ap,Wn,ftype,s)切比雪夫II型：b,a=cheby2(N,As,Wn,ftype,s)椭圆型： b,a=