自动控制实验报告二-二阶系统阶跃响应相关数据及其图像

自动控制实验报告二-二阶系统阶跃响应相关数据及其图像
自动控制实验报告二-二阶系统阶跃响应相关数据及其图像

二阶系统阶跃响应报告书-数据及图像

一、实验目的

1.研究二阶系统的特征参数,阻尼比ζ和无阻尼自然频率ωn对系统动态性能的影响。定量分析ζ和ωn与最大超调量Mp和调节时间t S之间的关系。

2.进一步学习实验系统的使用方法

3.学会根据系统阶跃响应曲线确定传递函数。

二、实验仪器

1.EL-AT-III型自动控制系统实验箱一台

2.计算机一台

三、实验原理

1.模拟实验的基本原理:

控制系统模拟实验采用复合网络法来模拟各种典型环节,即利用运算放大器不同的输入网络和反馈网络模拟各种典型环节,然后按照给定系统的结构图将这些模拟环节连接起来,便得到了相应的模拟系统。再将输入信号加到模拟系统的输入端,并利用计算机等测量仪器,测量系统的输出,便可得到系统的动态响应曲线及性能指标。若改变系统的参数,还可进一步分析研究参数对系统性能的影响。

2. 域性能指标的测量方法:

超调量ó%:

1)启动计算机,在桌面双击图标 [自动控制实验系统] 运行软件。

2) 检查USB线是否连接好,在实验项目下拉框中选中实验,点击按钮,出

现参数设置对话框设置好参数,按确定按钮,此时如无警告对话框出现表示通信正常,如出现警告表示通信不正常,找出原因使通信正常后才可以继续进行实验。

3)连接被测量典型环节的模拟电路。电路的输入U1接A/D、D/A卡的DA1输

出,电路的输出U2接A/D、D/A卡的AD1输入,将两个积分电容连在模拟开关上。

检查无误后接通电源。

4)在实验项目的下拉列表中选择实验二[二阶系统阶跃响应] 。

5)鼠标单击按钮,弹出实验课题参数设置对话框。在参数设置对话框中设置相

应的实验参数后鼠标单击确认等待屏幕的显示区显示实验结果

6)利用软件上的游标测量响应曲线上的最大值和稳态值,代入下式算出超调量:

Y MAX - Y∞

ó%=——————×100%

Y∞

T P与T P:

利用软件的游标测量水平方向上从零到达最大值与从零到达95%稳态值所需的时间值,便可得到T P与T P。

四、实验内容

典型二阶系统的闭环传递函数为

ω2n

?(S)= (1)

s2+2ζωn s+ω2n

其中ζ和ωn对系统的动态品质有决定的影响。

构成图2-1典型二阶系统的模拟电路,并测量其阶跃响应:

图2-1 二阶系统模拟电路图

电路的结构图如图2-2:

图2-2 二阶系统结构图

系统闭环传递函数为

(2)式中 T=RC,K=R2/R1。

比较(1)、(2)二式,可得

ωn=1/T=1/RC

ζ=K/2=R2/2R1 (3)

由(3)式可知,改变比值R2/R1,可以改变二阶系统的阻尼比。改变RC值可以改变无阻尼自然频率ωn。

今取R1=200K,R2=100KΩ和200KΩ,可得实验所需的阻尼比。电阻R取100KΩ,电容C 分别取1μf和0.1μf,可得两个无阻尼自然频率ωn。

五、实验步骤

1.连接被测量典型环节的模拟电路。电路的输入U1接A/D、D/A卡的DA1输出,电路的

输出U2接A/D、D/A卡的AD1输入,将两个积分电容得两端连在模拟开关上。检查无误后接通电源。

2.启动计算机,在桌面双击图标 [自动控制实验系统] 运行软件。

3.测查USB线是否连接好,在实验项目下拉框中选中任实验,点击按钮,出

现参数设置对话框设置好参数按确定按钮,此时如无警告对话框出现表示通信

正常,如出现警告表示通信不正常,找出原因使通信正常后才可以继续进行实验。

4.在实验项目的下拉列表中选择实验二[二阶系统阶跃响应], 鼠标单击按钮,弹出

实验课题参数设置对话框。在参数设置对话框中设置相应的实验参数后鼠标单击确认等待屏幕的显示区显示实验结果

5.取ωn=10rad/s, 即令R=100KΩ,C=1μf;分别取ζ=0.5、1、2,即取R1=100KΩ,R2分

别等于100KΩ、200KΩ、400KΩ。输入阶跃信号,测量不同的ζ时系统的阶跃响应,并由显示的波形记录最大超调量Mp和调节时间Ts的数值和响应动态曲线,并与理论值比较。

6.取ζ=0.5。即电阻R2取R1=R2=100KΩ;ωn=100rad/s, 即取R=100KΩ,改变电路中的电

容C=0.1μf(注意:二个电容值同时改变)。输入阶跃信号测量系统阶跃响应,并由显示的波形记录最大超调量σp和调节时间Tn。

7.取R=100KΩ;改变电路中的电容C=1μf,R1=100KΩ,调节电阻R2=50KΩ。输入阶跃信号

测量系统阶跃响应,记录响应曲线,特别要记录Tp和σp的数值。

8.测量二阶系统的阶跃响应并记入表中

六、实验报告

1.画出二阶系统的模拟电路图,讨论典型二阶系统性能指标与ζ,ωn的关系。

延迟时间td:增大无阻尼自然振荡频率或减小阻尼比,都可以减少延迟时间。即,当阻尼比不变时,闭环极点距s平面的坐标原点越远,系统的延迟时间越短;而当无阻尼自然频率不变时,闭环极点距s平面虚轴越近,系统的延迟时间越短。

上升时间tr:要减小上升时间,当阻尼比一定时,需增大无阻尼自然振荡频率ωn ;当ω

一定时,需减小ζ.

n

峰值时间T p:T p=π/ω,和阻尼振荡频率成反比。

最大超调量Mp:Mp=

可知最大超调量仅和阻尼比有关系,与无阻尼自然振荡频率无关。随着阻尼比的增大,最大超调量单调的减小。

调节时间T s:=

可知调节时间和闭环极点的实部数值成反比,闭环极点的实部数值越大,即极点离虚轴的距离越远,系统的调节时间越短。

2.

