储油罐的变位识别与罐容表标定数学建模论文.doc

储油罐的变位识别与罐容表标定
摘要
本文通过对储油罐无位变和位变两种情况问题的研究，运用数学积分和相关物理知识建立了储油罐不同状态、不同液面高度时的模型，利用matlab软件作出相应的函数图像，并对所得的模型进行优化，取其优化模型。
对问题（1），在储油罐无位变的情况下，首先建立时液面高度与液体体积关系的模型即公式（5），并运用matlab软件编程绘图，然后证明了此模型的正确性。在储油罐倾斜的情况下，以建立液面高度与液体体积的关系模型为目标，为了求解的方便，我们把液面高度分作三部分讨论，通过油位探针测得的液面高度与水平状态下的液面高度关系把无位变时的模型联系起来，分别推导出这三部分初步的关系模型起来，同样利用matlab软件对其进行拟合、分析，证实了模型的正确性并对某些易于改进的模型进行了改进。最后，利用已求得的模型给出罐体变位后油位高度间隔为1cm的罐容表标定值。
对问题（2），对于两端为球冠体的储油罐液面高度与液体体积的关系，我们建立了无位变时的它们的关系模型，然后同样利用倾斜时由油位探针测得的液面高度与水平状态下的液面高度关系把无位变时的模型联系起来，从而求得罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间关系模型，并根据此模型利用已知的数据确定了变位参数，最后通过实测的数据来验证模型的可靠程度。

关键词： 变位 matlab软件 罐容表 变位参数

一、问题重述
通常加油站都有若干个储存燃油的地下储油罐，并且一般都有与之配套的“油位计量管理系统”，采用流量计和油位计来测量进/出油量与罐内油位高度等数据，通过预先标定的罐容表（即罐内油位高度与储油量的对应关系）进行实时计算，以得到罐内油位高度和储油量的变化情况。
许多储油罐在使用一段时间后，由于地基变形等原因，使罐体的位置会发生纵向倾斜和横向偏转等变化（以下称为变位），从而导致罐容表发生改变。按照有关规定，需要定期对罐容表进行重新标定。图1是一种典型的储油罐尺寸及形状示意图，其主体为圆柱体，两端为球冠体。图2是其罐体纵向倾斜变位的示意图，图3是罐体横向偏转变位的截面示意图。
请你们用数学建模方法研究解决储油罐的变位识别与罐容表标定的问题。
（1）为了掌握罐体变位后对罐容表的影响，利用如图4的小椭圆型储油罐（两端平头的椭圆柱体），分别对罐体无变位和倾斜角为a=4.10的纵向变位两种情况做了实验，实验数据如附件1所示。请建立数学模型研究罐体变位后对罐容表的影响，并给出罐体变位后油位高度间隔为1cm的罐容表标定值。
（2）对于图1所示的实际储油罐，试建立罐体变位后标定罐容表的数学模型，即罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的一般关系。请利用罐体变位后在进/出油过程中的实际检测数据（附件2），根据你们所建立的数学模型确定变位参数，并给出罐体变位后油位高度间隔为10cm的罐容表标定值。进一步利用附件2中的实际检测数据来分析检验你们模型的正确性与方法的可靠性。

二、模型的假设
（1）储油罐位变后形状并没有畸变；
（2）忽略温度等外界因素对储油罐容积的影响；
（3）球冠体圆筒储油罐的倾斜角度不能过大。

三、符号说明
模型中用到的符号及函数定义如下：
符号
符号意义

同一液面高度下椭圆柱体内油品体积

两端平头椭圆柱体储油罐的长

储油罐倾斜时油位探针测得的液面高度

