湖水温度变化模型
湖水温度变化模型
前言: 夏季湖水温度有明显的正温层现象水的下层温度较低,秋季因湖区多风而发生湖水搅动,使水温分层温度现象基本消失,冬季湖面结冰,湖水温度出现逆温层现象。特别近几年来全球变暖越来越严重,这对夏季时湖水的温度的变化也照成了一定的影响,使得湖水照成水文变化异常的现象,影响了河中生物的生存与繁衍。使得水层上下循环不畅,造成下层水域缺氧,导致水生鱼类的死亡。论文利用数学建模理念和MATLAB软件对水温的变化与分布进行了分析和评论MATLAB、Word
题型关键词:深度、温度、多项式拟合问题重述
湖水在夏天会出现分层现象,其特点为接近湖面的水温度较高,越往下温度变低。这种上热下冷的现象影响了水的对流和混合过程,使得下层水域缺氧,导致水生鱼类的死亡。下面是某个湖的观测数据水层之间互不影响。符号说明:
h: 湖水深度,单位为m;
T: 湖水温度,单位为°C;
T T h :湖水深度的函数;MATLAB基本语句及解释:常用清空命令:(1)clear :清除内存中所有变量;(2)clc:清屏
常用数学函数:
三角函数:正弦 sin x ,余弦 cos x ,双曲正弦 sinh x ,反正弦 asin x ,反双曲正弦 asinh x ,正切 tan x ,余弦cot x ,正割 sec x ,余割 csc x 等;
指数函数:exp x , log x , log10 x ,log2 x , 平方根 sqrt x 等;
整值函数:朝0方向取整 fix x ,朝负无穷方向取整 floor x ,朝正无穷方向取整 ceil x ,四舍五入到最近的整数 round x ,符号函数 sign x 等
数组操作函数:siza A , size A,1 , size A,2 , length A ,flipud A ,fliplr A 等
(5)‘nearest’:最近邻点插值,直接完成计算; ‘spline’:三次样条函数插值;
‘linear’:线性插值(缺省方式),直接完成计算; ‘cubic’:三次函数插值;
(6) x [a b c d e] 创建包含指定元素的行向量; y [a b c d e]’ 求该矩阵的转置; plot x,y,’s’ 画实线,s为线型; a polyfit x,y,n 返回多项式系数从最高次系数到最底次系数,n是多项式的阶数 polyfit [a b c d],[e f g h],n 求其拟合曲线函数方程系数 b regress Y,X MATLAB统计工具箱; [b,bint,r,rint,statas] regress Y,X,alpha
3.问题的分析与建立模型思路
遇到这种数据表格问题,如果我们仅凭眼睛观察,很难看到其中的规律,也就更难写出有效的数学表达式从而建立数学模型问题给出了有限实验数据,想到用插值和多项式拟合的方法来解。将所给数据作图,横轴代表湖水深度,纵轴代表湖水温度,用MATLAB数学软件画出散点图,操作的命令假设湖水深度是温度的连续函数,其中一组统计数据为:x:湖水深度单位为my:湖水温度单位为 C,它是湖水深度的函数:y f x 这里多项式拟合的方法MATLAB中实现编程,先求出湖水温度函数y然后求出的拟合函数就可以MATLAB数学软件画出散点图,操作的命令为: x [0 2.3 4.9 9.1 13.7 18.3 22.9 27.2]; y [22.8 22.8 22.8 20.6 13.9 11.7 11.1 11.1]; plot x,y,''
图4-1
得到如上散点图观察散点图的特点,由散点图可知,图形观察散点图的特点,并通过实验选取不同的基函数类进行实验发现用4次多项式拟合比较好,相应的命令为次拟合命令 x [0 2.3 4.9 9.1 13.7 18.3 22.9 27.2]; y [22.8 22.8 22.8 20.6 13.9 11.7 11.1 11.1]; polyfit [0 2.3 4.9 9.1 13.7 18.3 22.9 27.2],[22.8 22.8 22.8 20.6 13.9 11.7 11.1 11.1],4
ans -0.0001 0.0103 -0.2279 0.9977 22.3743
拟合曲线函数方程为y *x.^4+0.0103*x.^3-0.2279*x.^2+0.9977*x+22.3743;
显示拟合函数命令: x [0:0.1:27.2];
y *x.^4+0.0103*x.^3-0.2279*x.^2+0.9977*x+22.3743; plot x,y
图4-2
为求出湖水在10m处的温度,键入命令: a polyfit x,y,4 ; polyval