Matlab在小型温度检测系统中的应用

 

　　1 前言

　　防潮、防霉、防腐、防爆是仓库日常工作的重要内容，是衡量仓库管理质量的重要指标。它直接影响到储备物资的使用寿命和工作可靠性。为保证日常工作的顺利进行，首要问题是加强仓库内温度与湿度的监测工作。但传统的方法是用与湿度表、毛发湿度表、双金属式测量计和湿度试纸等测试器材，通过人工进行检测，对不符合温度和湿度要求的库房进行通风、去湿和降温等工作。这种人工测试方法费时费力、效率低，且测试的温度及湿度误差大，随机性大。因此我们需要一种造价低廉、使用方便且测量准确的温湿度测量仪。

　　本实验基于Matlab环境下设计了一个小型温度检测系统，下位机使用AT89S51单片机和DS18B20完成温度数据采集，上位机在Matlab环境下，调用设备控制箱serial类操作RS-232串口，用串行通信方式交换数据，进而借助Matlab对数据进行分析和处理，得到了温度随时间变化的函数解析式，同时介绍了基于Matlab环境下PC机与单片机串行通信的实时数据处理的实现方法。

　　2 系统总体设计

　　
图1  系统结构图

　　温度检测系统的整体结构如图1所示。PC机串口与单片机USART口通过MAX232电平转换芯片相连，构成一个主从式通信系统。系统工作时，单片机对串口和DS18B20初始化，在读取温度的同时等待中断。PC机通过调用Matlab设备控制工具箱中的serial类及相关函数来创建串口设备对象，并以读写文件的方式实现对PC机串行口的访问，PC机通过Matlab向串行口发送特殊指令从而触发单片机中断系统，单片机调用中断服务例程，读取即使温度并将采集的数据通过串行口回送给PC机。此时，Matlab通过查询的方式，实时接收单片机发送的数据，并完成对数据的分析处理及图形显示。

　　MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

　　3 下位机部分

　　下位机部分由AT89S51单片机和DS18B20温度传感器构成，主要负责温度数据的采集工作，并通过串行通信实时地将数据传送到上位机进行处理，PC机与MUC串口通信技术相对而言已经比较成熟。

　　3.1串行通信协议

　　串口通信协议SPCP（Serial Port Communication Protocol）设计思想是基于帧传输方式，在本实验中，设定字符格式为1个起始位，8个数据位和一个停止位，无奇偶校验，中间8位即为有效数据，波特率设置为9600,为保证数据可靠传输，在传送数据前通过握手建立连接，软件握手协议规定如下：

　　上位机发送握手信号0xff给下位机，下位机如果接受到上位机的信号为握手信号，则回送数据包给上位机，其中个数据为握手信号，以二个数据为温度传感器采集到的温度数据，此时，上位机如果接受到的个数据不是握手信号，则丢弃该数据包，若是，则表示握手成功，直接存储第二个数据。

　　3.2温度数据采集（DS18B20）

　　本系统中采用DALLAS生产的"一线总线"可编程数字化温度传感器DS18B20,与微处理器连接时仅需要一条口线即可实现微处理器与DS18B20的双向通讯，在使用中不需要任何外围元件，设计可用数据线供电，简化系统的硬件，同时支持多点组网功能，多个DS18B20可以并联在惟一的三线上，实现多点测温，使用起来非常方便。尽管如此，DS18B20是以牺牲软件资源换取硬件资源的，由于采用单总线数据出数方式，DS18B20的数据I/O均由同一根线完成，因此，对读写的操作时序要求非常严格。

　　根据DTASHEET,对DS18B20的编程主要注意以下几个方面：

　　1,延时问题[2]:为了保证DS18B20的严格时序，可以将延时分为2种：10us以下的短延时和10us以上的长延时。短延时可以使用C51提供的内部函数_nop_（）来实现，一个nop（）函数相当于一条DJNZ汇编指令，约2us;长延时主要有15us,90us,270us,540us等，这些延时均为15us的整数倍，可以使用nop（）函数编写一个延时15us的函数delay15（n）。

　　2,基本操作：DS18B20的一线工作协议流程是：初始化→ROM操作指令→存储器操作指令→数据传输。对DS18B20 进行所有的读写操作都是从初始化开始的，主要分为初始化操作，读操作和写操作。

　　基于以上分析，对DS18B20的编程源码如下：

　　Init_DS18B20（void）  //初始化函数

　　{DQ = 1;    //DQ复位

　　Delay15（1）；   //稍做延时

　　DQ = 0;    //单片机将DQ拉低

　　Delay15（32）；  //延时 大于 480us

　　DQ = 1;    //拉高总线

　　Delay15（6）； //延时90us

　　x=DQ;      //读存在脉冲

　　delay15（20）； //延时约270us

　　}

　　Write_DS18B20（unsigned char dat）//写一个字节

