- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab插值与曲线拟合实验报告 湖南大学
湖南大学电气与信息工程学院
《数值计算》课程 上机实验报告
姓名: 班级: 学号: 日期: 指导老师: 本次实验题号:第 次实验
实验目的:
用MATLAB实现拉格朗日插值和分段线性插值。
0.46 0.47 0.48 0.49 Y 0.4846555 0.4937542 0.5027498 0.5116683
曲线拟合方法的实际应用
用电压V=10V的电池给电容器充电,电容器上t时刻的电压v(t)=V-(V-V0)e^(-t/T),其中V0是电容器的初始电压,T是充电常数。实验测量了一组数据如下,请根据数据表确定V0和T的大小。
t 0.5 1 2 3 4 5 7 9 V(t) 6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63
三. 算法介绍或方法基础
1.1 拉格朗日插值法
1.1.1 算法理论
对某个多项式函数,已知有给定的k+1个取值点:
其中对应着自变量的位置,而对应着函数在这个位置的取值。假设任意两个不同的xj都互不相同,那么应用拉格朗日插值公式所得到的拉格朗日插值多项式为:
其中每个为拉格朗日基本多项式(或称插值基函数),其表达式为:
拉格朗日基本多项式 的特点是在 上取值为1,在其它的点 上取值为0。对于给定的 个点,拉格朗日插值法的思路是找到一个在一点取值为1,而在其他点取值都是0的多项式。这样,多项式在点取值为,而在其他点取值都是0。而多项式
就可以满足
在其它点取值为0的多项式容易找到,例如:
它在点取值为:。由于已经假定两两互不相同,因此上面的取值不等于0。于是,将多项式除以这个取值,就得到一个满足“在取值为1,而在其他点取值都是0的多项式”:
这就是拉格朗日基本多项式。
拉格朗日插值法的公式结构整齐紧凑,在理论分析中十分方便,然而在计算中,当插值点增加或减少一个时,所对应的基本多项式就需要全部重新计算,于是整个公式都会变化,非常繁琐。此外,当插值点比较多的时候,拉格朗日插值多项式的次数可能会很高,因此具有数值不稳定的特点,也就是说尽管在已知的几个点取到给定的数值,但在附近却会和“实际上”的值之间有很大的偏差。这类现象也被称为龙格现象,解决的办法是分段用较低次数的插值多项式。
对于已给定的点 和待估计的点的横坐标x,如上述理论,将其值代入
计算出插值基函数的值,然后根据公式:
计算出纵坐标的估计值,由此完成对该点的插值过程,其中k为该点插值的阶数。
首先介绍线型插值,假设我们已知坐标 与 要得到 区间内某一位置?x?在直线上的值。根据图中所示,我们得到
由于?x?值已知,所以可以从公式得到 y 的值
已知求的过程与以上过程相同,只是与要进行交换。
线性插值经常用于已知函数?在两点的值要近似获得其它点数值的方法,这种近似方法的误差定义为:
其中?p?表示上面定义的线性插值多项式
根据罗尔定理,我们可以证明:如果?f?有二阶连续导数,那么误差范围是
正如所看到的,函数上两点之间的近似随着所近似的函数的二阶导数的增大而逐渐变差。从直观上来看也是这样:函数的曲率越大,简单线性插值近似的误差也越大。
分为k-1段,每段利用其端点进行线型插值。
1.2.2 算法描述
利用已给定的点 对插值区间分为段,将每段的端点与
在所构成的区间进行线性插值。
1.3曲线拟合
1.3.1 算法理论
令待求的未知量为,它们可由个直接测量通过下列函数关系求得:
若为真值,由上述已知函数求出真值,若其测量值为,则对应的误差为 最小二乘法可定量表示为:
将拟合函数取线性函数是一种简单的数据拟合方法,利用数据点
确定线性拟合函数
称为对数据的线性拟合。对于线性拟合问题,需要求函数
的最小值点,该问题的几何背景是寻求一条直线,使该直线与数据表所确定的平面散点的纵向距离的平方和最小。
由函数对两个变量求导得:
其余等于零,得正规方程组
也可将其矩阵形式写出来即:
解得的值,将其代入 即可得到拟合线性函数。
为了确定数据拟合问题,选用幂函数作为函数类,则
这就是多项式拟合函数.
为了确定拟合函数的系数,需要求解正规方程组
也可以用矩阵形式表示为
解得最优解即可,将其代入即可得到拟合多项式。若用指数函数对数据进行拟合,则使用线性化方
您可能关注的文档
- 2013篮球协协暑期社会实践(全民健身).doc
- 2013考研数学高效复习三大指南[厦门大学分数线].doc
- 沙冰奶昔制作方法2012年2月10日文档.doc
- 织物整理剂项目商业计划书(2013年融资成功案例范文)-专家免费咨询.pdf
- 特岗习题4.docx.doc
- 电气石超细粉体对废水中Zn^2+离子吸附的研究.pdf
- 基于有限元的磁流变弹性体设备磁致特性分析.pdf
- 绥化中医院妇产科介绍.doc
- (最新)国有企业审计.pdf
- 2013年晚会主持词修改 (最后定稿).doc
- 2024年学校党总支巡察整改专题民主生活会个人对照检查材料3.docx
- 2025年民主生活会个人对照检查发言材料(四个带头).docx
- 县委常委班子2025年专题生活会带头严守政治纪律和政治规矩,维护党的团结统一等“四个带头方面”对照检查材料四个带头:.docx
- 巡察整改专题民主生活会个人对照检查材料5.docx
- 2024年度围绕带头增强党性、严守纪律、砥砺作风方面等“四个方面”自我对照(问题、措施)7.docx
- 2025年度民主生活会领导班子对照检查材料(“四个带头”).docx
- 国企党委书记2025年度民主生活会个人对照检查材料(五个带头).docx
- 带头严守政治纪律和政治规矩,维护党的团结统一等(四个方面)存在的问题整改发言提纲.docx
- 党委书记党组书记2025年带头增强党性、严守纪律、砥砺作风方面等“四个带头”个人对照检查发言材料.docx
- 2025年巡视巡察专题民主生活会对照检查材料.docx
最近下载
- 《ISO31000:2024风险管理指南》指导手册(雷泽佳译2024-04).pptx VIP
- 2025届高考语文复习:《到橘子林去》高考真题说题课件.pptx VIP
- 广州市番禺区2015-2016学年九年级上学期政治期末统考试卷(含答案)详解.doc
- 《瑜伽》课程说课备课讲稿.ppt
- 2024年度医院整形美容外科学述职报告课件.pptx
- ISO31022-2020风险管理-法律风险管理指南(雷泽佳译).pdf VIP
- GBT 16180-2014《劳动能力鉴定 职工工伤与职业病致残等级》.doc
- TASCAM DR-40中文用户手册(说明书).pdf
- 体检中心一月质量控制重点工作计划和质控小结.docx VIP
- 2025最新字帖3年级下册_楷体.pdf
文档评论(0)