第七章函数逼近与曲线(面)拟合.docVIP

下载本文档

5
0
约1.95万字
约 12页
2017-08-17 发布于重庆
举报
版权申诉

第七章函数逼近与曲线(面)拟合.doc

1、本文档共12页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
8、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

第七章函数逼近与曲线(面)拟合.doc

7.1 曲线拟合、误差及其MATLAB程序例7.1.1 已知函数和一组数据列入表7–1中，比较最大误差，平均误差，均方根误差和误差平方和. 表7–1 例7.1.1的一组数据 xi -2.5 -1.7 -1.1 -0.8 0 0.1 0.5 3.6 yi -43.50 5.69 11.34 14.16 0 1.02 -6.37 185.84 解由给定的函数和数据，在MATLAB工作窗口输入 x=[-2.5,-1.7,-1.1,-0.8,0,0.1,0.5,3.6]; n=length(x); y=[-43.50 5.69 11.34 14.16 0 1.02 -6.37 185.84]; f=5.*x.^3-14.*x+7.*(sin(2*pi*x)).^2; fy=abs(f-y); fy2=fy.^2; [x,y,f,fy,fy2], Ew=max(fy), E1=sum(fy)/n, E2=sqrt((sum(fy2))/n), E=sum(fy2) 运行后屏幕显示如下 x y f fy fy2 -2.5000 -43.5000 -43.1250 0.3750 0.1406 -1.7000 5.6900 5.5666 0.1234 0.0152 -1.1000 11.3400 11.1634 0.1766 0.0312 -0.8000 14.1600 14.9716 0.8116 0.6586 0 0 0 0 0 0.1000 1.0200 1.0234 0.0034 0.0000 0.5000 -6.3700 -6.3750 0.0050 0.0000 3.6000 185.8400 185.2984 0.5416 0.2933 Ew = E1 = E2 = E = 0.8116 0.2546 0.3773 1.1390 7.2 曲线拟合的线性最小二乘法及其MATLAB程序例7.2.1 给出一组数据点列入表7–2中，试用线性最小二乘法求拟合曲线，并用（7.2），（7.3）和（7.4）式估计其误差，作出拟合曲线. 表7–2 例7.2.1的一组数据 xi -2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6 yi -192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04 解（1）在MATLAB工作窗口输入程序 x=[-2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6]; y=[-192.9 -85.50 -36.15 -26.52 -9.10 -8.43 -13.12 6.50 68.04]; plot(x,y,r*), legend(实验数据(xi,yi)) xlabel(x), ylabel(y), title(例7.2.1的数据点(xi,yi)的散点图) 在处的函数值，即输入程序 syms a1 a2 a3 a4 x=[-2.5 -1.7 -1.1 -0.8 0 0.1 1.5 2.7 3.6]; fi=a1.*x.^3+ a2.*x.^2+ a3.*x+ a4 运行后屏幕显示关于a1,a2, a3和a4的线性方程组 fi =[ -125/8*a1+25/4*a2-5/2*a3+a4, -4913/1000*a1+289/100*a2-17/10*a3+a4, -1331/1000*a1+121/100*a2-11/10*a3+a4, -64/125*a1+16/25*a2-4/5*a3+a4, a4, 1/1000*a1+1/100*a2+1/10*a3+a4, 27/8*a1+9/4*a2+3/2*a3+a4, 19683/1000*a1+729/100*a2+27/10*a3+a4, 5832/125*a1+324/25*a2+18/5*a3+a