[理学]matlab-曲线拟合工具箱讲义.ppt

曲线拟合工具箱 设有实验数据 ,寻找函数 使得函数在点 处的函数值与观测数据偏差的平方和达到最小.即求满足如下条件的函数 使得 曲线拟合定义 解决此类问题有以下几个步骤： （1）首先作出散点图，确定函数的类别； （2）根据已知数据确定待定参数的初始值，利用 Matlab软件计算最佳参数； （3）根据可决系数，比较拟合效果，计算可决系数的公式为 曲线拟合可分为： （1）参数拟合 ---- 最小二乘法 （2）非参数拟合 ---- 插值法 一、数据预处理 在曲线拟合之前必须对数据进行与处理，去 除界外值、不定值和重复值，以减少认为误 差，提高拟合的精度。 数据预处理包括： （1）数据输入与查看 （2）数据的预处理 传输数据通过数据GUI来实现，查看数据点 通过曲线拟合工具的散点图来实现。 1.输入和查看数据集 （1）打开曲线拟合工具界面 通过cftool命令打开曲线拟合工具界面 5个命令按钮 Data按钮：可输出、查看和平滑数据； Fitting按钮：可拟合数据、比较拟合曲线和数据集； Exclude按钮：可以从拟合曲线中排除特殊的数据点； Ploting按钮：在选定区间后，单击按钮，可以显示拟合曲线和数据集； Analysis按钮：可以做内插法、外推法、微分或积分拟合。 （2）输入数据集 在输入数据之前，数据变量必须存在于 matlab的工作区间。可以通过load命令输 入变量。单击曲线拟合工具界面中的Data 按钮，打开Data对话框，在对话框中进行设 置，可以输入数据。 Data对话框 包括两个选项卡：Data Sets 和 Smooth. Data Sets选项卡： .Import workspace vectors 把向量输 入工作区，主要以变量必须具有相同的维数 ，无穷大的值和不定值被忽略。 X data 用于选择观测数据 Y data 用于选择X的响应数据 Weight 用于选择权重，与响应数据相联系的向量，如果没选择，默认值为1. .Preview 对所选向量进行图形化预览 .Data set name 设置数据集的名称。工具箱可以随即产生唯一的文件名，但用户可以重命名。 .Data sets 选项以列表的形式显示所有拟合的数据集。当选择一个数据集时，可以对它做如下操作： .View 查看数据集，以图标形式和列表形式，可以选择方法排除异常值； .Rename 重命名 .Delete 删去数据组 例：输入数据，采用matlab自带的文件census census 有两个变量：cdate和pop。 cdate是一个年向量，包括1790-1990 年，间隔为10年； pop是对应年份的美国人口。 whos -file census Name Size Bytes Class Attributes cdate 21x1 168 double pop 21x1 168 double load census cftool(cdate,pop) 散点图 单击Data按钮 在X data和Y data两个下拉式列表框中选 择变量名，将在Data对话框中显示散点图的 预览效果： 当选择Data sets列表框中的数据集时，单 击View按钮，打开View Data Set对话框 工作表方式 2.数据的预处理 在曲线拟合工具箱中，数据的预处理主要包 括平滑法、排除法和区间排除法等。 (1)平滑数据 打开拟合工具箱，单击Data按钮，打开 Data对话框，选择Smooth选项卡 Smooth选项卡各选项的功能： .Original data set 用于挑选需要拟合的数据集； .Smoothed data set平滑数据的名称； .Method用于选择平滑数据的方法，每一个相应数据用通过特殊的曲线平滑方法所计算的结果来取代。平滑数据的方法包括： （ⅰ）Moving average 用移动平均值进行替换； （ⅱ）Lowess局部加权散点图平滑数据，采用线性最小二乘法和一阶多项式拟合得到的数据进行替换； （ⅲ）Loess局部加权散点图平滑数据，采用线性最小二乘法和二阶多项式拟合得到的数据进行交换； （ⅳ）Savitzky-Golay 采用未加权的线性最小二乘法过滤数据，利用指定阶数的多项式得到的数据进行替换； （ⅴ）Span用于进行平滑计算的数据点的数