1第一章数据建模常规方法的MATLAB实现精选.ppt

第一章 数据建模常规方法的MATLAB实现 1.1 数据的读入与写出 1.2 数据拟合方法 1.3 数据拟合应用实例 1.4 数据的可视化 1.1 数据的读入与写出 数学建模不可避免地要用到大量的数据，最简单的方法是复制、粘贴，但是不方便。另一种方法是与Excel和记事本（*.dat或者*.txt的文件）进行交互。 （算 法） 程 序 数据读入 输出 图形 数据 1.1.1 Excel与MATLAB的数据交互 首先要安装Excel和MATLAB。 第一，打开Excel的工具→宏→安全性→安全级（中） 第二，打开Excel的工具→加载宏→浏览→ 安装MATLAB的目录→toolbox →exlink →excllink.xla →确定，得到如下的工具条（即可使用）： 1.1.2 记事本与MATLAB的数据交互 %从记事本t.txt中读取数据 [name,type, x,y,answer]=textread(t.txt,%s Type %n %f %n %s,2) %将Matlab数据写入记事本 fid= fopen(tp.txt,wt); %文件扩展名可以为：*.dat或*.xls fprintf(fid,This is the database of class 1.\n); name=Sally;types=1;x=3.1;y=45;answer=Yes; fprintf(fid,% s Type %u %f %u %s \n,name,types,x,y,answer); name=Tom;types=2;x=2.5;y=20;answer=No; fprintf(fid,% s Type %u %f %u %s \n,name,types,x,y,answer); fclose(fid); 1.2 数据拟合方法 1.2.1 多项式拟合 一般情况下，数据点较少的用插值法。数据点较多的，只需要考察数据的总体变化趋势的，用拟合法。最常用的拟合方法是最小二乘拟合法。 1.多项式拟合指令 x=[1 2 3 4 5 6 7 8 9]; y=[9 7 6 3 -1 2 5 7 20]; P=polyfit(x,y,3); xi=0:0.2:10; yi=polyval(P,xi); plot(xi,yi,x,y,r*); 2.图形窗口的多项式拟合 自学（我的教材缺页！） 例子：d122 1.2.2 指定函数拟合 MATLAB主窗口左下角： start→toolboxes→CurveFitting→Curve Fitting Tool (cftool) 1.2.3 曲线拟合工具箱 1.3 数据拟合应用实例 1.3.1 人口预测模型 某地区1971-2000年人口数据，给出该地区人口增长的数学模型。 根据所给数据作出散点图。人口随时间的变化是非线性的，存在一条与x轴平行的渐近线，因此用Logistic曲线模型进行拟合。 程序见d131 1.3.2 薄膜渗透率的测定 程序见d132 背景 年 1625 1830 1930 1960 1974 1987 1999 人口(亿) 5 10 20 30 40 50 60 世界人口增长概况 中国人口增长概况 年 1908 1933 1953 1964 1982 1990 1995 2000 人口(亿) 3.0 4.7 6.0 7.2 10.3 11.3 12.0 13.0 研究人口变化规律 控制人口过快增长 1.3.3 如何预报人口的增长 指数增长模型——马尔萨斯提出 (1798) 常用的计算公式 x(t) ~时刻t的人口 基本假设 : 人口(相对)增长率 r 是常数 今年人口 x0, 年增长率 r k年后人口 随着时间增加，人口按指数规律无限增长. 与常用公式的一致 rt e x t x 0 ) ( = ? 指数增长模型的应用及局限性 与19世纪以前欧洲一些地区人口统计数据吻合. 适用于19世纪后迁往加拿大的欧洲移民后代. 可用于短期人口增长预测. 不符合19世纪后多数地区人口增长规律. 不能预测较长期的人口增长过程. 19世纪后人口数据 人口增长率r不是常数(逐渐下降) 阻滞增长模型(Logistic模型) 人口增长到一定数量后，增长率下降