　　{

　　unsigned char i="0";//定义循环变量

　　for （i=0; i<8; i++）

　　{

　　DQ = 0; //复位

　　DQ = dat&0x01;//取数据的第i位并送出

　　Delay15（1）；//延时

　　DQ = 1; //停止

　　dat》=1;//右移

　　}

　　}

　　Read_DS18B20（void）//读一个字节

　　{

　　unsigned char i="0";

　　unsigned char dat = 0;

　　for （i=8;i>0;i--）

　　{

　　DQ = 0; // 复位

　　dat》=1;

　　DQ = 1; // 给脉冲信号

　　if（DQ）

　　dat|=0x80;//取位脉冲并存入dat

　　delay15（1）；//延时

　　}

　　return（dat）；

　　}

　　4 上位机部分（PC）

　　上位机通过串口向下位机发送命令实现对下位机的控制，并实时地接受下位机传送过来的数据，对其分析处理，将结果用图形显示并存储，完成人机交互过程。

　　Matlab并不具备直接访问硬件的能力，但是支持面向对象技术，通过调用Instrument Control Toolbox中的serial类函数来创建串口对象，对串口对象操作就是对串口操作，使用起非常方便。同时，Matlab封装的串口对象支持对串口的异步读写操作，使得计算机在读写串口时能同时进行其他处理工作，因而能大大提高计算机执行效率。Matlab用多线程技术实现这种异步操作，通过异步读写设置，计算机在执行读写串口函数时能立即返回不必等待串口把数据传输完毕， 当指定的数据传输结束时就触发事件，执行事件回调函数，可以在事件回调函数中编程，进行数据处理，这样就不会造成因等待串口传输数据引起的时间浪费。

　　4.1 Matlab下串口编程

　　MATLAB的Instrument Control Toolbox提供了 MATLAB与仪器仪表通信的功能 ,它支持 GPIB 通用接口总线 、VISA、TCP/ IP、UDP、RS2232等多个协议 ,具有同步和异步读写功能以及事件处理和回调操作功能，可读写和记录二进制和ASCII文本数据。与串口有关的主要函数如下[3]:

　　（1）建立串口对象函数：obj=seril（'port','property name',propertyvalue……），其中主要的属性有：baudrate（波特率），databits（数据位），parity（校验方式），stopbits（终止位）等，可以在初始化时进行赋值或者使用set函数。

　　（2）打开串口设备对象：fopen（obj）

　　（3） 串口读写操作：当matlab通信数据采用二进制格式时，读写串口设备的命令为fread（）和fwrite（）；当通行数据采用文本（ASCII）格式时，读写串口设备的命令为fscanf（）和fprintf（）。

　　（4）关闭并清除设备对象：

　　Fclose（obj）；%关闭串口设备对象

　　Delete（obj）；%删除内存中的串口设备对象

　　基于本系统串口通信协议，对串口对象的读写部分程序如下：

　　Obj=serial（'com1','baudrate',9600,'parity','none','databits',8,'stopbits',1）；%初始化串口

　　Fopen（obj）；%打开串口对象

　　Fwrite（obj,256）；%向串口发送握手信号0xff

　　TMP=fread（obj,3,'unit8'）；%从串口读取3字节数据，后2个即是16bit温度数据

　　If TMP（1）= =256   %判断个字节受否是握手信号

　　For i = 1:3

　　Dat（i）=TMP（i+1）；%剔除个握手信号字节

　　End

　　End

　　Fclose（obj）；%关闭串口设备对象

　　Delete（obj）；%删除内存中的串口设备对象

　　4.2 Matlab数据分析

　　单片机一般能处理简单的8位无符号数的四则运算，而DS18B20可以程序设定9~12位的分辨率，可达±0.5℃，温度以16bit带符号位扩展的二进制补码形式读出，如果使用单片机进行快速的实时处理则比较费力，同时单片机还要与DS18B20及上位机通信，系统资源也比较紧张。因而可以将读取的16bit温度数据直接送往PC机，由上位机来完成。Matlab强大的计算能力和绘图功能给数据分析带来了极大的方便，这不仅可以合理利用系统资源，也使得系统的通信过程更流畅。

　　图2所示的是用DS18B20测得的连续10个温度数据的变化曲线图。使用max（）和min（）函数可以求出温度的极大极小值，调用polyfit（）函数还可以进行二乘回归分析与曲线拟合，进而求出温度变化的解析式，本实验中使用二阶拟合后，得到的温度随时间变化的函数解析式为：F（T）= 0.0125*T2 - 0.0145T+27.3083.

　　5 结语

　　MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特定类型的问题。

　　本文作者创新观点：使用DS18B20传感器简化了硬件结构，采用Matlab编程，充分利用其数据分析能力，合理利用系统资源，提高了开发效率，扩展了Matlab的使用范围。