理论力学matlab大作业

兰州理工大学2015级线性系统理论大作业

线性系统理论Matlab 实验报告

1、在造纸流程中，投料箱应该把纸浆流变成2cm的射流，并均匀喷洒在网状传送带上。为此，要精确控制喷射速度和传送速度之间的比例关系。投料箱内的压力是需要控制的主要变量，它决定了纸浆的喷射速度。投料箱内的总压力是纸浆液压和另外灌注的气压之和。由压力控制的投料箱是个耦合系统，因此，我们很难用手工方法保证纸张的质量。

在特定的工作点上，将投料箱线性化，可以得到下面的状态空间模型：

.??0.8?0.02??0.051? x???x??0.0010?u ?0.020???? y??x1,x2?

其中，系统的状态变量x1=液面高度，x2=压力，系统的控制变量u1=纸浆流量u2=气压阀门的开启量。在上述条件下，试设计合适的状态变量反馈控制器，使系统具有实特征根，且有一个根大于5

解：本题目是在已知状态空间描述的情况下要求设计一个状态反馈控制器，从而使得系统具有实数特征根，并要求要有一个根的模值要大于5，而特征根是正数时系统不

梧州学院

MATLAB与系统仿真

专 业：机械设计制造及其自动化班 级：学 号：姓 名：韦邮 箱：指导老师:

09机械5班 0901901012 靖

wei_jing1988@163.com 姚 金 2012年6月27日

1

一、春天是鲜花的季节，水仙花就是其中最迷人的代表，数学上有个水仙花数，它是这样定义的：“水仙花数”是指一个三位数，它的各位数字的立方和等于其本身，比如：153=1^3+5^3+3^3。要求编m文件输出所有存在的水仙花数。（本题10分） 解：

题意中“水仙花数”是指一个三位数说明100—999之间的数，在MATLAB编程M文件，运行可得。 程序如下： for i=100:999; a=floor(i/100); b=floor(i/10)-a*10; c=i-(a*100+b*10); if i==a^3+b^3+c^3; i end end 运行结果：

本程序源程序附件名为：

2

二、带饱和输出特性PID控制系统的被控对象为G0(s)=控制器为Gc(s)=Kp(1?1e?60s，其PID80s?1110.048?TDs)=K

基于MATLAB识别图片中的文字 知识点：图像识别 图像处理 计算机视觉

实验目的：以MATLAB为工具，识别图片中的文字。 问题描述：

交通是一个十分重要的问题，红绿灯上的摄像头可以辨别车主的身份。现在，我们是否能够用matlab设计一个图像识别的系统，通过车牌号的来识别相应的数字。 问题分析及模型建立： 第一步：确定车牌号的区域 对于一张图像

在matlab中，在处理图像元素时用（x,y,z）三坐标的形式表达的灰度集时，其中，x和y是空间坐标，z是f 在坐标（x,y）处的值，就可以表达图像在该坐标轴上的点。再将，将f(x,y)的数值简单地显示为一个矩阵，就可以定量地表达了一幅数值图像。

矩阵中的每个元素称为像素。

所以，如果我们想要截取一段车牌号的图像，我们就可以扫描图像每一个点的像素。然后我们进行恢复处理，算子边缘检测，腐蚀，闭操作，删除小对象得到

我们就可以利用RGB的值找到白色区域的边界，然后我们利用边界，截取我们的原图像就可以得到我们所需要的车牌号的区域了。

第二步：进行文字的分割

在进行分割之前，就需要将我们所获得的车牌号的区域的图像进行进一步的处理，突出我们需要的文字部分。 灰度处理

彩色的图片占用的空间比较大，处理会花费很长

基于MATLAB识别图片中的文字 知识点：图像识别 图像处理 计算机视觉

实验目的：以MATLAB为工具，识别图片中的文字。 问题描述：

交通是一个十分重要的问题，红绿灯上的摄像头可以辨别车主的身份。现在，我们是否能够用matlab设计一个图像识别的系统，通过车牌号的来识别相应的数字。 问题分析及模型建立： 第一步：确定车牌号的区域 对于一张图像

在matlab中，在处理图像元素时用（x,y,z）三坐标的形式表达的灰度集时，其中，x和y是空间坐标，z是f 在坐标（x,y）处的值，就可以表达图像在该坐标轴上的点。再将，将f(x,y)的数值简单地显示为一个矩阵，就可以定量地表达了一幅数值图像。

矩阵中的每个元素称为像素。

所以，如果我们想要截取一段车牌号的图像，我们就可以扫描图像每一个点的像素。然后我们进行恢复处理，算子边缘检测，腐蚀，闭操作，删除小对象得到

我们就可以利用RGB的值找到白色区域的边界，然后我们利用边界，截取我们的原图像就可以得到我们所需要的车牌号的区域了。

第二步：进行文字的分割

在进行分割之前，就需要将我们所获得的车牌号的区域的图像进行进一步的处理，突出我们需要的文字部分。 灰度处理

彩色的图片占用的空间比较大，处理会花费很长

基于MATLAB的通信系统仿真

专业：电子与通信工程

学号： 姓名： 成绩：

日期：

目录

基于Matlab的通信系统仿真 ······································································································ 3 系统综述 ······················································································································ 3

