信号与系统实验八
- 格式:doc
- 大小:176.00 KB
- 文档页数:5
实验八 信号的无失真传输一、实验目的1.了解信号的无失真传输的基本原理;2.熟悉信号无失真传输系统的结构与特性。
二、实验设备1. 信号与系统实验(一)2.虚拟示波器三、实验内容1.设计一个无源(或有源)的无失真传输系统;2.令幅值固定、频率可变化的正弦信号作为系统的输入信号,测量系统输出信号的幅值和相位(用李沙育图形法)。
四、实验原理1.信号的无失真传输是指通过系统后输出信号的波形与输入信号的波形完全相同,只允许有幅值上的差异和产生一定的延迟时间,具有这种特性的系统称为无失真传输系统。
令输入信号为X (t ),则系统的输出为)t -Y (t )=k x (t 0 式中k,t 0为常量,对上式取付氏变换,则有0-j ωω)e Y(j ω)=kx(j(ω)j -j ωt -e |H |==ke X(j ω)Y(j ω)H(j ω)=0ϕ |H|=k k 为常数 ωt (ω)=0-ϕ t 0>02.实验电路系统图8-1无失真传输的电路图其中R 1=R 2=20k ,C 1=C 2=1uF它的频率特性为=K R R R =C j ωR 1R C j ωR 1R C j ωR 1R =(j ω)U (j ω)U H(j ω)=122222111222i o +++++ 五、实验步骤1.分析信号无失真传输系统的模拟电路,如图8-1所示。
2.在模拟电路的输入端输入一个正弦信号,并改变其频率,用示波器观察输出信号的幅值和相位。
六、实验报告1.画出信号无失真传输系统的模拟电路。
2.分析无失真传输系统的结构特点,如果R2R1≠、C2C1≠,则系统的)(ωϕ和H(j ω)会产生什么变化?七、实验思考题1.为什么输出信号波形与输入信号波形相同?2.信号的无失真传输系统与全通滤波器有何不同?。
信号与系统实验指导手册沈阳工业大学信息科学与工程学院2005年10月前言“信号与系统”是电子工程、通信工程、信息工程、微电子技术、自动化、计算机等电类相关专业的一门重要的专业基础课,为国内、外各高等院校相关专业的主要课程。
由于本课程的理论性、系统性较强,为使抽象的概念和理论形象化、具体化,使学生能够比较深入的理解《信号与系统》课程的基本理论和分析方法,并提高学生分析问题和解决问题的能力。
为此,开设了基于本课程的实验。
《信号与系统》实验指导手册,将《信号与系统》课程的理论知识与“信号与系统”的实验系统设备结合,从内容上对教材起到了一定的补充作用,为学生具体实验进行了指导。
鉴于时间仓促,可能会存在一些不足与错误之处,欢迎大家批评指正,使之完善。
编者2005年10月目录实验一系统的特性测试 (1)实验二信号的采样与恢复 (8)实验三模拟滤波器分析 (14)实验四模拟滤波器的设计 (26)实验一系统的特性测试一、实验目的1、学会利用运算单元,搭建一些简单的实验系统。
2、学会测试系统的频率响应的方法。
3、了解二阶系统的阶跃响应特性。
4、学会对其零状态响应和零输入响应进行分析。
二、实验内容1、根据要求搭建一阶、二阶实验系统。
2、测试一阶、二阶系统的频响特性和阶跃响应。
三、预备知识学习使用波特图测试系统频响的方法。
四、实验仪器1、信号与系统实验箱一台(主板)。
2、线性系统综合设计性模块一块。
3、20MHz双踪示波器一台。
五、实验原理1、基本运算单元(1)比例放大1)反相数乘器由:2211R U R U -= 则有:1122R U R U = 2)同相数乘器 由:54443R R U R U +=则有:()45434R R R U U += (2) 积分微分器1)积分器:由:21211//1R SC U R U -= 则有:()1212121C SR R R U U +-= 2)微分器 由:14131R USC U -= 则有:S C R U U 1134-= (3) 加法器1)反相加法器有:⎪⎪⎭⎫⎝⎛+-=211032R U R U R U2)同相加法器由:()⎪⎪⎩⎪⎪⎨⎧+=⎪⎪⎭⎫ ⎝⎛++=+-+87576434433111R R U R U U R R R R U R U 令643*////R R R R = 则有:()⎪⎪⎭⎫⎝⎛++=4433787*5R U R U R R R R U 2、N 阶系统()()()()()()()()t e E t e dt dE t e dtd E te dt d E t r C t r dt dC t r dtd C t r dt d C m m m m m m n n n n n n ++∧++=++∧++------1111011110根据零状态响应(起始状态为零) ,则对其进行拉氏变换有:()()()()()()()()S e E S Se E S e S E S e S E S r C S Sr C S r S C S r S C m m m m n n n n ++∧++=++∧++----11101110则其传递函数表达式为:()()()nn n n mm m m C S C S C S C E S E S E S E S e S r S H ++∧++++∧++==----111011103、作为一阶系统,一般表达式为:()1010C S C E S E S H ++=一阶系统是构成复杂系统的基本单元,学习一阶的特点有助于对一般系统特性的了解。
信号与系统实验教程信号与系统实验是电子信息类专业中一门重要的实验课程。
在这门实验中,学生将学习如何利用实验仪器和软件工具来分析和处理信号,并理解信号在系统中的作用和相互之间的关系。
以下是一些常见的信号与系统实验教程:1. 实验一:信号的采集与表示- 学习使用信号采集仪器(例如信号发生器、示波器等)。
- 了解采样原理和采样频率对信号的影响。
- 学习如何将模拟信号转换为数字信号。
- 使用编程语言或工具对信号进行采样和表示。
2. 实验二:信号的变换与处理- 学习傅里叶变换和信号频谱分析的原理。
- 使用傅里叶变换工具(例如FFT算法)对信号进行频谱分析。
- 学习信号的时域和频域表示之间的转换关系。
- 学习数字滤波器的原理和应用。
3. 实验三:线性时不变系统的特性分析- 学习线性时不变系统的定义和性质。
- 了解系统的单位冲激响应和冲激响应与输入信号的卷积关系。
- 利用实验仪器测量系统的冲激响应。
- 使用软件工具对系统进行时域和频域特性分析。
4. 实验四:信号采样与重构- 学习信号采样和重构的理论基础。
- 利用实验仪器对信号进行采样和重构。
- 学习采样定理的应用和限制。
- 学习插值和抽取技术对信号进行采样和重构。
5. 实验五:系统的频率响应与稳定性- 学习系统的频率响应和稳定性分析。
- 使用频率响应仪器(例如频谱分析仪)对系统进行测量和分析。
- 学习系统的振荡和稳定条件。
- 学习系统的幅频特性和相频特性之间的关系。
以上是信号与系统实验教程的一些基本内容,具体的实验内容和教程可以根据教学大纲和教材进行更详细的设计和安排。
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2.符号运算表示法若一个连续时间信号可用一个符号表达式来表示,则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现(在命令窗口中输入,每行结束按回车键)。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下:图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮,进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分:1) 两个轴组成的坐标平面(横轴是时间,纵轴是信号值);2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮,点击各波形选择按钮可选择波形,点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V,频率为0.5Hz,相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条,改变滑块的位置,滑块上方实时显示滑块位置代表的数值,对应正弦波的三个参数:幅度、频率、相位;坐标平面内实时地显示随参数变化后的波形。
在七种信号中,除抽样函数信号外,对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5~图1-10所示。
图1-5 峰值为8V,频率为1Hz,占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V,滞后时间为5秒的阶跃信号图1-9 峰值为8V,频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中,通过对滑动块的控制修改信号的幅度、频率、相位,观察波形的变化。
电气学科大类2012 级《信号与控制综合实验》课程实验报告(基本实验一:信号与系统基本实验)姓名丁玮学号U201216149 专业班号水电1204 同组者1 余冬晴学号U201216150 专业班号水电1204 同组者2 学号专业班号指导教师日期实验成绩评阅人实验评分表基本实验实验编号名称/内容实验分值评分实验一常用信号的观察实验二零输入响应、零状态相应及完全响应实验五无源滤波器与有源滤波器实验六LPF、HPF、BPF、BEF间的变换实验七信号的采样与恢复实验八调制与解调设计性实验实验名称/内容实验分值评分创新性实验实验名称/内容实验分值评分教师评价意见总分目录1.实验一常用信号的观察 (1)2.实验二零输入响应、零状态响应及完全响应 (4)3.实验五无源滤波器与有源滤波器 (7)4.实验六 LPF、HPF、BPF、BEF间的转换 (14)5.实验七信号的采样与恢复 (19)6.实验八调制与解调 (29)7.实验心得与自我评价 (33)8.参考文献 (34)实验一常用信号的观察一.任务与目标1.了解常见信号的波形和特点;2.了解常见信号有关参数的测量,学会观察常见信号组合函数的波形;3.学会使用函数发生器和示波器,了解所用仪器原理与所观察信号的关系;4.掌握基本的误差观察与分析方法。
二.总体方案设计1.实验原理描述信号的方法有许多种,可以用数学表达式(时间的函数),也可以使用函数图形(信号的波形)。
信号可以分为周期信号和非周期信号两种。
普通示波器可以观察周期信号,具有暂态拍摄功能的示波器可以观察到非周期信号的波形。
目前,常用的数字示波器可以方便地观察周期信号及非周期信号的波形。
2.总体设计⑴观察常用的正弦波、方波、三角波、锯齿波等信号及一些组合函数的波形,如y=sin(nx)+cos(mx)。
⑵用示波器测量信号,读取信号的幅值与频率。
三.方案实现与具体设计1.用函数发生器产生正弦波,并且设定波形的峰值及频率,用示波器观察并记录波形,测量和读取信号的幅值与频率;2.用函数发生器产生方波,并且设定波形的峰值及频率,用示波器观察并记录波形,测量和读取信号的幅值与频率;3.用函数发生器产生三角波,并且设定波形的峰值及频率,用示波器观察并记录波形,测量和读取信号的幅值与频率;4.用函数发生器产生锯齿波,并且设定波形的峰值及频率,用示波器观察并记录波形,测量和读取信号的幅值与频率;5.用函数发生器产生两个不同频率的正弦波,分别设定波形的峰值及频率,用示波器叠加波形,并观察组合函数的波形。
信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解,掌握傅里叶系数的计算方法;(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法;(3) 熟悉傅里叶变换的性质,并能应用其性质实现信号的幅度调制;(4) 理解连续时间系统的频域分析原理和方法,掌握连续系统的频率响应求解方法,并画出相应的幅频、相频响应曲线。
二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ,它的周期为T ,角频率22fT,且满足狄里赫利条件,则该周期信号可以展开成傅里叶级数,即可表示为一系列不同频率的正弦或复指数信号之和。
傅里叶级数有三角形式和指数形式两种。
1)三角形式的傅里叶级数:01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ,n b 称为傅里叶系数,可由下式求得:222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2)指数形式的傅里叶级数:()jn tn nf t F e式中系数n F 称为傅里叶复系数,可由下式求得:221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时,本质上是对信号进行数值积分运算。
Matlab中进行数值积分运算的函数有quad函数和int函数。
其中int函数主要用于符号运算,而quad函数(包括quad8,quadl)可以直接对信号进行积分运算。
因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算,也可以采用符号运算方法。
quadl函数(quad系)的调用形式为:y=quadl(‘func’,a,b)或y=quadl(@myfun,a,b)。
其中func是一个字符串,表示被积函数的.m文件名(函数名);a、b分别表示定积分的下限和上限。
信号与系统实验报告
实验名称:信号与系统实验
一、实验目的:
1.了解信号与系统的基本概念
2.掌握信号的时域和频域表示方法
3.熟悉常见信号的特性及其对系统的影响
二、实验内容:
1.利用函数发生器产生不同频率的正弦信号,并通过示波器观察其时域和频域表示。
2.通过软件工具绘制不同信号的时域和频域图像。
3.利用滤波器对正弦信号进行滤波操作,并通过示波器观察滤波前后信号的变化。
三、实验结果分析:
1.通过实验仪器观察正弦信号的时域表示,可以看出信号的振幅、频率和相位信息。
2.通过实验仪器观察正弦信号的频域表示,可以看出信号的频率成分和幅度。
3.利用软件工具绘制信号的时域和频域图像,可以更直观地分析信号的特性。
4.经过滤波器处理的信号,可以通过示波器观察到滤波前后的信号波形和频谱的差异。
四、实验总结:
通过本次实验,我对信号与系统的概念有了更深入的理解,掌
握了信号的时域和频域表示方法。
通过观察实验仪器和绘制图像,我能够分析信号的特性及其对系统的影响。
此外,通过滤波器的处理,我也了解了滤波对信号的影响。
通过实验,我对信号与系统的理论知识有了更加直观的了解和应用。
实验一常用信号分类与观察一、实验目的1、观察常用信号的波形,了解其特点及产生方法。
2、学会用示波器测量常用波形的基本参数,了解信号及信号的特性。
二、实验内容1、了解几种常用典型信号的解析式及时域波形。
2、观察这些信号的波形,思考可以从那几个角度观察分析这些信号的参数。
三、实验仪器1、信号与系统实验箱一台(主板)。
2、20MHz双踪示波器一台。
四、实验原理对于一个系统特性的研究,其中重要的一个方面是研究它的输入输出关系,即在一特定的输入信号下,系统对应的输出响应信号。
因而对信号的研究是对系统研究的出发点,是对系统特性观察的基本手段与方法。
在本实验中,将对常用信号和特性进行分析、研究。
信号可以表示为一个或多个变量的函数,在这里仅对一维信号进行研究,自变量为时间。
常用信号有:指数信号、正弦信号、指数衰减正弦信号、抽样信号、钟形信号、脉冲信号等。
1、指数信号:指数信号可表示为()atf t Ke。
对于不同的a取值,其波形表现为不同的形式,如图1-1所示:武汉凌特电子技术有限公司 信号与系统综合实验指导书第 2 页 共 128 页图1-1 指数信号2、指数衰减正弦信号:其表达式为(0)()sin()(0)att f t Ke t t ω-<⎧=⎨>⎩,其波形如图1-2所示:图1-2 指数衰减正弦信号3、抽样信号:其表达式为:sin ()a t S t t=。
()a S t 是一个偶函数,t =±π,±2π,…,±n π时,函数值为零。
该函数在很多应用场合具有独特的运用。
其信号如图1-3所示:武汉凌特电子技术有限公司 信号与系统综合实验指导书第 3 页 共 128 页图1-3 抽样信号4、钟形信号(高斯函数):其表达式为:()2t f t Eeτ⎛⎫- ⎪⎝⎭=,其信号如图1-4所示:图1-4钟形信号5、脉冲信号:其表达式为)()()(T t u t u t f --=,其中)(t u 为单位阶跃函数。
《信号与系统》实验报告湖南工业大学电气与信息工程学院实验一用同时分析法观测50Hz非正弦周期信号的分解与合成一、实验目的1、用同时分析法观测50Hz非正弦周期信号的频谱,并与傅立叶级数各项的频率与系数作比较。
2、观测基波和其谐波的合成。
二、实验设备1、信号与系统实验箱:TKSS -A型或TKSS -B 型TKSS -C 型;2、双踪示波器三、实验原理1、 一个非正弦周期函数可以用一系列频率成整数倍的正弦函数来表示,其中与非正弦具有相同频率的成分称为基波或一次谐波,其他成分则根据其频率为基波频率的2、3、4、…、n 等倍数分别称为二次、三次、四次、…、n 次谐波,其幅度将随着谐波次数的增加而减小,直至无穷小。
2、 不同频率的谐波可以合成一个非正弦周期波,反过来,一个非正弦周期波也可以分解为无限个不同频率的谐波成分,3、 一个非正弦周期函数可以用傅立叶级数来表示,级数各项系数之间的关系可用一个频谱来表示,不同的非正弦周期函数具有不同的频谱图,各种不同波形及其傅氏级数表达式见表2-1,方波频谱图如图2-1表示Um1351/91/51/71/3790ωωωωωω图1-1 方波频谱图表2-1 各种不同波形的傅立叶级数表达式UmtTU 2τ方波Um0TU 2τ正弦整流全波UmTU 2τ三角波Um0T2τ正弦整流半波t tUm0tT U 2τ矩形波U1、方波 ())7s i n 715s i n 513s i n 31(s i n 4 ++++=t t t t u t u mωωωωπ 2、三角波())5s i n 2513sin 91(sin 82++-=t t t u t u mωωωπ3、半波())4c o s 1512cos 31sin 421(2 +--+=t t t u t u m ωωωππ 4、全波 ())6c o s 3514cos 1512cos 3121(4 +---=t t t u t u m ωωωπ5、 矩形波())3cos 3sin 312cos 2sin 21cos (sin 2 ++++=t T t T t T U T U t u m m ωτπωτπωτππτ实验装置的结构如图1-2所示DC20f f f f f f 3456图1-2信号分解于合成实验装置结构框图图中LPF 为低通滤波器,可分解出非正弦周期函数的直流分量。
信号与系统实验报告目录1. 内容概要 (2)1.1 研究背景 (3)1.2 研究目的 (4)1.3 研究意义 (4)2. 实验原理 (5)2.1 信号与系统基本概念 (7)2.2 信号的分类与表示 (8)2.3 系统的分类与表示 (9)2.4 信号与系统的运算法则 (11)3. 实验内容及步骤 (12)3.1 实验一 (13)3.1.1 实验目的 (14)3.1.2 实验仪器和设备 (15)3.1.4 实验数据记录与分析 (16)3.2 实验二 (16)3.2.1 实验目的 (17)3.2.2 实验仪器和设备 (18)3.2.3 实验步骤 (19)3.2.4 实验数据记录与分析 (19)3.3 实验三 (20)3.3.1 实验目的 (21)3.3.2 实验仪器和设备 (22)3.3.3 实验步骤 (23)3.3.4 实验数据记录与分析 (24)3.4 实验四 (26)3.4.1 实验目的 (27)3.4.2 实验仪器和设备 (27)3.4.4 实验数据记录与分析 (29)4. 结果与讨论 (29)4.1 实验结果汇总 (31)4.2 结果分析与讨论 (32)4.3 结果与理论知识的对比与验证 (33)1. 内容概要本实验报告旨在总结和回顾在信号与系统课程中所进行的实验内容,通过实践操作加深对理论知识的理解和应用能力。
实验涵盖了信号分析、信号处理方法以及系统响应等多个方面。
实验一:信号的基本特性与运算。
学生掌握了信号的表示方法,包括连续时间信号和离散时间信号,以及信号的基本运算规则,如加法、减法、乘法和除法。
实验二:信号的时间域分析。
在本实验中,学生学习了信号的波形变换、信号的卷积以及信号的频谱分析等基本概念和方法,利用MATLAB工具进行了实际的信号处理。
实验三:系统的时域分析。
学生了解了线性时不变系统的动态响应特性,包括零状态响应、阶跃响应以及脉冲响应,并学会了利用MATLAB进行系统响应的计算和分析。
信号与系统实验
实验八:系统的复频域分析
小组成员:
黄涛13084220
胡焰焰13084219
洪燕东13084217
一、实验目的
1、掌握系统的复频域分析方法。
2、掌握测试系统的频率响应的方法。
二、预习内容
1、系统频响的方法。
(见第四章波特图的介绍)
三、实验原理
1. N 阶系统系统的传递函数
用微分方程描述的N 阶系统为:
)()()()()()()()(1111011110t e E t e dt d E t e dt d E t e dt d E t r C t r dt d C t r dt d C t r dt d C m m m m m m n n n n n n ++++=++++------ 根据零状态响应(起始状态为零),则对其进行拉氏变换有:
)()()()()()()()(11101110S e E S Se E S e S E S e S E S r C S Sr C S r S C S r S C n n m m n n n n ++++=++++---- 则系统传递函数可表达为:
n n n n n n m m C S C S C S C E S E S E S E S e S r S H ++++++++==----11101110)()()(
用差分方程描述的N 阶系统为:
0101[][1]...[][][1]...[]N M a y n a y n a y n N b x n b x n b x n M +-++-=+-++-
根据零状态响应(起始状态为零),则对其进行拉氏变换有:
110101()()...()()()...()N M N M a Y z a z Y z a z Y z b X z b z X z b z X z ----+++=+++
则系统传递函数可表达为:
101101...()()()...M M N
N b b z b z Y z H z X z a a z a z ----++==+++ 2.根据系统传递函数的零极点图分析系统
零点:传递函数分子多项式的根。
极点:传递函数分母多项式的根。
根据零极点图的不同分布分析系统。
3.涉及到的Matlab 函数
(1)freqz 函数:实验六中出现过,可用来求单位圆上的有理z 变换的值。
调用格式:同实验六
(2)zplane 函数:得到有理z 变换的零极点图。
调用格式:zplane(num,den)
其中,num 和den 是按1
z -的升幂排列的、z 变换分子分母多项式系数的行向量。
(3)roots 函数:求多项式的根。
调用格式:r=roots(c), c 为多项式系数向量;r 为根向量。
四、实验内容
1.系统零极点的求解 (1)求解系统2
321)(232+++-=s s s s s H 和321223211)(----+++-=z z z z z H 的零极点,验
证下面程序的运行结果,根据系统零极点图分析系统性质。
b=[1,0,-1]; a=[1,2,3,2]; zr=roots(b); pr=roots(a);
plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2);
grid; legend('零点','极点');
figure; zplane(b,a);
图7-1 系统零极点图 图7-2 由zplane 函数直接绘制系统零极点图
(2)参考上述程序,绘制系统231
)(2++-=s s s s H 和43214
321245535952)(--------++++++++=z z z z z z z z z H 的零极点图,并分析系统性质。
与用zplane 函数直接绘制系统零极点图(注:圆心的圆圈并非系统的零点)做比较。
clc;
b=[1,-1];
a=[1,3,2];
zr=roots(b);
pr=roots(a);
subplot(2,1,1);
plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2);
grid;
legend('零点','极点');
subplot(2,1,2);
zplane(b,a);
legend('零点','极点');
figure;
zplane(b,a);
b=[2,5,9,5,3];
a=[5,45,2,1,1];
zr=roots(b);
pr=roots(a);
subplot(2,1,1);
plot(real(zr),imag(zr),'go',real(pr),imag(pr),'mx','markersize',12,'linewidth',2);
grid;
%legend('零点','极点');
subplot(2,1,2);
zplane(b,a);
2. 求解z 变换
(1)对上题中的系统)(z H 在单位圆上求z 变换。
w=-4*pi:8*pi/511:4*pi;
num=[2 5 9 5 3];
den=[5 45 2 1 1];
h=freqz(num,den,w);
subplot(2,1,1);
plot(w/pi,abs(h));
grid;title('幅度谱');
xlabel('omega/\pi');ylabel('振幅');
subplot(2,1,2);
plot(w/pi,angle(h));
grid;title('相位谱');
xlabel('omega/\pi');
ylabel('以弧度为单位的相位');。