实验一 Matlab软件的使用..
Matlab实验第一次实验答案
实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找sqrt（开方）函数的使用方法；解：sqrtSquare rootSyntaxB = sqrt(X)DescriptionB = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results.RemarksSee sqrtm for the matrix square root.Examplessqrt((-2:2)')ans =0 + 1.4142i0 + 1.0000i1.00001.41422、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B解：A=[1 2;3 4 ];B=[5 5;7 8 ];A^2*B（2）矩阵除法已知A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/B解：A=[1 2 3;4 5 6;7 8 9 ];B=[1 0 0;0 2 0;0 0 3 ];A\B,A/B（3）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'解：A=[5+1i,2-1i,1;6*1i,4,9-1i ];A1=A.',A2=A'（4）使用冒号表达式选出指定元素已知：A=[1 2 3;4 5 6;7 8 9];求A 中第3列前2个元素；A 中所有列第2，3行的元素；方括号[]解：A=[1 2 3;4 5 6;7 8 9 ];B=A([1,2],[3]),C=A(2:end, : )用magic 函数生成一个4阶魔术矩阵，删除该矩阵的第四列3、多项式（1）求多项式42)(3--=x x x p 的根解：A=[1 0 -2 -4];B=roots(A)（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征多项式;解：A=[1.2 3 5 .9 ; 5 1.7 5 6 ;3 9 0 1 ;1 2 3 4];A=poly(A);A=poly2sym(A)把矩阵A作为未知数代入到多项式中；4、基本绘图命令（1）绘制余弦曲线y=cos(t)，t∈[0，2π]解：t=0:.1:2*pi;y=cos(t);plot(t,y),grid（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]t=0:.1:2*pi;y1=cos(t-0.25);y2=sin(t-0.25);plot(t,y1，t,y2)grid5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；t=0:.01:4*pi;y=10*sin(t);plot(t,y,'-.',t,y,'r')grid6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000时的最大n值；解法1：s=0;i=0;while(s<2000) i=i+1;s=s+i; ends=s-i,i=i-1解法2：s=0;for i=1:1000; s=s+i;if(s>2000) ,break;endends=s-i,i=i-1（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
实验1Matlab的基本使用方法
7
path(path,’d:\mydir’)。用这种方法设置的目录信息只在当前环境中有效，一旦 Matlab 重新启 动，以前的设置就无效。
等于）
1
2） 逻辑运算：&（与）、|（或）、~（非）（非零元素都表示为真） 关系函数和逻辑函数：
函数
功能
函数
功能
all(x)
检查 x 是否全为 1
any(x)
检查 x 是否有不为零的元素
find(x)
找出非零元素的位置标识
isempty(x) 检查 x 是否为空阵
isinf(x)
检查 x 是否为无穷大
Your variables are:
RXYZtxyz
>> whos
Name
Size
Bytes Class
R
31x31
7688 double array
X
31x31
7688 double array
Y
31x31
7688 double array
Z
31x31
7688 double array
t
1x501
(11) 帮助.
3、实验指导 Matlab 是美国 Mathworks 公司推出的一套高性能的数值分析和计算软件，它将矩阵运
算、数值分析、图形处理、编程技术结合在一起，为用户提供了一个强有力的科学及工程问 题分析计算和程序设计的工具。
Matlab 语言可以被认为是一种解释性语言，用户可以在 Matlab 的工作空间中输入一个 指令，也可以在编辑器中编写应用程序，应用程序执行时，Matlab 软件对其中的命令和函 数进行翻译，然后在 Matlab 环境中对它进行处理，最后返回结果。 （1）数值、变量和表达式
实验1_Matlab操作入门
实验一 Matlab操作入门实验目的：1．了解Matlab的工作环境。
掌握help命令的使用。
了解general函数库中有哪些常用的函数，掌握who, whos, clear, save, load 等命令的使用。
2．掌握矩阵的输入方法；掌握利用elmat函数库中一些函数，如linspace, logspace, zeros, ones, eye, rand, randn, magic等创建矩阵的方法。
3．了解elmat函数库中的预定义变量。
4．掌握elmat函数库中diag, blkdiag, triu, tril, fliplr, flipud, rot90等函数的用法。
5．掌握elmat函数库中ndims, size, numel, length,disp等函数的使用。
6．掌握datafun函数库中的一些基本统计函数max, min, mean, median, sum, prod 等的基本用法。
7．掌握matfun函数库中的一些基本矩阵操作函数的使用：inv(A)—求矩阵A的逆，det(A)—求矩阵A的行列式值， rank(A)-求矩阵的秩，trace(A)-求矩阵Ａ的迹（即主对角线元素之和），norm(A)-求向量或矩阵的范数，cond(A)-求矩阵Ａ的条件数，etc。
实验内容：一、Matlab环境1．启动Matlab，了解Matlab的主要窗口：命令窗口，当前目录窗口，工作空间窗口，历史命令窗口。
2．输入help, 了解你的计算机上安装的Matlab中，有哪些帮助主题（函数库、工具箱）。
3．输入 help general（通用命令函数库）, 了解Matlab中有哪些通用的命令。
4．输入 help syntax，阅读具体的内容，了解Matlab语句或命令的格式。
5．输入 help format, 了解format命令的使用格式。
6．输入help timefun（时间、日期函数库）, 了解有哪些时间和日期函数。
MATLAB实验报告
实验一、Matlab 软件的安装与基本使用一、实验目的和要求：学会MA TLAB 软件的安装方法，熟悉MA TLAB 的工作环境及其基本功能的设置，以及基本指令的使用，了解MA TLAB 数值数组及其矩阵操作的使用、复数的输入方法。
二、实验内容和步骤：1、完成MA TLAB 软件的安装。
2、安装完成后启动MA TLAB 软件，熟悉MA TLAB 的工作环境。
打开MA TLAB 7，默认打开的窗口包括： (1)命令窗口(Command Window)； (2)命令历史窗口(Command History)； (3)工作间管理窗口(Workspace)；(4)当前路径窗口(Current Directory)此外，还有捷径键、编译窗口、图形窗口和帮助窗口等其他种类的窗口，结合课本实例，熟悉每个窗口的功能和基本使用。
3、（1）求23)]47(212[÷-⨯+的运算结果，注意命令后有无分号的作用。
命令及结果如下：t=[12+2*(7-4)]/3^2t = 2（2）求20012x x v t at =++的值，其中t=5s 、010x m =、015/v m s =、29.81/sec a m =-。
命令及结果如下：syms t x0 v0 a; t=5,x0=10,v0=15,a=-9.18; x=x0+v0*t+1/2*a*t^2 t = 5 x0 = 10 v0 = 15 x = -29.7500(3) 设3x =, 4v =，用MA TLAB 命令求下列各式的值：1）22log()x v + 命令及结果如下：syms x v x=3,v=4; log(x^2+v^2)x = 3 ans =3.21892)2(2)x v - 命令及结果如下：syms x vx=3,v=4;f=(sqrt(x-3))/(x-2*v)^2 x = 3 f = 0 3) |sin 2|vx e命令及结果如下： syms x v x=3,v=4;t=abs(2*x)*exp(v) x = 3t = 327.58894、用三种不同的方法输入矩阵1234246836912a ⎛⎫ ⎪= ⎪ ⎪⎝⎭，运行并观察实验结果，矩阵的三种基本输入方法如下：a=[1,2,3,4;3,4,6,8;3,6,9,12] a =1 2 3 4 3 4 6 8 3 6 9 12>> a=[1,2,3,43,4,6,8 3,6,9,12]a =1 2 3 4 3 4 6 8 3 6 9 12>> >>利用M 文件创建 1 2 3 4 3 4 6 8 3 6 9125、输入复数矩阵⎥⎦⎤⎢⎣⎡-+++=i ii i24356223B ，并求其实部、虚部、模和幅角（分别用弧度和角度表示），了解MA TLAB 向量化运算的特点。
信息论编码实验一MATLAB的基本操作
实验一 MATLAB 的基本操作一、实验目的1、掌握Matlab 软件使用的基本方法；2、熟悉Matlab 的数据表示、基本运算方法；3、熟悉Matlab 绘图命令及基本绘图控制。
二、实验仪器与软件1、PC 机 1台2、MATLAB7.0环境三、实验原理MATLAB 环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。
MATLAB 有3种窗口，即：命令窗口（The Command Window ）、m-文件编辑窗口（The Edit Window ）和图形窗口（The Figure Window ），而Simulink 另外又有Simulink 模型编辑窗口。
1、命令窗口（The Command Window ）当MATLAB 启动后，出现的最大的窗口就是命令窗口。
用户可以在提示符“>>”后面输入交互的命令，这些命令就立即被执行。
在MATLAB 中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。
在命令窗口中输入该文件名，这一连串命令就被执行了。
因为这样的文件都是以“.m ”为后缀，所以称为m-文件。
2、m-文件编辑窗口（The Edit Window ）我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m-文件。
在MATLAB 主界面上选择菜单“File/New/M-file ”就打开了一个新的m-文件编辑窗口；选择菜单“File/Open ”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。
如，计算3[(12)34)]2+⨯-÷，只需在提示符“>>”后输入“((1+2)*3-4)/2^3”，然后按Enter 键。
该命令行涉及加、减、乘、除及幂运算符，MATLAB 运算的执行次序遵循的优先规则为：从左到右执行；幂运算具有最高的优先级，乘法和除法具有相同的次优先级，加法和减法有相同的最低优先级；使用括号可以改变前述优先次序，并由最内层括号向外执行。
MATLAB实验代码
实验一MATLAB软件使用初步一、实验目的∙学习MATLAB软件的安装过程，熟悉MATLAB软件界面的组成及基本使用方法。
∙理解数组（Array）的分类，及标量（scalar）、矢量（vector）和矩阵（Matrix）的区别，熟悉数组与矩阵的构造方法，掌握数组与矩阵的基本运算法则。
二、实验要求1、掌握MATLAB软件的启动与退出的方法。
2、掌握MATLAB软件的桌面环境。
3、掌握MATLAB软件的基本使用方法。
4、掌握M文件编辑器的使用。
5、学会使用MATLAB的帮助系统。
6、掌握标量、矢量和矩阵的构造方法。
7、掌握数组的四则运算。
8、掌握矩阵的四则运算。
9、掌握基本的数据输入和输出函数。
三、实验内容1．通过桌面快捷图标或“开始—程序”菜单运行MATLAB，熟悉MATLAB的桌面环境，如图1所示。
图1 MATLAB桌面环境2．浏览各个菜单项的内容，试着改变各个窗口的字体与字号。
3．熟悉matlab 中的特殊变量：i 、j 、pi 、ans 等等。
4．在实验报告中回答下面哪些赋值语句是正确的，如果不正确给出理由。
a=1正确a20=1.2;正确 年薪=100000错误不能是汉字_fee=20.0错误开头不能下划线coordinate.x=1正确5．已知，，用MATLAB 分别执行下列语句。
并在实验报告中记录语句和1u =3v =结果。
a. b. c. d. 43u v ()222v u v -+333v v u -243v πu=1,v=3 (4*u)/(3*v) (2*v^-2)/(u+v)^2 (v^3)/(u^3-v^3)(4**pi*v^2)/3ans=0.4444 0.2222 -1.0385 37.69916. 运行以下命令并将运行结果与工作空间浏览器中的状态进行比较。
>> a=sin(pi); b=[1, 2]; c=[1; 2];>> who 所有可用变量 >> whos 显示变量所占的矩阵大小和所占的字节在实验报告中回答Size 下的1x1、1x2、2x1的含义>> whos Name Size Bytes Classa 1x1 8 double arrayb 1x2 16 double arrayc 16 double array矩阵大小7．通过工作目录浏览器（Current Directory Browser ）或者“File”菜单中相应菜单项改变当前工作目录。
matlab实验
实验一 MATLAB基本操作一、实验目的1、了解MATLAB应用程序环境2、掌握MATLAB语言程序的书写格式和MATLAB语言程序的结构。
3、掌握在MATLAB应用环境下编写程序4、掌握MATALB调试过程，帮助文件5、掌握MATLAB语言上机步骤，了解运行一个MATLAB程序的方法。
6、本实验可在学习完教材第一章后进行。
二、主要仪器及耗材PC电脑，MATLAB6.5软件三、实验内容和步骤1、MATLAB语言上机步骤：（1）、进入系统在C盘或其他盘上找到MATLAB或MATLAB6.5，然后双击其图标打开文件夹。
然后进行编辑源程序->编译->连接->执行程序->显示结果（2）、常用命令编辑切换（F6），编译（F9），运行（CTRL+F9），显示结果（ALT+F5）其它常用命令见“附录一”。
2、有下面的MATLAB程序。
（1）数值计算功能：如，求方程 3x4+7x3 +9x2-23=0的全部根p=[3,7,9,0,-23]; %建立多项式系数向量x=roots(p) %求根（2）绘图功能：如，绘制正弦曲线和余弦曲线x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));（3）仿真功能：如，请调试上述程序。
3、熟悉MATLAB环境下的编辑命令，具体见附录一。
三、实验步骤１、静态地检查上述程序，改正程序中的错误。
２、在编辑状态下照原样键入上述程序。
３、编译并运行上述程序，记下所给出的出错信息。
４、按照事先静态检查后所改正的情况，进行纠错。
５、再编译执行纠错后的程序。
如还有错误，再编辑改正，直到不出现语法错误为止。
四、实验注意事项１、记下在调试过程中所发现的错误、系统给出的出错信息和对策。
分析讨论对策成功或失败的原因。
2、总结MATLAB程序的结构和书写规则。
五、思考题1、matlab到底有多少功能？2、MATLAB的搜索路径3、掌握使用MATLAB帮助文件实验二 MATLAB 矩阵及其运算一、 实验目的1、了解矩阵的操作，包括矩阵的建立、矩阵的拆分、矩阵分析等2、了解MATLAB 运算，包括算术运算、关系运算、逻辑运算等3、掌握字符串的操作，了解结构数据和单元数据。
实验1 Matlab软件环境的基本使用
实验1 Matlab 的基本运算（1）实验目的：1、 了解Matlab 的发展和主要功能；2、 熟悉Matlab 工作环境的各个窗口；3、 掌握建立矩阵的方法；4、 掌握Matlab 各种表达式的书写规则以及常用函数的使用。
实验内容：1. 在命令窗口依次输入下列命令，根据执行结果分析其功能：helplookfor invhelp invwhich inv2. 建立自己的工作目录，再将自己的工作目录设置到Matlab 搜索路径下。
用help 命令查询自己的工作目录；3. 利用Matlab 的帮助功能，分别查询inv 、plot 、max 、round 等函数的功能及用法；4. 先求下列表达式的值，然后显示Matlab 工作空间的使用情况并保存全部变量。
（1）sin z e =+122851 （2）ln(z x =+212，其中.i x +⎡⎤=⎢⎥-⎣⎦2120455 （3）...sin(.)ln ,.,.,.,...,.,.,.a a e e a z a a --+=++=---03033030330292828293022提示：利用冒号表达式生成a 向量，求各点的函数值时用点乘运算。
（4）,,,:.:.,t t z t t t t t t ⎧≤<⎪=-≤<=⎨⎪-+≤<⎩224201112005252123其中提示：用逻辑表达式求分段函数值。
5. 已知：A -⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦12344347873657，B -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦131203327 求下列表达式的值。
（1）*A B A B I +-+6和，其中I 为单位矩阵；（2）*A B 和.A B *（3）^A 3和.^A 3（4）/A B 和\A B（5）[],A B 和()^[,],:;A B ⎡⎤⎣⎦1326. 将课件中的相关例题进行演示实验。
思考题：1. 与其他高级语言相比，Matlab 有哪些显著的特点？2. Matlab 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立的窗口？又怎样恢复？ 实验报告提交：将所有实验内容进行总结（实验内容6不必写），要求表述清晰，结果明确。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一 Matlab软件的使用1、实验目的和要求㈠掌握Matlab软件的使用方法；㈡常用离散时间信号的产生、显示和运算。
2、实验内容㈠Matlab软件常用命令和工具箱的使用，学会简单的编程；㈡编程实现常用离散时间信号；㈢离散时间信号的叠加、移位、线性卷积等基本运算。
3、实验原理利用软件生成数字信号处理系统中所涉及的信号及信号的加、移位、卷积等运算。
4、主要仪器设备计算机，Matlab软件。
实验步骤：程序1：n=-5:30;x=[zeros(1,5),1,zeros(1,30)];stem(n,x,'fill');grid on程序2：n=-5:30;x=[zeros(1,5),ones(1,8),zeros(1,23)];stem(n,x,'fill');grid on；程序3：n=-15:20;x=[zeros(1,15),1,2,3,4,5,zeros(1,16)];stem(n,x,'fill');grid on;程序4：N=10;M=10;L=N+M-1;x=[1,2,3,4,5,zeros(1,5)];h=[1,2,1,2,zeros(1,6)];y=conv(x,h);n=0:L-1;stem(n,y,'*k');grid on ;实验结果：1、2、3、4、实验二 DTFT和Z变换的Matlab实现1、实验目的和要求㈠掌握DTFT和Z变换的Matlab实现的方法；㈡进一步了解离散时域系统的时频特性。
2、实验内容㈠ DTFT和z变换的计算和基本性质；㈡离散时域系统的仿真；㈢利用函数实现离散时域系统的频率响应，分析零、极点的分布，进一步了解离散时域系统的时域特性。
3、实验原理利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统频率响应等的图像并于笔算结果进行比较，找出异同。
编译合适程序能计算取值范围不同的离散卷积。
4、主要仪器设备计算机，Matlab软件。
实验步骤：程序1：n=0:30;x=zeros(1,length(n));h=zeros(1,length(n));x([find((n>=0)&(n<=4))])=1;h([find((n>=0)&(n<=8))])=0.5; figure(1)subplot(3,1,1); stem(n,x);axis([0,30,0,2]); title('输入序列'); xlabel('n'); ylabel('x(n)'); subplot(3,1,2); stem(n,h);axis([0,30,0,2]);title('冲激响应序列'); xlabel('n');ylabel('h(n)');y=conv(x,h);subplot(3,1,3);n=0:length(y)-1; stem(n,y);title('输出响应'); xlabel('n');ylabel('y(n)');程序2：b=[0.0181,0.0543,0.0543,0.0181];a=[1.000,-1.76,1.1829,-0.2781];w=pi*freqspace(500); H=freqz(b,a,w);MH=abs(H);AH=angle(H);subplot(2,1,1);plot(w/pi,MH);grid;axis([0,1,0,1]);xlabel('w(pi)');ylabel('|H|');title('幅度、相位响应'); subplot(2,1,2);plot(w/pi,AH);grid;xlabel('w(pi)'); ylabel('angle(H)');实验结果：1、2、实验三 DFS和DFT的Matlab实现1、实验目的和要求㈠加深DFS和DFT的算法原理和基本性质的理解；㈡熟悉利用FFT对典型数字信号进行频谱分析。
2、实验内容㈠使用Matlab实现程序周期序列的DFS；㈡利用FFT和IFFT编程实现两个有限长序列的循环卷积；3、实验原理对有限长序列使用离散Fouier变换(DFT)可以很好的反映序列的频谱特性，而且易于用快速算法在计算机上实现。
4、主要仪器设备计算机，Matlab软件。
实验步骤：程序1：function shiyan22()F=50;N=64;T=0.000625;n=1:N;x=cos(2*pi*F*n*T);subplot(2,1,1);plot(n,x);hold ontitle('x(n)');xlabel('n');X=fft(x);subplot(2,1,2);plot(n,X);gridtitle('FFT|X|');xlabel('f(pi)');程序2：n=0:30;%输入x(n)和冲激响应h(n)x=zeros(1,length(n));h=zeros(1,length(n));x([find((n>=0)&(n<=4))])=1;h([find((n>=0)&(n<=8))]) =0.5;subplot(3,1,1);stem(x);title('x(n)');axis([0,30,0,2]);subplot(3,1,2);stem(h);title('h(n)');axis([0,30,0,2]);X=fft(x);H=fft(h);Y=X.*H;y=ifft(Y);subplot(3,1,3);实验结果：1、2、实验四 IIR数字滤波器的设计和实现1、实验目的和要求㈠掌握冲击响应不变法和双线形变换法设计IIR数字滤波器的原理和方法；㈡通过观察对实际信号的滤波，获得对数字滤波的感性认识。
2、实验内容㈠巴特沃思滤波器设计实现模拟低通滤波器；㈡利用冲击响应不变法和双线形变换法对模拟滤波器进行A/D转换；3、实验原理用数字滤波器的单位脉冲响应序列)(nh模仿模拟滤波器的冲激响应)(tha,让)(nh正好等于)(tha的采样值，即)()(nThnha=。
4、主要仪器设备计算机，Matlab软件。
实验步骤：设计一个巴特沃思数字低通滤波器，设计指标如下：通带内πωω2.0=<p幅度衰减不大于1dB；阻带πωω35.0=>s幅度衰减不小于15dB。
编制实现该数字滤波器程序并且实现数字滤波。
function shiyan3()fp=200*pi;fs=300*pi;wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1;Fs=1000;[N,fc]=buttord(fp,fs,Rp,As,'s'); [b,a]=butter(N,fc,'s');w=[0:1000*2*pi];[hf,w]=freqs(b,a,1000);[d,c]=impinvar(b,a,Fs);wd=[0:512]*pi/512;hw1=freqz(d,c,wd);[f,e]=bilinear(b,a,Fs);hw2=freqz(f,e,wd);OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);ep=sqrt(10^(Rp/10)-1);Ripple=sqrt(1/(1+ep.^2));Attn=1/10^(As/20);subplot(3,1,1);plot(wd/pi,abs(hw1)/abs(hw1(1)));title('幅度响应(冲击响应不变法)');xlabel('w(pi)');ylabel('H');ylabel('H');axis([0,1,0,1.1]);set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[0,Attn,Ripple,1]); grid;subplot(3,1,2);plot(wd/pi,20*log10((abs(hw2)+eps)/abs(hw2(1))));title('幅度响应(双线性变换法(dB))');xlabel('w(pi)');ylabel('H');axis([0,1,-40,5]);set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[-50,-15,-1,0]); grid;subplot(3,1,3);plot(wd/pi,-angle(hw2));title('相位响应');xlabel('w(pi)');ylabel('pi unit');%axis([0,1,0,1.1]);set(gca,'XTickmode','manual','XTick',[0,0.2,0.35,1.1]); set(gca,'YTickmode','manual','YTick',[-1,0,1]);grid;实验结果：实验五 FIR数字滤波器的设计和实现1、实验目的和要求㈠掌握用窗函数法设计FIR数字滤波器的原理和方法；㈡熟悉线性相位FIR滤波器的幅频特性和相频特性；㈢了解不同窗函数对滤波器性能的影响。
2、实验内容㈠用Hanning窗设计线性相位带通滤波器；㈡分别改用矩形窗和Blackman窗设计线性带通滤波器；3、实验原理确定数字滤波器的性能要求，根据性能要求，合理选择单位脉冲响应)(nh的奇偶对称性，从而确定理想频率响应)(ωjdeH的幅频特性和相频特性；求理想单位脉冲响应)(nhd，选择适当的窗函数)(nw，根据)()()(nwnhnhd=求所需设计的FIR滤波器单位脉冲响应；4、主要仪器设备计算机，Matlab软件。