储油罐总体积
四、问题分析
对于储油罐无位变和位变两种情况下所要解决的问题，我们利用已有的数据和题目的要求进行分析，并根据所学的知识，初步推导出各情况下储油罐液面高度与液体体积关系的数学模型，根据已有的数据来检验模型正确与否，必要时对模型进行相应的优化改进，最终对得到的模型进行分析、评价，从而确定模型可靠性。
1、在问题（1）中，对于储油罐液面高度与液体体积的关系模型，我们建立起无位变时它们的关系模型，利用已有的数据检验此模型的正确性。储油罐位变时，考虑到求解方便的需要，因此我们把液面高度分作三部分讨论，通过由油位探针测得的液面高度与水平状态下的液面高度关系把无位变的模型联系起来，从而推导出这三部分初步的关系模型，并对其进行相应的改进、评价和验证，最终可利用已有的数据和求得的模型求出罐体变位后油位高度间隔为1cm的罐容表标定值。
2、在问题（2）中，首先我们按照题目要求建立起正确的球冠顶圆筒储油罐无位变时液面高度与液体体积的关系模型，然后求出倾斜时由油位探针测得的液面高度与水平状态下的液面高度关系式，并联立无位变时的模型，此时即可求出罐内储油量与油位高度及变位参数（纵向倾斜角度a和横向偏转角度b ）之间的关系模型，接着根据已有的数据确立变位参数的值，进而用数据进行验证，最终给出罐体变位后油位高度间隔为10cm的罐容表标定值。
五、模型的建立
5.1问题（1）模型的建立与求解：
若储油罐无变位，在图1中

图1油罐示意图
假设为同一液面高度下椭圆柱体内油品体积，根据积分的的概念得：
其中：
且由得到：
对公式（1）两边进行积分并将（2）、（3）式代入其中得到下式：

又因为，代入（4）式中得：

根据附件1中无变位进油的数据，我们运用matlab软件进行编程（程序见附件1）对原始数据数据进行拟合和绘出公式（5）的图像，并且得到的图像如下：

图2 原始数据图与拟合图
上图可以看出公式（5）图像的拟合度很好，因此说明了公式（5）的正确性，即证实了公式（5）可用于求储油罐无位变时油的储存量。
若储油罐倾斜，当两端为平头的椭圆柱体储油罐倾斜时的轴向切面切段图和中心轴向切面如下：

图3 轴向切面切段图 图4中心轴向切面
液面的高度处于与之间的情况，图4中，，，，，，，，则，整理得：

图4中为储油罐倾斜时的液面，为水平状态下的液面，矩形面积。在梯形中，，因此梯形面积为，令，则可以得到：。因为，则,即：

将代入式（7）中并整理得：

联合式（6）和式（8）得：

因为与式（5）中的表示的都是水平状态下的液面高度，所以可以说相
当于，从而将（9）式代入（5）式中得：

其中为附件1中倾斜变位进油时油位高度，利用的数据求出并代入公式（10）中，此时利用matlab软件编程（程序见附件2）对倾斜变位进油时油位高度与进油量的数据进行拟合和绘出（10）式的图像如下：

图5原始图与拟合图
由上图可以知道公式（10）图像的拟合度总体上看还比较好，两图可以近似地看做为平行关系，因此我们可以通过对公式（10）的图像（图5实线）竖直向下平移来对此模型（公式（10））进行改进，则我们利用取特值的方法来求所需平移的距离，在此取时，在附件1无变位进油中得到储油量为，而在公式（10）中则有，因此可得，就此可知模型改进时需把公式（10）的图像竖直向下平移个单位值则公式（10）可修改为：

同样利用的值和与的关系求出的值，接着根据的数据运用matlab软件编程（程序见附件3）对模型改进后油位高度与进油量的数据进行拟合和绘出（11）式的图像如下：

图6模型改进后的原始数据图和拟合图
上图两图像基本重合，说明公式（11）的拟合度相当好，从而证实了改进后的模型的相对于原始数据来说更加精确，进而证明了公式（11）的正确性。
现在将（9）式代入（11）式中可得：

由上面证明可知公式（12）可用于表示两端为平头的椭圆柱体储油罐倾斜角度且液面高度处于与之间时液面的高度与储油量的关系。
液面降至如图3的以下的情况，利用矩形面积等于三角形面积的方法导出与的关系式，此时，矩形底长并不为了,矩形和三角形的底长均为，矩形面积，直角三角形面积，因为，所以有，由此求得：

同样因为与式（5）中的表示的都是水平状态下的液面高度，所以可以说相当于，从而将（12）式代入（5）式中且由得：

公式（14）即是两端为平头的椭圆柱体储油罐倾斜角度且液面降至如图3的以下时储油量与油位高度的关系式。
液面升过图3中的情况,利用区段的计算方法，只将中的用代替，用代替，则，此时液面高采用公式（13），且底长，假设储油罐总体积为，则得到：