由计算结果与理论值比较发现,测得结果中,超调量较理论值偏小,调节时间较理论值偏大,但很接近理论值,在误差允许范围内。产生偏差的原因可能是实际测量中有能量的损耗。

3.画出系统响应曲线,再由ts和Mp计算出传递函数,并与由模拟电路计算的传递函数相比较。

当C=1uf,ζ=0,ωn=10rad/s时,响应曲线:

计算得:

ζ=0.0157,ωn=12.62rad/s,

所以实际系统传递函数是

实际值与理论值存在误差,但是在可允许的范围内。

当C=1uf,ζ=0.25,ωn=10rad/s时,响应曲线是:

计算得:

ζ=0.262,ωn=12.47rad/s

理论传递函数是:Array

实际值与理论值存在误差,在误差允许范围内。

当C=1uf,ζ=0.5,ωn=10rad/s时,响应曲线是:

计算得:

ζ=0.510,ωn=11.00rad/s

实际传递函数是

理论传递函数是

实际值与理论值有误差,在误差范围内允许。

当C=1uf,ζ=1,ωn=10rad/s时,响应曲线:

计算得:

ζ=1,ωn=10rad/

实际值与理论值基本不存在差异。

n

计算得:

ζ=0.522,ωn=129.85rad/s

实际值与理论值存在误差,在允许范围内。

n

计算得:

ζ=1,ωn=100rad/s

实际值与理论值基本不存在误差。

实验总结:

通过这次实验,我了解了ωn,ζ对二阶系统的影响,更深入的了解二阶系统的响应特点。同时,实际值和理论值存在着一定的误差,可能是系统内部的能量损耗导致的,这在以后的自动控制实践中,是一定要考虑的。

二阶系统的阶跃响应及频率特性

实验二二阶系统的阶跃响应及频率特性 实验简介:通过本实验学生能够学习二阶系统的频率响应和幅频特性的测试方法,对实验装置和仪器的调试操作,具备对实验数据、结果的 处理及其与理论计算分析比较的能力。 适用课程:控制工程基础 实验目的:A 学习运算放大器在控制工程中的应用及传递函数的求取。 B 学习二阶系统阶跃响应曲线的实验测试方法。 C 研究二阶系统的两个重要参数ζ、ω n 对阶跃瞬态响应 指标的影响。 D 学习频率特性的实验测试方法。 E 掌握根据频率响应实验结果绘制Bode图的方法。 F 根据实验结果所绘制的Bode图,分析二阶系统的主要 动态特性(M P ,t s )。 面向专业:机械类 实验性质:综合性/必做 知 识 点:A《模拟电子技术》课程中运算放大器的相关知识; B《数字电子技术》课程中采样及采样定理的相关知识; C《机械工程控制基础》课程中,传递函数,时域响应, 频率响应三章的内容。 学 时 数:2 设备仪器:XMN-2自动控制原理学习机,CAE-98型微机接口卡,计算机辅助实验系统2.0软件,万用表。 材料消耗:运算放大器,电阻,电容,插接线。 要 求:实验前认真预习实验指导书的实验内容,完成下述项目, 做实验时交于指导教师检查并与实验报告一起记入实验成绩。 B推导图2所示积分放大器的输出输入时域关系和传递函数。

C 推导图3所示加法和积分放大器的输出输入时域关系(两输入单输出) 和S <1>.写出op1,op2,op9,0p6对应的微分方程组(4个方程)。 <2>.画出系统方框图。 <3>.用方框图化简或方程组联立消元的方法求取实验电路所示系统的 传递函数,写出求解过程。 和ζ。 <4>.求取该系统的ω n 实验地点:教一楼327室 实验照片:实验装置及仪器

自动控制原理实验报告

《自动控制原理》 实验报告 姓名: 学号: 专业: 班级: 时段: 成绩: 工学院自动化系

