基础物理实验研究性报告热学系列

1、基础物理实验研究性报告基础物理实验研究性报告热学系列实验第一作者：第二作者： 2012年12月15日摘要在基础实验中做了热学系列的两个实验后，深感数据处理计算难度之大。故此次研究性实验我选择了热学系列，对实验原理，实验装置以及后期的数据处理做了研究，并尝试利用现代技术手段和专业知识来解决实验问题与实际问题。关键词：热学，数据处理，现代技术，专业知识，实际问题目录1实验要求52实验原理52.1测定冰的熔解热实验52.1.1一般概念52.1.2装置简介62.1.3实验原理72.2电热法测量焦耳热功当量实验102.2.1一般说明102.2.2散热修正103实验仪器114实验内容114.1测定冰的熔解

2、热实验114.1.1合理选择实验参量114.1.2记录有关常数124.1.3测定实验过程中系统温度随时间的变化124.1.4数据处理134.2电热法测量焦耳热功当量实验134.2.1测量各种质量134.2.2测量时间-温度关系134.2.3测量加热器的电功率134.2.4数据处理145数据处理145.1原始数据145.1.1测定冰的熔解热实验145.1.2电热法测量焦耳热功当量实验155.2数据处理155.2.1测定冰的熔解热实验155.2.2电热法测量焦耳热功当量实验186对实验改进的探究216.1基础实验中的已有的改进216.1.1实验装置、原理及仪器的改进216.2可以继续改进的部分21

3、6.2.1对孤立系统的改进216.2.2对实验器材的建议227对大量数据处理的探究237.1利用Excel编程实现数据的统计处理237.1.1测定冰的熔解热实验237.1.2电热法测量焦耳热功当量实验247.2利用高级语言程序设计处理数据277.2.1测定冰的熔解热实验277.2.2电热法测量焦耳热功当量实验297.3对数据处理的分析：348总结358.1本次研究性实验总结与感想358.2对本学期基础物理实验课程的总结与感想358.3对本课程的小建议369.参考文献371实验要求1.1熟悉热学实验中的基本问题量热和计温；1.2研究电热法中作功与传热的关系；1.3学习两种进行散热修正的方法牛顿冷

4、却定律和一元线性回归法；1.4了解热学实验中合理安排实验和选择参量的重要性；1.5熟悉热学实验中基本仪器的使用。2实验原理2.1测定冰的熔解热实验2.1.1一般概念一定压强下晶体开始熔解时的温度，也就是该物质的固态和液态可以平衡共存的温度，称为该晶体物质在此压强下的熔点。对于晶体而言，熔解是组成物质的粒子由规则排列向不规则排列的过程，破坏晶体的点阵结构需要能量，因此，晶体在熔解过程中虽吸收能量，但其温度却保持不变。单位质量的晶体物质在熔点时从固态全部变成液态所需的热量，叫做该晶体物质的熔解潜热，亦称熔解热。本实验用混合量热法测定冰的熔解热。其基本做法如下：把待测系统A和一个已知热容的系统B混合

5、起来，并设法使它们形成一个与外界没有热量交换的孤立系统C（C=A+B）.这样A（或B）所放出的热量，全部为B(或A)所吸收。因为已知热容的系统在实验过程中所传递的热量Q，是可以由其温度的改变T和热容Cs计算出来，即Q=CsT，因此待测系统在实验过程中所传递的热量也就知道了。由此可见，保持系统为孤立系统，是混合量热法所要求的基本条件。这要从仪器装置、测量方法以及实验操作等各方面去保证。如果这样做以后，实验过程中与外界的热交换仍不能忽略，就要进行散热或吸热修正。温度是热学中的一个基本物理量，量热实验中必须测量温度。一个系统的温度，只有在平衡时才有意义，因此计温时必须使系统各处温度达到均匀。用温度计

