线性系统的频率响应分析实验

实验三 线性系统的频率响应分析

一、实验目的：掌握实验方法测量系统的波特图和极坐标图。

二、实验设备：PC 机一台，TD-ACC 系列教学实验系统一套。

三．实验对象的结构框图及模拟电路图。

图1

开环传函为： 1()0.1(0.11)

G S S S =+ 闭环传函： 221100()0.010.1110100

S S S S S Φ==++++ 得转折频率ω=10(rad/s) 阻尼比ξ＝0.5。

四、实验步骤:

此次实验，采用直接测量方法测量对象的闭环波特图及奈奎斯特图。将信号源单元的“ST”插针分别与“S”插针和“+5V”插针断开，运放的锁零控制端“ST”此时接至控制计算机单元的“DOUT0”插针处，锁零端受“DOUT0”来控制。将数模转换单元的“/CS”接至控制计算机的“/IOY1”，数模转换单元的“OUT1”，接至图1中的信号输入端.

1．实验接线：按模拟电路图图1接线，检查无误后方可开启设备电源。

2．直接测量方法 (测对象的闭环波特图)

(1)“CH1”路表笔插至图1中的4＃运放的输出端。

(2) 打开集成软件中的频率特性测量界面，弹出时域窗口，点击按钮，在弹出的窗口中根据需要设置好几组正弦波信号的角频率和幅值，选择测量方式为“直接”测量，每

组参数应选择合适的波形比例系数，具体如下图所示

(3) 确认设置的各项参数后，点击按钮，发送一组参数，待测试完毕，显示时域波形，此时需要自行移动游标，将两路游标同时放置在两路信号的相邻的波峰（波谷）处，或零点处，来确定两路信号的相位移。两路信号的幅值系统将自动读出。重复操作(3)，直到所有参数测量完毕。

(4) 待所有参数测量完毕后，点击按钮，弹出波特图窗口，观察所测得的波特图，该图由若干点构成，幅频和相频上同一角频率下两个点对应一组参数下的测量结果。将波特图绘制或保存下来。

点击极坐标图按钮，可以得到对象的闭环极坐标。将极坐标图绘制或保存下来。

(5) 根据所测图形可适当修改正弦波信号的角频率和幅值重新测量，达到满意的效果。

五、注意：

(1) 测量过程中要去除运放本身的反相的作用，即保持两路测量点的相位关系与运放无关，所以在测量过程中可能要适当加入反相器，滤除由运放所导致的相位问题。

(2) 测量过程中，可能会由于所测信号幅值衰减太大，信号很难读出，须放大，若放大的比例系数不合适，会导致测量误差较大。所以要适当地调整误差或反馈比例系数。

离散系统频率响应和零极点分布实验报告

实验报告 课程名称数字信号处理实验 实验名称离散系统频率响应和零极点分布 学生姓名 学生学号 学生班级 实验日期 实验目的：通过matlab仿真简单的离散时间系统，研究其的时频域特性，加深对离散系统的频率响应分析和零、极点分布的概念理解。 实验原理：离散系统的时域方程为