实验一 典型环节的MATLAB 仿真 一、实验目的 1.熟悉MATLAB 桌面和命令窗口,初步了解SIMULINK 功能模块的使用方法。 2.通过观察典型环节在单位阶跃信号作用下的动态特性,加深对各典型环节响应曲线的理解。 3.定性了解各参数变化对典型环节动态特性的影响。 二、实验原理 1.比例环节的传递函数为 K R K R R R Z Z s G 200,1002)(211 212==-=-=- = 其对应的模拟电路及SIMULINK 图形如图1-3所示。 三、实验内容 按下列各典型环节的传递函数,建立相应的SIMULINK 仿真模型,观察并记录其单位阶跃响应波形。 ① 比例环节1)(1=s G 和2)(1=s G ; ② 惯性环节11)(1+= s s G 和1 5.01 )(2+=s s G ③ 积分环节s s G 1)(1= ④ 微分环节s s G =)(1 ⑤ 比例+微分环节(PD )2)(1+=s s G 和1)(2+=s s G ⑥ 比例+积分环节(PI )s s G 11)(1+=和s s G 211)(2+= 四、实验结果及分析 图1-3 比例环节的模拟电路及SIMULINK 图形

① 仿真模型及波形图1)(1=s G 和2)(1=s G ② 仿真模型及波形图11)(1+= s s G 和1 5.01)(2+=s s G 11)(1+= s s G 1 5.01 )(2+=s s G ③ 积分环节s s G 1)(1= ④ 微分环节

二阶系统阶跃响应实验报告

实验一 二阶系统阶跃响应 一、实验目的 (1)研究二阶系统的两个重要参数:阻尼比ξ和无阻尼自振角频率ωn 对系统动 态性能的影响。 (2)学会根据模拟电路,确定系统传递函数。 二、实验内容 二阶系统模拟电路图如图2-1 所示。 系统特征方程为T 2s 2+KTs+1=0,其中T=RC ,K=R0/R1。根据二阶系统的标准 形式可知,ξ=K/2,通过调整K 可使ξ获得期望值。 三、预习要求 (1) 分别计算出T=0.5,ξ= 0.25,0.5,0.75 时,系统阶跃响应的超调量σP 和过渡过 程时间tS 。 ) 1( p 2 e ζζπσ--=, ζ T 3t s ≈

代入公式得: T=0.5,ξ= 0.25,σp=44.43% ,t s=6s; T=0.5,ξ= 0.5,σp=16.3% ,t s=3s; T=0.5,ξ= 0.75,σp=2.84% ,t s=2s; (2)分别计算出ξ= 0.25,T=0.2,0.5,1.0 时,系统阶跃响应的超调量σP 和过渡过程时间tS。 ξ= 0.25,T=0.2,σp=44.43% ,t s=2.4s; ξ= 0.25,T=0.5,σp=44.43% ,t s=6s; ξ= 0.25,T=1.0,σp=44.43% ,t s=12s; 四、实验步骤 (1)通过改变K,使ξ获得0,0.25,0.5,0.75,1.0 等值,在输入端加同样幅值的阶跃信号,观察过渡过程曲线,记下超调量σP 和过渡过程时间tS,将实验值和理论值进行比较。 (2)当ξ=0.25 时,令T=0.2 秒,0.5 秒,1.0 秒(T=RC,改变两个C),分别测出超调量σP 和过渡过程tS,比较三条阶跃响应曲线的异同。 五、实验数据记录与处理: 阶跃响应曲线图见后面附图。 原始数据记录: (1)T=0.5,通过改变R0的大小改变K值

数字图像处理实验报告

数字图像处理实验报告 实验一数字图像基本操作及灰度调整 一、实验目的 1)掌握读、写图像的基本方法。 2)掌握MATLAB语言中图像数据与信息的读取方法。 3)理解图像灰度变换处理在图像增强的作用。 4)掌握绘制灰度直方图的方法,理解灰度直方图的灰度变换及均衡化的方 法。 二、实验内容与要求 1.熟悉MATLAB语言中对图像数据读取,显示等基本函数 特别需要熟悉下列命令:熟悉imread()函数、imwrite()函数、size()函数、Subplot()函数、Figure()函数。 1)将MATLAB目录下work文件夹中的forest.tif图像文件读出.用到imread, imfinfo 等文件,观察一下图像数据,了解一下数字图像在MATLAB中的处理就是处理一个矩阵。将这个图像显示出来(用imshow)。尝试修改map颜色矩阵的值,再将图像显示出来,观察图像颜色的变化。 2)将MATLAB目录下work文件夹中的b747.jpg图像文件读出,用rgb2gray() 将其 转化为灰度图像,记为变量B。 2.图像灰度变换处理在图像增强的作用 读入不同情况的图像,请自己编程和调用Matlab函数用常用灰度变换函数对输入图像进行灰度变换,比较相应的处理效果。 3.绘制图像灰度直方图的方法,对图像进行均衡化处理 请自己编程和调用Matlab函数完成如下实验。 1)显示B的图像及灰度直方图,可以发现其灰度值集中在一段区域,用 imadjust函 数将它的灰度值调整到[0,1]之间,并观察调整后的图像与原图像的差别,调整后的灰

度直方图与原灰度直方图的区别。 2) 对B 进行直方图均衡化处理,试比较与源图的异同。 3) 对B 进行如图所示的分段线形变换处理,试比较与直方图均衡化处理的异同。 图1.1 分段线性变换函数 三、实验原理与算法分析 1. 灰度变换 灰度变换是图像增强的一种重要手段,它常用于改变图象的灰度范围及分布,是图象数字化及图象显示的重要工具。 1) 图像反转 灰度级范围为[0, L-1]的图像反转可由下式获得 r L s --=1 2) 对数运算:有时原图的动态范围太大,超出某些显示设备的允许动态范围, 如直接使用原图,则一部分细节可能丢失。解决的方法是对原图进行灰度压缩,如对数变换: s = c log(1 + r ),c 为常数,r ≥ 0 3) 幂次变换: 0,0,≥≥=γγc cr s 4) 对比拉伸:在实际应用中,为了突出图像中感兴趣的研究对象,常常要求 局部扩展拉伸某一范围的灰度值,或对不同范围的灰度值进行不同的拉伸处理,即分段线性拉伸: 其对应的数学表达式为:

图像处理实验报告

重庆交通大学 学生实验报告 实验课程名称数字图像处理 开课实验室数学实验室 学院理学院年级信息与计算科学专业 2 班学生姓名李伟凯学号631122020203 开课时间2014 至2015 学年第 1 学期

实验(一)图像处理基础 ?实验目的 学习Matlab软件的图像处理工具箱,掌握常用的一些图像处理命令;通过编程实现几种简单的图像增强算法,加强对图像增强的理解。 ?实验内容 题目A.打开Matlab软件帮助,学习了解Matlab中图像处理工具箱的基本功能;题目B.掌握以下常见图像处理函数的使用: imread( ) imageinfo( ) imwrite( ) imopen( ) imclose( ) imshow( ) impixel( ) imresize( ) imadjust( ) imnoise( ) imrotate( ) im2bw( ) rgb2gray( ) 题目C.编程实现对图像的线性灰度拉伸y = ax + b,函数形式为:imstrech(I, a, b); 题目D.编程实现对图像进行直方图均衡化处理,并将实验结果与Matab中imhist 命令结果比较。 三、实验结果 1).基本图像处理函数的使用: I=imread('rice.png'); se = strel('disk',1); I_opened = imopen(I,se); %对边缘进行平滑 subplot(1,2,1), imshow(I), title('原始图像') subplot(1,2,2), imshow(I_opened), title('平滑图像') 原始图像平滑图像

数字图像处理实验报告 (2)

目录 实验一:数字图像的基本处理操作 (2) 1.1:实验目的 (2) 1.2:实验任务和要求 (2) 1.3:实验步骤和结果 (2) 1.4:结果分析 (6) 实验二:图像的灰度变换和直方图变换 (7) 2.1:实验目的 (7) 2.2:实验任务和要求 (7) 2.3:实验步骤和结果 (7) 2.4:结果分析 (11) 实验三:图像的平滑处理 (11) 3.1:实验目的 (11) 3.2:实验任务和要求 (11) 3.3:实验步骤和结果 (12) 3.4:结果分析 (15) 实验四:图像的锐化处理 (16) 4.1:实验目的 (16) 4.2:实验任务和要求 (16) 4.3:实验步骤和结果 (16) 4.4:结果分析 (18)