结构框图 ······················································································ 错误！未定义书签。 系统实现 ··········································································

2008 级

《MATAB语言与控制系统仿真》课程

大作业

姓 名 学 号 所在院系 班 级 日 期 同组人员

作业评分 评阅人

设计报告评分表

项 目 方案论证9分 设计过程应包括的主要内容或考核要点 性能指标分析；控制方法及实现方案 常见问题 未分析指标 没有说明方案理由 未进行方案比较 有分析指标，未用程序分析 有方案说明，未用程序说明 有设计过程，未运用程序 有简单程序应用，无算法 有程序，没有程序说明 未说明设计理由 扣分 -5 -1 -1 -1 -1 -20 -5 -5 -2 -9 -3 -2 -2 -2 -2 -2 -1 -1 -1 -1 -1 -9 -4 得分 控制器设计与参数计算 基本35分 部分 结果63分 分析10分 格式规范9分 第1没有误差分析 对设计

计算机实习

姓名：

班级： 学号： 分数：

题目一：X： 1 2 3 4 5 6 7 8 9；y：9 7 6 3 -1 2 5 7 20。对XY进行曲线拟合。 分析：

对于已给数据先用matlab画出散点图，根据散点图假想xy之间的函数关系，在这里就不写出散点图的程序，由散点图可预测为三次函数。所以调用matlab Ployfit函数（三次） 拟合所写程序： x=[1 2 3 4 5 6 7 8 9]; y=[9 7 6 3 -1 2 5 7 20]; P=polyfit(x,y,3); xi=0:.2:10; yi=polyval(P,xi); plot(xi,yi,x,y,'r*'); 总结：

根据上述程序我们能得到有关XY的拟合曲线，也可进一步求出函数中的相关系数。

? 拟合的曲线呢？

题目二：编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 分析：

编写一个M文件写出此算法，在解决问题时直接输入数字，然后调用此M文件即可解决问题。

M文件程序： function c=collatz(n) % collatz

《MATLAB》期末大作业

学院土木工程与建筑学院

专业

班级

姓名

指导教师李琳

2018 年 5 月16 日

明

作业内容题目2:问题描述：在[0 2π]范围内绘制二维曲线图y=cos(5x)*sin(x)

（1）问题分析

这是一个二维绘图问题，先划定x的范围与间距，再列出y的表达式，利用plot函数绘制二维曲线。

（2）软件说明及源代码

>> x = 0:pi/10:2*pi;

>>y = cos(5*x).*sin(x);

>>plot(x,y)

（3）实验结果

题目4:问题描述：创建符号函数并求解，要求写出步骤和运行结果

(1)创建符号函数f=ax2+bx+c

(2)求f=0的解

（1）问题分析

这是一个符号函数显示以及符号函数的求解问题，第一问先定义常量与变量，在写出f表达式，利用pretty函数显示f。第二问利用solve函数求解f=0时的解。

（2）软件说明及源代码

第一问

>> syms a b c x;

>> f=a*x^2+b*x+c;

>> pretty(f)

第二问

>>syms a b c x;

>>f=a*x^2+b*x+c;

>> solve(f)

（3）实

目 录

1 前言： .......................................................... 1 2 利用自然伽马测井数据简单分析地层岩性 ............................ 1 2．1 估算地层泥质含量 ........................................... 1 2．2 估算地层biot系数 .......................................... 2 3 利用测井数据计算分析地层的弹性模量、泊松比。 .................... 2 3．1纵、横声波速度.............................................. 2 3．2岩石的动态弹性参数确定...................................... 3 3．3岩石的动态弹性参数与静态弹性参数的转化...................... 3 4 利用测井数据计算粘聚力、内摩擦角与地层抗拉强度的连续剖面 ........ 3 4．1岩石单轴抗压强度的确定.....................

一、作业题目

结合所学的《岩石力学》课程及相关知识，利用给出的测井数据，对地层力学参数、孔隙压力、地应力、地层坍塌压力与破裂压力进行分析计算，分析储层出砂可能性，作出地层力学参数、地层主应力、地层坍塌、破裂压力剖面，分析井壁坍塌原因；研究储层段的出砂可能性，形成结课作业报告。

1. 已知条件

1) A井测井数据，分析孔隙压力，建议采用Eaton法，Eaton指数3.0。. 2) B井对地层力学参数地应力、地层坍塌压力与破裂压力进行分析计算，

结合实用泥浆密度分析井壁失稳原因，并提出合理化建议，分析储层出砂可能性，推荐合理的完井方式。 3) 已知：

? 地层孔隙压力当量密度为1.03g/cm3，

? 地层岩性：3000米以前为典型砂泥岩地层，3000米为砂泥岩，夹薄

层煤。

? 储层段：2800-3000米砂岩层。

? 地应力实测值：在3690m处实测水平最大主应力大约70MPa，水平

最小主应力大约63MPa；

? 测井过程中钻井液密度为1.25g/cm3； 4) 岩心抗压强度试验结果：

井深（m） 3030 围压值（MPa） 0 25 0 30 破坏强度（MPa） 16.1 106.5 22.6 127.1 3110 2. 要求

1） 编写程序读