∑∑==-=-M k k N k k k n x p k n y d ) ()( 其变换域分析方法如下： 频域 ) ()()(][][][][][ωωωj j j m e H e X e Y m n h m x n h n x n y =?-= *=∑∞ -∞ = 系统的频率响应为 ω ωω ωωωω jN N j jM M j j j j e d e d d e p e p p e D e p e H ----++++++==......)()()(1010 Z 域 ) ()()(][][][][][z H z X z Y m n h m x n h n x n y m =?-= *=∑∞ -∞ = 系统的转移函数为 N N M M z d z d d z p z p p z D z p z H ----++++++==......)()()(110110 分解因式 ∏-∏-=∑∑= =-=-=-=-N i i M i i N i i k M i i k z z K z d z p z H 11 11 )1()1()(λξ ，其中 i ξ和i λ称为零、极 点。 在MATLAB 中，可以用函数[z ，p ，K]=tf2zp （num ，den ）求得有理分式形式的系统转移函数的零、极点，用函数zplane （z ，p ）绘出零、极点分布图；也可以用函数zplane （num ，den ）直接绘出有理分式形式的系统转移函数的零、极点分布图。 另外，在MATLAB 中，可以用函数 [r ，p ，k]=residuez （num ，den ）完成部分分式展开计算；可以用函数sos=zp2sos （z ，p ，K ）完成将高阶系统分解为2阶系统的串联。 实验内容：一个LTI 离散时间系统的输入输出差分方程为 y(n)-1.6y(n-1)+1.28y(n-2) =0.5x(n)+0.1x(n-1) 实验要求：(1)编程求出此系统的单位冲激响应序列，并画出其波形。 （2）若输入序列x(n)=δ(n)+2δ(n-1)+3δ(n-2)+4δ(n-3)+5δ (n-4),編程求此系統輸出序列y(n),并画出其波形。 （3）编程得到系统频响的幅度响应和相位响应，并画图。 （4）编程得到系统的零极点分布图，分析系统的因果性和稳定性。 实验过程： （1）编程求此系统的单位冲激响应序列，并画出其波形 程序

视觉分辨率及空间频率响应测试实验报告

视觉分辨率及空间频率响应（SFR）测试实验报告 班级：学号：姓名： 一、实验目的： 1、理解数码相机视觉分辨率的定义及其度量单位。 2、了解数码相机分辨率测试标准ISO12233以及GB/T 19953-2005《数码相机分辨率的测量》，熟悉测试标板构成，掌握其使用方法。 3、掌握数码相机视觉分辨率测试方法，能够通过目视判别数码相机的分辨率特性。 4、了解数码相机空间频率响应（SFR）的测试原理，理解空间频率响应（SFR）曲线的含义。 5、掌握数码相机空间频率响应（SFR）的测试方法，能够通过SFR曲线判别数码相机的分辨率特性。 二、实验要求： 1、使用数码相机拍摄ISO12233标准分辨率靶板,要求连续拍摄三幅图。 2、目视判别数码相机的视觉分辨率，需分别判别水平、垂直、和斜45度方向的视觉分辨率(注意：若拍摄的靶板有效区域高度仅占据相机幅面高度的一部分，需将目视判别结果乘以修正系数以得到真实的测量结果。修正系数=以像素为单位的相机幅面高度/以像素为单位的靶板有效区域高度)。 3、使用Imatest软件测量数码相机空间频率响应（SFR）曲线，需分别测量水平及垂直方向的SFR，并取MTF50、MTF20作为测量结果，与视觉分辨率测试结果进行比较。 4、独立完成实验报告，需明确相机型号、相机基本设置、并包含所拍摄图案以及判别结果和相应说明。 三、实验过程 在光学测量实验室使用手机（iPhone6s）连续拍摄三张ISO12233标准分辨率靶板。拍摄过程中使手机上下屏幕边缘尽量与靶板上下边缘对齐，以减小修正系数。其中使用的相机参数如下：

拍摄的照片如下： 照片一（修正系数为）

实验二实验报告

PAM和PCM编译码器系统 一、实验目的 1.观察了解PAM信号形成的过程；验证抽样定理；了解混叠效应形 成的原因； 2.验证PCM编译码原理；熟悉PCM抽样时钟、编码数据和输入/输出 时钟之间的关系；了解PCM专用大规模集成电路的工作原理和应用。 二、实验内容和步骤 1.PAM编译码器系统 1.1自然抽样脉冲序列测量 （1）准备工作； （2）PAM脉冲抽样序列观察； （3）PAM脉冲抽样序列重建信号观测。 1.2平顶抽样脉冲序列测量 （1）准备工作； （2）PAM平顶抽样序列观察； （3）平顶抽样重建信号观测。 1.3信号混叠观测 （1）准备工作 （2）用示波器观测重建信号输出的波形。 2.PCM编译码器系统 2.1PCM串行接口时序观察 （1）输出时钟和帧同步时隙信号的观察； （2）抽样时钟信号与PCM编码数据测量； 2.2用示波器同时观察抽样时钟信号和编码输出数据信号端口 （TP502），观测时以TP504同步，分析掌握PCM编码输数据和抽样时钟信号（同步沿、脉冲宽度）及输出时钟的对应关系； 2.3PCM译码器输出模拟信号观测，定性观测解码信号与输入信号的 关系：质量，电平，延时。 2.4PCM频率响应测量：调整测试信号频率，定性观察解码恢复出的 模拟信号电平，观测输出信号电平相对变化随输入信号频率变化的相对关系；

2.5PCM动态范围测量：将测试信号频率固定在1000Hz，改变测试信 号电平，定性观测解码恢复出的模拟信号的质量。 三、实验数据处理与分析 1.PAM编译码器系统 （1）观察得到的抽样脉冲序列和正弦波输入信号如下所示： 上图中上方波形为输入的正弦波信号，下方为得到的抽样脉冲序列，可见抽样序列和正弦波信号基本同步。 （2）观测得到的重建信号和正弦波输入信号如下所示：

信号与系统实验报告1

学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟

20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求

实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。

实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;

线性控制系统的频率响应分析

一．实验目的 1．了解和掌握对数幅频曲线和相频曲线（波德图）、幅相曲线（奈奎斯特图）的构造及绘制方法。 2．二阶开环系统中的相位裕度和幅值穿越频率的计算。 二．实验内容及要求 1．一阶惯性环节的频率特性曲线测试。 2．二阶开环系统的频率特性测试，研究表征系统稳定程度的相位裕度和 幅值穿越频率对系统的影响。 三、实验主要仪器设备和材料 1．labACT自控/计控原理实验机一台 2．数字存储示波器一台 四、实验方法、步骤及结果测试 1．一阶惯性环节的频率特性曲线 惯性环节的频率特性测试模拟电路见图4-1。 图4-1 惯性环节的频率特性测试模拟电路 实验步骤：注：‘S ST'不能用“短路套”短接！ （1）将数/模转换器（B2）输出OUT2作为被测系统的输入。 （2）按图4-1安置短路套及测孔联线。 （3）运行、观察、记录： ①运行LABACT程序，选择自动控制菜单下的线性控制系统的频率响应分析-实验项目，选择一阶系统，再选择开始实验，点击开始，实验机将自动产生0.5Hz~64Hz多个频率信号，测试被测系统的频率特性，等待将近十分钟，测试结束。 ②测试结束后，可点击界面下方的“频率特性”选择框中的任意一项进行切换，将显示被测系统的对数幅频、相频特性曲线（伯德图）和幅相曲线（奈 奎斯特图），同时在界面上方将显示点取的频率点的L、、Im、Re等相关数

据。如点击停止，将停止示波器运行，不能再测量数据。 ③分别改变惯性环节开环增益与时间常数，观察被测系统的开环对数幅频曲线、相频曲线及幅相曲线，在幅频曲线或相频曲线上点取相同的频率点，测量、记录数据于实验数据表中。 实验数据表1：改变惯性环节开环增益，（T=0.05，C=1u，R2=50K） 实验数据表2： 改变惯性环节时间常数， K=1（R1=50K、R2=50K） 2．二阶开环系统的频率特性曲线 二阶系统模拟电路图的构成如图4-2所示。

数字信号处理实验报告 -频率响应与系统稳定性

专业：电子信息工程班级：N11级-1F 姓名： 学号：

实验项目：系统响应及系统稳定性 实验台号：同组者： 1、实验目的 （1）掌握求系统响应的方法 （2）掌握时域离散系统的时域特性 （3）分析、观察及判断系统的稳定性 2、实验原理与方法 描述系统特性有多种方式，时域描述有差分方程和单位脉冲响应，频域描述有系统函数和频率响应。已知输入信号可以由差分方程、单位脉冲响应、系统函数或频率响应来求系统的输出信号。 （1）求系统响应：本实验仅在时域求系统响应。在计算机上，已知差分方程可调用filter函数求系统响应；已知单位脉冲响应可调用conv函数计算系统响应。 （2）系统的时域特性：系统时域特性是指系统的线性、时不变性、因果性和稳定性。本实验重点分析系统的稳定性，包括观察系统的暂态响应和稳态响应。 （3）系统的稳定性判断：系统的稳定性是指对任意有外接信号输入，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和条件。实际中，检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否是有界输出,或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳 定的。

（4）系统的稳态响应 系统的稳态输出是指当∞→n 时，系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n 的加大，幅度趋于稳定，达到稳态输出。注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤 （1）已知差分方程求系统响应 设输入信号 )()(81n R n x =，) ()(2n u n x =。已知低通滤波器的差分方程为 )1(9.0)1(05.0)(05.0)(-+-+=n y n x n x n y 。 试求系统的单位冲响应，及系统对)()(81n R n x =和)()(2n u n x =的输出信号，画出输出波形。 051015 20253035404550 n h n 系统的单位脉冲响应 5 10 15 20 253035 40 45 50 n y 1n 系统对R8(n)的响应 05101520 253035404550 n y 2n 系统对u(n)的响应 实验图（1） （2）已知单位脉冲响应求系统响应 设输入信号 )()(8n R n x =，已知系统的单位脉冲响应分别为)()(101n R n h =， )3()2(5.2)1(5.2)()(2-+-+-+=n n n n n h δδδδ，试用线性卷积法分别求出 各系统的输出响应，并画出波形。

实验二：频率响应测试

成绩 北京航空航天大学 自动控制原理实验报告 院（系）名称自动化科学与电气工程学院 专业名称自动化 学生学号13191006________ 学生________ 万赫__________ 指导老师_____ 王艳东 自动控制与测试教学实验中心

实验二频率响应测试 实验时间2015.11.13 实验编号30 同组同学无 一、实验目的 1、掌握频率特性的测试原理及方法 2、学习根据所测定出的系统的频率特性，确定系统传递函数的方法 目的。 二、实验容 1. 测定给定环节的频率特性。 2. 系统模拟电路图如下图： 系统结构图如下图：

系统的传递函数： 取R=100KΩ，则G(s)=错误!未找到引用源。 取R=200KΩ，则G(s)=错误!未找到引用源。 取R=500KΩ，则G(s)=错误!未找到引用源。 若正弦输入信号为Ui(t)=A1Sin(ωt),则当输出达到稳态时，其输出信号为 Uo(t)=A2Sin(ωt+ψ)。改变输入信号频率f=错误!未找到引用源。值，便可测得二组A1/A2和ψ随f(或ω)变化的数值，这个变化规律就是系统的幅频特性和相频特性。 三、实验原理 1. 幅频特性即测量输入与输出信号幅值A1及A2，然后计算其比值A2/A1。 2. 实验采用“沙育图形”法进行相频特性的测试。 设有两个正弦信号: X(ωt)=XmSin(ωt) ，Y(ωt)=YmSin(ωt+ψ) 若以X(t)为横轴，Y(t)为纵轴，而以ω作为参变量，则随着ωt的变化，X(t)和Y(t)所确定的点的轨迹，将在X-Y平面上描绘出一条封闭的曲线。这个图形就是物理学上成称

实验三 二阶系统频率响应

实验三 二阶系统频率响应 一、实验目的 （1）学习系统频率特性响应的实验测试方法。 （2）了解二阶闭环系统中的对数幅频特性和相频特性的计算。 （3）掌握根据频率响应实验结果绘制波特图的方法。 （4）掌握欠阻尼二阶闭环系统中的自然频率、阻尼比对谐振频率、谐振峰值和带宽的影响及对应的计算。 二、实验设备 （1）XMN-2型学习机； （2）CAE-USE 辅助实验系统 （3）万用表 （4）计算机 三、实验内容 本实验用于观察和分析二阶系统瞬态响应的稳定性。 二阶闭环系统模拟电路如图3-1所示，它由两个积分环节（OP1和OP2）及其反馈回路构成。 图3-1 二阶闭环系统模拟电路图 OP1和OP2为两个积分环节，传递函数为s T s G i 1 )(-=（时间常数RC T i =）。二阶闭环系统等效结构图如图3-2所示。 图3-2 二阶闭环系统等效结构图 该二阶系统的自然振荡角频率为RC T n 11==ω，阻尼为i f R R K 22= =ζ。 四、实验步骤 （1）调整Rf=40K ，使K=0.4（即ζ=0.2）；取R=1M ，C=1μ，使T=1秒（ωn=1/1）。 （2）输入信号位)sin(t X ω=，改变角频率使ω分别为 0.2,0.6,0.8,0.9,1.0,1.2,1.6,2.0,3.0rad/s 。稳态时，记录下输出响应)sin(φω+=t Y y 五、数据采集及处理 输出信号幅值Y 输出信号初相φ L(ω) φ(ω) ω(rad/s) T 0.2 0.6 0.8 0.9 1.0 1.2

1.6 2.0 3.0 六、实验报告 1、绘制系统结构图，并求出系统传递函数，写出其频率特性表达式。 2、用坐标纸画出二阶闭环系统的对数幅频、相频曲线（波特图）。 3、其波特图上分别标示出谐振峰值（Mr）、谐振频率（ωr）和带宽频率（ωb）。 4、观察和分析曲线中的谐振频率（ωr）、谐振峰值（Mr）和带宽（ωb），并与理论计算值作对比。

信号检测实验报告

Harbin Institute of Technology 匹配滤波器实验报告 课程名称：信号检测理论 院系：电子与信息工程学院 姓名：高亚豪 学号：14SD05003 授课教师：郑薇 哈尔滨工业大学

1. 实验目的 通过Matlab 编程实现对白噪声条件下的匹配滤波器的仿真，从而加深对匹配滤波器及其实现过程的理解。通过观察输入输出信号波形及频谱图，对匹配处理有一个更加直观的理解，同时验证匹配滤波器具有时间上的适应性。 2. 实验原理 对于一个观测信号()r t ，已知它或是干扰与噪声之和，或是单纯的干扰， 即 0()()()()a u t n t r t n t +?=?? 这里()r t ，()u t ，()n t 都是复包络，其中0a 是信号的复幅度，()u t 是确知的归一化信号的复包络，它们满足如下条件。 2|()|d 1u t t +∞ -∞=? 201||2 a E = 其中E 为信号的能量。()n t 是干扰的均值为0，方差为0N 的白噪声干扰。 使该信号通过一个线性滤波系统，有效地滤除干扰，使输出信号的信噪比在某一时刻0t 达到最大，以便判断信号的有无。该线性系统即为匹配滤波器。 以()h t 代表系统的脉冲响应，则在信号存在的条件下，滤波器的输出为 0000()()()d ()()d ()()d y t r t h a u t h n t h τττττττττ+∞+∞+∞ =-=-+-???

右边的第一项和第二项分别为滤波器输出的信号成分和噪声成分，即 00()()()d x t a u t h τττ+∞ =-? 0 ()()()d t n t h ?τττ+∞ =-? 则输出噪声成分的平均功率（统计平均）为 2 20E[|()|]=E[|()()d |]t n t h ?τττ+∞ -? **00*000200 =E[()(')]()(')d d '=2()(')(')d d ' 2|()|d n t n t h h N h h N h ττττττδττττττττ+∞+∞+∞+∞+∞ ---=?? ?? ? 而信号成分在0t 时刻的峰值功率为 22 20000|()||||()()d |x t a u t h τττ+∞ =-? 输出信号在0t 时刻的总功率为 22000E[|()|]E[|()()|]y t x t t ?=+ 22**0000002200E[|()||()|()()()()] |()|E[|()|] x t t x t t t x t x t t ????=+++=+ 上式中输出噪声成分的期望值为0，即0E[()]0t ?=，因此输出信号的功率 成分中只包含信号功率和噪声功率。 则该滤波器的输出信噪比为 222000022000|||()()d ||()|E[|()|]2|()|d a u t h x t t N h τττρ?ττ+∞ +∞-==?? 根据Schwartz 不等式有

《测试信号分析与处理》实验报告

测控1005班齐伟0121004931725 （18号）实验一差分方程、卷积、z变换 一、实验目的 通过该实验熟悉 matlab软件的基本操作指令，掌握matlab软件的使用方法，掌握数字信号处理中的基本原理、方法以及matlab函数的调用。 二、实验设备 1、微型计算机1台； 2、matlab软件1套 三、实验原理 Matlab 软件是由mathworks公司于1984年推出的一套科学计算软件，分为总包和若干个工具箱，其中包含用于信号分析与处理的sptool工具箱和用于滤波器设计的fdatool工具箱。它具有强大的矩阵计算和数据可视化能力，是广泛应用于信号分析与处理中的功能强大且使用简单方便的成熟软件。Matlab软件中已有大量的关于数字信号处理的运算函数可供调用，本实验主要是针对数字信号处理中的差分方程、卷积、z变换等基本运算的matlab函数的熟悉和应用。 差分方程（difference equation)可用来描述线性时不变、因果数字滤波器。用x表示滤波器的输入，用y表示滤波器的输出。 a0y[n]+a1y[n-1]+…+a N y[n-N]=b0x[n]+b1x[n-1]+…+b M x[n-M] (1) ak,bk 为权系数，称为滤波器系数。 N为所需过去输出的个数，M 为所需输入的个数卷积是滤波器另一种实现方法。 y[n]= ∑x[k] h[n-k] = x[n]*h[n] (2) 等式定义了数字卷积，*是卷积运算符。输出y[n] 取决于输入x[n] 和系统的脉冲响应h[n]。 传输函数H(z)是滤波器的第三种实现方法。 H(z)=输出/输入= Y(z)/X(z) （3）即分别对滤波器的输入和输出信号求z变换，二者的比值就是数字滤波器的传输函数。 序列x[n]的z变换定义为 X (z)=∑x[n]z-n (4) 把序列x[n] 的z 变换记为Z{x[n]} = X(z)。

信号与系统实验报告实验三 连续时间LTI系统的频域分析

实验三 连续时间LTI 系统的频域分析 一、实验目的 1、掌握系统频率响应特性的概念及其物理意义； 2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用； 3、学习和掌握幅度特性、相位特性以及群延时的物理意义； 4、掌握用MA TLAB 语言进行系统频响特性分析的方法。 基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理解LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。 二、实验原理及方法 1 连续时间LTI 系统的频率响应 所谓频率特性，也称为频率响应特性，简称频率响应（Frequency response ），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。 上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：)(*)()(t h t x t y =，由傅里叶变换的时域卷积定理可得到： )()()(ωωωj H j X j Y = 3.1 或者： ) () ()(ωωωj X j Y j H = 3.2 )(ωj H 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即 ? ∞ ∞ --= dt e t h j H t j ωω)()( 3.3

由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数，因此，也可以表示成复数的不同表达形式。在研究系统的频率响应时，更多的是把它表示成极坐标形式： ) ()()(ω?ωωj e j H j H = 3.4 上式中，)j (ωH 称为幅度频率相应（Magnitude response ），反映信号经过系统之后，信号各频率分量的幅度发生变化的情况，)(ω?称为相位特性（Phase response ），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。)(ωj H 和)(ω?都是频率ω的函数。 对于一个系统，其频率响应为H(j ω)，其幅度响应和相位响应分别为)(ωj H 和)(ω?，如果作用于系统的信号为t j e t x 0)(ω=，则其响应信号为 t j e j H t y 0)()(0ωω= t j j e e j H 00)(0)(ωω?ω=))((000)(ω?ωω+=t j e j H 3.5 若输入信号为正弦信号，即x(t) = sin(ω0t )，则系统响应为 ))(sin(|)(|)sin()()(00000ω?ωωωω+==t j H t j H t y 3.6 可见，系统对某一频率分量的影响表现为两个方面，一是信号的幅度要被)(ωj H 加权，二是信号的相位要被)(ω?移相。 由于)(ωj H 和)(ω?都是频率ω的函数，所以，系统对不同频率的频率分量造成的幅度和相位上的影响是不同的。 2 LTI 系统的群延时 从信号频谱的观点看，信号是由无穷多个不同频率的正弦信号的加权和（Weighted sum ）所组成。正如刚才所述，信号经过LTI 系统传输与处理时，系统将会对信号中的所有频率分量造成幅度和相位上的不同影响。从相位上来看，系统对各个频率分量造成一定的相位移（Phase shifting ），相位移实际上就是延时（Time delay ）。群延时（Group delay ）的概念能够较好地反

振动测试技术模态实验报告

研究生课程论文(2016-2017学年第二学期) 振动测试技术 研究生：

模态试验大作业 0 模态试验概述 模态试验（modal test）又称试验模态分析。为确定线性振动系统的模态参数所进行的振动试验。模态参数是在频率域中对振动系统固有特性的一种描述，一般指的是系统的固有频率、阻尼比、振型和模态质量等。 模态试验中通过对给定激励的系统进行测量，得到响应信号，再应用模态参数辨识方法得到系统的模态参数。由于振动在机械中的应用非常普遍。振动信号中包含着机械及结构的内在特性和运行状况的信息。振动的性质体现着机械运行的品质，如车辆、航空航天设备等运载工具的安全性与舒适性；也反映出诸如桥梁、水坝以及其它大型结构的承载情况、寿命等。同时，振动信号的发生和提取也相对容易因此，振动测试与分析已成为最常用、最基本的试验手段之一。 模态分析及参数识别是研究复杂机械和工程结构振动的重要方法，通常需要通过模态实验获得结构的模态参数即固有频率、阻尼比和振型。模态实验的方法可以分为两大类：一类是经典的纯模态实验方法，该方法是通过多个激振器对结构进行激励，当激振频率等于结构的某阶固有频率，激振力抵消机构内部阻尼力时，结构处于共振状态，这是一种物理分离模态的方法。这种技术要求配备复杂昂贵的仪器设备，测试周期也比较长；另一类是数学上分离模态的方法，最常见的方法是对结构施加激励，测量系统频率响应函数矩阵，然后再进行模态参数的识别。 为获得系统动态特性，常需要测量系统频响函数。目前频响函数测试技术可以分为单点激励单点测量( SISO)、单点激励多点测量( SIMO) 、多点激励多点测量( MIMO)等。单点激励一般适用于较小结构的频响函数测量，多点激励适用于大型复杂机构，如机体、船体或大型车辆机构等。按激励力性质的不同，频响函数测试分为稳态正弦激励、随机激励及瞬态激励三类，其中随机激励又有纯随机、伪随机、周期随机之分。瞬态激励则有快速正弦扫描激励、脉冲激励和阶跃激励等几种方式。按激励力性质的不同，频响函数测试分为稳态正弦激励、随机激励及瞬态激励三类，其中随机激励又有纯随机、伪随机、周期随机之分，瞬态激励则有快速正弦扫描激励、脉冲激励和阶跃激励等几种方式。 振动信号的分析和处理技术一般可分为时域分析、频域分析、时频域分析和时间序列建模分析等。这些分析处理技术从不同的角度对信号进行观察和分析，为提取与设备运行状态有关的特征信息提供了不同的手段。信号的时域分析包括时域统计分析、时域波形分析和时域相关分析。对评价设备运行状态和

振动实验报告

振动力学实验报告 学院：___________________ 班级：___________________ 学号：___________________ 姓名：___________________ 山东科技大学

单自由度系统振动实验报告 实验者姓名：________ 院系：_______系_______专业_______班_______组实验日期：________年________月________日 自由振动法测量单自由度系统的参数 一、实验目的 二、实验对象和装置 三、实验步骤 四、实验数据记录和整理 1、无阻尼单自由度自由振动系统实验测量：

计算单自由度振动的振动频率、周期、固有频率、衰减系数、相对阻尼系数周期、频率和阻尼系数： 2、有阻尼单自由度自由振动系统实验测量： 计算单自由度振动的振动频率、周期、固有频率、阻尼系数、相对阻尼系数： 五、简答 1、上述无阻尼自由振动实验中，为什么振动曲线呈现衰减状态？ 2、简述阻尼对于自由振动周期、频率的影响。

用冲击激励法测量系统的频率响应函数 实验者姓名：________ 院系：_______系_______专业_______班_______组实验日期：________年________月________日 一、实验目的 二、实验对象和装置 三．实验步骤

四、实验数据记录和整理 1、无阻尼单自由度自由振动系统实验测量： 2、有阻尼单自由度自由振动系统实验测量： 五、简答 1、力锤施加力的大小是否影响单自由度系统的振动频率和阻尼，为什么？ 2、实验过程中，力锤敲击质量块时应注意什么？

第三章 系统频率特性

第三章 系统频率特性 系统的时域分析是分析系统的直接方法，比较直观，但离开计算机仿真，分析高阶系统是困难的。系统频域分析是工程广为应用的系统分析和综合的间接方法。频率分析不仅可以了解系统频率特性，如截止频率、谐振频率等，而且可以间接了解系统时域特性，如快速性，稳定性等，为分析和设计系统提供更简便更可靠的方法。 本章首先阐明频率响应的特点，给出计算频率响应的方法，接着介绍Nyquist 图和Bode 图的绘制方法、系统的稳定裕度及系统时域性能指标计算。 3.1 频率响应和频率特性 3.1.1 一般概念 频率响应是指系统对正弦输入的稳态响应。考虑传递函数为G(s)的线性系统，若输入正弦信号 t X t x i i ωsin )(= (3.1-1) 根据微分方程解的理论，系统的稳态输出仍然为与输入信号同频率的正弦信号，只是其幅值和相位发生了变化。输出幅值正比于输入的幅值i X ，而且是输入正弦频率ω的函数。输出的相位与i X 无关，只与输入信号产生一个相位差?，且也是输入信号频率ω的函数。即线性系统的稳态输出为 )](sin[)()(00ω?ωω+=t X t x (3.1-2)

由此可知，输出信号与输入信号的幅值比是ω的函数，称为系统的幅频特性，记为)(ωA 。输出信号与输入信号相位差也是ω的函数，称为系统的相频特性，记为)(ω?。 幅频特性： )()()(0ωωωi X X A = (3.1-3) 相频特性： )()()(0ω?ω?ω?i -= (3.1-4) 频率特性是指系统在正弦信号作用下，稳态输出与输入之比对频率的关系特性，可表示为： )()()(0ωωωj X j X j G i = (3.1-5) 频率特性)(ωj G 是传递函数)(s G 的一种特殊形式。任何线性连续时间系统的频率特性都可由系统传递函数中的s 以ωj 代替而求得。 )(ωj G 有三种表示方法： )()()(ω?ωωj e A j G = (3.1-6) )()()(ωωωjV U j G += (3.1-7) )(sin )()cos()()(ω?ωωωωjA A j G += (3.1-8) 式中，实频特性： )(cos )()(ω?ωωA U = 虚频特性：

实验报告三_频率特性测量

实验报告 课程名称： 自动控制理论实验 指导老师： 吴越 成绩： 实验名称： 频率特性测量 实验类型： 同组学生姓名： 鲍婷婷 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得 一、实验目的 1. 掌握用超低频信号发生器和示波器测定系统或环节频率特性的方法； 2. 了解用TD4010型频率响应分析测试仪测定系统或环节的频率特性方法。 二、主要仪器设备 1.超低频信号发生器 2.电子模拟实验装置 3.超低频慢扫描示波器 三、实验步骤 1．测量微分积分环节的频率特性； (1)相频特性 相频特性的测试线路如图4-3-1所示，其中R 1=10k Ω、C 1=1uF 、R 2=2k Ω、C 2=50uF 。测量时，示波器的扫描旋钮指向X-Y 档。把超低频信号发生器的正弦信号同时送入被测系统和X 轴，被测系统的输出信号送入示波器Y 轴，此时在示波器上可得到一李沙育图形。 然后将椭圆移至示波器屏幕中间，椭圆与X 轴两交点的间的距离即为2X 0，将 Y 输入接地，此时得到的延X 轴光线长度 即为2X m ，因此求得θ=sin -1 (2X 0/2X m )，变化输入信号频率ω(rad/s)，即可得到一 组θ(ω)。测量时必须注意椭圆光点的转动方向，以判别相频特性是超前还是迟后。当系统或环节的相频特性是迟后时，光点为逆时针转动；反之超前时，光点为顺时针转动。测试时，ω取值应匀称，否则会影响曲线的准确度。 (2) 幅频特性：示波器选择停止扫描档，超低频信号发生的正弦信号同时送入X 轴和被测系统；被测环节的输出信号仍送入Y 轴；分别将X 通道和Y 通道接地，示波器上出现的两条光线对应的两条光线长度为2X m 、2Y m ，改变频率ω，则可得一组L(ω)。

频率响应介绍_频率响应概念

频率响应介绍_频率响应概念 频率响应是指将一个以恒电压输出的音频信号与系统相连接时，音箱产生的声压随频率的变化而发生增大或衰减、相位随频率而发生变化的现象，这种声压和相位与频率的相关联的变化关系称为频率响应。也是指在振幅允许的范围内音响系统能够重放的频率范围，以及在此范围内信号的变化量称为频率响应，也叫频率特性。在额定的频率范围内，输出电压幅度的最大值与最小值之比，以分贝数（dB）来表示其不均匀度。频率响应在电能质量概念中通常是指系统或计量传感器的阻抗随频率的变化。 频率响应确定方法分析法基于物理机理的理论计算方法，只适用于系统结构组成易于确定的情况。在系统的结构组成给定后，运用相应的物理定律，通过推导和计算即可定出系统的频率响应。分析的正确程度取决于对系统结构了解的精确程度。对于复杂系统，分析法的计算工作量很大。 实验法频率响应图册采用仪表直接量测的方法，可用于系统结构难以确定的情况。常用的实验方式是以正弦信号作为试验信号，在所考察的频率范围内选择若干个频率值，分别测量各个频率下输入和稳态输出正弦信号的振幅和相角值。输出与输入的振幅比值随频率的变化特性是幅频特性，输出与输入的相角差值随频率的变化特性是相频特性。 频率响应性能系统的过渡过程与频率响应有着确定的关系，可用数学方法来求出。但是除一阶和二阶系统外，这样做常需要很多时间，而且在很多情况下实际意义不大。常用的方法是根据频率响应的特征量来直接估计系统过渡过程的性能。频率响应的主要特征量有：增益裕量和相角裕量、谐振峰值和谐振频率、带宽和截止频率。 增益裕量和相角裕量它可提供控制系统是否稳定和具有多大稳定裕量的信息。 谐振峰值Mr和谐振频率rMr和r规定为幅频特性|G（j）|的最大值和相应的频率值。对于具有一对共轭复数主导极点（见根轨迹法）的高阶线性定常系统，当Mr值在（1.0～1.4）M0范围内时，可获得比较满意的过渡过程性能。其中M0是=0时频率响应的幅值。r的大小表征过渡过程的快速性：r值越大，系统在单位阶跃作用下输出响应的快速性越好。带宽和截止频率截止频率c规定为幅频特性|G（j）|达到0.7M0并继续下降时的临界频率。

语音信号处理实验报告

通信与信息工程学院 信息处理综合实验报告 通信与信息工程学院 二〇一八年

实验题目：语音信号分析与处理 一、实验内容 1. 设计内容 利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分 （1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。 （4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分 （5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。 （7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。

二、实验原理 1.设计原理分析 本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。 首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。 对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析 在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。 上述程序是在Matlab软件中画出语音信号的时域波形(图1),回放录入声音。从图中可以看出在时域环境下,信号呈现出3个不规则的信号峰值。然后对语音信号进行频谱分析,在Matlab中可以利用函数fft对信号行快速傅里叶变换,得到信号的频谱图如图1所示。 3．MATLAB中的窗函数 实际应用的窗函数，可分为以下主要类型： a) 幂窗--采用时间变量某种幂次的函数，如矩形、三角形、梯形或其它时间（t）的高次幂； b)三角函数窗--应用三角函数，即正弦或余弦函数等组合成复合函数，例如汉宁窗、海明窗等； c)指数窗--采用指数时间函数，如形式，例如高斯窗等。 d)平顶窗--平顶窗在频域时的表现就象它的名称一样有非常小的通带波动。