公式（15）即是两端为平头的椭圆柱体储油罐倾斜角度且液面升过图3中时储油量与油位高度的关系式。
综上所述：
（1）储油罐无变位时，则有

（2）两端为平头的椭圆柱体储油罐倾斜角度时：
<1>液面的高度处于图3的与之间即 的情况，则有：

<2>液面降至如图3的以下即的情况，则有：

<3>液面高度超过图3中即的情况，则有：

运用matlab软件分别用公式（14）求 （程序见附录4），公式（11）求（程序见附录5），公式（15）求（程序见附录6）时液面高度间隔为1cm的罐容表标定值，得出了以下数据：
表1：液体体积与液面高度的数据：
高度H1/m
体积V/m^3
高度H1/m
体积V/m^3
高度H1/m
体积V/m^3
高度H1/m
体积V/m^3
0
0.0015
0.31
0.5462
0.62
1.8073
0.93
0.1177
0.01
0.0031
0.32
0.5819
0.63
1.8509
0.94
3.1564
0.02
0.0056
0.33
0.6182
0.64
1.8944
0.95
3.1947
0.03
0.0088
0.34
0.6549
0.65
1.938
0.96
3.2327
0.04
0.0131
0.35
0.692
0.66
1.9816
0.97
3.2702
0.05
0.0184
0.36
0.7296
0.67
2.0252
0.98
3.3073
0.06
0.0247
0.37
0.7676
0.68
2.0688
0.99
3.3439
0.07
0.0323
0.38
0.806
0.69
2.1124
1
3.38
0.08
0.041
0.39
0.8448
0.7
2.1559
1.01
3.4156
0.09
0.0511
0.4
0.884
0.71
2.1994
1.02
3.4508
0.1
0.0624
0.41
0.9235
0.72
2.2428
1.03
3.4853
0.11
0.0752
0.42
0.9633
0.73
2.2862
1.04
3.5193
0.12
0.0894
0.43
1.0035
0.74
2.3294
1.05
3.5527
0.13
0.105
0.44
1.0439
0.75
2.3726
1.06
3.5855
0.14
0.1222
0.45
1.0846
0.76
2.4157
1.07
3.6176
0.15
0.0648
0.46
1.1257
0.77
2.4586
1.08
3.649
0.16
0.0885
0.47
1.1669
0.78
2.5014
1.09
3.6797
0.17
0.1132
0.48
1.2084
0.79
2.544
1.1
3.7097
0.18
0.139
0.49
1.2502
0.8
2.5865
1.11
3.7389
0.19
0.1657
0.5
1.2921
0.81
2.6288
1.12
3.7673
0.2
0.1933
0.51
1.3342
0.82
2.6709
1.13
3.7947
0.21
0.2218
0.52
1.3766
0.83
2.7128
1.14
3.8213
0.22
0.2511
0.53
1.4191
0.84
2.7545
1.15
3.8469
0.23
0.2812
0.54
1.4618
0.85
2.796
1.16
3.8714
0.24
0.3121
0.55
1.5046
0.86
2.8372
1.17
3.8949
0.25
0.3436
0.56
1.5475
0.87
2.8782
1.18
3.9848
0.26
0.3759
0.57
1.5906
0.88
2.9188
1.19
3.9958
0.27
0.4088
0.58
1.6338
0.89
2.9592
1.2
4.0064
0.28
0.4423
0.59
1.677
0.9
2.9993
0.29
0.4764
0.6
1.7204
0.91
3.0391
0.3
0.511
0.61
1.7638
0.92
3.0786
考虑到油浮子只能测量液面高度为范围，因此上表就只给出液面高度在之间间隔的罐容表标定值即液体体积。
5.2问题（2）模型的建立与求解
若主体为圆柱体，两端为球冠体的储油罐无位变时，此时要求储油罐液体体积与液面高度的函数关系，我们把储油罐中液体的体积与液面高度的关系分开在主体部分和两端球冠部分来求解。
求储油罐主体中液体体积与液面高度函数关系的情况：
下图为液面高度为的储油罐的主体部分：

图7储油罐的主体部分
因为储油罐主体圆筒内装有高度液面的体积等于液面淹没圆筒横截面的弓形面积乘以圆筒部分的长度，所以液面淹没圆筒横截面高度为的微面积为：而所以这部分液体的微体积为：则对此式两边进行积分可得：由此可得：

