华南农业大学期末考试试卷（A 卷） 2012-2013学年第 二 学期 考试科目：数学建模 考试类型：（闭卷）考试 考试时间： 120 分钟 学号 姓名 年级专业 一篮白菜从河岸一边带到河岸对面，由于船的限制，一次只能带 一样东西过河，绝不能在无人看守的情况下将狼和羊放在一起；羊和白菜放在一起，怎样才能将它们安全的带到河对岸去? 建立多步决策模型,将人、狼、羊、白菜分别记为i = 1，2，3，4，当i 在此岸时记x i = 1，否则为0；此岸的状态下用s =（x 1，x 2，x 3，x 4）表示。该问题中决策为乘船方案，记为d = (u 1, u 2, u 3, u 4)，当i 在船上时记u i = 1，否则记u i = 0。 (1) 写出该问题的所有允许状态集合；（3分） (2) 写出该问题的所有允许决策集合；（3分） (3) 写出该问题的状态转移率。（3分） (4) 利用图解法给出渡河方案. （3分） 解：(1) S={(1,1,1,1), (1,1,1,0), (1,1,0,1), (1,0,1,1), (1,0,1,0)} 及他们的5个反状（3分） (2) D = {(1,1,0,0), (1,0,1,0), (1,0,0,1), (1,0,0,0)} （6分） (3) s k+1 = s k + (-1) k d k （9分） (4)方法：人先带羊，然后回来，带狼过河，然后把羊带回来，放下羊，带白菜过去，然后再回来把羊带过去。 ?或: 人先带羊过河，然后自己回来，带白菜过去，放下白菜，带着羊回来，然后放下羊，把狼带过去，最后再回转来，带羊过去。 （12分） 1、 二、(满分12分) 在举重比赛中，运动员在高度和体重方面差别很大，请就下面两种假设，建立一个举重能力和体重之间关系的模型： （1） 假设肌肉的强度和其横截面的面积成比例。6分 （2） 假定体重中有一部分是与成年人的尺寸无关，请给出一个改进模型。6分 解：设体重w （千克）与举重成绩y （千克） （1） 由于肌肉强度(I)与其横截面积(S)成比例，所以 y ?I ?S 设h 为个人身高，又横截面积正比于身高的平方，则S ? h 2 再体重正比于身高的三次方，则w ? h 3 （6分） ( 12分) 14分) 某学校规定，运筹学专业的学生毕业时必须至少学
温度控制系统研究背景与现状 1 研究背景 (1) 2 国内外现状 (1) 定值开关温度控制法 (1) PID线性温度控制法 (2) 智能温度控制法 (3) 国内外实例 (4) 1 研究背景 温度是生活及生产中最基本的物理量，它表征的是物体的冷热程度。自然界中任何物理、化学过程都紧密地与温度相联系。在很多生产过程中，温度的测量和控制都直接和安全生产、提高生产效率、保证产品质量、节约能源等重大技术经济指标相联系。自18世纪工业革命以来，工业过程离不开温度控制。温度控制广泛应用于社会生活的各个领域，如家电、汽车、材料、电力电子等。温度控制的精度以及不同控制对象的控制方法选择都起着至关重要的作用，温度是锅炉生产质量的重要指标之一，也是保证锅炉设备安全的重要参数。同时，温度是影响锅炉传热过程和设备效率的主要因素。基于此，运用反馈控制理论对锅炉进行温度控制，满足了工业生产的需求，提高了生产力。 2 国内外现状 温度控制技术按照控制目标的不同可分为两类：动态温度跟踪与恒值温度控制。动态温度跟踪实现的控制目标是使被控对象的温度值按预先设定好的曲线进行变化。在工业生产中很多场合需要实现这一控制目标，如在发酵过程控制，化工生产中的化学反应温度控制，冶金工厂中燃烧炉中的温度控制等。恒值温度控制的目的是使被控对象的温度恒定在某一数值上，且要求其波动幅度（即稳态误差）不能超过某一给定值。从工业温度控制器的发展过程来看，温度控制技术大致可分以下几种： 定值开关温度控制法 所谓定值开关控温法，就是通过硬件电路或软件计算判别当前温度值与设定目标温度值之间的关系，进而对系统加热源（或冷却装置）进行通断控制。若当前温度值比设定温度值高，则关断加热器，或者开动制冷装置；若当前温度值比设定温度值低，则开启加热器并同时关断制冷器。这种开关控温方法比较简单，在没有计算机参与的情况下，用很简单的模拟电路就能够实现。目前，采用这种控制方法的温度控制器在我国许多工厂的老式工业电炉中仍被使用。由于这种控制方式是当系统温度上升至设定点时关断电源，当系统温度下降至设定点时开通
上机练习题一 班级： 姓名： 学号： 1.建立起始值=3,增量值=5.5,终止值=44的一维数组x 答案： x=(3:5.5:44) 2.写出计算 Sin(30o )的程序语句. 答案： sin(pi*30/180) 或 sin(pi/6) 3.矩阵??????????=187624323A ,矩阵???? ??????=333222111B ;分别求出B A ?及A 与B 中对应元素之间的乘积. 答案：A = [3,2,3; 4,2,6; 7,8,1] B = [1,1,1; 2,2,2; 3,3,3] A*B ；A.*B 4计算行列式的值1 876243 23=A 。答案：det(A) 5对矩阵 ???? ??????=187624323A 进行下述操作。 (1)求秩。答案：rank(A) (2)求转置。答案：A' (3) 对矩阵求逆，求伪逆。答案：inv(A) ，pinv(A) (4) 左右反转,上下反转。答案：fliplr(A)，flipud(A) (5) 求矩阵的特征值. 答案：[u,v]=eig(A) (6) 取出上三角和下三角. 答案：triu(A) tril(A) (7)以A 为分块作一个3行2列的分块矩阵。答案：repmat(a) 6 计算矩阵??????????897473535与???? ??????638976242之和。 >> a=[5 3 5;3 7 4;7 9 8]; >> b=[2 4 2;6 7 9;8 3 6]; >> a+b 7 计算??????=572396a 与?? ????=864142b 的数组乘积。 >> a=[6 9 3;2 7 5]; >> b=[2 4 1;4 6 8];
数学建模部分课后习题解答 中国地质大学 能源学院 华文静 1.在稳定的椅子问题中，如设椅子的四脚连线呈长方形，结论如何？ 解： 模型假设 （1） 椅子四条腿一样长，椅脚与地面接触处视为一点，四脚的连线呈长方形 （2） 地面高度是连续变化的，沿任何方向都不会出现间断（没有像台阶那样的情况）， 即从数学角度来看，地面是连续曲面。这个假设相当于给出了椅子能放稳的必要条件 （3） 椅子在任何位置至少有三只脚同时着地。为了保证这一点，要求对于椅脚的间 距和椅腿的长度而言，地面是相对平坦的。因为在地面上椅脚间距和椅腿长度的尺寸大小相当的范围内，如果出现深沟或凸峰（即使是连续变化的），此时三只脚是无法同时着地的。 模型建立 在上述假设下，解决问题的关键在于选择合适的变量，把椅子四只脚同时着地表示出来。首先，引入合适的变量来表示椅子位置的挪动。生活经验告诉我们，要把椅子通过挪动放稳，通常有拖动或转动椅子两种办法，也就是数学上所说的平移与旋转变换。然而，平移椅子后问题的条件没有发生本质变化，所以用平移的办法是不能解决问题的。于是可尝试将椅子就地旋转，并试图在旋转过程中找到一种椅子能放稳的情形。 注意到椅脚连线呈长方形，长方形是中心对称图形，绕它的对称中心旋转180度后，椅子仍在原地。把长方形绕它的对称中心旋转，这可以表示椅子位置的改变。于是，旋转角度θ这一变量就表示了椅子的位置。为此，在平面上建立直角坐标系来解决问题。 设椅脚连线为长方形ABCD,以对角线AC 所在的直线为x 轴，对称中心O 为原点，建立平面直角坐标系。椅子绕O 点沿逆时针方向旋转角度θ后，长方形ABCD 转至A1B1C1D1的位置，这样就可以用旋转角）0（πθθ≤≤表示出椅子绕点O 旋转θ后的位置。 其次，把椅脚是否着地用数学形式表示出来。当椅脚与地面的竖直距离为零时，椅脚就着地了，而当这个距离大于零时，椅脚不着地。由于椅子在不同的位置是θ的函数，因此，椅脚与地面的竖直距离也是θ的函数。 由于椅子有四只脚，因而椅脚与地面的竖直距离有四个，它们都是θ的函数，而由假设（3）可知，椅子在任何位置至少有三只脚同时着地，即这四个函数对于任意的θ，其函数值至少有三个同时为0。因此，只需引入两个距离函数即可。考虑到长方形ABCD 是对称中心图形，绕其对称中心O 沿逆时针方向旋转180度后，长方形位置不变，但A,C 和B,D 对换了。因此，记A ，B 两脚与地面竖直距离之和为)(θf ，C,D 两脚之和为 )(θg ，其中[]πθ，0∈，使得)()(00θθg f =成立。 模型求解 如果0)0()0(== g f ，那么结论成立。
过程控制工程课程设计 课题名称空调温度控制系统的建模与仿真 学院 专业 班级 学生 学号 时间 6 月13日至 6月19日 指导教师(签字) 2011 年 6 月 19 日
目录 第一章设计题目及要求 (1) 1.1设计背景 (1) 1.2设计任务 (1) 1.3主要参数 (2) 1.3.1恒温室： (2) 1.3.2热水加热器ⅠSR、ⅡSR： (2) 1.3.3电动调节阀： (2) 1.3.4温度测量环节： (2) 1.3.5调节器： (2) 第二章空调温度控制系统的数学模型 (3) 2.1恒温室的微分方程 (3) 2.1.1微分方程的列写 (3) 2.1.2 增量微分方程式的列写 (5) 2.2 热水加热器对象的微分方程 (5) 2.3敏感元件及变送器的特性 (6) 2.3.1敏感元件的微分方程 (7) 2.3.2变送器的特性 (7) 2．3．3敏感元件及变送器特性 (8) 2.4 执行器的特性 (8) 第三章控制系统方案设计 (9) 3.1系统分析 (9) 3.2 单回路控制系统设计 (10) 3.2.1单回路控制系统原理 (10) 3.2.2单回路系统框图 (10) 3.3串级控制系统的设计 (11) 3.3.1串级控制系统原理 (11) 3.3.2串级控制系统框图 (12) 第四章单回路系统调节器参数整定 (13) 5.1.1、PI控制仿真 (16) 5.1.2 PID控制仿真 (17) 5.1.3、PI与PID控制方式比较 (17) 第六章设计小结 (18) 参考文献 (18)
第一章设计题目及要求 1.1设计背景 设计背景为一个集中式空调系统的冬季温度控制环节，简化系统图如附图所示。 系统由空调房间、送风道、送风机、加热设备及调节阀门等组成。为了节约能量，利用一部分室循环风与室外新风混合，二者的比例由空调工艺决定，并假定在整个冬季保持不变。用两个蒸汽盘管加热器1SR、2SR对混合后的空气进行加热，加热后的空气通过送风机送入空调房间。本设计中假设送风量保持不变。 1.2设计任务 设计主要任务是根据所选定的控制方案，建立起控制系统的数学模型，然后用MATLAB对控制系统进行仿真，通过对仿真结果的分析、比较，总结不同的控
实验一 控制系统的数学模型 一 实验目的 1、学习用MATLAB 创建各种控制系统模型。 2、掌握传递函数模型、零-极点增益模型以及连续系统模型与离散系统模型之间的转化，模型的简化。 二 相关理论 1传递函数描述 (1)连续系统的传递函数模型 连续系统的传递函数如下： ? 对线性定常系统，式中s 的系数均为常数，且a1不等于零，这时系统在MATLAB 中 可以方便地由分子和分母系数构成的两个向量唯一地确定出来，这两个向量分别用num 和den 表示。 num=[b1,b2,…,bm,bm+1] den=[a1,a2,…,an,an+1] 注意：它们都是按s 的降幂进行排列的。 tf （）函数可以表示传递函数模型：G=tf(num, den) 举例： num=[12,24,0,20];den=[2 4 6 2 2]; G=tf(num, den) (2)零极点增益模型 ? 零极点模型实际上是传递函数模型的另一种表现形式，其原理是分别对原系统传递 函数的分子、分母进行分解因式处理，以获得系统的零点和极点的表示形式。 K 为系统增益，zi 为零点，pj 为极点 在MATLAB 中零极点增益模型用[z,p,K]矢量组表示。即： z=[z1,z2,…,zm] p=[p1,p2,...,pn] K=[k] zpk （）函数可以表示零极点增益模型：G=zpk(z,p,k) (3)部分分式展开 ? 控制系统常用到并联系统，这时就要对系统函数进行分解，使其表现为一些基本控 制单元的和的形式。 ? 函数[r,p,k]=residue(b,a)对两个多项式的比进行部分展开，以及把传函分解为微 分单元的形式。 ? 向量b 和a 是按s 的降幂排列的多项式系数。部分分式展开后，余数返回到向量r ， 极点返回到列向量p ，常数项返回到k 。 ? [b,a]=residue(r,p,k)可以将部分分式转化为多项式比p(s)/q(s)。 11 211121......)()()(+-+-++++++++==n n n n m n m m a s a s a s a b s b s b s b s R s C s G ))...()(())...()(()(2121n m p s p s p s z s z s z s K s G ------=22642202412)(23423++++++=s s s s s s s G
第二章自动控制系统的数学模型 教学目的： （1）建立动态模拟的概念，能编写系统的微分方程。 （2）掌握传递函数的概念及求法。 （3）通过本课学习掌握电路或系统动态结构图的求法，并能应用各环节的传递函数，求系统的动态结构图。 （4）通过本课学习掌握电路或自动控制系统动态结构图的求法，并对系统结构图进行变换。 （5）掌握信号流图的概念，会用梅逊公式求系统闭环传递函数。 （6）通过本次课学习，使学生加深对以前所学的知识的理解，培养学生分析问题的能力 教学要求： （1）正确理解数学模型的特点； （2）了解动态微分方程建立的一般步骤和方法； （3）牢固掌握传递函数的定义和性质，掌握典型环节及传递函数； （4）掌握系统结构图的建立、等效变换及其系统开环、闭环传递函数的求取，并对重要的传递函数如：控制输入下的闭环传递函数、扰动输入 下的闭环传递函数、误差传递函数，能够熟练的掌握； （5）掌握运用梅逊公式求闭环传递函数的方法； （6）掌握结构图和信号流图的定义和组成方法，熟练掌握等效变换代数法则，简化图形结构，掌握从其它不同形式的数学模型求取系统传递函 数的方法。 教学重点： 有源网络和无源网络微分方程的编写；有源网络和无源网络求传递函数；传递函数的概念及求法；由各环节的传递函数，求系统的动态结构图；由各环节的传递函数对系统的动态结构图进行变换；梅逊增益公式的应用。 教学难点：举典型例题说明微分方程建立的方法；求高阶系统响应；求复杂系统的动态结构图；对复杂系统的动态结构图进行变换；求第K条前向通道特记式 的余子式 。 k 教学方法：讲授 本章学时：10学时 主要内容： 2.0 引言 2.1 动态微分方程的建立 2.2 线性系统的传递函数 2.3 典型环节及其传递函数 2.4系统的结构图 2.5 信号流图及梅逊公式
《数学模型》作业解答 第七章（ 2008 年 12 月 4 日） 1．对于节蛛网模型讨论下列问题： （ 1）因为一个时段上市的商品不能立即售完，其数量也会影响到下一时段的价格，所以第 k 1时段的价格y k 1由第k 1 和第 k 时段的数量x k 1和x k决定，如果仍设x k 1仍只取
决于 y k ，给出稳定平衡的条件，并与节的结果进行比较 . （ 2）若除了 y k 1 由 x k 1 和 x k 决定之外， x k 1 也由前两个时段的价格 析稳定平衡的条件是否还会放宽 . 解：（ 1）由题设条件可得需求函数、供应函数分别为： y k 1 f x k 1 x k ) ( 2 x k 1 h( y k ) 在 P 0 (x 0 , y 0 ) 点附近用直线来近似曲线 f , h ，得到 y k 1 y 0 ( x k 1 x k x 0 ), 2 x k 1 x 0 ( y k y 0 ) , 由（ 2）得 x k 2 x 0 ( y k 1 y 0 ) （ 1）代入（ 3）得 x k 2 x 0 ( x k 1x k x 0 ) 2 2x k 2 x k 1 x k 2x 0 2 x 0 对应齐次方程的特征方程为 2 2 ( ) 2 8 特征根为 1, 2 4 y k 和 y k 1 确定 . 试分 (1) ( 2) (3) 当 8 时，则有特征根在单位圆外，设 8 ，则
1,2 ( ) 2 ( ) 2 8 42 2 4 1,2 1 2 即平衡稳定的条件为 2与 P 207 的结果一致 . （ 2）此时需求函数、供应函数在 P 0 (x 0 , y 0 ) 处附近的直线近似表达式分别为： y k 1 y 0 ( x k 1 x k x 0 ), ( 4) 2 x k 1 x 0 ( y k y k 1 y 0 ) , ( 5) 2 由（ 5）得， (x x 0 ) β(y y y k 1 y 0 ) ( 6 ) 2 k 3 k 2 将（ 4）代入（ 6），得 2( x k 3 x 0 ) ( x k 2 x k 1 x 0 ) ( x k 1 x k x 0 ) 2 2 4 x k 3x k 2 2 x k 1 x k 4 x 0 4 x 0 对应齐次方程的特征方程为 4 3 2 2 0 (7) 代数方程（ 7 ）无正实根，且 αβ , , 2 4 不是（ 7）的根 . 设（ 7）的三个非零根分 别为 1, 2, 3，则 1 2 3 4 1 2 2 3 3 1 2 1 2 3 4 对（ 7）作变换： , 则 12 3 q 0, p 其中 p 1 (2 2 2 ), q 1(833 2 2 ) 4 12 4 123 6
第一部分课后习题 1.学校共1000名学生，235人住在A宿舍，333人住在B宿舍，432人住在C宿舍。学生 们要组织一个10人的委员会，试用下列办法分配各宿舍的委员数： （1）按比例分配取整数的名额后，剩下的名额按惯例分给小数部分较大者。 （2）2.1节中的Q值方法。 （3）d’Hondt方法：将A，B，C各宿舍的人数用正整数n=1，2，3，…相除，其商数如 将所得商数从大到小取前10个（10为席位数），在数字下标以横线，表中A，B，C行有横线的数分别为2，3，5，这就是3个宿舍分配的席位。你能解释这种方法的道理吗。 如果委员会从10人增至15人，用以上3种方法再分配名额。将3种方法两次分配的结果列表比较。 （4）你能提出其他的方法吗。用你的方法分配上面的名额。 2.在超市购物时你注意到大包装商品比小包装商品便宜这种现象了吗。比如洁银牙膏50g 装的每支1.50元，120g装的3.00元，二者单位重量的价格比是1.2：1。试用比例方法构造模型解释这个现象。 （1）分析商品价格C与商品重量w的关系。价格由生产成本、包装成本和其他成本等决定，这些成本中有的与重量w成正比，有的与表面积成正比，还有与w无关的因素。 （2）给出单位重量价格c与w的关系，画出它的简图，说明w越大c越小，但是随着w 的增加c减少的程度变小。解释实际意义是什么。 3.一垂钓俱乐部鼓励垂钓者将调上的鱼放生，打算按照放生的鱼的重量给予奖励，俱乐部 只准备了一把软尺用于测量，请你设计按照测量的长度估计鱼的重量的方法。假定鱼池中只有一种鲈鱼，并且得到8条鱼的如下数据（胸围指鱼身的最大周长）： 先用机理分析建立模型，再用数据确定参数 4.用宽w的布条缠绕直径d的圆形管道，要求布条不重叠，问布条与管道轴线的夹角 应 多大（如图）。若知道管道长度，需用多长布条（可考虑两端的影响）。如果管道是其他形状呢。
2 发酵罐温度控制系统的数学模型 发酵罐温度控制系统实验平台是以一个7L 发酵罐为主体，罐壁设置有冷却套,相应的设立测温点和调节阀,通过阀门调节冷却套内冷却液的流量来实现对发酵罐内温度的控制，发酵罐示意图如图1所示。 图1 发酵罐示意图 在白酒发酵的过程中，发酵罐内由于酵母的作用，在发酵过程中会产生生化反应热，热量的逐渐释放导致发酵温度逐渐上升。在整个发酵过程中，发酵温度必须根据具体的生产工艺进行严格控制，罐内温度通过控制冷却夹套内的冷却水的流量进行降温，整套系统没有外部加热措施。罐内发酵反应热有一部分使罐内温度升高，一部分热量散失到罐壁和冷媒中，在此不考虑发酵体与罐壁之间的热量传递，罐内的热平衡方程为： ? =-Tdt mC Q Q 21 （2-1） 式中 1Q ：发酵过程产生的热量；2Q :发酵过程散失的热量；m ：反应物质量 C ：发酵罐内反应物的比热容；T 发酵罐温度。 公式1-1可以写成： ? =?Tdt MC Q （2-2） 式中 21Q Q Q -=? 对公式1-2求拉普拉斯变换得： s m C T Q S S )()(=? （2-3） 即可由罐内的热平衡方程式可以得到发酵罐内的传递函数为: m C s Q T G S S S 1 ) ()()(= ?= （2-4） 考虑到在实际的过程中的干扰因素,所以被控对象的数学模型中添加一个滞后环节。因此,用一阶惯性加纯滞后环节来表示,其传递函数为 mCs e Q T G s S S S τ-= ?= ) ()()( （2-5）
3 模糊预测控制器的设计及仿真结果 针对发酵罐中发酵对象大时滞、大时变、严格的非线性、多变量耦合等特点。采用了将模糊控制与预测控制结合的方法，利用模糊建模方法建立对象预测模型。将设定值与预测输入值之间的预测误差值及预测误差值的变化率作为模糊控制器的输入，模糊控制器再根据模糊规则来推理得到控制量，通过执行机构控制被控对象。其结构图如图2所示。 图2模糊控制系统结构图 3.1预测控制部分 预测控制算法与动态矩阵控制算法类似, 主要通过预测模型,利用系统的输入输出数据预测未来时刻系统输出,作为糊控制器的输入。 3.1.1预测模型 假设被控对象基于阶跃响应的预测模型向量为T N a a a a ],...,,[21=，N 为建模时域。则在k 时刻对系统施加一个控制增量Δu(k)时,即可算出在其作用下未来时刻N 个输出值的向量形式： )()()(k u a k y k y po m ??+= （3-1） 式中)(k y po 为k 时刻未加Δu(k)时的初始预测值,)(k y m 为k 时刻在Δu(k)作用下的模型预测值。 3.1.2在线校正 当k 时刻对系统施加控制u(k)时，利用预测模型即可得出未来时刻的输出预测值 )(k y m 。但是，由于实际存在的模型时变、非线性、环境干扰等因素的影响，预测值会偏离 实际值，故在k+l 时刻要利用系统的实际输出y (k+1)进行在线校正： )]|1()1([)()(k k y k y h k y k y m m p +-++= （3-2） 式中h 为N 维误差校正向量,这里取0.11=h ，9.0=i h ，i=2,3...,N 。)(k y p 为校正后的预测值，经过移位后即可作为k+1时刻的初始预测值，用向量形式可表示为： )()1(k y S k y p po ?=+ （3-3） 式中S 为位移阵。
基于单片机的温度控制 系统设计文献综述 前言 随着现代工业的发展，人们需要对工业生产中有关温度系统进行控制，如钢铁冶炼过程需要对刚出炉的钢铁进行热处理，塑料的定型及各种加热炉、热处理炉、反应炉和锅炉中温度进行实时监测和精确控制温度是日常生活、工业、医学、环境保护、化工、石油等领域最常遇到的一个物理量。而且，很多领域的温度可能较高或较低，现场也会较复杂，有时人无法靠近或现场无需人力来监控。如加热炉大都采用简单的温控仪表和温控电路进行控制, 存在控制精度低、超调量大等缺点, 很难达到生产工艺要求。且在很多热处理行业都存在类似的问题，所以，设计一个较为通用的温度控制系统具有重要意义。这时我们可以采用单片机控制，这些控制技术会大大提高控制精度，不但使控制简捷，降低了产品的成本，还可以和计算机通讯，提高了生产效率. 单片机是指芯片本身，而单片机系统是为实现某一个控制应用需要由用户设计的，是一个围绕单片机芯片而组建的计算机应用系统，这是单片机应用系统。单片机自问世以来，性能不断提高和完善，其资源又能满足很多应用场合的需要，加之单片机具
有集成度高、功能强、速度快、体积小、功耗低、使用方便、价格低廉等特点，因此，应用日益广泛，并且正在逐步取代现有的多片微机应用系统。 1.陈岩《基于ARM 的远程控制温控系统的设计》一个基于ARM的远程控制系统的设计.该系统以无线寻呼网络接收POCSAG编码的控制命令字,同时利用DIMF信号发送器将要反馈的数据通过公用电话网络以DTMF编码传送回去,从而实现了一个功能完整的远程控制系统,弥补了以往远程控制系统的不足同。 2．金凯鹏胡即明《基于模糊PID 算法远程温度控制系统的实现》针对实时温度控制对象,算法远程温度控制系统是一套远程控制系统,并结合了模糊PID控制算法,利用其电路组成和设计原理,实现了对远程温度系统的监视和控制功能.采集端主要实现温度采集、数码显示、温度设定、无线编码发射、加热开关控制等功能;监控部分主要实现无线解码接收、温度显示、报警等功能模块.本系统实现了实时控制与无线传输结合. 3.王晓员《基于单片机多点温度控制的硬件构建设计》针对目前许多塑料反应炉温度控制不准确的现状,进行了基于MCS-51系列单片机多点温度控制的硬件构建的设计.采用数字化温度传感器DS18820,TLC2543型号的12位开关电容运次逼近模数A/D转换器.成本低、可靠性高 4.王芳《利用单片机实现温度智能控制》温度控制系统是
第一章 课后习题6. 利用1.5节药物中毒施救模型确定对于孩子及成人服用氨茶碱能引起严重中毒和致命的最小剂量。 解：假设病人服用氨茶碱的总剂量为a ，由书中已建立的模型和假设得出肠胃中的药量为： )()0(mg M x = 由于肠胃中药物向血液系统的转移率与药量)(t x 成正比，比例系数0>λ，得到微分方程 M x x dt dx =-=)0(,λ（1） 原模型已假设0=t 时血液中药量无药物，则0)0(=y ，)(t y 的增长速度为x λ。由于治疗而减少的速度与)(t y 本身成正比，比例系数0>μ，所以得到方程： 0)0(,=-=y y x dt dy μλ（2） 方程（1）可转换为：t Me t x λ-=)( 带入方程（2）可得：)()(t t e e M t y λμμ λλ ----= 将01386=λ和1155.0=μ带入以上两方程，得： t Me t x 1386.0)(-= )(6)(13866.01155.0---=e e M t y t 针对孩子求解，得： 严重中毒时间及服用最小剂量：h t 876.7=，mg M 87.494=； 致命中毒时间及服用最小剂量：h t 876.7=，mg M 8.4694= 针对成人求解： 严重中毒时间及服用最小剂量：h t 876.7=，mg M 83.945= 致命时间及服用最小剂量：h t 876.7=，mg M 74.1987= 课后习题7. 对于1.5节的模型，如果采用的是体外血液透析的办法，求解药物中毒施救模型的血液用药量的变化并作图。
解：已知血液透析法是自身排除率的6倍，所以639.06==μu t e t x λ-=1100)(，x 为胃肠道中的药量，1386.0=λ )(6600)(t t e e t y λμ---= 1386.0,639.0,5.236)2(,1100,2,====≥-=-λλλu z e x t uz x dt dz t 解得：()2,274.112275693.01386.0≥+=--t e e t z t t 用matlab 画图: 图中绿色线条代表采用体外血液透析血液中药物浓度的变化情况。 从图中可以看出，采取血液透析时血液中药物浓度就开始下降。T=2时，血液中药物浓度最高，为236.5；当z=200时，t=2.8731，血液透析0.8731小时后就开始解毒。 第二章 1.用 2.4节实物交换模型中介绍的无差别曲线的概念，讨论以下的雇员和雇主之间的关系： 1）以雇员一天的工作时间和工资分别为横坐标和纵坐标，画出雇员无差别曲线族的示意图，解释曲线为什么是那种形状； 2）如果雇主付计时费，对不同的工资率画出计时工资线族，根据雇员的无差别曲线族和雇主的计时工资线族，讨论双方将在怎样的一条曲线上达成协议； 3）雇员和雇主已经达成了协议，如果雇主想使用雇员的工作时间增加到t 2，他有两种
数学模型第三版课后习 题答案 TPMK standardization office【 TPMK5AB- TPMK08- TPMK2C- TPMK18】
《数学模型》作业解答 第二章(1)（2008年9月16日） 1． 学校共1000名学生，235人住在A 宿舍，333人住在B 宿舍，432人住 在C 宿舍.学生们要组织一个10人的委员会，试用下列办法分配各宿舍的委员数： （1）. 按比例分配取整数的名额后,剩下的名额按惯例分给小数部分较大者; （2）. §1中的Q 值方法； （3）.d ’Hondt 方法：将A 、B 、C 各宿舍的人数用正整数n=1,2,3,……相除，其商数如下表： 将所得商数从大到小取前10个（10为席位数），在数字下标以横线，表中A 、B 、C 行有横线的数分别为2，3，5，这就是3个宿舍分配的席位.你能解释这种方法的道理吗 如果委员会从10个人增至15人，用以上3种方法再分配名额，将3种方法两次分配的结果列表比较. 解：先考虑N=10的分配方案， 方法一（按比例分配） 分配结果为： 4 ,3 ,3321===n n n 方法二（Q 值方法） 9个席位的分配结果（可用按比例分配）为：
第10个席位：计算Q 值为 3Q 最大，第10个席位应给C.分配结果为 5 ,3 ,2321===n n n 方法三（d ’Hondt 方法） 此方法的分配结果为：5 ,3 ,2321===n n n 此方法的道理是：记i p 和i n 为各宿舍的人数和席位（i=1,2,3代表A 、B 、C 宿舍）. i i n p 是每席位代表的人数，取,,2,1 =i n 从而得到的i i n p 中选较大者，可使对所有的, i i i n p 尽量接近. 再考虑15=N 的分配方案，类似地可得名额分配结果.现将3种方法两次分配的结果列表如下： 2． 试用微积分方法，建立录像带记数器读数n 与转过时间的数学模型. 解： 设录像带记数器读数为n 时，录像带转过时间为t.其模型的假设见课本. 考虑t 到t t ?+时间内录像带缠绕在右轮盘上的长度，可得 ,2)(kdn wkn r vdt π+=两边积分，得 ??+=n t dn wkn r k vdt 0 )(2π 第二章(2)（2008年10月9日） 15．速度为v 的风吹在迎风面积为s 的风车上，空气密度是ρ ，用量纲分析方法确定风车获得的功率P 与v 、S 、ρ的关系.
数学建模试题（带答案） 第一章 4.在1.3节“椅子能在不平的地面上放稳吗”的假设条件中，将四脚的连线呈正方形改为长方形，其余不变。试构造模型并求解。 答：相邻两椅脚与地面距离之和分别定义为)()(a g a f 和。f 和g 都是连续函数。椅子在任何位置至少有三只脚着地，所以对于任意的a ，)()(a g a f 和中至少有一个不为零。不妨设0)0(,0)0(g >=f 。当椅子旋转90°后，对角线互换， 0π/2)(,0)π/2(>=g f 。这样，改变椅子的位置使四只脚同时着地。就归结为证 明如下的数学命题： 已知a a g a f 是和)()(的连续函数，对任意0)π/2()0(,0)()(,===?f g a g a f a 且， 0)π/2(,0)0(>>g f 。证明存在0a ，使0)()(00==a g a f 证：令0)π/2(0)0(),()()(<>-=h h a g a f a h 和则， 由g f 和的连续性知h 也是连续函数。 根据连续函数的基本性质， 必存在0a （0＜0a ＜π/2）使0)(0=a h ，即0)()(00==a g a f 因为0)()(00=?a g a f ，所以0)()(00==a g a f
8
第二章 7. 10.用已知尺寸的矩形板材加工半径一定的圆盘，给出几种简便有效的排列方法，使加工出尽可能多的圆盘。
第三章 5.根据最优定价模型 考虑成本随着销售量的增加而减少，则设 kx q x q -=0)( (1)k 是产量增加一个单位时成本的降低 ， 销售量x 与价格p 呈线性关系0,,>-=b a bp a x (2) 收入等于销售量乘以价格p :px x f =)( (3) 利润)()()(x q x f x r -= (4) 将（1）（2）（3）代入（4）求出 ka q kbp pa bp x r --++-=02)( 当k q b a ,,,0给定后容易求出使利润达到最大的定价*p 为 b a kb ka q p 2220*+--= 6.根据最优定价模型 px x f =)( x 是销售量 p 是价格，成本q 随着时间增长，ββ,0t q q +=为增长率，0q 为边际成本（单位成本）。销售量与价格二者呈线性关系0,,>-=b a bp a x . 利润)()()(x q x f x u -=.假设前一半销售量的销售价格为1p ，后一半销售量的销售价格为2p 。 前期利润 dt bp a t q p p u T ))](([)(12 /011--=? 后期利润 dt bp a t q p p u T T ))](([)(22/22--=? 总利润 )()(21p u p u U += 由 0,02 1=??=??p U p U 可得到最优价格: )]4([2101T q b a b p β++= )]4 3([2102T q b a b P β++=
： 远程与继续教育学院 本科毕业论文（设计） 题目：温控系统的设计及仿真(MATLAB) 、 学习中心： 学号： 姓名： 专业：机械设计制造及自动化 指导教师： " 2013 年 2 月 28 日
) 摘要 温度是工业对象中一个主要的被控参数，它是一种常见的过程变量，因为它直接影响燃烧、化学反应、发酵、烘烤、煅烧、蒸馏、浓度、挤压成形，结晶以及空气流动等物理和化学过程。温度控制不好就可能引起生产安全，产品质量和产量等一系列问题。温度控制是许多设备的重要的构成部分，它的功能是将温度控制在所需要的温度范围内，以利于进行工件的加工与处理。 一直以来，人们采用了各种方法来进行温度控制，都没有取得很好的控制效果。如今，随着以微机为核心的温度控制技术不断发展，用微机取代常规控制已成必然，因为它确保了生产过程的正常进行，提高了产品的数量与质量，减轻了工人的劳动强度以及节约了能源，并且能够使加热对象的温度按照某种指定规律变化。 实践证明，用于工业生产中的炉温控制的微机控制系统具有高精度、功能强、经济性好的特点，无论在提高产品质量还是产品数量，节约能源，还是改善劳动条件等方面都显示出无比的优越性。 本设计以89C51单片机为核心控制器件，以ADC0809作为A/D转换器件，采用闭环直接数字控制算法，通过控制可控硅来控制热电阻，进而控制电炉温度，最终设计了一个满足要求的电阻炉微型计算机温度控制系统。 关键词：1、单片机；2、PLC；3、MATLAB &
（ 目录 1单片机在炉温控制系统中的运用 (6) 1、1系统的基本工作原理 (6) 2温控系统控制算法设计 (7) 温度控制算法的比较 (7) 数字PID算法 (11) 、 3 结论 (21) 致谢 (22) 参考文献 (23) [
。 例１差分方程—-资金的时间价值 问题1：抵押贷款买房——从一则广告谈起 每家人家都希望有一套(甚至一栋)属于自己的住房，但又没有足够的资金一次买下，这就产生了贷款买房的问题。先看一下下面的广告（这是1991年1月1日某大城市晚报上登的一则广告）,任何人看了这则广告都会产生许多疑问,且不谈广告中没有谈住房面积、设施等等,人们关心的是：如果一次付款买这栋房要多少钱呢?银行贷款的利息是多少呢？为什么每个月要付12０0元呢？是怎样算出来的?因为人们都知道,若知道了房价(一次付款买房的价格），如果自己只能支付一部分款,那就要把其余的款项通过借贷方式来解决,只要知道利息，就应该可以算出五年还清每月要付多少钱才能按时还清贷款了,从而也就可以对是否要去买该广告中所说的房子作出决策了。现在我们来进行数学建模。由于本问题比较简单无需太多的抽象和简化。 a。明确变量、参数，显然下面的量是要考虑的： 需要借多少钱，用记; 月利率(贷款通常按复利计）用R记； 每月还多少钱用x记; 借期记为N个月。 b.建立变量之间的明确的数学关系。若用记第k个月时尚欠的款数,则一个月后（加上利息后)欠款 , 不过我们又还了x元所以总的欠款为 ｋ=0，1，２，３， 而一开始的借款为.所以我们的数学模型可表述如下 (1) c. （1)的求解。由
（2）这就是之间的显式关系。 d.针对广告中的情形我们来看（１）和（2）中哪些量是已知的。Ｎ=5年=60个月，已知;每月还款x＝１20０元，已知Ａ.即一次性付款购买价减去７０000元后剩下的要另外去借的款，并没有告诉你，此外银行贷款利率Ｒ也没告诉你,这造成了我们决策的困难.然而,由（2）可知6０个月后还清,即,从而得 (3） A和x之间的关系式，如果我们已经知道银（３)表示Ｎ=6０,x=1２00给定时0 A。例如，若Ｒ=0．01,则由(3)可算得行的贷款利息R，就可以算出0 5３94６元。如果该房地产公司说一次性付款的房价大于70０0０十53946＝12３946元的话，你就应自己去银行借款。事实上，利用图形计算器或Maｔheｍａtiｃa这样的 数学软件可把（3）的图形画出来,从而可以进行估算决策。以下我们进一步考虑下面两个问题。 注1问题1标题中“抵押贷款”的意思无非是银行伯你借了钱不还，因而要你用某种不动产（包括房子的产权）作抵押，即万一你还不出钱了，就没收你的不动产。 例题１某高校一对年青夫妇为买房要用银行贷款６000０元，月利率0.01，贷款期25年=300月，这对夫妇希望知道每月要还多少钱，25年就可还清。假设这对
基于模糊PID算法的电阻炉温度控制系统设计引言 电加热炉是典型工业过程控制对象，其温度控制具有升温单向性，大惯性，纯滞后，时变性等特点，很难用数学方法建立精确的模型和确定参数。而PID控制因其成熟，容易实现，并具有可消除稳态误差的优点，在大多数情况下可以满足系统性能要求，但其性能取决于参数的整定情况。且快速性和超调量之间存在矛盾，使其不一定满足快速升温、超调小的技术要求。模糊控制在快速性和保持较小的超调量方面有着自身的优势，但其理论并不完善，算法复杂，控制过程会存在稳态误差。 将模糊控制算法引入传统的加热炉控制系统构成智能模糊控制系统，利用模糊控制规则自适应在线修改PID参数，构成模糊自整定：PID控制系统，借此提高其控制效果。基于PID控制算法，以ADUC845单片机为主体，构成一个能处理较复杂数据和控制功能的智能控制器，使其既可作为独立的单片机控制系统，又可与微机配合构成两级控制系统。该控制器控制精度高，具有较高的灵活性和可靠性。 2 温度控制系统硬件设计 该系统设计的硬件设计主要由单片机主控、前向通道、后向通道、人机接口和接口扩展等模块组成，如图l所示。由图1可见，以内含C52兼容单片机的ADUC845为控制核心．配有640 KB的非易失RAM数据存储器、外扩键盘输人、320x240点阵的图形液晶显示器进行汉字、图形、曲线和数据显示，超温报警装置等外围电路；预留微型打印机接口，可以现场打印输出结果；预留RS232接口，能和PC机联机，将现场检测的数据传输至PC机来进一步处理、显示、打印和存档。
电阻炉的温度先由热电偶温度传感器检测并转换成微弱的电压信号，温度变送器将此弱信号进行非线性校正及电压放大后，由单片机内部A／D转换器将其转换成数字量。此数字量经数字滤波、误差校正、标度变换、线性拟合、查表等处理后。一方面将炉窑温度经人机面板上的LCD显示：另一方面将该温度值与被控制值(由键盘输入的设定温度值)比较，根据其偏差值的大小，提供给控制算法进行运算，最后输出移相控制脉冲，放大后触发可控硅导通(即控制电阻炉平均功率)。达到控制电炉温度的目的。如果实际测得的温度值超过了该系统所要求的温度范围，单片机就向报警装置发出指令，系统进行报警。 2．1 系统主控模块 系统主控模块电路如图2所示，它主要由CPU及数据存储器，晶体振荡器、复位电路、图形液晶显示器(LCD)及控制电路、微型打印机接口控制电路、实时日历时钟，热电偶信号处理电路等构成。这里，该系统设计可测量3点温度。传感器选择K型(镍铬-镍硅)热电偶，可用于从室温到1 200°C的温度测量，测量范围宽，精度高。在温度测量范围内K型热电偶的输出热电势只有0～45．119 mV，为了和ADUC845的A ／D转换器相匹配，采用ACl226和1B51作为信号调理电路，由AC1226、1B51构成热电偶冷端温度补偿及信号调理器电路。当热端距测温仪表较远时，需利用热电偶匹配导线将冷端延长。CD4051为多路模拟开关，由ABC控制接通，当5～3接通时，输入接地，UO输出UOmin，用于零点校准；当4～3接通时，单片机1．25 V稳定参考电压Uref，再经电阻R1、R2分压，得到毫伏级参考输入电压，UO输出UOmax，用于增益校准；当2～3、1～3、12～3分别分时接通时，依次输入3个热电偶正常测温所得变换电压，UO从而输出3个温度点所对应的电压UOA，UOB，UOC。在HI端与+UISO端之间串上一只220 MΩ上拉电阻，一旦热电偶开路，HI端即被偏置为+UISO，迫使1B51的输出电压超量程，由此判定热电偶已开路。多路模拟开关和测量数据采集过程在单片机协调下工作，每次数据采集都进行自动判断和校准阁。 2．2 控制输出驱动电路 对温度的控制是通过可控硅调功器电路实现，如图3所示。双向可控硅管和硅碳棒串接在交流220 V、50 Hz交流市电回路中，图3中只给出了A相。移相触发脉冲由ADUC845用软件在P1．3引脚上产生的，零同步脉冲同步后，经光耦合管和驱动器输出送到可控硅的控制极。过零同步脉冲由过零触发电路产生，利用同步变压器和电压比较器LM311组成正弦交流电的正半波过零检测电路，它在交流电每一个正半周的