实验一:数字图像的基本处理操作 1.1:实验目的 1、熟悉并掌握MATLAB、PHOTOSHOP等工具的使用; 2、实现图像的读取、显示、代数运算和简单变换。 3、熟悉及掌握图像的傅里叶变换原理及性质,实现图像的傅里叶变换。 1.2:实验任务和要求 1.读入一幅RGB图像,变换为灰度图像和二值图像,并在同一个窗口内分 成三个子窗口来分别显示RGB图像和灰度图像,注上文字标题。 2.对两幅不同图像执行加、减、乘、除操作,在同一个窗口内分成五个子窗口来分 别显示,注上文字标题。 3.对一幅图像进行平移,显示原始图像与处理后图像,分别对其进行傅里叶变换, 显示变换后结果,分析原图的傅里叶谱与平移后傅里叶频谱的对应关系。 4.对一幅图像进行旋转,显示原始图像与处理后图像,分别对其进行傅里 叶变换,显示变换后结果,分析原图的傅里叶谱与旋转后傅里叶频谱的 对应关系。 1.3:实验步骤和结果 1.对实验任务1的实现代码如下: a=imread('d:\tp.jpg'); i=rgb2gray(a); I=im2bw(a,0.5); subplot(1,3,1);imshow(a);title('原图像'); subplot(1,3,2);imshow(i);title('灰度图像'); subplot(1,3,3);imshow(I);title('二值图像'); subplot(1,3,1);imshow(a);title('原图像'); 结果如图1.1 所示:

自动控制实验报告1

东南大学自动控制实验室 实验报告 课程名称:自动控制原理 实验名称:闭环电压控制系统研究 院(系):仪器科学与工程专业:测控技术与仪器姓名:学号: 实验室:常州楼五楼实验组别:/ 同组人员:实验时间:2018/10/17 评定成绩:审阅教师: 实验三闭环电压控制系统研究

一、实验目的: (1)通过实例展示,认识自动控制系统的组成、功能。 (2)会正确实现闭环负反馈。 (3)通过开、闭环实验数据说明闭环控制效果。 二、实验原理: (1)利用各种实际物理装置(如电子装置、机械装置、化工装置等)在数学上的“相似性”,将各种实际物理装置从感兴趣的角度经过简化、并抽象成相同的数学形式。我们在设计控制系统时,不必研究每一种实际装置,而用几种“等价”的数学形式来表达、研究和设计。又由于人本身的自然属性,人对数学而言,不能直接感受它的自然物理属性,这给我们分析和设计带来了困难。所以,我们又用替代、模拟、仿真的形式把数学形式再变成“模拟实物”来研究。这样,就可以“秀才不出门,遍知天下事”。实际上,在后面的课程里,不同专业的学生将面对不同的实际物理对象,而“模拟实物”的实验方式可以做到举一反三,我们就是用下列“模拟实物”——电路系统,替代各种实际物理对象。 (2)自动控制的根本是闭环,尽管有的系统不能直接感受到它的闭环形式,如步进电机控制,专家系统等,从大局看,还是闭环。闭环控制可以带来想象不到的好处,本实验就是用开环和闭环在负载扰动下的实验数据,说明闭环控制效果。自动控制系统性能的优劣,其原因之一就是取决调节器的结构和算法的设计(本课程主要用串联调节、状态反馈),本实验为了简洁,采用单闭环、比例调节器K。通过实验证明:不同的K,对系性能产生不同的影响,以说明正确设计调节器算法的重要性。 (3)为了使实验有代表性,本实验采用三阶(高阶)系统。这样,当调节器K值过大时,控制系统会产生典型的现象——振荡。本实验也可以认为是一个真实的电压控制系统。 三、实验设备: THBDC-1实验平台 四、实验线路图: 五、实验步骤:

matlab图像处理实验报告

图像处理实验报告 姓名:陈琼暖 班级:07计科一班 学号:20070810104

目录: 实验一:灰度图像处理 (3) 实验二:灰度图像增强 (5) 实验三:二值图像处理 (8) 实验四:图像变换 (13) 大实验:车牌检测 (15)

实验一:灰度图像处理题目:直方图与灰度均衡 基本要求: (1) BMP灰度图像读取、显示、保存; (2)编程实现得出灰度图像的直方图; (3)实现灰度均衡算法. 实验过程: 1、BMP灰度图像读取、显示、保存; ?图像的读写与显示操作:用imread( )读取图像。 ?图像显示于屏幕:imshow( ) 。 ?

2、编程实现得出灰度图像的直方图; 3、实现灰度均衡算法; ?直方图均衡化可用histeq( )函数实现。 ?imhist(I) 显示直方图。直方图中bin的数目有图像的类型决定。如果I是个灰度图像,imhist将 使用默认值256个bins。如果I是一个二值图像,imhist使用两bins。 实验总结: Matlab 语言是一种简洁,可读性较强的高效率编程软件,通过运用图像处理工具箱中的有关函数,就可以对原图像进行简单的处理。 通过比较灰度原图和经均衡化后的图形可见图像变得清晰,均衡化后的直方图形状比原直方图的形状更理想。

实验二:灰度图像增强 题目:图像平滑与锐化 基本要求: (1)使用邻域平均法实现平滑运算; (2)使用中值滤波实现平滑运算; (3)使用拉普拉斯算子实现锐化运算. 实验过程: 1、 使用邻域平均法实现平滑运算; 步骤:对图像添加噪声,对带噪声的图像数据进行平滑处理; ? 对图像添加噪声 J = imnoise(I,type,parameters)

图像处理实验报告

实验报告 实验课程名称:数字图像处理 班级:学号:姓名: 注:1、每个实验中各项成绩按照10分制评定,每个实验成绩为两项总和20分。 2、平均成绩取三个实验平均成绩。 2016年 4 月18日

