连续系统复频域分析报告附MATLAB实现信号与系统实验报告.docxVIP

计算机与信息工程学院设计性实验报告 专业： 通信工程 年级 / 班级： 2011 级 第二学年 第二学期 课程名称 信号与系统 指导教师 谢自梅 学号姓名 实验地点 计算机学院 111 实验时间 项目名称 连续系统的复频域分析及 实验类型 设计性 MATLAB实现 一、实验目的 1.掌握用 matlab 分析系统时间响应的方法 2.掌握用 matlab 分析系统频率响应的方法 3.掌握系统零、极点分布与系统稳定性关系 二、实验原理 1.系统函数 H(s) 系统函数 : 系统零状态响应的拉氏变换与激励的拉氏变换之比 . H(s)=R(s)/E(s) matlab中可采用多种方法描述系统 ,本文采用传递函数 ( 系统函数 ) 描述法 . matlab中 , 传递函数描述法是通过传递函数分子和分母关于 s降幂排列的多项式系数来表示的 . 例如 , 某系统传递函数如下 H (s) s 1 (1) s2 1.3s 0.8 则可用如下二个向量 num和 den来表示 : num=[1,1] ； den=[1,1.3,0.8] 2.用matlab 分析系统时间响应 脉冲响应 y=impulse(num,den,T) 为等间隔的时间向量 , 指明要计算响应的时间点 . 2) 阶跃响应 y=setp(num,den,T) T同上 . 对任意输入的响应 y=lsim(num,den,U,T) U:任意输入信号 . T 同上 . 用matlab 分析系统频率响应特性 频响特性 : 系统在正弦激励下稳态响应随信号频率变化的特性 . j ( ) H ( j ) H (s) s j H ( j ) e 1 / 8 ( ): 相频响应特性 ( 或相移特性 ). Matlab 求系统频响特性函数 freqs 的调用格式 : h=freqs(num,den, ) 为等间隔的角频率向量 , 指明要计算响应的频率点 . 系统零、极点分布与系统稳定性关系 系统函数 H(s) 集中表现了系统的性能 , 研究 H(s) 在S平面中极点分布的位置 , 可很方面地判断系统稳定性 . 稳定系统 : H(s) 全部极点落于 S左半平面 ( 不包括虚轴 ), 则可以满足 lim [ h(t)] 0 t 系统是稳定的 . 不稳定系统 : H(s) 极点落于 S右半平面 , 或在虚轴上具有二阶以上极点 , 则在足够长时间后 ,h(t) 仍继续增长 , 系统是不稳定的 . 临界稳定系统 : H(s) 极点落于 S平面虚轴上 , 且只有一阶 , 则在足够长时间后 ,h(t) 趋于一个非零数值或形成一个等幅振荡 . 系统函数 H(s) 的零、极点可用 matlab 的多项式求根函数 roots() 求得 . 极点 :p=roots(den) 零点 :z=roots(num) 根据 p和z用 plot() 命令即可画出系统零、极点分布图 , 进而分析判断系统稳定性 . 三、实验内容 s 设 H (s) ( s p1)( s p2) 设① p1=-2,p2=-30; ②p1=-2,p2=3 针对极点参数①② , 画出系统零、极点分布图 , 判断该系统稳定性 . 针对极点参数①② ,绘出系统的脉冲响应曲线 , 并观察 t →∞时 , 脉冲响应变化趋势 . 针对极点参数① , 绘出系统的频响曲线 . 四、实验要求 1．预习实验原理； 2．对实验内容编写程序 (M文件 ), 上机运行； 3．绘出实验内容的各相应曲线或图。 五、实验设备 1.装 MATLAB 软件的计算机 1 台 2 / 8 六、实验步骤及结果 1、针对极点参数①② , 画出系统零、极点分布图 , 判断该系统稳定性 . 已知： H ( s)  s ( s p1)(s p2) p1=-2,p2=-30 s 则 H ( s) (s+2)( s+30) 则 den1 =[1 32 60]; num1 =[1 0]; p1=-2,p2=3 s 则 H (s) (s+2)( s+3) 则 den2 =[1 5 6]; num2 =[1 0]; 编写程序如下： num=[1 0]; den=conv([1,2],[1,30]); sys=tf(num,den); poles=roots(den) figure(1); pzmap(sys); grid on num=[1 0]; den=conv([1,2],[1,-3]); sys=tf(num,den); poles=roots(den) figure(2); pzmap(sys); grid on 系统零、极点分布图如下图所示： 3 / 8 ) 1 d n o c e s ( si x A y r a ni g a m I  Pole-Zero Map