MATLAB在数学建模中的应用 随手笔记(二)----- 数据拟合方法
' I* N7 O u! c一、多项式拟合' x! Z( d/ l( {% C3 y/ u3 s
将数据点按多项式的形式进行拟合,使用最小二乘法,可以确定多项式的系数。多项式拟合有指令语句和图形窗口两种方法:
0 `8 {' d" ~, P% @8 x1、多项式拟合指令7 \% l9 a: N& v" |
polyfit(x,y,n) :多项式拟合,返回降幂排列的多项式系数。
- J0 z4 u. f. f( h) Mpolyval(p,xi) :计算多项式的值。: v p! f5 E, C8 ^% a5 F
其中,x,y是数据点的值;n是拟合的最高次幂;p是返回的多项式系数;xi是要求的点的横坐标。, \ @, @, E9 A' D
2、图像窗口的多项式拟合" X7 i$ r5 A4 ^! |8 w: }8 v
在图形窗口中可以用菜单的方式对数据进行简单、快速、高效的拟合。
( a1 s8 ~3 j% j; \具体步骤(2017b):
1 R# V8 K- p1 Z) ]% t . G/ J6 L9 L V3 \3 V6 p& ?
& S5 [0 O( g9 ]/ z* {; _
2 j8 q' d( v8 I
![]()
$ p( l/ l: r9 ?' Y& O, F ) f: u0 a& N0 e% J1 B' s3 `
二、指定函数拟合
8 P: g1 [4 ]' F9 O K* h) ^1 b在MATLAB中也可以用用户自定义的函数进行拟合,通过下面的例子读者可以了解指定函数进行数据拟合的基本方法。4 p' b9 f3 m7 O) L
对该数据进行指定拟合:4 m+ o$ C; m/ r! B/ t1 x8 `' R9 \
x=[0;0.4;1.2;2;2.8;3.6;4.4;5.2;6;7.2;8;9.2;10.4;11.6;12.4;13.6;14.4;15];%列向量% ~$ l. O. |' }5 o+ W, t& f
y=[1;0.85;0.29;-0.27;-0.53;-0.4;-0.12;0.17;0.28;0.15;-0.03;-0.15;-0.071;0.059;0.08;0.032;-0.015;-0.02];
h' t5 |8 Q. ?在MATLAB中可以用如下命令画出上述点的散点图( N7 J7 r+ c& {6 t; t
/ ^, d6 r; i6 X2 l% G% S
, I% n. S6 O4 U2 f- J1 I3 f2 _& p
知道其对应的函数形式为f(t)=acos(kt)e^wt,则可用MATLAB进行拟合。编写如下M文件:/ x# S5 A; L0 ]% x2 t. F
& w: Q0 x4 I* K; \
运行此程序,结果如下:& a3 j7 r1 F5 S9 D
8 c; M% |1 J# t
, l1 S: w: Q+ E* X
程序中,fittype函数是自定义拟合函数;cfun=fit(x,y,f)是根据自定义的拟合函数f来拟合数据x,y。注意:此处数据必须为列向量的形式。fittype函数和fit函数的用法和参数规则可参考Help( {# i) p( ?7 X
从结果可以看出,拟合的曲线为f(x)=0.9987cos(1.001t)e^-0.2066t。拟合曲线给出了数据的大致趋势,效果很好,并给出了各参数的置信区间。( o. @0 q$ u. ~+ v
注意:command window里给出了warning,是由a,k,w三个参数的初始值未给出导致的,因此如果拟合结果不理想,可以多运行几次。7 N8 [7 l9 y& \" T5 A8 R6 f
三、曲线拟合工具箱. |" B7 B+ s5 a* c, _: ~$ V" S
MATLAB的曲线拟合工具箱功能非常的强大,使用也很方便。 详细步骤:![]()
6 y' n+ P6 n8 L' g6 \8 V 3 n7 o2 c4 d4 v% J2 y& j4 k
( s& ?$ y% g( V2 v# P! p6 @
界面中有五个按钮,功能是:
: l/ s6 T; u/ Z( _% Y1 UData:输出、查看和平滑数据;
2 q7 c6 M8 P8 U3 r! SFitting:拟合数据、比较拟合曲线和数据集; ^5 ^* q" y& Z- p F
Exclude:可从拟合曲线中排除特殊的数据点;: o$ e6 R! y) H0 T
Plotting:选定区间后,单击按钮可选择原始数据和拟合数据作图;
7 _8 W, n; v& d$ M" AAnalysis:对拟合进行满意度、偏差等分析。
5 C+ s Z, e1 m曲线拟合工具箱中包含了各种常用的数据拟合方法,可以对各种函数进行拟合,具体的操作方法请参考Help。6 h0 }! D2 l5 y; A# x: u
3 v2 M8 k# V7 t. X* o$ ]# @3 z
. U9 T N) u, |* O0 U5 U) T) c |