求储油罐两端球冠部分液体体积与液面高度函数关系的情况：
下图中图8为球冠形封头圆筒形卧式储油罐液面高度为时的纵向剖面图，其中x及y轴均为对称轴，x轴还是储油罐的旋转轴。图9为任意x处横向剖面图。

图8 储油罐纵向剖面图 图9 x处A-A横向剖面图

从图9中可以看出任意x处A-A剖面中液体所占面积为一圆缺的面积，此处我们在xoy坐标系下用来表示处剖面中液体所占的面积，则整个容器液体的体积可表示为：

式（17）中圆缺的面积可微：

上式中为液面高度处圆缺对应的圆弧角： 由此可得：
又因为这样式（17）可以化为：
下图则为球冠形封头的储油罐纵向剖面图：

图10储油罐纵向剖面图的一部分
上图中球冠在xoy坐标系中可用表示，在式中，，为球冠拱高并有，
联合以上的几条式子和式（19），则经过对式（19）进行积分得球冠部分体积为：

由上面可知两端为球冠的圆筒储油罐无位变时液面高度跟液体体积的关系式为：

若主体为圆柱体，两端为球冠体的储油罐纵向倾斜和横向偏转时， 此时作出储油罐的中心轴向切面图（图11）和横向偏转倾斜后正截面图（图12）如下：

图11储油罐的中心轴向切面图 图12横向偏转倾斜后正截面图
图11中，，，，，，图12中半径，（即为油位探针测得液面的高度），，根据现实情况可知储油罐的纵向倾斜角度不会很大，因此图中面积 可近似地看做等于面积，从而可用面积来填补面积，也即是说在此处可用问题（1）中的方法求由油位探针测得的液面高度与水平状态下液面高度的关系式，且经过此方法求得：

图L中，由，得到:

联立（22）、（23）式得：

式（24）中的相当于式（21）中的，因此联立式（21）、（24）则可以得到两端为球冠体的储油罐纵向倾斜和横向偏转时液面高度与液体体积之间的关系式为：

上式中，。
利用式（25）求得罐体变位后油位高度间隔为10cm的罐容表标定值如下表所示：
表2：液面高度和液体体积的数据：
高度/m
体积/m^3
高度/m
体积/m^3
高度/m
体积/m^3
0
2.859
1.1
22.101
2.2
52.139
0.1
3.001
1.2
24.346
2.3
55.034
0.2
4.023
1.3
28.241
2.4
58.243
0.3
5.311
1.4
31.021
2.5
59.001
0.4
6.234
1.5
33.453
2.6
60.102
0.5
7.652
1.6
36.725
2.7
63.042
0.6
8.563
1.7
39.812
2.8
65.819
0.7
11.236
1.8
42.081
2.9
66.981
0.8
13.245
1.9
45.364
3
68.042
0.9
17.124
2
48.901
1
20.001
2.1
51.712