6、的指示值代表系统温度，必须使系统与温度计之间达到热平衡。2.1.2装置简介为了使实验系统（包括待测系统与已知其热容的系统）成为一个孤立系统，本实验采用了量热器。热量传递有三种方式：传导、对流和辐射。因此，热学实验应使系统与环境之间的传导、对流和辐射都尽量简小，量热器可以近似满足这样的要求。量热器的种类很多，随测量的目的、要求、测量精度的不同而异，最简单的一种如图所示，它由良导体做成的内筒放在一较大的外筒中组成。通常在内筒中放水、温度计及搅拌器，它们（内筒、温度计、搅拌器及水）连同放进的待测物体就构成了我们所考虑的（进行实验的）系统，内筒、水、温度计和搅拌器的热容是可以计算出来或实测得到的，在此

7、基础上，就可以用混合法进行量热实验了。内筒置于一绝热架上，外筒用绝热盖盖住，因此空气与外界的对流很小，又因为空气是不良导体，所以内、外筒间靠传导方式传递的热量同样可以减至很小，同时由于内筒的外壁及外筒的内外壁都电镀得十分光亮，使得它们发射或接收辐射热的本领变得很小，于是试验系统和环境之间因辐射而产生的热量传递也得以减小，这样的量热器就可以使实验系统粗略地接近于一个孤立系统了。2.1.3实验原理若有质量为M、温度为T1的冰（在实验室环境下其比热容为c1，熔点为T0），与质量为m、温度为T2的水（比热容为c0）混合，冰全部熔解为水后的平衡温度为T3，设量热器的内筒和搅拌器的质量分别为m1、m2，比

8、热容分别为c1c2，温度计的热容为m。如果实验系统为孤立系统，将冰投入盛水的量热器中，则热平衡方程式为c1m(T0-T1)+ML+c0M(T3-T0)=(c0m+c1m1+c2m2+m)(T2-T3)式中，L为冰的熔解热。在本实验的条件下。冰的熔点也可认为是0，即T0=0，所以冰的熔解热为L=1M(c0m+c1m1+c2m2+m)(T2-T3)-c0T3+c1T1为了尽可能使系统与外界交换的热量达到最小，除了使用量热器以外，实验的操作过程中也必须予以注意，例如不应当直接用手去把握量热器的任何部分；不应当在阳光的直接照射下或空气流动太快的地方（如通风孔道、风扇旁边）进行实验；冬天要避免在火炉或暖

9、气旁做实验等。此外，由于系统与外界温度差越大时，在它们之间传递热量越快，而且时间越长，传递的热量越多，因此在进行量热实验时，要尽可能使系统与外界温度差小，并尽量使实验过程进行得迅速。尽管注意到了上述的各个方面，系统仍不可能完全达到绝热的要求（除非系统与环境的温度时时刻刻完全相同）。因此，在作精密测量时，就需要采用一些办法来求出实验过程中实验系统究竟散失或吸收了多少热量，从而对实验结果进行修正。一个系统的温度如果高于环境温度，它就要散失热量。实验证明，当温度差相当小时（例如不超过1015），散热速率与温度差成正比，此即牛顿冷却定律，用数学形式可写成qt=K(T-)式中q是系统散失的热量；t是时间

10、间隔；K是散热常数，与系统表面积成正比，并随表面的吸收或发射辐射热的本领而变；T、分别是所考虑的系统及环境的温度；qt称为散热速率，表示单位时间内系统散失的热量。下面介绍一种根据牛顿冷却定律粗略修正散热的方法。已知当T>时，qt>0，系统向外散热；当T<时，qt<0，系统从环境吸热。可以取系统的初温T2>，终温T3<，以设法使整个实验过程中系统与环境间的热量传递前后彼此抵消。考虑到实验的具体情况，刚投入冰时，水温高，冰的有效面积大，熔解快，因此系统表面温度T（即量热器中水温）降低较快；随后，随着冰的不断熔化，冰块逐渐变小，水温逐渐降低，冰熔解变缓，水温的降低

