插值与拟合演示文稿.pptVIP

To MATLAB hd1 返回 4.作出水深小于5的海域范围,即z=5的等高线. 2.在矩形区域(75,200)×(-50,150)进行插值。 1.输入插值基点数据 3. 作海底曲面图 当前第31页\共有59页\编于星期三\1点 %程序一：插值并作海底曲面图 x =[129.0 140.0 103.5 88.0 185.5 195.0 105.5 157.5 107.5 77.0 81.0 162.0 162.0 117.5 ]; y =[ 7.5 141.5 23.0 147.0 22.5 137.5 85.5 -6.5 -81 3.0 56.5 -66.5 84.0 -33.5 ]; z =[ 4 8 6 8 6 8 8 9 9 8 8 9 4 9 ]; x1=75:1:200; y1=-50:1:150; [x1,y1]=meshgrid(x1,y1); z1=griddata(x,y,z,x1,y1,v4); meshc(x1,y1,z1) 当前第32页\共有59页\编于星期三\1点 海底曲面图 当前第33页\共有59页\编于星期三\1点 %程序二：插值并作出水深小于5的海域范围。 x1=75:1:200; y1=-50:1:150; [x1,y1]=meshgrid(x1,y1); z1=griddata(x,y,z,x1,y1,v4); %插值 z1(z1=5)=nan; %将水深大于5的置为nan，这样绘图就不会显示出来 meshc(x1,y1,z1) 当前第34页\共有59页\编于星期三\1点 水深小于5的海域范围 当前第35页\共有59页\编于星期三\1点 实验作业1 山区地貌：在某山区测得一些地点的高程如下表：(平面区域1200≤x ≤4000,1200≤y ≤3600)，试作出该山区的地貌图和等高线图，并对几种插值方法进行比较． 当前第36页\共有59页\编于星期三\1点 5.1 引言 对于情况较复杂的实际问题（因素不易化简，作用机理不详）可直接使用数据组建模，寻找简单的因果变量之间的数量关系， 从而对未知的情形作预报。这样组建的模型为拟合模型。 拟合模型的组建主要是处理好观测数据的误差，使用数学表达式从数量上近似因果变量之间的关系。拟合模型的组建是通过对有关变量的观测数据的观察、分析和选择恰当的数学表达方式得到的。 五、拟合的使用及求解 当前第37页\共有59页\编于星期三\1点 5.2 拟合模型的分类 5.2.1 直线拟合 5.2.2 曲线拟合 5.2.3 观察数据修匀 对于已给一批实测数据，由于实测方法、实验环境等一些外界因素的影响，不可避免地会产生随机干扰和误差。我们自然希望根据数据分布的总趋势去剔除观察数据中的偶然误差，这就是所谓的数据修匀（或称数据平滑）问题。 当前第38页\共有59页\编于星期三\1点 直 线 拟 合 问 题 引 例 1 温度t(oC) 20.5 32.7 51.0 73.0 95.7 电阻R(?) 765 826 873 942 1032 已知热敏电阻数据： 求60oC时的电阻R． 设 R=at+b a,b为待定系数 当前第39页\共有59页\编于星期三\1点 曲 线 拟 合 问 题 引 例 2 t (h) 0.25 0.5 1 1.5 2 3 4 6 8 c (?g/ml) 19.21 18.15 15.36 14.10 12.89 9.32 7.45 5.24 3.01 已知一室模型快速静脉注射下的血药浓度数据(t=0注射300mg) 求血药浓度随时间的变化规律c(t). 在直角坐标系下作图如下(plot) MATLAB(aa1) 当前第40页\共有59页\编于星期三\1点 曲 线 拟 合 问 题 的 提 法 已知一组（二维）数据，即平面上 n个点（xi,yi) i=1,…,n, 寻求一个函数（曲线）y=f(x), 使 f(x) 在某种准则下与所有数据点最为接近，即曲线拟合得最好． + + + + + + + + + x y y=f(x) (xi,yi) ?i ?i 为点（xi,yi) 与曲线 y=f(x) 的距离 当前第41页\共有59页\编于星期三\1点 曲线拟合问题最常用的解法——线性最小二乘法的基本思路 第一步:先选定一组函数 r1(x), r2(x), …,rm(x), mn, 令