实验一 图像的二维离散傅立叶变换 一、实验目的 掌握图像的二维离散傅立叶变换以及性质 二、实验要求 1) 建立输入图像,在64?64的黑色图像矩阵的中心建立16?16的白色矩形图像点阵, 形成图像文件。对输入图像进行二维傅立叶变换,将原始图像及变换图像(三维、中心化)都显示于屏幕上。 2) 调整输入图像中白色矩形的位置,再进行变换,将原始图像及变换图像(三维、中 心化)都显示于屏幕上,比较变换结果。 3) 调整输入图像中白色矩形的尺寸(40?40,4?4),再进行变换,将原始图像及变 换图像(三维、中心化)都显示于屏幕上,比较变换结果。 三、实验仪器设备及软件 HP D538、MATLAB 四、实验原理 傅里叶变换作为分析数字图像的有利工具,因其可分离性、平移性、周期性和共轭对称性可以定量地方分析数字化系统,并且变换后的图像使得时间域和频域间的联系能够方便直观地解决许多问题。实验通过MATLAB 实验该项技能。 设),(y x f 是在空间域上等间隔采样得到的M ×N 的二维离散信号,x 和y 是离散实变量,u 和v 为离散频率变量,则二维离散傅里叶变换对一般地定义为 ∑∑ -=-=+-= 101 )],( 2ex p[),(1 ),(M x N y N yu M xu j y x f MN v u F π,1,0=u …,M-1;y=0,1,…N-1 ∑∑-=-=+=101 )],( 2ex p[),(),(M x N y N uy M ux j v u F y x f π ,1,0=x …,M-1;y=0,1,…N-1 在图像处理中,有事为了讨论上的方便,取M=N ,这样二维离散傅里叶变换对就定义为 ,]) (2ex p[),(1 ),(101 ∑∑ -=-=+- = N x N y N yu xu j y x f N v u F π 1,0,=v u …,N-1 ,]) (2ex p[ ),(1 ),(101 ∑∑-=-=+= N u N v N vy ux j v u F N y x f π 1,0,=y x ,…,N-1 其中,]/)(2exp[N yv xu j +-π是正变换核,]/)(2exp[N vy ux j +π是反变换核。将二维离散傅里叶变换的频谱的平方定义为),(y x f 的功率谱,记为 ),(),(|),(|),(222v u I v u R v u F v u P +== 功率谱反映了二维离散信号的能量在空间频率域上的分布情况。 五、实验步骤、程序及结果: 1、实验步骤: (1)、编写程序建立输入图像; (2)、对上述图像进行二维傅立叶变换,观察其频谱 (3)、改变输入图像中白框的位置,在进行二维傅里叶变换,观察频谱;

自动控制原理实验报告2

自动控制原理课程实验 2010-2011学年第一学期 02020801班 张驰2008300566

? 课本实验内容 6-26 热轧厂的主要工序是将炽热的钢坯轧成具有预定厚度和尺度的钢板,所得到的最终产品之一是宽为3300mm 、厚为180mm 的标准板材。他有两台主要的辊轧台:1号台与2号台。辊轧台上装有直径为508mm 的大型辊轧台,由4470km 大功率电机驱动,并通过大型液压缸来调节轧制宽度和力度。 热轧机的典型工作流程是:钢坯首先在熔炉中加热,加热后的钢坯通过1号台,被辊轧机轧制成具有预期宽度的钢坯,然后通过2号台,由辊轧机轧制成具有与其厚度的钢板,最后再由热整平设备加以整平成型。 热轧机系统控制的关键技术是通过调整热轧机的间隙来控制钢板的厚度。热轧机控制系统框图如下: 扰动)(s N )(s R (1)已知)54(/)(20++=s s s s s G ,而)(s G c 为具有两个相同实零点的PID 控制器。要求:选择PID 控制器的零点和增益,使闭环系统有两对对等的特征根; (2)考察(1)中得到的闭环系统,给出不考虑前置滤波器)(s G P 与配置适当)(s G P 时,系统的单位阶跃响应; (3)当)(s R =0,)(s N =1/s 时,计算系统对单位阶跃扰动的响应。 ? 求解过程 解:(1)已知 )54(/)(20++=s s s s s G )(s G P )(s G C )(0s G

选择 s z s K s G c /)()(2+= 当取K=4,Z=1.25时,有 s s s s s G c 4/25.610/)25.1(4)(2++=+= 系统开环传递函数 )54(/)25.1(4)()(2220+++=s s s s s G s G c 闭环传递函数:)25.61094/()5625.15.2(4))()(1/()()()(2 34200++++++=+=s s s s s s s G s G s G s G s c c φ (2) 当不考虑前置滤波器时,单位阶跃输入作用下的系统输出 )25.61094(/)5625.15.2(4)()()(2342++++++==s s s s s s s s R s s C φ 系统单位阶跃响应如图1中(1)中实线所示。 当考虑前置滤波器时,选 2)25.1/(5625.1)(+=s s G p 则系统在单位阶跃输入作用下的系统输出 )25.61094(/25.6)()()()(234++++==s s s s s s R s s G s C p φ 系统单位阶跃曲线如图1中(1)虚线所示。 (3)当)(s R =0,)(s N =1/s 时,扰动作用下的闭环传递函数 )25.61094/())()(1/()()(23400++++-=+-=s s s s s s G s G s G s c n φ 系统输出 )25.61094/(1)()()(2 34++++-==s s s s s N s s C n n φ 单位阶跃响应曲线如图1中(2)所示。 MATLAB 程序代码: MA TLAB 程序:exe626.m K=4;z=1.25; G0=tf(1,conv([1,0],[1,4,5])); Gc=tf(K*conv([1,z],[1,z]),[1,0]); Gp=tf(1.5625,conv([1,z],[1,z])); G1=feedback(Gc*G0,1); G2=series(Gp,G1); G3=-feedback(G0,Gc); t=0:0.01:10; [x,y]=step(G1,t);[x1,y1]=step(G2,t); figure(1);plot(t,x,'-',t,x1,':');grid

数字图像处理实验报告92184

数字图像处理试验报告 实验二:数字图像的空间滤波和频域滤波 姓名:XX学号:2XXXXXXX 实验日期:2017 年4 月26 日1.实验目的 1. 掌握图像滤波的基本定义及目的。 2. 理解空间域滤波的基本原理及方法。 3. 掌握进行图像的空域滤波的方法。 4. 掌握傅立叶变换及逆变换的基本原理方法。 5. 理解频域滤波的基本原理及方法。 6. 掌握进行图像的频域滤波的方法。 2.实验内容与要求 1. 平滑空间滤波: 1) 读出一幅图像,给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同 一图像窗口中。 2) 对加入噪声图像选用不同的平滑(低通)模板做运算,对比不同模板所形成的效果, 要求在同一窗口中显示。 3) 使用函数 imfilter 时,分别采用不同的填充方法(或边界选项,如零填 充、’replicate’、’symmetric’、’circular’)进行低通滤波,显示处理后 的图像。 4) 运用for 循环,将加有椒盐噪声的图像进行10 次,20 次均值滤波,查看其特点, 显示均值处理后的图像(提示:利用fspecial 函数的’average’类型生成均值滤 波器)。 5) 对加入椒盐噪声的图像分别采用均值滤波法,和中值滤波法对有噪声的图像做处理, 要求在同一窗口中显示结果。 6) 自己设计平滑空间滤波器,并将其对噪声图像进行处理,显示处理后的图像。 2. 锐化空间滤波 1) 读出一幅图像,采用3×3 的拉普拉斯算子 w = [ 1, 1, 1; 1 – 8 1; 1, 1, 1] 对其进行滤波。 2) 编写函数w = genlaplacian(n),自动产生任一奇数尺寸n 的拉普拉斯算子,如5 ×5的拉普拉斯算子 w = [ 1 1 1 1 1 1 1 1 1 1 1 1 -24 1 1 1 1 1 1 1 1 1 1 1 1] 3) 分别采用5×5,9×9,15×15和25×25大小的拉普拉斯算子对blurry_moon.tif

