基于Matlab实现热敏电阻温度特性实验的数据处理

官华夏

(东莞市第六高级中学 广东 东莞 523000)

1 引言

用非平衡电桥研究热敏电阻的温度特性曲线是普通物理实验中一个比较重要的电磁学实验，学生在进行实验操作后，对于实验数据的处理与分析，是一个比较困难的过程．实验数据曲线拟合及分析的传统方法，是采用最小二乘法和坐标纸作图法．这种手工做法，计算量大，精度不高．本文根据最小二乘法原理并结合数值分析方法，将曲线化直进行分析，利用Matlab软件编程计算对实验数据进行线性拟合，从而得出热敏电阻的温度特性．

Matlab软件是由Mathworks公司发布的主要面对科学计算、数据可视化、系统仿真及交互程序设计的高科技计算环境．由于其功能强大，而且简单易学，Matlab软件已经成为高校教师、科研人员和工程技术人员的必学软件，它能够极大地提高工作效率和质量．我们知道，Matlab一词是矩阵实验室(Matrix Laboratory)的简称，与其他高级语言相比，Matlab程序编写简单，计算高效，且提供了大量的专业工具箱，便于专业应用．特别是工具箱中的Matlab数学函数库包括了大量的计算方法，让笔者体现到了其强大的数学计算功能[1,2]．

2 实验方法与原理

利用作图法研究半导体负温度系数热敏电阻随温度变化的规律，负温度系数热敏电阻Rt随热力学温度T的变化成指数函数关系

(1)

式中A和B为常数．为了能够更精确的测量，采用了将曲线化直的方法，对上式两边取对数得

采用非平衡电桥测量热敏电阻，电路如图1所示．

图1 非平衡电桥测量热敏电阻电路图

当电桥不平衡时，电流计有电流Ig流过，我们用支路电流法求出Ig与热敏电阻Rt的关系，根据基尔霍夫定律列式，最后可得

Rt=

(2)

从式(2)以及(1)式可以看出，Ig与Rt以及Rt与T都是一一对应的，也就是Ig与T有着确定的关系．因此我们可以通过测出Ig与T的关系来推出Rt与T的关系，即测得热敏电阻的温度特性[3,4]．

原始实验数据如表1所示.

表1 用非平衡电桥测热敏电阻实验数据

3 利用Matlab软件进行实验数据处理

通过实验原理的学习，可以知道Rt与T成指数关系

利用Matlab编程求出系数B和lnA的最小二乘估计，并拟合出曲线．

3.1 利用Matlab进行编程

打开Matlab7.0界面，在命令窗口输入命令语句：

t=input('t=');

ig=input('Ig=');%以数组的方式输入相应实验数据，即[...]的形式

接着对“t”和“ig”的数组大小进行定义，获取行数和列数，调用size函数：

m=size(t);

n=size(ig);

以下是主要的程序内容：

for i=1:n(2) %从第一个数据到最后一个连续运算

r(i)=(5000*1950*0.9-ig(i)*(5000*

1950*1950+2156*5000*1950+2156*5000*

1950))/(ig(i)*(5000*1950+1950*1950+1950*5000+2156*1950+2156*1950)+1950*0.9); %根据公式和相对应系数得出电阻值

end

for i=1:m(2)

tf(i)=1/t(i);

end

for i=1:n(2)

lnrt(i)=log(r(i));

end

f=@(x)a(1)*x+a(2);%其中a(1)是B，a(2)是lnA

g=@(y)exp(a(1)/y+a(2));

subplot(2,2,1);%将显示框分为两行两列，共两个图，创建第一个子图

fplot(f,[0.0026,0.0034],'k');%调用fplot函数绘制图形

hold on;

hold on;

subplot(2,2,2);%创建第二个子图

fplot(g,[300,380],'k');

hold on;

plot(t,r,'*k');% 绘制T-Rt关系曲线

hold off;

r1=corrcoef(tf,lnrt);%拟合直线的相关系数

3.2 实验数据处理

保存之后按enter运行，分别以数组的形式输入对应的实验数据，可以得到

B=3.348 3×103

lnA=-0.002 7×103

r1=0.999 7

即拟合后的直线是

拟合直线的相关系数是

r1=0.999 7

并得到

Rt=0.5419 0.6005 0.6972 0.7940 0.8982

1.0151 1.1820 1.3521 1.5544 1.9096

2.1286 2.4728 3.0046 3.6178 4.3226

(×103Ω)

0.0029 0.0029 0.0030 0.0030 0.0030

0.0031 0.0031 0.0032 0.0032 0.0033

(K-1)

lnRt=6.2951 6.3977 6.5471 6.6771 6.8004

6.9228 7.0750 7.2094 7.3489 7.5546

7.6632 7.8131 8.0079 8.1936 8.3716

由此也可以得出该热敏电阻随温度的变化规律为

进一步可以计算出该热敏电阻的激活能

E=4.62×10-20J

图2 Matlab拟合的ln Rt-直线图

图3 Matlab拟合的Rt-T曲线图

图中的“*”表示实验数据点，实线表示拟合直线，从图2、图3中可以看出，所测量的实验数据没有粗大误差，实验数据点较合理地分布在拟合直线近旁，直线拟合相关系数r1约为1，因此用线性回归函数比较合理．Matlab软件编程过程简单明了，学生容易上手并学会使用．

3 结束语

随着智能化时代的来临，大量计算软件应运而生．Matlab摆脱了传统非交互式程序语言(如C，Fortran)的编辑模式，将数值分析、矩阵计算、数据可视化等诸多强大功能集成在一个易于使用的视窗环境中，已成为教学和科研中最基础的软件之一．

利用Matlab软件进行实验数据处理，过程简单方便、准确度高、直观可视，让学生从繁重的数据处理中寻找到了帮助，这也大大减小了学生的学业负担．学生在使用计算机软件解决实际问题的过程中在，既锻炼了能力，也激发了对物理实验的兴趣．这不仅培养了学生对所学知识的综合运用，也为后续的科研训练打下良好的基础．当然，对实验数据的处理也要分层次，对简单的数据分析，笔者还是建议学生独立对数据进行分析处理、作图计算，对于复杂的实验可以采用Matlab软件进行处理．

Matlab软件，作为一种通用软件，使用者需要掌握一定的计算机编程知识和矩阵知识，也需要熟悉其中的函数和命令，是一种高级的基于矩阵、数组的语言，使用者既可以用它来快速编写简单的程序，也可以编写庞大复杂的应用程序．编程语言作为人类和计算机交流的通用语言必将会是未来人类主要的语言之一，这也让笔者认为介绍Matlab软件是相当有必要的．