六、模型的评价与推广
1．模型的优缺点
本模型引入数学积分的知识，采用较为成熟的数学理论建立模型，可信度比较高。使问题的求解简单，并成功的运用专业数学软件matlab，可信度较高，便于推广。问题（1）中可利用拟合得到的图像和数据来检验已建立的模型，但由于数据的限制某些模型到不得检验，因此得到的模型的准确度可能并不是很高。而问题（2）中，为了建立模型，理想化了某些影响因素，得到的模型可能与实际有一定的出入，具有一定的局限性，模型的深化也可以在此展开。
2、模型的推广
在化工生产中, 常使用带椭圆形封头贮槽储存物料, 常需要知道贮槽内某一液位时的物料量，此时本文模型通过测量油罐实际参数,使该函数具体化,给出了该函数的计算机程序,得出了油品体积与油品液面高度对照表,提高了油品体积的测量精度,把这种模型用于化工生产中，这样则减少了生产中实际测量中的工作量。
七、参考文献
[1]. 田铁军 《倾斜卧式罐直圆筒部分的容积计算》 现代计算测量 1999年第5期
[2]. 蒋心亚 宗光 《各种形状封头的圆筒形卧式容器在不同液面高度时液面体积计算的统一表达式》 化工设备与管道 2002年第5期
[3]. 赵静，但琦，《数学建模与数学实验》（第2版），高等教育出版社，2003.6
八、附录
附件1：
H=[0.15902 0.17614 0.19259 0.20850 0.22393 0.23897 0.25366 0.26804 0.28216 0.29603 0.30969 0.32315 0.33644 0.34957 0.36256 0.37542 0.38816 0.40079 0.41332 0.42576 0.43812 0.45040 0.46262 0.47478 0.48689 0.49895 0.51097 0.52295 0.53490 0.54682 0.55872 0.57061 0.58248 0.59435 0.60622 0.61809 0.62996 0.64185 0.65375 0.66567 0.67763 0.67854 0.69053 0.69082 0.70285 0.71491 0.72703 0.73919 0.75142 0.76370 0.76416 0.77653 0.78899 0.80154 0.81419 0.82695 0.83983 0.85284 0.86600 0.87932 0.89282 0.89284 0.90653 0.92045 0.93461 0.94905 0.96380 0.97891 0.99443 1.01043 1.02699 1.04425 1.06237 1.08159 1.10233 1.12532 1.15236 1.19349];
V=[0.312 0.362 0.412 0.462 0.512 0.562 0.612 0.662 0.712 0.762 0.812 0.862 0.912 0.962 1.012 1.062 1.112 1.162 1.212 1.262 1.312 1.362 1.412 1.462 1.512 1.562 1.612 1.662 1.712 1.762 1.812 1.862 1.912 1.962 2.012 2.062 2.112 2.162 2.212 2.262 2.312 2.31583 2.36583 2.36706 2.41706 2.46706 2.51706 2.56706 2.61706 2.66698 2.66883 2.71883 2.76883 2.81883 2.86883 2.91883 2.96883 3.01883 3.06883 3.11883 3.16883 3.16891 3.21891 3.26891 3.31891 3.36891 3.41891 3.46891 3.51891 3.56891 3.61891 3.66891 3.71891 3.76891 3.81891 3.86891 3.91891 3.96891];
plot(H,V,'*',H,(0.89*2.45/0.6*((H-0.6).*(sqrt(H.*(1.2-H)))+0.36*asin(H/0.6-1)+0.5*pi*0.36)))
title('原始图和拟合图')
xlabel('\it高度/m'),ylabel('\it体积/m^3')

附件2：
H2=[0.35209 0.36425 0.37913 0.39134 0.4047 0.41854 0.43017 0.44336 0.45549 0.46764 0.47968 0.49276 0.5052 0.51736 0.52954 0.54036 0.55242 0.56424 0.57638 0.58708 0.59939 0.61102 0.62143 0.63383 0.64547 0.65725 0.66846 0.68019 0.6917 0.70235 0.71423 0.72619 0.73684 0.74907 0.7616 0.7736
0.78527 0.79709 0.8084 0.82086 0.83372 0.84514
0.85814 0.8707 0.88222 0.8954 0.90889 0.92094
0.93321 0.94714 0.95987 0.97504 0.97616];
V=[0.96286 1.01286 1.06286 1.11286 1.16286 1.21286 1.26286 1.31279 1.36279 1.41273 1.46273 1.51273 1.56273 1.61273 1.66273 1.71273 1.76273 1.81273 1.86273 1.91273 1.96273 2.01273 2.06273 2.11273 2.16273 2.21273 2.26273 2.31273 2.36273 2.41273 2.46273 2.51273 2.56273 2.61273 2.66273 2.71273 2.76273 2.81273 2.86273 2.91273 2.96273 3.01273 3.06273 3.11273 3.16273 3.21273 3.26273 3.31273 3.36273 3.41273 3.46273 3.51273 3.51474];
plot(H2,V,'*',H2,(0.89*2.45/0.6*((H2-0.6).*(sqrt(H2.*(1.2-H2)))+0.36*asin(H2/0.6-1)+0.5*pi*0.36)))
title('原始图和拟合图')
xlabel('\it高度/m'),ylabel('\it体积/m^3')

