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

《数学实验》报告 实验名称?? Matlab大作业 学 院 材料科学与工程 专业班级?????材料080?? 姓??? 名????????? ?? 2010年6月 【实验目的】 运用数学知识与MATLAB相结合，运用数学方法，建立数学模型，用MATLAB软件辅助求解模型，解决实际问题。 【实验任务】 海水温度随着深度的变化而变化，海面温度较高，随着深度的增加温度越来越低，这样也就影响了海水的对流和混合，使得深层海水中的氧气越来越少，这是潜水员必须考虑的问题，同时根据这个规律也可以对海水鱼层作一划分。现在通过实验观测得一组海水深度h与温度的t的数据如下： t 23.5 22.9 20.1 19.1 15.4 11.5 9.5 8.2 h 0 1.5 2.5 4.6 8.2 12.5 16.5 26.5 要求： 找出温度t与深度h之间的一个近似函数关系； 找出温度变化最快的深度位置（实际上该位置就是潜水员在潜水时，随着下潜深度的不同，需要更换呼入气体种类的位置，也是不同种类鱼层的分界位置）。 【实验分析及求解】 模型准备 模型 因为所取数据点是在0~27m之间，同时由拟合曲线可以发现在其他区间内，确实与实际情况不符，故取区间为（0~30）。 （2）利用拟合曲线求温度随深度变化最快的位置： 即求 的最大值，进而可以转化为求二阶导数的极值即，对应最大的极值点即为所求的位置。 令解得的极值点，对应最大 的点即为所求点。 用MATLAB解得的结果为： 由得： H= 1.6197811681049637254053846911074 4.7013743475793350231729782711085 8.7169265460971614409156533926428 13.417167700861440255839235695813 21.667720502617871613364782131819 代入得： -3.1156124188526912076380861099649, 0.5570442013360668953956637191788, -3.216335574505053725682056116640, 4.47086321730533459443013039762, -122.709127151164101883805173491 显然，H=21.667720502617871613364782131819时，对应的 最大，所以H=21.67即为所求点。 即温度随深度变化最快的位置为：H=21.67m 3..模型检验： 针对求出的温度随深度变化的近似函数关系，用所给数据进行检验： 测量值t 23.5 22.9 20.1 19.1 15.4 11.5 9.5 8.2 测量值h 0 1.5 2.5 4.6 8.2 12.5 16.5 26.5 计算值t 23.5000 22.9000 20.1001 19.1004 15.4026 11.5148 9.5530 8.6817 可以发现，数据基本吻合，可见近似关系符合实际. 【实验程序及结果】 Matlab程序： %曲线拟合 clf t=[23.5 22.9 20.1 19.1 15.4 11.5 9.5 8.2]; h=[0 1.5 2.5 4.6 8.2 12.5 16.5 26.5]; plot(h,t,p),hold on y=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 = 1.674e-005 H^7 - 0.0011747 H^6 + 0.030973 H^5 - 0.38897 H^4 + 2.4095 H^3 - 6.7843 H^2 + 5.5198 H + 23.5 %最值求取 syms H T =1.674e-005* H^7 - 0.0011747 *H^6 + 0.030973 *H^5 - 0.38897* H^4 +