数据拟合的最小二乘.doc

  1. 1、本文档共5页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数据拟合的最小二乘

实验七 数据拟合实验 信息与计算科学金融 崔振威 201002034031 实验目的: 1、掌握数据拟合的操作步骤及其原理 实验内容: 1、p213-5.2.8-2(a)、p213-3、p215-5.2.9-1 实验要求: 根据给定数据画出散布图并画出最终的拟合曲线图 用列表的形式给出中间所使用向量的各分量 写出相应的拟合曲线函数 P213-5.2.8-2 2、对下列数据集,求解最小二乘抛物线f(x)=Ax^2+Bx+C (a) xk yk -2 -1 0 1 2 -5.8 1.1 3.8 3.3 -1.5 在matlab窗口中输入: x=[-2 -1 0 1 2];y=[-5.8;1.1;3.8;3.3;-1.5];A=[ones(size(x)) x x.^2];c=A\y;c 结果输出: ans = 3.9800 1.0800 -1.9000 z=[-2:0.01:2];t=[ones(size(z)) z z.^2];f=t*c;plot(z,f,x,y,r) 结果输出,如下图所示: p213-3 (a) 解: 在matlab窗口中输入: a=[1 2 3 4 5];b=[1 1 1 1 1];y=[log(0.6) log(1.9) log(4.3) log(7.6) log(12.6)] 可以得出下面结果: y = -0.5108 0.6419 1.4586 2.0281 2.5337 在matlab窗口中输入: A=[5,15;15,55];B=[6.1515;25.9298];T=inv(A)*B T = -1.0123 0.7475 由上面输出结果以及根据公式Ax=B求出x,即得出: 画图: 在matlab窗口中输入: t=[1:0.01:5];f=exp(-1.0123+0.7475*t); plot(t,f,g);hold on x=[1 2 3 4 5]; z=[0.6 1.9 4.3 7.6 12.6];plot(x,z,o) 可以得出下图: 由公式:即可得出4.7913 (b) 解:在matlab窗口中输入: a=[log(1) log(2) log(3) log(4) log(5)];b=[1 1 1 1 1];sum(b.*b) 可得: ans = 5 在matlab窗口中输入: sum(b.*a) 得出: ans = 4.7875 在matlab窗口中输入并得出 sum(a.*a) ans = 6.1995 在matlab窗口中输入并得出 y=[log(0.6) log(1.9) log(4.3) log(7.6) log(12.6)];sum(b.*y) ans = 6.1515 sum(a.*y) ans = 8.9368 再在窗口中输入可得出结果: A=[5,4.7875;4.7875,6.1995];b=[6.1515;8.9368];y=inv(A)*b y = -0.5755 1.8860 由上述结果以及根据可以求出: 在matlab中画出图: 由公式 即可求出0.402 分析:由上面求出的结果(a)、(b)中可以看出,(b)的平方根误差小于(a),平方根误差越小,拟合得效果越好,所以最佳拟合曲线应选择:

文档评论(0)

zilaiye + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档