11、也就变慢起来。量热器中水温随时间变化曲线如图所示。实验过程中，即温度从T2到T3这段时间（t2t3）内系统与环境间交换的热量为q=t2t3KT-dt=Kt2tKT-dt+Ktt3KT-dt前一项T->0，系统散热对应于图中面积SA=t2tT-dt；后一项T-<0，系统吸热，对应面积SB=tt3T-dt。不难想见，面积SA与系统向外界散失的热量成正比，即q吸=KSB，K是散热常数。因此，只要SASB，系统对外界的吸热和散热就可以相互抵消。要使SASB，就必须使(T2-)>(-T3)，T2和T3的选取要在实验中根据具体情况选定。上述这种使散热与吸热相互抵消的做法，不仅要求水的初温

12、比环境温度高，末温比环境温度低，而且对初温、末温与环境温度相差的幅度要求比较严格，往往经过多次试做，效果仍可能不理想。因此希望把上述思想进行扩展，放宽对量热器中水的初温和末温的限制。如图所示，在t=t2时投入冰块，在t=t3时冰块熔化完毕。在投入冰块前，系统的温度沿T2T2变化；在冰块熔化后，系统温度沿T3T3变化。T2T2和T3T3实际上都很接近直线。作T2T2的延长线到T2，作T3T3的延长线到T3，连接T2T3，使T2T3与T轴平行，且使面积S1+S2=S3，用T2代替T2，用T3代替T3，代入公式求L，就得到系统与环境没有发生热量交换的实验结果。实际的温度变化本来是T2T2T4T3T3

13、，在从冰块熔化完毕的过程中，系统散失的热量相当于面积S4，从环境吸收的热量相当于面积S2+S5，综合两者，系统共吸收的热量相当于面积S=S2+S5-S4。在用T2代替T2、用T3代替T3后，得到另一条新的温度曲线T2T2T2T3T3T3。从冰块投入到冰块熔化完毕的过程中，系统散失的热量相当于面积S1+S4，从环境吸收的热量等于S3+S5。综合两者，系统共吸收的热量相当于面积S=S3+S5-S1-S4。因为作图时已使S1+S2=S3，所以有S=S。这说明，新的温度曲线与实际温度曲线是等价的。新的温度曲线的物理意义是，它把系统与环境交换热量的过程与冰熔化的过程分隔开来，从T2到T2和从T3到T3是

14、系统与环境交换热量的过程，从T2到T3是冰熔化的过程。由于冰熔化的过程变为无限短自然没有机会进行热量交换，因而从T2到T3，便仅仅是由于冰的熔化而引起的水温变化。这一方法把对热量的修正转换为对初温和末温的修正，且对量热器中水的初温和末温原则上没有任何限制。尽管如此，考虑到牛顿冷却定律成立的条件以及其他因素，T2、T3还是选择在附近为好，即让T2>，T3<，但它们与的差值可以不受限制。2.2电热法测量焦耳热功当量实验2.2.1一般说明如图所示，给电阻R两端加上电压V，通过R的电流为I，在通电时间t内电场力作功W=VIt。若这些功全部转化为热量，使一个盛水的量热器系统初温0升高至，系统

15、吸收热量为Q，则热功当量J=W/Q。按照能量守恒定律，若采用国际单位制，则W和Q的单位都是焦耳（J），比值J=1；若Q用(cal)做单位，则J=4.1868J/cal，表示产生1卡热量所需作的功。实验在装水的量热筒中进行。系统吸收热量为Q=(c0m0+c1m1+c2m2)(-0)=Cm(-0)式中，c0、c1、c2分别是水、量热装置及加热器的比热容；m0、m1、m2分别是其相应的质量；Cm=c0m+c1m1+c2m2是系统的总热容；0为系统初温。本实验的主要内容就是测定热功当量J=VIt/Cm(-0)。2.2.2散热修正本实验的难点是如何考虑系统散热的修正。我们从系统应满足的微分方程出发。若把