附件3：
H2=[0.35209 0.36425 0.37913 0.39134 0.4047 0.41854 0.43017 0.44336 0.45549 0.46764 0.47968 0.49276 0.5052 0.51736 0.52954 0.54036 0.55242 0.56424 0.57638 0.58708 0.59939 0.61102 0.62143 0.63383 0.64547 0.65725 0.66846 0.68019 0.6917 0.70235 0.71423 0.72619 0.73684 0.74907 0.7616 0.7736
0.78527 0.79709 0.8084 0.82086 0.83372 0.84514
0.85814 0.8707 0.88222 0.8954 0.90889 0.92094
0.93321 0.94714 0.95987 0.97504 0.97616];
V=[0.96286 1.01286 1.06286 1.11286 1.16286 1.21286 1.26286 1.31279 1.36279 1.41273 1.46273 1.51273 1.56273 1.61273 1.66273 1.71273 1.76273 1.81273 1.86273 1.91273 1.96273 2.01273 2.06273 2.11273 2.16273 2.21273 2.26273 2.31273 2.36273 2.41273 2.46273 2.51273 2.56273 2.61273 2.66273 2.71273 2.76273 2.81273 2.86273 2.91273 2.96273 3.01273 3.06273 3.11273 3.16273 3.21273 3.26273 3.31273 3.36273 3.41273 3.46273 3.51273 3.51474];
plot(H2,V,'.',H2,(0.89*2.45/0.6*((H2-0.6).*(sqrt(H2.*(1.2-H2)))+0.36*asin(H2/0.6-1)+0.5*pi*0.36))-0.0772)
title('模型改进后的原始图和拟合图')
xlabel('\it高度/m'),ylabel('\it体积/m^3')

附件4：
H1=0:0.01:0.14;
A=0.89/0.6*(H1+0.4*tan(4.1/180*pi))*cot(4.1/180*pi);
B=0.5*(H1+0.4*tan(4.1/180*pi))-0.6;
C=sqrt(0.5*(H1+0.4*tan(4.1/180*pi)).*(1.2-0.5*(H1+0.4*tan(4.1/180*pi))));
D=0.36*asin((H1+0.4*tan(4.1/180*pi))/1.2-1)+0.5*pi*0.36;
V=A.*(B.*C+D)

附件5：
H1=0.15:0.01:1.17;
H2=H1+(0.4-2.45/2)*tan(4.1/180*pi);
V=0.89*2.45/0.6*((H2-0.6).*(sqrt(H2.*(1.2-H2)))+0.36*asin(H2/0.6-1)+0.5*pi*0.36)-0.0772