图像处理 实验报告

摘要: 图像处理,用计算机对图像进行分析,以达到所需结果的技术。又称影像处理。基本内容图像处理一般指数字图像处理。数字图像是指用数字摄像机、扫描仪等设备经过采样和数字化得到的一个大的二维数组,该数组的元素称为像素,其值为一整数,称为灰度值。图像处理技术的主要内容包括图像压缩,增强和复原,匹配、描述和识别3个部分。图像处理一般指数字图像处理。 数字图像处理的目的是改善图像的质量,它以人为对象,以改善人的视觉效果为目的。目前,图像处理演示系统应用领域广泛医学、军事、科研、商业等领域。因为数字图像处理技术易于实现非线性处理,处理程序和处理参数可变,故是一项通用性强,精度高,处理方法灵活,信息保存、传送可靠的图像处理技术。本图像处理演示系统以数字图像处理理论为基础,对某些常用功能进行界面化设计,便于初级用户的操作。 设计要求 可视化界面,采用多幅不同形式图像验证系统的正确性; 合理选择不同形式图像,反应各功能模块的效果及验证系统的正确性 对图像进行灰度级映射,对比分析变换前后的直方图变化; 1.课题目的与要求 目的: 基本功能:彩色图像转灰度图像 图像的几何空间变换:平移,旋转,剪切,缩放 图像的算术处理:加、减、乘 图像的灰度拉伸方法(包含参数设置); 直方图的统计和绘制;直方图均衡化和规定化; 要求: 1、熟悉图像点运算、代数运算、几何运算的基本定

义和常见方法; 2、掌握在MTLAB中对图像进行点运算、代数运算、几何运算的方法 3、掌握在MATLAB中进行插值的方法 4、运用MATLAB语言进行图像的插值缩放和插值旋转等 5、学会运用图像的灰度拉伸方法 6、学会运用图像的直方图设计和绘制;以及均衡化和规定化 7、进一步熟悉了解MATLAB语言的应用,将数字图像处理更好的应用于实际2.课题设计内容描述 1>彩色图像转化灰度图像: 大部分图像都是RGB格式。RGB是指红,绿,蓝三色。通常是每一色都是256个级。相当于过去摄影里提到了8级灰阶。 真彩色图像通常是就是指RGB。通常是三个8位,合起来是24位。不过每一个颜色并不一定是8位。比如有些显卡可以显示16位,或者是32位。所以就有16位真彩和32位真彩。 在一些特殊环境下需要将真彩色转换成灰度图像。 1单独处理每一个颜色分量。 2.处理图像的“灰度“,有时候又称为“高度”。边缘加强,平滑,去噪,加 锐度等。 3.当用黑白打印机打印照片时,通常也需要将彩色转成灰白,处理后再打印 4.摄影里,通过黑白照片体现“型体”与“线条”,“光线”。 2>图像的几何空间变化: 图像平移是将图像进行上下左右的等比例变化,不改变图像的特征,只改变位置。 图像比例缩放是指将给定的图像在x轴方向按比例缩放fx倍,在y轴按比例缩放fy倍,从而获得一幅新的图像。如果fx=fy,即在x轴方向和y轴方向缩放的比率相同,称这样的比例缩放为图像的全比例缩放。如果fx≠fy,图像的比例缩放会改变原始图象的像素间的相对位置,产生几何畸变。 旋转。一般图像的旋转是以图像的中心为原点,旋转一定的角度,也就是将图像上的所有像素都旋转一个相同的角度。旋转后图像的的大小一般会改变,即可以把转出显示区域的图像截去,或者扩大图像范围来显示所有的图像。图像的旋转变换也可以用矩阵变换来表示。

MATLAB下二阶系统的单位阶跃响应

二阶系统在不同参数下对单位阶跃信号的响应 一、二阶系统 所谓二阶系统就是其输入信号、输出信号的关系可用二阶微分方程来表征的系统。比如常见的RLC电路(图a)、单自由度振动系统等。 图a 图b 二阶系统传递函数的标准形式为 2 22 () 2 n n n H s s s ω ξωω = ++ 二、二阶系统的Bode图(nω=1) MATLAB程序为 >> clear >> num=[1]; >> den=[1 0.2 1]; >> bode(num,den); grid on hold on den=[1 0.4 1]; bode(num,den); >> den=[1 0.6 1]; >> bode(num,den); >> den=[1 0.8 1]; >> bode(num,den); >> den=[1 1.4 1]; >> bode(num,den); >> den=[1 2 1]; >> bode(num,den); >> legend('0.1','0.2','0.3','0.4','0.7','1.0')

运行结果为 三、二阶系统对单位阶跃信号的响应( =1) n MATLAB程序为 >> clear >> num=[1]; >> den=[1 0 1]; >> t=0:0.01:25; >> step(num,den,t) >> grid on >> hold on >> den=[1 0.2 1]; >> step(num,den,t) >> den=[1 0.4 1]; >> step(num,den,t) >> den=[1 0.6 1]; >> step(num,den,t) >> den=[1 0.8 1]; >> step(num,den,t) >> den=[1 1.0 1]; >> step(num,den,t)

数字图像处理实验报告实验三

中南大学 数字图像处理实验报告实验三数学形态学及其应用

