数学建模插值法与曲线拟合讲课.pptVIP

第30页，共62页，星期日，2025年，2月5日二维插值的MATLAB实现2、另一个是离散数据节点的插值命令：z=griddata(x0,y0,z0,x,y,’method’)其中，z：被插值点处的函数值；x0,y0,z0：插值节点，x0,y0,z0均为向量；x,y：被插值点；method：插值方法，包括：linear——线性插值；cubic——三次插值；第31页，共62页，星期日，2025年，2月5日船在该海域会搁浅吗？在某海域测得一些点(x,y)处的水深z由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）*（-50，150）里的哪些地方船要避免进入.第32页，共62页，星期日，2025年，2月5日解决问题的步骤：1.作出测量点的分布图；2.求出矩形区域（75，200）*（-50，150）的细分网格节点之横、纵坐标向量；3.利用MATLAB中的散点插值函数求网格节点的水深；4.作出海底曲面图形和等高线图；5.作出水深小于5的海域范围.第33页，共62页，星期日，2025年，2月5日程序clearx0=[129140103.588185.5195105157.5107.57781162162117.5];y0=[7.5141.52314722.5137.585.5-6.5-81356.5-66.584-33.5];z0=[-4-8-6-8-6-8-8-9-9-8-8-9-4-9];subplot(2,2,1)plot(x0,y0,+);%作出测量点的分布图；x=75:1:200;%加密y=-50:1:150;[x,y]=meshgrid(x,y);z=griddata(x0,y0,z0,x,y,cubic);subplot(2,2,2)mesh(x,y,z),%用插值方法求出网格节点处的z坐标矩阵,绘制出三维图形subplot(2,2,3)meshc(x,y,z),%绘制等高线subplot(2,2,4)contour(x,y,z,[-5-5]);%水深5英尺处海底曲面的等高线gridon第34页，共62页，星期日，2025年，2月5日第35页，共62页，星期日，2025年，2月5日拟合的标准（1）用各点误差绝对值的和表示（2）用各点误差按绝对值的最大值表示（3）用各点误差的平方和表示第36页，共62页，星期日，2025年，2月5日最小二乘拟合式中R2称为均方误差。由于计算均方误差的最小值的原则容易实现而被广泛采用。按均方误差达到极小构造拟合曲线的方法称为最小二乘法。第37页，共62页，星期日，2025年，2月5日++++++++++++++++++++++++++++++p=a1+a2xp=a1+a2x+a3x2p=a1+a2x+a3x2p=a1+a2/xp=aebxp=ae-bx将数据(xi,yi)i=1,…,n作图，通过直观判断确定p(x)：第38页，共62页，星期日，2025年，2月5日MATLAB---曲线拟合工具箱Matlab有一个功能强大的曲线拟合工具箱（CurveFittingToolbox）cftool，使用方便，能实现多种类型的线性、非线性曲线拟合。调用：cftool界面如下所示第39页，共62页，星期日，2025年，2月5日第40页，共62页，星期日，2025年，2月5日“Data”按钮数据的选取点击“Data”按钮，弹出“Data”窗口；利用Xdata和Ydata的下拉菜单读入数据x,y，可修改数据集名“Datasetname”，然后点击“Createdataset”按钮，退出“Data”窗口，返回工具箱界面，这时会自动画出数据集的曲线图；

第41页，共62页，星期日，2025年，2月5日第42页，共62页，星期日，2025年，2月5日“Fitting”按钮曲线拟合点击“Fitting”按钮，弹出“Fitting”窗口；点击“Newfit”按钮，可修改拟合项目名称“Fitname”，通过“Dataset”下拉菜单选择数据集，然后通过下拉菜单“Type