- 1、本文档共3页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
MATLAB对实验数据进行曲线拟合画图
MATLAB软件提供了基本的曲线拟合函数的命令 Page 1 of 3
您正在学习:第五章 数据拟合 第四节 曲线拟合的MATLAB实现
§5.4 曲线拟合的MATLAB实现
MATLAB软件提供了基本的曲线拟合函数的命令.
多项式函数拟合:a=polyfit(xdata,ydata,n)
其中n表示多项式的最高阶数,xdata ,ydata为将要拟合的数据,它是用数组的方式输入.输出参数a为
y a x n a x a a [a , , a , a ]
拟合多项式 1 n n1 的系数 1 n n1
多项式在x处的值y可用下面程序计算.
y=polyval(a,x)
一般的曲线拟合:p=curvefit(‘Fun’,p0,xdata,ydata)
其中Fun表示函数Fun(p,data) 的M 函数文件,p0表示函数的初值.curvefit()命令的求解问题形式是
min sum (Fun(p ,xdata) ydata ).^2
p
若要求解点x处的函数值可用程序f=Fun(p,x)计算.
y aebx cedx (x , y ), i 1,2, , n,
例如已知函数形式 ,并且已知数据点 i i 要确定四个未知参数a,
b ,c,d .
xdata [x ,x , , x ];ydata [y , y , , y ]
使用curvefit命令,数据输入 1 2 n 1 2 n ;初值输
p 0 [a ,b ,c , d ] y aebx cedx
0 0 0 0 ;并且建立函数 的M文件(Fun .m ).若定义
p 1 a, p 2 b, p 3 c, p 4 d ,则输出p [p 1 , p 2 , p 3 , p 4 ]
又如引例的求解,MATLAB程序:
t=[l :16]; %数据输人
y=[ 4 6 .4 8 8 .4 9 .28 9 .5 9 .7 9 .86 10.2 10 .32 10 .42 10 .5 10 .55 10 .58 10 .6] ;
plot(t,y,’o’) %画散点图
p=polyfit(t,y,2) (二次多项式拟合)
计算结果:
p=-0.0445 1.0711 4.3252 %二次多项式的系数
由此得到某化合物的浓度y与时间t 的拟合函数
y 4.3252 1.0711t 0.0445t 2
对函数的精度如何检测呢?仍然以图形来检测,将散点与拟合曲线画在一个画面上.参见图5.3 .
xi linspace(0,16,160)
yi polyval (p , xi )
plot (x , y , o , xi , yi )
/eduwest/web_courseware/maths/0092/5/5-4.htm 2009-11-26
MATLAB软件提供了基本的曲线拟合函数的命令 Page 2 of 3
由此看见上述曲线拟合是比较吻合的。
文档评论(0)