实验三 数学形态学及其应用 一.实验目的 1.了解二值形态学的基本运算 2.掌握基本形态学运算的实现 3.了解形态操作的应用 二.实验基本原理 腐蚀和膨胀是数学形态学最基本的变换,数学形态学的应用几乎覆盖了图像处理的所有领域,给出利用数学形态学对二值图像处理的一些运算。 膨胀就是把连接成分的边界扩大一层的处理。而收缩则是把连接成分的边界点去掉从而缩小一层的处理。 二值形态学 I(x,y), T(i,j)为 0/1图像Θ 腐蚀:[]),(&),(),)((),(0,j i T j y i x I AND y x T I y x E m j i ++=Θ== 膨胀:[]),(&),(),)((),(0 ,j i T j y i x I OR y x T I y x D m j i ++=⊕== 灰度形态学T(i,j)可取10以外的值 腐蚀: []),(),(min ),)((),(1 ,0j i T j y i x I y x T I y x E m j i -++=Θ=-≤≤ 膨胀: []),(),(max ),)((),(1 ,0j i T j y i x I y x T I y x D m j i +++=⊕=-≤≤ 1.腐蚀Erosion: {}x B x B X x ?=Θ: 1B 删两边 2B 删右上 图5-1 剥去一层(皮) 2.膨胀Dilation: {}X B x B X x ↑⊕:= 1B 补两边 2B 补左下 图5-2 添上一层(漆) 3.开运算open :

B B X ⊕Θ=)(X B 4.闭close :∨ Θ⊕=B B X X B )( 5.HMT(Hit-Miss Transform:击中——击不中变换) 条件严格的模板匹配 ),(21T T T =模板由两部分组成。1T :物体,2T :背景。 {} C x x i X T X T X T X ??=?21, 图5-3 击不中变换示意图 性质: (1)φ=2T 时,1T X T X Θ=? (2))()()(21T X T X T X C Θ?Θ=? C T X T X )()(21Θ?Θ= )/()(21T X T X ΘΘ= 6.细化/粗化 (1)细化(Thin ) C T X X T X XoT )(/??=?= 去掉满足匹配条件的点。 图5-4 细化示意图 系统细化{}n B oB XoB T Xo ))(((21=, i B 是1-i B 旋转的结果(90?,180?,270?)共8种情况 适于细化的结构元素 1111000d d I = d d d L 10110 0= (2)粗化(Thick ) )(T X X T X ??=? 用(){}0,01=T (){}0,12=T 时,X X X T X =?=? X 21 1 1 2 3 T ? XoT X ? X X ?T X ΘT T ⊕

图形图像处理实验报告

第四次实验报告 实验课程:图像图像处理实验人:尹丽(200921020047) 实验时间:2012年4月19日实验地点:5-602 指导老师:夏倩老师成绩: 一、实验内容: ⑴图像的锐化:使用Sobel,Laplacian 算子分别对图像进行运算,观察并体会运算结果。 ⑵综合练习:对需要进行处理的图像分析,正确运用所学的知识,采用正确的步骤,对图像进行各类处理,以得到令人满意的图像效果。 二、实验目的: 学会用Matlab中的下列函数对输入图像按实验内容进行运算;感受各种不同的图像处理方法对最终图像效果的影响。(imfilter;fspecial;) 三、实验步骤:

1、仔细阅读Matlab 帮助文件中有关以上函数的使用说明,能充分理解其使用方法并能运用它们完成实验内容。 2、将Fig3.41(c).jpg 图像文件读入Matlab ,使用filter2函数分别采用不同的算子对其作锐化运算,显示运算前后的图像。 3、算子的输入可采用直接输入法。其中Sobel ,Laplacian ,也可用fspecial 函数产生。 4、各类算子如下: ???? ??????---121000121 ??????????-111181111 5、将Fig3.46(a).jpg 图像文件读入Matlab ,按照以下步骤对其进行处理: (1)用带对角线的Laplacian 对其处理,以增强边缘。 (2)用imadd 函数叠加原始图像。可以看出噪声增强了,应想法降低。 (3)获取Sobel 模板并用filter2对其进行5×5邻域平均,以减少噪声。 5(1)实验代码如图: 对角线Laplacian Sobel 垂直梯度

自动控制原理实验报告73809

-150-100 -50 50 实验一 典型环节的模拟研究及阶跃响应分析 1、比例环节 可知比例环节的传递函数为一个常数: 当Kp 分别为0.5,1,2时,输入幅值为1.84的正向阶跃信号,理论上依次输出幅值为0.92,1.84,3.68的反向阶跃信号。实验中,输出信号依次为幅值为0.94,1.88,3.70的反向阶跃信号, 相对误差分别为1.8%,2.2%,0.2%. 在误差允许范围内可认为实际输出满足理论值。 2、 积分环节 积分环节传递函数为: (1)T=0.1(0.033)时,C=1μf (0.33μf ),利用MATLAB ,模拟阶跃信号输入下的输出信号如图: T=0.1 T=0.033 与实验测得波形比较可知,实际与理论值较为吻合,理论上T=0.033时的波形斜率近似为T=0.1时的三倍,实际上为8/2.6=3.08,在误差允许范围内可认为满足理论条件。 3、 惯性环节 i f i o R R U U -=TS 1 CS R 1Z Z U U i i f i 0-=-=-=15 20

惯性环节传递函数为: K = R f /R 1,T = R f C, (1) 保持K = R f /R 1 = 1不变,观测T = 0.1秒,0.01秒(既R 1 = 100K,C = 1μf , 0.1μf )时的输出波形。利用matlab 仿真得到理论波形如下: T=0.1时 t s (5%)理论值为300ms,实际测得t s =400ms 相对误差为:(400-300)/300=33.3%,读数误差较大。 K 理论值为1,实验值2.12/2.28, 相对误差为(2.28-2.12)/2.28=7%与理论值 较为接近。 T=0.01时 t s (5%)理论值为30ms,实际测得t s =40ms 相对误差为:(40-30)/30=33.3% 由于ts 较小,所以读数时误差较大。 K 理论值为1,实验值2.12/2.28, 相对误差为(2.28-2.12)/2.28=7%与理论值较为接近 (2) 保持T = R f C = 0.1s 不变,分别观测K = 1,2时的输出波形。 K=1时波形即为(1)中T0.1时波形 K=2时,利用matlab 仿真得到如下结果: t s (5%)理论值为300ms,实际测得t s =400ms 相对误差为:(400-300)/300=33.3% 读数误差较大 K 理论值为2,实验值4.30/2.28, 1 TS K )s (R )s (C +-=

2. 实验二 二阶系统阶跃响应

实验二二阶系统阶跃响应 一、实验目的 1. 研究二阶系统的特征参数,阻尼比ζ和无阻尼自然频率ωn对系统动态性能的影响,定量分析ζ和ωn与最大超调量σp和调节时间ts之间的关系。 2. 进一步学习实验系统的使用。 3. 学会根据系统的阶跃响应曲线确定传递函数。 4. 学习用MATLAB仿真软件对实验内容中的电路进行仿真。 二、实验原理 典型二阶闭环系统的单位阶跃响应分为四种情况: 1)欠阻尼二阶系统 如图1所示,由稳态和瞬态两部分组成:稳态部分等于1,瞬态部分是振荡衰减的过程,振荡角频率为阻尼振荡角频率,其值由阻尼比ζ和自然振荡角频率ωn决定。 (1)性能指标: : 单位阶跃响应C(t)进人±5%(有时也取±2%)误差带,并且不再超出该误差带的调节时间t S 最小时间。 超调量σ% ;单位阶跃响应中最大超出量与稳态值之比。 单位阶跃响应C(t)超过稳态值达到第一个峰值所需要的时间。 峰值时间t P : 结构参数ξ:直接影响单位阶跃响应性能。 (2)平稳性:阻尼比ξ越小,平稳性越差 长,ξ过大时,系统响应迟钝,(3)快速性:ξ过小时因振荡强烈,衰减缓慢,调节时间t S 也长,快速性差。ξ=0.7调节时间最短,快速性最好。ξ=0.7时超调量σ%<5%,调节时间t S 平稳性也好,故称ξ=0.7为最佳阻尼比。 2)临界阻尼二阶系统(即ξ=1) 系统有两个相同的负实根,临界阻尼二阶系统单位阶跃响应是无超调的,无振荡单调上升的,不存在稳态误差。

3)无阻尼二阶系统(ξ=0时)此时系统有两个纯虚根。 4)过阻尼二阶系统(ξ>1)时 此时系统有两个不相等的负实根,过阻尼二阶系统的单位阶跃响应无振荡无超调无稳态误差,上升速度由小加大有一拐点。 三、实验内容 1. 搭建模拟电路 典型二阶系统的闭环传递函数为: 其中,ζ 和ωn对系统的动态品质有决定的影响。 搭建典型二阶系统的模拟电路,并测量其阶跃响应: 二阶系统模拟电路图其结构图为: 系统闭环传递函数为: 式中, T=RC,K=R2/R1。 比较上面二式,可得:ωn=1/T=1/RC ζ=K/2=R2/2R1。 2 2 2 2 ) ( ) ( ) ( n n n w s w s w s R s C S + + = = ξ φ

东北大学图像处理实验报告

计算机图像处理实验报告 哈哈哈哈哈哈实验台31 1.应用MATLAB语言编写显示一幅灰度图像、二值图像、索引图像及 彩色图像的程序,并进行相互之间的转换 1)彩色图像转换为灰度图像、索引图像、二值图像 A=imread('F:\colorful.jpg'); subplot(221);imshow(A);title('彩色图像'); I1=rgb2gray(A); subplot(222);imshow(I1);title('灰度图像'); [X1,map]=rgb2ind(A,256); subplot(223);imshow(X1);title('索引图像'); BW=im2bw(A); subplot(224);imshow(BW);title('二值图像'); 彩色图像灰度图像 索引图像二值图像

