MATLAB数学实验报告-海水温度随深度变化

1、数学实验报告实验名称Matlab大作业学院材料科学与工程专业班级姓名2010年6月【实验目的】运用数学知识与MATLAB相结合，运用数学方法，建立数学模型，用MATLAB软件辅助求解模型，解决实际问题。【实验任务】海水温度随着深度的变化而变化，海面温度较高，随着深度的增加温度越来越低，这样也就影响了海水的对流和混合，使得深层海水中的氧气越来越少，这是潜水员必须考虑的问题，同时根据这个规律也可以对海水鱼层作一划分。现在通过实验观测得一组海水深度h与温度的t的数据如下：th0要求：找出温度t与深度h之间的一个近似函数关系；找出温度变化最快的深度位置（实际上该位置就是潜水员在潜水时，随着下潜深度的不

2、同，需要更换呼入气体种类的位置，也是不同种类鱼层的分界位置）。【实验分析及求解】模型准备：该问题是要运用已测数据，研究海水温度随深度的变化关系。找出温度随深度变化的函数关系，并利用此近似关系，找到温度变化最快的深度位置，进而为潜水员潜水时换气和鱼类分层的分界位置提供理论依据和实践指导。模型分析及求解：（1）建立海水温度与深度变化的函数关系：实际上，海水的温度变化不仅仅是由深度决定的，可能还与海水的盐度、密度、比热等因素有关，但深度对海水的温度起决定性作用，因此可以假设海水的温度和深度之间存在简单的函数关系。题目中给出了一组海水深度h与温度的t的数据，因此可对该组数据进行多项式的拟合，找出合适的

3、近似关系，得出海水温度随深度变化的规律。th0为了找出符合深度h和温度t的合适多项式，于是对其进行了3、5、7阶的拟合，以下为拟合后的图形:深度H05101520253035404550深度H当然，拟合的次数越高，在一定区间内吻合程度就会越高，但对计算也会造成一定的难度。对比以上图形，可以发现，只有在030m深度范围内，拟合出的曲线才于实际情况相符合。同时，七阶曲线的拟合程度已经基本可以表达二者之间的关系：因此，选用7阶函数最为二者之间的近似函数关系.拟合的函数关系式如下:T二0.167404x10-4H7-0.117471x10-2H6+0.309730x10-1H5-0.388967H4+

4、2.40949H3-6.78428H2+5.51977H+23.5000因为所取数据点是在027m之间，同时由拟合曲线可以发现在其他区间内，确实与实际情况不符，故取区间为（030）。（2）利用拟合曲线求温度随深度变化最快的位置：dTd2TdT即求的最大值，进而可以转化为求二阶导数的极值即=0，对应最大的极值点即为所求的位置。dHdH2dH鴛=0-1171828x10-3H6-0-704826x10-2H5+0-1548650H4-1-555868H3+7-22847H2-13-56856H+5-51977=0.7030968x10-3H5-0.3524130x10-1H4+0.6194600H

5、3-4.667604H2+14.45694H-13.56856dH2令dH=0解得dH的极值点，对应最大dH的点即为所求点。dH2dHdH用MATLAB解得的结果为：由dT=0得:dH2H=dT代入dH得:dTdH4.,dT显然，日=时，对应的二122.709127151164101883805173491dH最大，所以H=即为所求点。即温度随深度变化最快的位置为：H=3.模型检验：针对求出的温度随深度变化的近似函数关系，用所给数据进行检验：测量值t测量值h0计算值t可以发现，数据基本吻合，可见近似关系符合实际.【实验程序及结果】Matlab程序：%曲线拟合clft=;h=0;plot(h,t

6、,p),holdony=polyfit(h,t,7);T=poly2str(y,H)h=0:1:50;y1=polyval(y,h);plot(h,y1)xlabel(深度H),ylabel(温度T)legend(数据点,拟合曲线)title(7阶拟合曲线)运行结果：T=HA7-HA6+HA5-HA4+HA3-HA2+H+%最值求取symsHT=*HA7-*HA6+*HA5-*HA4+*HA3-*HA2+*H+;dT=diff(T,H)dT2=diff(T,H,2)Tj=solve(dT2,H)Tj1=subs(dT,H,Tj)运行结果：dT=(8281*HA6)/-(9959*HA5)/052+(*HA4)/0-(38897*HA3)/25000+(14457*HA2)/2000-(67843*H)/5000+27599/5000dT2=(5843*HA5)/6464-(8795*HA4)/052+(*HA3)/9007-(116691*HA2)/25000+(14457*H)/1000-67843/5000Tj=Tj1=-122.【实验总结】1海水温度随深度变化的近似函数关系为：0,30)T=0.167404x1