一阶网络频响特性测量

一阶网络频响特性测量 信号与系统实验报告 实验名称: 一阶网络频响特性测量 姓名: 姚敏 学号: 110404212 班级: 通信(2)班 时间: 2013.6.7 南京理工大学紫金学院电光系 一、实验目的 1、掌握一阶网络的构成方法; 2、掌握一阶网络的系统响应特性; 3、了解一阶网络频响特性图的测量方法; 二、实验基本原理 系统响应特性是指系统在正弦信号激励下，稳态响应随信号频率变化而变化的特性，称为系统的频率响应特性(frequency response)简称频响特性。 一阶系统是构成复杂系统的基本单元。学习一阶系统的特点有助于对一般系统特性的了解。一阶系统的系统函数为H(s)，表达式可以写成: 1H(s),k, k为一常数 (3-1) s，, 激励信号x(t)为: xtEt()sin(),, (3-2) m0 按照系统频响特性的定义可求得该一阶系统的稳态响应为: ytEHt()sin(),,，,(3-3) ssm000 j,0H(s)|,H(j,),H(j,)eH,H(j,)其中，。，，，，sj,,00000

可见，当改变系统输入信号的频率时，稳态响应的幅度和相位也随之而改变。 1,k,,,,0因果系统是稳定的要求:，不失一般性可设。该系统的频响,特性为: 1H(j,), (3-4) j,,，1 1从其频响函数中可以看出系统响应呈低通方式，其3dB带宽点。系统的频,响特性图如下图: , 图1 一阶网络频响特性图 一阶低通系统的单位冲击响应与单位阶跃响应如下图: 图2 一阶网络单位冲击响应与单位阶跃响应图 三、实验内容及结果 1、填写表1: 输入信号频输出信号幅度(mV) 相对幅度(dB) 相位差() ,,,21率f0 10Hz 1.96V -0.175 -181.4 1kHz 1.68V -1.514 -212.6 2kHz 1.27V -3.945 -229.7 3kHz 952mV -6.448 118.0 4kHz 760mV -8.404 114.0 5kHz 628mV - 10.061 108.9 6kHz 524mV -11.634 105.7 7kHz 456mV -12.841 -255.7 8kHz 400mV -13.979 102.2 9kHz 360mV -14.895 101.9 10kHz 324mV -15.810 100.7