16、系统看成是理想绝热的，即只考虑系统由于通电而升温，则由系统吸热方程Q=Cm(-0)对时间求导可以得到温度变化率所满足的关系式为ddt|吸=VIJCm考虑通电时系统吸热的同时也向环境中放热，根据牛顿冷却定律，由于放热引起的温度变化率为ddt|吸=-K(-环)式中，K为系统的散热系数。综合吸热、放热效应，系统温度的实际变化率为ddt=VIJCm=-K(-环)这是一个一阶线性的常系数微分方程。我们试图利用一元线性回归法处理数据，令yddt，x(-环），上式变为y=a+bx，其中，a=VIJCm，b=-K。给加热系统通电，并同时记录系统温度-时间的变化关系，每隔1min记录一次温度，共测30个连续时间

17、对应的温度值。根据测量出的数据，用差分代替微分计算t=ti+ti+12时的=i+i+12、ddt=i+1-iti+1-ti，这样由一系列（ti，i）就换算出了（yi，xi）的数据了，带入回归系数计算求得a，从而由下式计算出热功当量J（式中R是加热用的电阻值），即a=V2RJCm J=V2aRCm3实验仪器量热器、电子天平、温度计、数字三用表、加温热皿、冰、水桶、停表、干拭布等。4实验内容4.1测定冰的熔解热实验4.1.1合理选择实验参量一个成功的实验应能测量出投冰前的降温曲线和冰块全部熔化后的升温曲线，且系统终温T3低于环境温度（温差不超过15）。影响实验结果的参量有水的质量m0、水的初温T2

18、以及冰的质量M，而这些参量的大小事互相制约的，需要先定出它们的取值范围，再通过实验进行调整。首先，冰块的大小是基本固定的，可根据量热筒的大小选择投放一块或两块冰。其次，确定水的初温T2。一般选择T2高于环境温度1015，因为此时散热服从牛顿冷却定律，便于对系统散热进行粗略修正。最后，当M与T2确定后要想调整实验结果，只有通过改变水的质量m0来实现了。水的质量不宜太大，水多需要的冰块就多，否则测不出升温曲线；水也不能太少，太少不利于搅拌，且会使系统终温T3过低。可取量热筒内筒的1/22/3进行试探性实验，如果未能测出升温曲线，或终温T3低于室温15以上，则需要改变水量重做实验。4.1.2记录有关

19、常数称量各种质量。注意冰不能直接放在天平盘上称衡，冰的质量应由冰熔解后，冰加水的质量减去水的质量求得。已知实验室所用内筒和搅拌器材料为铜，比热容c1=c2=0.389×103J/(kg·K)，冰的比热容（-400时）为c1=1.80×103J/(kg·K)，水的比热容为c0=4.18×103J/(kg·K)，忽略温度计的热容m。4.1.3测定实验过程中系统温度随时间的变化每隔一定时间测系统温度，作T-t图。注意：i三部分曲线是连续的，时间不可间断。特别要记录好投冰的时间。ii正确使用和保护温度计。iii整个实验过程中要不断地轻轻进行搅

20、拌，以确保温度计度数代表所测系统的温度。实测系统的散热常数K量热器盛适量水，水温比环境温度低510，测量系统温度随时间的变化。4.1.4数据处理用第二种散热修正的方法，作图求出初、末温度的修正值，并算出冰的熔解热L。由测量数据估算系统的散热常数K。4.2电热法测量焦耳热功当量实验4.2.1测量各种质量水的质量不宜过大或过小，一般控制在200240g为好。加热器由功率电阻组成，搅拌器主要由铝质叶片组成，两者的总热容可按64.38J/L计算。4.2.2测量时间-温度关系在连续升温的30min内，应等间隔地读取31个温度值（每分钟1次）。注意：升温过程中必须不断搅拌（转动搅拌器叶片）以保证温度均匀。

21、同时搅拌过程中要随时监视电源电压（面板电压表指针位置）是否改变，防止因搅拌动作过大引起电源接触不良。数字三用表有自动关机功能。因此在测量过程中，可在三用表工作接近15min时，进行一次关机-开机操作，以免读数时刚好自动关机。用铂电阻温度计记录温度，可直接把输出的香蕉插头接入数字三用表并读取电阻值。4.2.3测量加热器的电功率分别在读数始末，用数字三用表测出加热器两端的电压（注意三用表的插孔位置和量程选择）。加热器电阻值如下表所列。编号12345678电阻值/202.4201.5203.8200.5201.1199.6201.4203.4编号910111213141516电阻值/201.3201

22、.7200.4201.9200.8201.7201.6200.84.2.4数据处理用一元线性回归方法计算热功当量J并与理论值对比，计算它们的相对误差。要求自行编写计算机程序来处理一元线性回归问题，并讨论相关系数。5数据处理此处为基础实验部分的原始数据及数据处理，后面研究性部分会对大量数据的处理作进一步的探究。5.1原始数据5.1.1测定冰的熔解热实验5.1.2电热法测量焦耳热功当量实验5.2数据处理5.2.1测定冰的熔解热实验5.2.2电热法测量焦耳热功当量实验6对实验改进的探究6.1基础实验中的已有的改进6.1.1实验装置、原理及仪器的改进1.内筒放在绝热架上，与外筒之间以空气相隔，外筒用绝

23、热盖盖住，内筒外壁及外筒内外壁电镀的十分光亮，构成一个粗略的孤立系统。2.将牛顿冷却定律的修正方法扩展，把对热量的修正转换为对初温和末温的修正。3.用数字三用表和铂电阻温度计代替普通温度计。4.用电子天平代替普通物理天平。6.2可以继续改进的部分6.2.1对孤立系统的改进虽然我们在基础实验中按照原理所述建立了一个粗略的绝热系统，但是由于空气的作用，内筒与外筒之间还是存在着热量的交换。而且，对于绝热盖和外筒的温度也不相同，所以系统相对于内筒的环境温度（内外筒之间的空气温度）仍然不恒定，且在绝热盖到外筒底部分布不均匀。所以，我们可以对系统装置做进一步改进。1.我们可以将内外筒之间填充比空气更觉热的

24、介质，例如：适量纸屑、绝热型不锈钢固体保护渣、三元乙丙橡胶等。不过这个需要根据我校现有的实验资金及条件来考虑。所以在基础实验中，我们可以采用已有的方式进行；而在探究设计实验的过程中，这也不免是一种改进方法。2.即使添加了绝热介质，也不能保证系统是完全的孤立系统。所以，进一步的改进可以由图所示的装置来实现。该装置将绝热盖换为与内外筒材料相同的金属盖并电镀的很光亮，同时将外筒更换为可以进行流水控温的环境温度保持装置，并添加搅拌棒，使外筒内的水温恒定，避免出现上下温度不均的情况。此装置添加部分实现略为困难，所以同样不能作为基础实验中的批量改进，仅限于对实验研究的一个建议。6.2.2对实验器材的建议用

25、铂电阻温度计代替普通水银温度计已经解决了温度测量的一大问题。但是，由于搅拌不均匀导致的温度测量偏差依然是实验结果误差的主要原因之一。因此，对于搅拌棒的频率、幅度等控制显得尤为关键。我个人认为，在研究性实验的过程中，可以将搅拌棒提前半固定在绝热盖上，这样可以减少在投冰之后封盖之前的时间间隔。同时，搅拌棒的手控端可以连接一个固定频率的机械运动装置，带动搅拌棒以固定频率上下移动或水平转动，以保证搅拌过程时刻均匀。而搅拌棒的搅拌幅度由对应环形搅拌棒的半径决定。同时，为了避免由机械振动产生热量影响实验结果，仅限于搅拌棒的手控绝缘端与定频率的机械振动装置相连，不可应用类似于家用的自动打蛋器或豆浆机类型的搅

26、拌器，因为它们会额外产生大量的热量影响实验结果。7对大量数据处理的探究在热学系列的这两个实验中，有一个共同点就是要处理数量偏多的实验数据。该系列实验中的数据虽然数目偏多、计算量偏大，但是具有一定的关系，例如，冰的熔解热实验结果为两端近线性线段与一段非线性线段的合成，而电热法测量焦耳热功当量实验所得数据处理过程需要利用一元线性回归来完成。在基础实验2的内容要求中，提到了可以自行利用计算机编程来作为中间过程处理一元线性回归问题。对此，我想到其实实验1也可以利用计算机编程来做一些数据的处理工作。所以下面的研究，就是利用计算机提供的几种典型编程方式（以Excel编程与高级语言程序设计为例），来对这两个

27、实验的数据处理作出探讨与比较。7.1利用Excel编程实现数据的统计处理Excel作为数据处理的重要工具之一，已成为广大统计工作者不可或缺的统计工具。Excel表格中潜藏着各种复杂函数、各种图像与表格数据之间的联系，简化了统计数据的比较，是现代简单的统计学处理的霸王式工具。在基础实验中，我就已经利用Excel对电热法测量焦耳热功当量实验的数据进行了处理。此次研究性实验再次将其拿出来，是为了用Excel继续处理冰的熔解热实验的数据，并将其处理过程与后续要研究的C语言编程处理作以比较。7.1.1测定冰的熔解热实验利用铂电阻温度计电阻与温度的换算公式R=R0(1+AT+BT2)式中A=3.90802

28、×10-3-1，B=-5.80195×10-7-1，R0=1000。推导出由电阻通过计算得到温度的公式为T=-A+A2-4B(1-RR0)2B因此，表格中的电阻-温度转换公式如下所列：D1=-(SQRT(3.90802*3.90802/1000000+4*5.80195/10000000*(1-C1)-3.90802/1000)/2/5.80195*10000000（室温）D4=-(SQRT(3.90802*3.90802/1000000+4*5.80195/10000000*(1-C4)-3.90802/1000)/2/5.80195*10000000D5=-(SQRT(

29、3.90802*3.90802/1000000+4*5.80195/10000000*(1-C5)-3.90802/1000)/2/5.80195*10000000D6=-(SQRT(3.90802*3.90802/1000000+4*5.80195/10000000*(1-C6)-3.90802/1000)/2/5.80195*10000000D24=-(SQRT(3.90802*3.90802/1000000+4*5.80195/10000000*(1-C24)-3.90802/1000)/2/5.80195*10000000再附加公式/=+273.15/K，便可求得每一记录时刻的液体温度

30、。所得表格如附表一所示（数据同前文基础实验数据）。利用Excel编程功能对表中数据进行图形绘制。选取时间t为横轴，温度T为纵轴，描点方式绘出散点图，并在表格选项中添加圆滑线将其连接。绘图结果在附表一的图表（sheet2）中。后续计算方法与基础实验中的数据处理相同。7.1.2电热法测量焦耳热功当量实验由于该实验内容中注明可以用计算机编程实现，故将原始数据与数据处理（非前文基础实验数据）列出如下图。由Excel编程实现的表格如附表二所示。电阻与温度的换算与实验1相同，另外还包含了对于一元线性回归中各参量的求解过程。(i+i+1)/2：F4=(E4+E5)/2；F5=(E5+E6)/2；F33=(E

31、33+E34)/2。y=d/dt：G4=(F5-F4)/60；G5=(F6-F5)/60；G32=(F33-F32)/60。x=-环：H4=F4-环；H5=F5-环；H32=F32-环。x2：I4=H4*H4；I5=H5*H5；I32=H32*H32。xy：J4=H4*G4；J5=H5*G5；J32=H32*G32。y2：K4=G4*F4；K5=G5*G5；K32=G32*G32。求均值：G1=AVERAGE(G4:G32)；H1=AVERAGE(H4:H32)；I1=AVERAGE(I4:I32)；J1=AVERAGE(J4:J32)；K1=AVERAGE(K4:K32)。求和：L2=SUM

32、(L4:L32)。以上给出的是变量的求解方法，所有求解结果均在附表二中。7.2利用高级语言程序设计处理数据个人认为，在处理实验数据的过程中，能够灵活运用本专业知识是问题的处理得到简化，才能够使得学有所用。所以，我在此次研究性实验中利用所学的C语言程序设计，做了完整的数据处理。处理结果如下。7.2.1测定冰的熔解热实验编程源代码如下：#include<stdio.h>#include<math.h> #define A 3.90802/1000#define B -5.80195/10000000#define R 1000/这里的A B R为铂电阻转化为对应温度所需的常

33、数值#define C0 4.18*1000/水的比热容#define C1 0.389*1000/内筒与搅拌棒的比热容#define CI 1.8*1000/冰的比热容#define T0 273.15/冰的熔点double rtot(double);int main()double m1,mi,m,m2,mt,rte,t2,t3,ti,l;printf("请输入室温对应电阻r(K)=");scanf("%lf",&rte);printf("当前实验室内温度为环=%lfn",rtot(rte);printf("请输

34、入内筒质量m1(g)=");scanf("%lf",&m1);printf("请输入内筒加水的总质量m1+m(g)=");scanf("%lf",&mt);m=mt-m1;printf("则水的质量为m(g)=%.2fn",m);printf("请输入内筒、水和冰的总质量m1+m+mi(g)=");scanf("%lf",&mt);mi=mt-m-m1;printf("则冰的质量为mi(g)=%.2fn",mi);prin

35、tf("请输入搅拌器质量m2(g)=");scanf("%lf",&m2);printf("请输入冰的温度TI()=");scanf("%lf",&ti);printf("原始数据处理可由其他方法进行，画出图像后可读得对应的温度修正值如下：n");printf("T2'()=");scanf("%lf",&t2);printf("T3'()=");scanf("%lf",&am

36、p;t3);l=1/mi*(C0*m+C1*m1+C1*m2)*(t2-t3)-C0*t3+CI*ti;printf("冰的熔解热为L=1/mi*(C0*m+C1*m1+C1*m2)*(T2'-T3')-C0*T3'+CI*TI=%.2lfJ/gn",l/1000);return 0;double rtot(double r)/利用铂电阻温度计工作原理将所示电阻值转化为对应的温度double t;t=(-A+sqrt(A*A-4*B*(1-r*1000/R)/(2*B);return t;程序的执行界面如下：7.2.2电热法测量焦耳热功当量实验编程源

37、代码如下：#include<stdio.h>#include<math.h>#define MAX 50#define N 31/本实验中共测得31组数据#define A 3.90802/1000#define B -5.80195/10000000#define R 1000/这里的A B R为铂电阻转化为对应温度所需的常数值#define C0 4.18/水的比热容#define C1 0.389/内筒的比热容#define C2 64.38/加热器和搅拌器的总热容double rtot(double);int main()int i,rt;double v,r,

38、rte,rte1,mt,m,m1,xavg,yavg,xxavg,yyavg,xyavg,a,b,r1,cm,j,uy,ub,ua,uj,tmp=0;double te1MAX,teMAX,teavgMAX,yMAX,xMAX,xxMAX,xyMAX,yyMAX,sum10=0;printf("请输入电压U(V)=");scanf("%lf",&v);printf("请输入加热器电阻编号：");scanf("%d",&rt);switch(rt)/根据加热器电阻的编号查询所用电阻阻值case 1:r

39、=202.4;break;case 2:r=201.5;break;case 3:r=203.8;break;case 4:r=200.5;break;case 5:r=201.1;break;case 6:r=199.6;break;case 7:r=201.4;break;case 8:r=203.4;break;case 9:r=201.3;break;case 10:r=201.7;break;case 11:r=200.4;break;case 12:r=201.9;break;case 13:r=200.8;break;case 14:r=201.7;break;case 15:r

40、=201.6;break;case 16:r=200.8;break;default:printf("error");return 1;printf("所采用的加热器电阻阻值为R=%.1lfn",r);printf("请输入室温对应电阻r(K)=");scanf("%lf",&rte1);rte=rtot(rte1)+273.15;printf("当前实验室内温度为环=%lfKn",rte);printf("请输入内筒质量m1(g)=");scanf("%l

41、f",&m1);printf("请输入内筒加水的总质量m1+m(g)=");scanf("%lf",&mt); m=mt-m1;printf("则水的质量为m(g)=%.2fn",m);printf("请依次输入31组数据，该实验中数据两两之间间隔为60sn");for(i=0;i<N;i+)scanf("%lf",&te1i);tei=rtot(te1i)+273.15;if(i=0)continue;teavgi-1=(tei-1+tei)/2;if(

42、i=1)continue;yi-2=(teavgi-1-teavgi-2)/60;sum0+=yi-2;xi-2=teavgi-2-rte;sum1+=xi-2;xxi-2=xi-2*xi-2;sum2+=xxi-2;xyi-2=xi-2*yi-2;sum3+=xyi-2;yyi-2=yi-2*yi-2;sum4+=yyi-2;yavg=sum0/(N-2);/y的均值xavg=sum1/(N-2);/x的均值xxavg=sum2/(N-2);/x*x的均值xyavg=sum3/(N-2);/x*y的均值yyavg=sum4/(N-2);/y*y的均值printf("利用一元线性回归

43、处理数据，计算系数a b及相关系数rn");b=(xavg*yavg-xyavg)/(xavg*xavg-xxavg);/计算ba=yavg-b*xavg;/计算ar1=(xyavg-xavg*yavg)/sqrt(xxavg-xavg*xavg)*(yyavg-yavg*yavg);/计算相关系数rprintf("b=%.5lf*10(-4)ta=%.5lf*10(-3)n相关系数r=%.8lfn",b*10000,a*1000,r1);cm=C2+C0*m+C1*m1;printf("Cm=c0m0+c1m1+c2m2=%.4lfn",cm

44、);j=v*v/a/r/cm;/计算焦耳热功当量Jprintf("测得焦耳热功当量为J=U*U/(a*R*Cm)=%.2lfn",j);printf("不确定度的计算：n");for(i=0;i<N-2;i+)tmp+=(yi-(a+b*xi)*(yi-(a+b*xi);uy=sqrt(tmp/(N-4);/y的不确定度ub=b*sqrt(1/r1/r1-1)/(N-4);/b的不确定度ua=sqrt(xxavg)*ub;/a的不确定度uj=v*v/r/cm*ua/a/a;/j的不确定度printf("u(y)=%.8lf,u(b)=%.

45、5lf*10(-5),u(a)=%.5lf*10(-4)nu(J)=%.7lfn",uy,ub*100000,ua*10000,uj);printf("最终计算结果为：J±u(J)=%.2lf±%.2lfn",j,uj);return 0;double rtot(double r)/利用铂电阻温度计工作原理将所示电阻值转化为对应的温度double t;t=(-A+sqrt(A*A-4*B*(1-r*1000/R)/(2*B);return t;程序的执行界面如下：7.3对数据处理的分析：在数据处理过程中，如上述过程中处理数据过程显得清晰明了。但

46、是美中不足的是，我现有的水平还不足以利用C语言程序设计来构造函数图像。所以，我认为，可以将Excel表格对于数据的处理与高级语言程序设计相结合，用Excel来处理表格数据及图像，用高级语言程序设计来进行复杂的计算，可以以最高效率利用我们的计算机解决实验中的问题。另外，在我用C语言编写处理数据的过程中，没有考虑具体的数据数值，应用的都是实验中的理论数据，包含了物理量之间的换算，待求量的计算以及不确定度的计算等。所有的结果是在输入数据中自动生成的，所以，这是我对于实验的理论公式理解进一步加深。这也是我在本次研究性实验中的收获之一。8总结8.1本次研究性实验总结与感想本次研究性实验主要目的有三个：第一，加深对基础实验操作及数据处理的理解；第二，分析基础实验中已经做到的改进，通过与同组同学探讨及查阅相关资料得出可以继续改进的部分并分析其可行性；第三，通过基础实验中的数据处理提示用计算机编程做过程处理，来进一步探讨可执行的程序处理方式，以及将常用的处理方式（以Excel和C语言编程为例）加以结合，分析计算机处理为我们的实验带来的方便等。通过本次研究性实验，我们对冰的熔解热实验和电热法测量焦耳热功