2)灰度图像转换为索引图像、二值图像 clear A=imread('F:\colorful.jpg'); B=rgb2gray(A); subplot(131);imshow(B);title('灰度图像'); [X2,map]=gray2ind(B,128); subplot(132);imshow(X2);title('索引图像'); BW2=im2bw(B); subplot(133);imshow(BW2);title('二值图像'); 灰度图像索引图像二值图像 3)索引图像转为灰度图像、二值图像、彩色图像 clear A=imread('F:\colorful.jpg'); [X,map]=rgb2ind(A,256); subplot(221);imshow(X);title('索引图像'); I3=ind2gray(X,map); subplot(222);imshow(I3);title('灰度图像'); BW3=im2bw(X,map,0.5); subplot(223);imshow(BW3);title('二值图像'); RGB=ind2rgb(X,map); subplot(24);imshow(RGB);title('还原彩色图像'); 索引图像灰度图像 二值图像还原彩色图像

自动控制实验报告一-控制系统的稳定性分析

实验一控制系统的稳定性分析 一、实验目的 1.观察系统的不稳定现象。 2.研究系统开环增益和时间常数对稳定性的影响。 二、实验仪器 1.自动控制系统实验箱一台 2.计算机一台 三、实验内容 系统模拟电路图如图 系统模拟电路图 其开环传递函数为: G(s)=10K/s(0.1s+1)(Ts+1) 式中K1=R3/R2,R2=100K,R3=0~500K;T=RC,R=100K,C=1f或C=0.1f两种情况。 四、实验步骤

1.连接被测量典型环节的模拟电路。电路的输入U1接A/D、D/A卡的DA1输出,电 路的输出U2接A/D、D/A卡的AD1输入,将纯积分电容两端连在模拟开关上。检查无误后接通电源。 2.启动计算机,在桌面双击图标[自动控制实验系统] 运行软件。 3.在实验项目的下拉列表中选择实验三[控制系统的稳定性分析] 5.取R3的值为50K,100K,200K,此时相应的K=10,K1=5,10,20。观察 不同R3值时显示区内的输出波形(既U2的波形),找到系统输出产生增幅振荡时相应的R3及K值。再把电阻R3由大至小变化,即R3=200k,100k,50k,观察不同R3值时显示区内的输出波形, 找出系统输出产生等幅振荡变化的R3及K 值,并观察U2的输出波形。 五、实验数据 1模拟电路图

2.画出系统增幅或减幅振荡的波形图。C=1uf时: R3=50K K=5: R3=100K K=10

R3=200K K=20: 等幅振荡:R3=220k:

增幅振荡:R3=220k: R3=260k:

C=0.1uf时:R3=50k:

相关文档
最新文档