Matlab中曲线拟合的方法 投稿钱癔癕.doc

Matlab中曲线拟合的方法 投稿：钱癔癕 福建电脑 ２０１０年第７期 Ｍａｔｌａｂ中曲线拟合的方法 申红莲 （衡水学院教计学院河北衡水０５３０００） 【摘要】：本文从最小二乘法出发，介绍了在Ｍａｄａｂ中对数据进行曲线拟合的方法：命令行和图形用户界面方式，尤其是给出了大量数据的拟合方法，通过实例。给使用者以参考。 【关键词】：Ｍａｄａｂ；曲线拟舍；最小二乘法 １、引言 【ｐ６，ｓ６］＝ｐｏｌｙｔｉｔ（ｘＹ，６）ｐ６＝０．０４０３－０．０５４２－０．５３４７ Ｏ．６０４２ １．４９４４—２．５５００ Ｏ．００００ 曲线拟合（ｃｕｒｖｅｆｉｔｔｉｎｇ）是用连续曲线近似地刻画或比拟平面上离散点组所表示的坐标之问的函数关系的一种数据处理方法。在科学实验或社会活动中．通过实验或观测得到量茁与ｙ的一组数据对（ｚ。竹），ｉ＝ｌ，２，…，ｍ，其中各戈．是彼此不同的。人们希与，之问的依赖关系．即在一定意义下”最佳”地逼近或拟合已 ｙ６＝ｐｏｌｙｖａｌ（ｐ６，ｘ）； ｒｌｏｒｌｎｒ：２．８４４８ｅ一０１４ ｓ６＝Ｒ：【７ｘ７ｄｏｕｂｌｅ］ｄｆ：０ 望用一类与数据的规律相吻合的解析表达式，锨来反映量ｘ 知数据。似称作拟合函数，够的图像称作拟合曲线１１－３］。 进行曲线拟合可按如下步骤进行．第一步把所给的数据画在一个坐标图上，通过图表来判断其数学形式：第二步决定数学形式中的待定参数；第三步求得数学模型后，有时需要将实际测定的数据与用公式求出的理论值进行比较．判定其误差程度．这也是通常衡量拟合曲线优劣的标准．即选择合适的参数使得拟合函数与实际观测值的均方误差 ｐｌｏｔ（ｘ，ｙ，，ｍ，）ｐｌｏｔ（ｘ，ｙ２，，ｇ一，ｐｌｏｔ（ｘ，ｙ４，，ｍ．，）ｐｌｏｔ（ｘ，ｙ６，ｂｓ－＇）．ｘｌａｂｅｌ（＇ｘ，） ｙｌａｂｅｌ（＇ｙ３ ｌｅｇｅｎｄ（原始数据ｊ２次拟合，，４次多项式拟合，＇，６次多项式拟合１： 结果：以ｐ２，ｐ４，ｐ６为多项式的系数即可构造Ｙ关于ｘ的２次、４次和６次多项式函数。 Ｒ＝∑（钗ｔ）一Ｍ）２ ｆ＝ｌ 达到最小，此时称所求曲线为在最ｄｘ－乘法意义下的拟合曲线。 近年来利用数学软件Ｍａｆｌａｂ进行曲线拟合的情况较多．下边分别就几种情况进行说明。２、曲线拟合２．１命令行曲线拟合 【命令】ｐｏｌｙｆｉｔ 【调用格式】【ｐ，ｓ】＝ｐｏｌｙｔｉｔｘ，Ｙ，ｎ） 【功能】对于数据组ｘ，Ｙ进行多项式拟合，拟合的多项式的最高阶数为ｎ，其中Ｐ为多项式的系数矩阵．８为测差估计值的矩阵． 例ｌ给定下列数据，分别用２次，４次，６次多项式进行曲线拟合。 （一３，４），（－２，２），（－１，３），（Ｏ，Ｏ），（１，一１），（２，－２），（３，－５）解在ＭＡＴＬＡＢ命令窗口输入如下的程序代码：ｘ＝【一３—２－１ｙ＝【４ ｈｏｌｄ ２３ Ｏｎ 图１为拟合的多项式的曲线图像．从图中可以看出６次多项式曲线与原数据吻合情况较好。另外从ｓ６。ｓ４，ｓ２的第三项ｎｏｒｍｒ（均方误差）也可以看出，１．４４７６ｅ一０１４小于１．７０９４，１．７０９４又小于１．７５９３。由此可以看出在提高多项式的次数的情况下．有可能得到较好的拟合曲线，但是对于某些题目。也不是次数越高。精度越高，所以使用多项式拟合时关键是选择适中的次数。 ２．２图形用户界面的曲线拟合 为方便使用．在Ｍａｔｌａｂ中还提供了支持曲线拟合的图形用户界面嘲。它位于”Ｆｉｇｕｒｅ”窗１３的”Ｔｏｏｌｓ般使用步骤如下： （１）画图。 Ｂａｓｉｃ ０１２ ３】 ０－１－２－５】 ０． ㈣ ［ｐ２，ｓ２］＝ｐｏｌｙｆｉｔ（ｘ，ｙ，２）ｐ２＝－０．１３１０－１．３９２９ Ｆｉｔｔｉｎｇ”菜单中。一 ｙ２＝ｐｏｌｙｖＭ为； １１０１＂１１１１＂：．１．７５９３ ｓ２＝Ｒ：【３ｘ３ｄｏｕｂｌｅ］ｄｆ：４ ｐ４＝－０．００３８－０．０２７８－０． 在命令窗口中输入待拟合的数据，并利用Ｐｌｏｔ命令在“Ｆｉｇｕｒｅ”窗口中点击菜单栏中的”ＴｏｏｌｓｋＢａｓｉｃ 【ｐ４，叫＝ｐｏｌｙｆｉｔｘ，ｙ＇４） ０９４７－１．１９８４ （２） Ｆｉｔ－ ０．６２７７ ｒｉｎｇ＂得到”ＢａｓｉｃＦｉｔｔｉｎｇ”窗Ｉ：１。点击右下角的向右按钮，得”Ｂａｓｉｃ Ｙ４＝ｐｏｌｙｖａｌ（ｐ４，ｘ）； ｎｏｒｍＥＩ．７ＩＤ９４ ｓ４＝Ｒ：［５ｘ５ｄｏｕｂｌｅ］ｄｒ：２ Ｆｉｔｔｉｎｇ”窗口的全貌。 （３）在”Ｐｌｏｔｆｉｔｓ”复选框中选择”ｌｉｎｅａｒ＂”ｃｕｂｉｃ”、”５ｔｈ ｄｅｇｒｅｅ 鲞全项目：衡水学院２００９年度校扭青年专项课题（２００９０５６） 万方数据 ２０１０年第７期福建电脑 结果：Ｐ＝一３．６６；３６ ｄｒ：４８ ｎｏｒｍｒ：３４．０３８７ １０．８１０４ ７．０５３１ １１ ｐｏｌｙｎｏｍｉａｌ“等选项．即可进行线性．二次和多次多项式的拟合过程．并且可以观察