信号与系统连续时间系统的频率响应

实验报告 实验名称：连续时间系统的频率响应

一、实验目的： 1 加深对连续时间系统频率响应理解； 2 掌握借助计算机计算任意连续时间系统频率响应的方法。 二、实验原理： 连续时间系统的频率响应可以直接通过所得表达式计算，也可以通过零极点 图通过用几何的方法来计算，而且通过零极点图可以迅速地判断系统的滤波特 性。 根据系统函数H(s)在s平面的零、极点分布可以绘制频响特性曲线，包括幅 频特性 H(jw) 曲线和相频特性?(w)曲线。这种方法的原理如下： 假定，系统函数H(s)的表达式为 当收敛域含虚轴时，取s = jw，也即在s平面中，s沿虚轴从- j∞移动到+ j∞时， 得到 容易看出，频率特性取决于零、极点的分布，即取决于Zj 、Pi 的位置，而式中K 是系数，对于频率特性的研究无关紧要。分母中任一因子(jw- Pi )相当于由极点 p 引向虚轴上某点 jw的一个矢量；分子中任一因子(jw-Zj)相当于由零点Zj引至虚轴上某点 jw的一个矢量。 在右图示意画出由零点Zj和极点 Pi 与 jw点连接构成的两个矢量，图中Nj、Mi 分别表示矢量的模，ψj、θi 表示矢量的辐角（矢量与正实轴的夹角，逆时针为正）。对于任意零点Zj 、极点Pi ，相应的复数因子（矢量）都可表示为： 于是，系统函数可以改写为

当ω延虚轴移动时，各复数因子（矢量）的模和辐角都随之改变，于是得出幅频特性曲线和相频特性曲线。这种方法称为s 平面几何分析。通过零极点图进行计算的方法是： 1 在S 平面上标出系统的零、极点位置； 2 选择S 平面的坐标原点为起始点，沿虚轴向上移动，计算此时各极点和零点与该点的膜模和夹角； 3 将所有零点的模相乘，再除以各极点的模，得到对应频率处的幅频特性的值； 4 将所有零点的幅角相加，减去各极点的幅角，得到对应频率处的相角。 三、实验内容 用 C 语言编制相应的计算程序进行计算，要求程序具有零极点输入模块， 可以手工输入不同数目的零极点。 计算频率从0～5频段的频谱，计算步长为0.1，分别计算上面两个系统的幅频特性和相频特性，将所得结果用表格列出，并画出相应的幅频特性曲线和相频特性曲线。 判断所给系统的滤波特性，对于带通滤波器，计算出 3dB 带宽的起始频点和结束频点；对于低通或高通滤波器，计算出3dB 带宽的截止或开始的频率。 四、画出系统一和系统二的零极点图 系统一 系统二 五、程序流程图和程序代码 程序流程图如下：