下载文档原格式
实验七 连续信号与系统复频域分析的MATLAB 实现一、实验目的1. 掌握连续时间信号拉普拉斯变换的MATLAB 实现方法；2. 掌握连续系统复频域分析的MATLAB 实现方法。
二、实验原理1. 连续时间信号的拉普拉斯变换连续时间信号的拉普拉斯正变换和逆变换分别为：⎰∞∞--=dt e t f s F st )()(⎰∞+∞-=j j stds e s F j t f σσπ)(21)(Matlab 的符号数学工具箱（Symbolic Math Toolbox ）提供了能直接求解拉普拉斯变换和逆变换的符号运算函数laplace()和ilaplace ()。
下面举例说明两函数的调用方法。
（1）拉普拉斯变换例1.求以下函数的拉普拉斯变换。
)()()2()()()1(221t te t f t e t f t t εε--==解：输入如下M 文件：syms tf1=sym('exp(-2*t)*Heaviside(t)'); F1=laplace(f1) %求f1(t)的拉普拉斯变换 f2=sym('t*exp(-t)*Heaviside(t)'); F2=laplace(f2) 运行后，可得如下结果：F1 = 1/(s+2) F2 = 1/(s+1)^2 （2）拉普拉斯逆变换例2.若系统的系统函数为1]Re[,231)(2->++=s s s s H 。
求冲激响应)(t h 。
解：输入如下M 文件：H=sym('1/(s^2+3*s+2)');h=ilaplace(H) %求拉普拉斯逆变换运行后，可得如下结果：h=exp(-t)-exp(-2*t) 2. 连续系统的复频域分析 若描述系统的微分方程为∑∑===Mj j j Ni i i t f b t ya 0)(0)()()(则系统函数为)()()()()(00s A s B sa sb s F s Y s H Ni ii Mj jj===∑∑== 其中，∑∑====Mj j j Ni i i s b s B s a s A 0)(,)(。
《信号与系统》课程实验报告《信号与系统》课程实验报告一图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 抽样函数信号仿真途中，通过对滑动块的控制修改信号的幅度、频率、相位，观察波形的变化。
实验7 频带传输系统综合实验一、实验目的：1、将所做过的独立实验进行组合，构成可通话的频带传输系统。
掌握系统工作原理，了解信号流程并完成两人通话功能。
2、分析系统连接后所遇到的问题，且提出解决向题的方料3、了解载波相位误差，位同步相位误差对系统的性能的影响。
4、了解相位抖动对系统的影响。
二、实验内容：1、了解二人通话的2DPSK频带传输系统的原理，掌握实验模块的正确组合和连线。
2、了解决收端载波相位误差，位同步相位误差对系统造成影响的方法。
3、观察相位抖动对系统的影响。
4、分析系统组成后容易出现什么问题，提出解决问题的方法。
5、当发端发送固定数字信号时，将收、发端关键波形按相位关系绘画出来。
且加以说明。
三、预习要求1、复习前面各章有关的独立实验内容。
2、了解可通话2DPSK频带传输实验的原理。
3、自行设计实验方案及测试步骤。
四、实验仪器1、直流电源一台2、示波器一台3、数字信源模块、数字调制模块、载波时钟提取模块、数字解调模块、帧同步/终端模块、PCM编译码模块各一块。
五、实验原理图3-16-1 系统流程图1. 数字信源中的多路信号的复接原理数字信源模块的原理框图如附图所示。
1.1时序信号的产生本模块通过二进制分频器，得到１６kＨＺ和８kＨＺ方波信号，然后送入２／４译码器，得到反相的脉宽为八个时钟周期的四个脉冲信号，经反相器后得到正相的时序脉冲信号。
其波形及相位关系如图1所示。
图1 复接器中的时序脉冲信号1.2四路数据码的复接本信源模块中的四路独立的八位数码，在以上四路时序信号的控制下，依次选通模拟开关１、２、３、４，按顺序依次将四路数码接入同一通道，形成了一路串行码，完成四路数据码的复接。
在本基带传输实验中，将其中第一路数据码置成帧同步码（巴克码）Ｘ１１１００１０。
2. 数字终端模块的分接器原理原理框图如附图所示。
2.1时序脉冲产生电路由Ｕ７、Ｕ８、Ｕ９（７４ＬＳ１６４）八位移存器和Ｕ１２（７４ＬＳ０４）非门、Ｕ１０（７４ＬＳ７４）Ｄ触发器组成。
电气学科大类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、验证抽样定理。
二、实验设备1、信号与系统实验箱2、双踪示波器三、原理说明1、离散时间信号可以从离散信号源获得，也可以从连续时间信号抽样而得。
抽样信号f s(t)可以看成连续f(t)和一组开关函数s (t)的乘积。
s (t)是一组周期性窄脉冲，见实验图5-1，T s(t)称为抽样周期，其倒数f s(t)= 1/T s称为抽样频率。
图5-1 矩形抽样脉冲对抽样信号进行傅立叶分析可知，抽样信号的频率包括了原连续信号以及无限个经过平移的信号频率。
平移的频率等于抽样频率f s(t)及其谐波频率2f s、3f s》》》》》》。
当抽样信号是周期性窄脉冲时，平移后的频率幅度(sinx)/x规律衰减。
抽样信号的频谱是原信号频谱周期的延拓，它占有的频带要比原信号频谱宽得多。
2、正如测得了足够的实验数据以后，我们可以在坐标纸上把一系列数据点连起来，得到一条光滑的曲线一样，抽样信号在一定条件下也可以恢复到原信号。
只要用一截止频率等于原信号频谱中最高频率f n的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器输出可以得到恢复后的原信号。
3、但原信号得以恢复的条件是f s 2，其中f s为抽样频率，为原信号占有的频带宽度。
而f min=2 为最低抽样频率又称“柰奎斯特抽样率”。
当f s<2 时，抽样信号的频谱会发生混迭，从发生混迭后的频谱中我们无法用低通滤波器获得原信号频谱的全部内容。
在实际使用中，仅包含有限频率的信号是及少的，因此即使f s=2 ，恢复后的信号失真还是难免的。
图5-2画出了当抽样频率f s>2 （不混叠时）f s<2 （混叠时）两种情况下冲激抽样信号的频谱。
t f(t)0F()t 0m ωm ω-(a)连续信号的频谱Ts t 0f s (t)F()t0m ωm ω-s ω-s ω（）(b)高抽样频率时的抽样信号及频谱 不混叠图5-2 冲激抽样信号的频谱实验中f s >2 、f s =2 、f s <2 三种抽样频率对连续信号进行抽样，以验证抽样定理——要使信号采样后能不失真地还原，抽样频率f s 必须大于信号频率中最高频率的两倍。
信号与系统实验报告好啦，今天咱们来聊聊信号与系统实验报告。
这话题有点儿“高大上”，但咱们不妨来点轻松的，把它聊得有趣一些。
先说说信号是什么。
信号其实就是一种信息传递的方式，可能是声音，可能是光，甚至是你手机屏幕上刷过的每一条消息。
简单来说，信号就是承载着信息的载体。
你看，像咱们日常生活中，电台广播，手机接收到的短信，甚至你家电视里放的广告，它们都是信号的一种表现形式。
啊，听起来有点儿复杂吧？其实不难，就像你一收到朋友发来的微信，手机屏幕上跳出来的就是一个信号。
信号怎么才能“正常工作”呢？这就得说到“系统”了。
系统呢，说白了就是一套能够处理信号的工具。
你想啊，信号如果没有一个合适的“平台”去接收、传递和处理，那就变得一团乱麻了。
就像是你给朋友发了个短信，但他手机坏了，信号接收不进去，结果信息就白发了。
系统在这里就相当于是一个“修理工”，它能让信号顺利通过、准确无误地到达目的地。
接下来说说我们在实验中的“主角”——信号与系统。
你看，实验嘛，往往让我们有点“心慌慌”。
不过，信号与系统的实验其实有点像玩拼图。
你得先弄清楚信号的各种“形状”，然后用系统去“加工处理”，让它变得符合要求。
比如，咱们常用的模拟信号，它是一个连续的过程，类似于咱们生活中的声音一样，是没有间断的。
而数字信号呢，就像你手机屏幕上的数字，离散的，断断续续的。
每种信号都有自己独特的“脾气”，你得了解它们的特点，才能搭配合适的系统。
你要是觉得这些实验有点儿复杂，那就来点儿幽默的比喻吧。
信号就像是你的朋友说的话，而系统就是你听的耳朵。
朋友说话的声音，可能因为距离远近，语速快慢，甚至音量的大小而有所不同。
系统就得根据这些变化去处理，比如调节音量、清晰度，甚至过滤掉不必要的噪声。
你想想，假如你能在嘈杂的环境下清楚地听到朋友的声音，那就是系统给你提供的帮助。
信号与系统的实验，就是在这种“听”和“说”之间找到平衡点。
咱们得说说实验中的一些基本工具了。
信号与系统实验报告一、信号的时域基本运算1．连续时间信号的时域基本运算两实验之一实验分析：输出信号值就等于两输入信号相加（乘）。
由于b=2，故平移量为2时，实际是右移1，符合平移性质。
两实验之二心得体会：时域中的基本运算具有连续性，当输入信号为连续时，输出信号也为连续。
平移，伸缩变化都会导致输出结果相对应的平移伸缩。
2.离散时间信号的时域基本运算两实验之一实验分析：输出信号的值是对应输入信号在每个n值所对应的运算值，当进行拉伸变化后，n值数量不会变，但范围会拉伸所输入的拉伸系数。
两实验之二心得体会：离散时间信号可以看做对连续时间信号的采样，而得到的输出信号值，也可以看成是连续信号所得之后的采样值。
二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析：当两相互卷积函数为冲激函数时，所卷积得到的也是一个冲激函数，且该函数的冲激t值为函数x，函数y冲激t值之和。
两实验之二心得体会：连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数，并一直移动k值直至最后，最后累和出来的最终函数便是所得到的卷积函数。
3.RC电路时域积分两实验之一实验分析：全响应结果正好等于零状态响应与零输入响应之和。
两实验之二心得体会：具体学习了零状态，零输入，全响应过程的状态及变化，与之前所学的电路知识联系在一起了。
三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析：输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会：直观地观察到卷积和的产生，可以看成连续卷积的采样形式，从这个方面去想，更能深入地理解卷积以及采样的知识。
2.离散差分方程求解两实验之一实验分析：其零状态响应序列为0 0 4 5 7.5，零输入响应序列为2 4 5 5.5 5.75，全状态响应序列为2 4 9 10.5 13.25，即全状态=零输入+零状态。
两实验之二心得体会：求差分方程时，可以根据全状态响应是由零输入输入以及零状态相加所得，分开来求，同时也加深了自己对差分方程的求解问题的理解。
学号：14142500190信号与系统实验报告学生姓名张章班级电子14-2BF 院部物理与电子学院专业电子科学与技术任课老师王晓明指导老师王晓明二0一五——二0一六学年第二学期实验项目名称：离散线性时不变系统的分析实验成绩：实验日期： 2016. 6.14 实验室： 6404一、实验目的深刻理解离散时间系统的系统函数在分析离散系统的时域特性、频域特性以及稳定性中的重要作用及意义，熟练掌握利用MATLAB 分析离散系统的时域响应、频响特性和零极点的方法。
二、相关知识点1．离散时间线性时不变系统的时域响应（单位冲击响应、单位阶跃响应、零输入响应、零状态响应和全响应）。
2．离散时间线性时不变系统的频率特性（幅度谱和相位谱）。
3．离散时间线性时不变系统的零极点及稳定性。
三、实验原理、方法和手段1. 离散系统的时域响应离散时间LTI 系统可用如下的线性常系数差分方程来描述：已知输入信号x[k]以及系统初始状态y[-1]，y[-2] ...，就可以求出系统的响应。
MATLAB 提供了计算差分方程的数值解的函数，可以计算上述差分方程描述的离散系统的单位脉冲响应、零状态响应和完全响应等。
在调用MATLAB 函数时，需要利用描述该离散系统的系数函数。
对差分方程进行Z 变换即可得系统函数：在MATLAB 中可用向量a 和向量b 分别保存分母多项式和分子多项式的系数：这些系数均从z 0按z 的降幂排列。
（1） 离散系统的单位冲激响应h [k ]的计算[h,k] = impz(b, a)：计算系统的单位脉冲响应h[k]和相应的时间向量k ；也可简写为：h = impz(b, a)。
][)]1([]1[][11N k y a N k y a k y a k y N N -+--++-+- ][)]1([]1[][110M k x b M k x b k x b k x b M M -+--++-+=- 1(1)0111(1)11()()()()1()M M M M N N N N b b z b z b z Y z b z H z X z a z a z a z a z ----------++++===++++],,,,[110M M b b b b b -= ],,,,1[110-=N a a a a[h,k] = impz(b, a, n)：计算n 点单位脉冲响应h[k]； 也可简写为：h = impz(b, a, n)。
impz(b, a)：绘制单位脉冲响应h[k]的图形。
（2） 离散系统零状态响应y [k ]的计算y = filter(b, a, x)：计算系统在输入x 作用下的零状态响应y[k]； （3）离散系统全响应y [k ] 的计算y = filter(b, a, x, zi)：%计算系统在输入x 和初始状态作用下的完全响应y[k]。
zi 是由系统的初始状态经过filtic 函数转换而得到的初始条件：zi= filtic (b, a, Y0) ，Y0为系统的初始状态，Y0= [ y[-1]，y[-2]，y[-3]，...]。
已知系统函数为 ，计算：1) 离散系统的单位脉冲响应；2) 输入x[k]＝ u[k]，求系统的零状态响应y[k]；3) 输入x[k]＝ u[k]，初始条件y[-1]＝ 1， y[-2]＝ 2，求系统的完全响应y[k]。
1) 计算前40个点的单位脉冲响应N =40; a=[1,0.4,-0.12]; b=[1,2]; y=impz(b,a,N); stem(y)xlabel('k');title('h[k]') 2) 计算前100个点的零状态响应N=100; b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N); y=filter(b,a,x)3) 计算前100个时刻的完全响应filter(b,a,x,zi)中的初始值zi 不是y [-1]＝ 1， y [-2]＝ 2， 它可以由filtic 函数计算。
N=100;11212()10.40.12z H z z z ---+=+-b=[1,2]; a=[1,0.4,-0.12]; x=ones(1,N);zi=filtic(b,a,[1,2]); y=filter(b,a,x,zi); 2．离散系统的频率响应当离散因果LTI 系统的系统函数H (z )的极点全部位于z 平面单位圆内时，系统的频率响应可由H (z )求出，即[H, w]=freqz(b, a, n)；%计算系统的n 点频率响应H ，w 为频率点向量。
H=freqz(b, a, w)； %计算系统在指定频率点向量w 上的频响； freqz(b,a)；%绘制频率响应曲线。
abs(H)；%幅度响应函数，angle(H)；%相位响应。
其中：b 和a 分别为系统函数H (z )的分子分母系数矩阵；已知某离散系统的系统函数为： 分析系统的幅频特性。
b=[1,1]; a=[1,-1,0.5]; [H,w]=freqz(b,a); plot(w,abs(H))xlabel('Frequency(rad)'); ylabel('Magnitude'); title('Magnitude response'); 3．离散系统的系统函数零极点分析离散LTI 系统的系统函数H (z )可以表示为零极点形式：使用MATLAB 提供的roots 函数计算离散系统的零极点； 使用zplane 函数绘制离散系统的零极点分布图。
注意：在利用这些函数时，要求H(z)的分子多项式和分母多项式的系数的个数相等，若不等则需要补零。
1(1)011121(1)1112()()...()()()()1()()...()M M M M M N N N N N b b z b z b z z z z z z z Y z H z kX z a z a z a z z p z p z p ----------++++---===++++---)(j j e j e )e ()()e (j ΩΩ=Ω==ΩϕH z H H z 1121()10.5z H z z z ---+=-+已知系统函数为 计算该系统函数的零极点，并画出系统函数零极点分布图。
b=[1,2,0]; a=[1,0.4,-0.12]; z=roots(b) p=roots(a) zplane(b,a) 四、实验条件计算机一台，MATLAB 软件。
五、实验步骤及结果测试1．某离散线性时不变系统的差分方程如下： ()3(1)2(2)()y n y n y n x n +-+-=满足初始状态y(-1)=0，y(-2)=0.5，求系统输入为x(n)=2n u(n)时的零输入5911212()10.40.12z H z z z---+=+-2．已知离散系统的系统函数为：1231234 1.6 1.64()10.40.350.4z z z H z z z z --------+=++-，求该系统的零极点图，并绘出零极点图，通过零极点的分布判断系统的因果性和稳定性。
b=[4,-1.6,-1.6,4];a=[1,0.4,0.35,-0.4]; z=roots(b); p=roots(a);zplane(b,a);3．已知离散系统的系统函数为：2462460.10.20.20.1()10.60.40.04z z z H z z z z -------+-=+++，绘出系统的零极点分布图，系统在0~π频率范围内的绝对幅频响应、相对幅频响应、相位频率响应和群延迟，并观察零极点分布与系统幅频响应的关系。
b=[0.1,-0.2,0.2,-0.1]; a=[1,0.6,0.4,0.04]; z=roots(b); p=roots(a); zplane(b,a);b=[0.1,-0.2,0.2,-0.1]; a=[1,0.6,0.4,0.04]; [H,w]=freqz(b,a); plot(w,abs(H));xlaber('Frequency(rad)'); ylaber('Magnitude');title('Magnitude response');b=[0.1,-0.2,0.2,-0.1]; a=[1,0.6,0.4,0.04]; [H,w]=freqz(b,a); plot(w,angle(H));xlaber('Frequency(rad)'); ylaber('Magnitude'); title('Magnitude response');4．已知一阶离散系统的系统函数为：11()H z z p =- ， （1）假设系统函数的零点在原点，极点分别取0.2、0.5、0.8，比较它们的幅频响应曲线。
（2）假设系统的极点在原点，零点分别取0.2、0.5、0.8，比较它们的幅频响应曲线，从中总结零极点位置对幅频响应的影响。
b=[1,0]; a=[1,-0.2];[H,w]=freqz(b,a); plot(w,abs(H));xlaber('Frequency(rad)'); ylaber('Magnitude');title('Magnitude response');b=[1,0]; a=[1,-0.5];[H,w]=freqz(b,a); plot(w,abs(H));xlaber('Frequency(rad)'); ylaber('Magnitude');title('Magnitude response');b=[1,0]; a=[1,-0.8];[H,w]=freqz(b,a); plot(w,abs(H));xlaber('Frequency(rad)'); ylaber('Magnitude');title('Magnitude response');b=[1,-0.2];[H,w]=freqz(b,a);plot(w,abs(H));ylaber('Magnitude');a=[1,0];[H,w]=freqz(b,a);plot(w,abs(H));xlaber('Frequency(rad)');ylaber('Magnitude');a=[1,0];[H,w]=freqz(b,a);plot(w,abs(H));xlaber('Frequency(rad)');ylaber('Magnitude');。