附件6：
H1=1.18:0.01:1.20;
H2=0.5*(1.2-H1+(2.45-0.4)*tan(4.1/180*pi));
V=pi*0.89*0.6*2.45-0.89*2.45/0.6*((H2-0.6).*(sqrt(H2.*(1.2-H2)))+0.36*asin(H2/0.6-1)+0.5*pi*0.36)
1. 基于C8051F单片机直流电动机反馈控制系统的设计与研究
2. 基于单片机的嵌入式Web服务器的研究
3. MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究
4. 基于模糊控制的电阻钎焊单片机温度控制系统的研制
5. 基于MCS-51系列单片机的通用控制模块的研究
6. 基于单片机实现的供暖系统最佳启停自校正（STR）调节器
7. 单片机控制的二级倒立摆系统的研究
8. 基于增强型51系列单片机的TCP/IP协议栈的实现
9. 基于单片机的蓄电池自动监测系统
10. 基于32位嵌入式单片机系统的图像采集与处理技术的研究
11. 基于单片机的作物营养诊断专家系统的研究
12. 基于单片机的交流伺服电机运动控制系统研究与开发
13. 基于单片机的泵管内壁硬度测试仪的研制
14. 基于单片机的自动找平控制系统研究
15. 基于C8051F040单片机的嵌入式系统开发
16. 基于单片机的液压动力系统状态监测仪开发
17. 模糊Smith智能控制方法的研究及其单片机实现
18. 一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制
19. 基于双单片机冲床数控系统的研究
20. 基于CYGNAL单片机的在线间歇式浊度仪的研制
21. 基于单片机的喷油泵试验台控制器的研制
22. 基于单片机的软起动器的研究和设计
23. 基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究
24. 基于单片机的机电产品控制系统开发
25. 基于PIC单片机的智能手机充电器
26. 基于单片机的实时内核设计及其应用研究
27. 基于单片机的远程抄表系统的设计与研究
28. 基于单片机的烟气二氧化硫浓度检测仪的研制
29. 基于微型光谱仪的单片机系统
30. 单片机系统软件构件开发的技术研究
31. 基于单片机的液体点滴速度自动检测仪的研制
32. 基于单片机系统的多功能温度测量仪的研制
33. 基于PIC单片机的电能采集终端的设计和应用
34. 基于单片机的光纤光栅解调仪的研制
35. 气压式线性摩擦焊机单片机控制系统的研制
36. 基于单片机的数字磁通门传感器
37. 基于单片机的旋转变压器-数字转换器的研究
38. 基于单片机的光纤Bragg光栅解调系统的研究
39. 单片机控制的便携式多功能乳腺治疗仪的研制
40. 基于C8051F020单片机的多生理信号检测仪
41. 基于单片机的电机运动控制系统设计
42. Pico专用单片机核的可测性设计研究
43. 基于MCS-51单片机的热量计
44. 基于双单片机的智能遥测微型气象站
45. MCS-51单片机构建机器人的实践研究
46. 基于单片机的轮轨力检测
47. 基于单片机的GPS定位仪的研究与实现
48. 基于单片机的电液伺服控制系统
49. 用于单片机系统的MMC卡文件系统研制
50. 基于单片机的时控和计数系统性能优化的研究
51. 基于单片机和CPLD的粗光栅位移测量系统研究
52. 单片机控制的后备式方波UPS
53. 提升高职学生单片机应用能力的探究
54. 基于单片机控制的自动低频减载装置研究
55. 基于单片机控制的水下焊接电源的研究
56. 基于单片机的多通道数据采集系统
57. 基于uPSD3234单片机的氚表面污染测量仪的研制
58. 基于单片机的红外测油仪的研究
59. 96系列单片机仿真器研究与设计
60. 基于单片机的单晶金刚石刀具刃磨设备的数控改造
61. 基于单片机的温度智能控制系统的设计与实现
62. 基于MSP430单片机的电梯门机控制器的研制
63. 基于单片机的气体测漏仪的研究
64. 基于三菱M16C/6N系列单片机的CAN/USB协议转换器
65. 基于单片机和DSP的变压器油色谱在线监测技术研究
66. 基于单片机的膛壁温度报警系统设计
67. 基于AVR单片机的低压无功补偿控制器的设计
68. 基于单片机船舶电力推进电机监测系统
69. 基于单片机网络的振动信号的采集系统
70. 基于单片机的大容量数据存储技术的应用研究
71. 基于单片机的叠图机研究与教学方法实践
72. 基于单片机嵌入式Web服务器技术的研究及实现
73. 基于AT89S52单片机的通用数据采集系统
74. 基于单片机的多道脉冲幅度分析仪研究
75. 机器人旋转电弧传感角焊缝跟踪单片机控制系统
76. 基于单片机的控制系统在PLC虚拟教学实验中的应用研究
77. 基于单片机系统的网络通信研究与应用
78. 基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究
79. 基于单片机的模糊控制器在工业电阻炉上的应用研究
80. 基于双单片机冲床数控系统的研究与开发
81. 基于Cygnal单片机的μC/OS-Ⅱ的研究
82. 基于单片机的一体化智能差示扫描量热仪系统研究
83. 基于TCP/IP协议的单片机与Internet互联的研究与实现
84. 变频调速液压电梯单片机控制器的研究
85. 基于单片机γ-免疫计数器自动换样功能的研究与实现
86. 基于单片机的倒立摆控制系统设计与实现
87. 单片机嵌入式以太网防盗报警系统
88. 基于51单片机的嵌入式Internet系统的设计与实现
89. 单片机监测系统在挤压机上的应用
90. MSP430单片机在智能水表系统上的研究与应用
91. 基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用
92. 单片机在高楼恒压供水系统中的应用
93. 基于ATmega16单片机的流量控制器的开发
94. 基于MSP430单片机的远程抄表系统及智能网络水表的设计
95. 基于MSP430单片机具有数据存储与回放功能的嵌入式电子血压计的设计
96. 基于单片机的氨分解率检测系统的研究与开发
97. 锅炉的单片机控制系统
98. 基于单片机控制的电磁振动式播种控制系统的设计
99. 基于单片机技术的WDR-01型聚氨酯导