云师大数字信号处理实验07第二篇实验一连续系统分析..docVIP

本科学生实验报告 学号 124090380 姓名 赵玲 学院 物电学院 专业、班级 12电子 实验课程名称 连续系统法分析 教师及职称 和伟 开课学期 2014 至 2015 学年 下 学期 填报时间 2015 年 5 月 13 日 云南师范大学教务处编印 实验序号 07 实验名称 连续系统分析 实验时间 2015年5月13 实验室 同析3栋数字信号处理实验室 一．实验预习 实验目的 理解连续时间系统的系统函数在分析连续系统的时域特性、频率特性及稳定性中的重要作用及意义。 掌握利用MATLAB分析连续系统的时域响应、频率响应和零极点的基本方法。 掌握根据系统函数的零极点设计简单的滤波器的方法。 实验原理、实验流程或装置示意图 matlab提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频率响应等分析函数。 连续时间LTI系统的输入x(t)与输出y(t)可用如下的线性常系数微分方程系统函数，系统时域响应，系统频率响应等分析函数。 连续时间LTI系统的输入x（t）与输出y（t）可用如下的线性常系数微分方程来描述： . 如果已知系统的输入信号的表示式以及系统的初始状态就可以利用解析方法求出系统的响应，但对于高阶微分方程描述的连续系统，解析计算将会变得非常繁琐和困难。MATLAB提供了微分方程的数值计算函数，可以计算上述n阶微分方程描述的连续系统的响应，包括系统的单位冲激响应，单位阶跃响应，零输入响应，零状态响应和完全响应。在调用这些函数时，需要利用连续系统对应的系统函数对微分方程进行拉氏变换即可得系统函数: 在MATLAB中可使用向量分别保存分子多项式和分母多项式的系数，这些系数均按s的降幂直至s0排列。 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频率响应等分析函数。 1．系统的零极点分析 连续LTI系统的系统函数可以表示为部分分式，即p1,p2,……,pn为H(s)的极点（pole） z1,z2,……,zn为H(s)的零点（zero）k为增益） 通过拉氏反变换可求得系统的单位冲激响应h(n)。设m=n，且H(s)的极点pi全部为单极点，则： 对应可以看出，系统函数H(s)的极点pi决定了冲激响应h(t)的基本形式，而零点和极点共同确定了冲激响应h(t)的幅值ki。 在MATLAB中提供了roots函数计算系统的零极点，pzmap函数绘制连续系统的零极点分布图。 连续系统的频率响应: 若连续因果LTI系统函数H(s)的极点全部位于s左半平面（前提条件），则系统的频率响应|H(jw)|可由H(s)求出，即 在MATLAB中freqs函数可以分析连续系统的频率响应，其格式如下： H=freqs(b,a,w)计算系统在指定频率点向量w上的频率响应Ｈ,b=[bm,bm-1,…,b1,b0]和a=[am,am-1,…,a1,a0]分别保存分子多项式和分母多项式的系数。 [H,W]=freqs(b,a)自动选取200个频率点计算频率响应，W为频率点向量。 [H,W]=freqs(b,a,n)计算n个点的频率响应。freqs(b,a)自动绘制系统频率响应曲线。 impulse（sys）计算并画出连续系统的冲激响应，sys可由函数获得。其中b和a分别是系统函数的分子多项式和分母多项式的系数矩阵。 b= impulse（sys,t）计算并画出连续系统在向量t定义的时间范围内的冲激响应，向量h保存对时间的冲激响应的输出量。 3． 实验设备及材料 一台带matlab软件的计算机。 4．实验方法步骤及注意事项 实验方法步骤： 打开MATLAB软件 根据题目要求编写程序 运行程序 分析实验结果 关闭计算机 注意事项： (1)对于实验电脑要爱惜，遵守实验的规则。 (2)程序运行前要检查程序是否正确。 在使用matlab编程时，应该养成良好的编写习惯，新建一个flies编写。 一些快捷键的使用，能提高编程效率。 Help能查询到不懂使用的函数使用方法,比如用到的tf、zpk等函数。如下界面： 二、实验内容 1、已知描述连续系统的微风方程为dy(t)/dt+10y(t)=2x(t),输入x(t )=u(t),初始状态位y(0- )= 1,计算该系统的响应，并与理论结果比较，列出系统响应分析的步骤。 %1_1……该系统的单位冲激响应 clc,