已阅读5页,还剩1页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
西京学院数学软件实验任务书课程名称数学软件实验班级数0901学号0912020107姓名李亚强实验课题等距节点插值,Hermite插值,分段插值(线性,二次,三次)实验目的熟悉等距节点插值,Hermite插值,分段插值(线性,二次,三次)实验要求运用Matlab/C/C+/Java/Maple/Mathematica等其中一种语言完成实验内容等距节点插值,Hermite插值,分段插值(线性,二次,三次)成绩教师实验十六实验报告1、 实验名称:等距节点插值,Hermite插值,分段插值(线性,二次,三次)。2、 实验目的:进一步熟悉等距节点插值,Hermite插值,分段插值(线性,二次,三次)。3、 实验要求:运用Matlab/C/C+/Java/Maple/Mathematica等其中一种语言完成程序设计。4、 实验原理:1 等距节点插值:差分分为前向差分、后向差分和中心差分三种,它们的记法及定义如下所示:阶前向差分公式阶后向差分公式阶中心差分公式其中: -前向差分; -后向差分; -中心差分。假设,为了方便计算,构造差分表()。这里只说明前向牛顿插值,其多项式可表示为如下形式: 其中为步长,且的取值范围为。2 埃尔米特插值:埃尔米特插值法满足在节点上等于给定函数值,而且在节点上的导数值也等于给定的导数值,对于有高阶导数的情况,埃尔米特插值多项式比较复杂,在实际应用中,常常遇到的是函数值与一阶导数值给定的情况,在这种情况下,个节点的埃尔米特插值多项式的表达形式如下所示:其中 3 分段插值:给定插值节点、节点函数值及对应的导数值,则满足下面条件的分段埃尔米特插值函数的表达式如下所示:5、 实验内容:%等距节点插值function f,f0= dengjujiedian(x,y,x0)syms t;if(length(x) = length(y) n = length(x); c(1:n) = 0.0;else disp(x和y的维数不相等!); return;endf = y(1);y1 = 0;xx =linspace(x(1),x(n),(x(2)-x(1);if(xx = x) disp(节点之间不是等距的!); return;endfor(i=1:n-1) for(j=1:n-i) y1(j) = y(j+1)-y(j); end c(i) = y1(1); l = t; for(k=1:i-1) l = l*(t-k); end; f = f + c(i)*l/factorial(i); simplify(f); y = y1; endf0=subs(f,t,(x0-x(1)/(x(2)-x(1);%埃尔米特插值function f,f0= Hermite(x,y,y_1,x0)syms t;f = 0.0;if(length(x) = length(y) if(length(y) = length(y_1) n = length(x); else disp(y和y的导数的维数不相等!); return; endelse disp(x和y的维数不相等!); return;endfor i=1:n h = 1.0; a = 0.0; for j=1:n if( j = i) h = h*(t-x(j)2/(x(i)-x(j)2); a = a + 1/(x(i)-x(j); end end f = f + h*(x(i)-t)*(2*a*y(i)-y_1(i)+y(i);endf0=subs(f,t,x0);%分段差值function f,f0 = fenduan(x,y,y_1,x0)syms t;f = 0.0;f0 = 0.0;if(length(x) = length(y) if(length(y) = length(y_1) n = length(x); else disp(y和y的导数的维数不相等!); return; endelse disp(x和y的维数不相等!); return;end for i=1:n if(x(i)=x0) index = i; break; endend h = x(index+1) - x(index); fl = y(index)*(1+2*(t-x(index)/h)*(t-x(index+1)2/h/h + . y(index+1)*(1-2*(t-x(index+1)/h)*(t-x(index)2/h/h;fr = y_1(index)*(t-x(index)*(t-x(
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 24.3 正多边形和圆 人教版数学九年级上册堂堂练(含答案)
- 2024-2025学年陕西省渭南市合阳县多校联考四年级(上)第一次月考数学试卷
- 江苏省南京市六区重点名校2024届中考数学仿真试卷含解析
- 植物病虫害诊断与防治行业市场发展趋势及投资咨询报告
- 线上职业发展咨询平台行业调研及投资前景分析报告
- 2023年浙江数智交院科技股份有限公司招聘笔试真题
- 2023年内蒙古自治区科学技术厅所属事业单位招聘考试试题及答案
- 2023年陇南市宕昌县招聘司法协理员笔试真题
- 2023年福建省晋江水务集团有限公司招聘笔试真题
- 2024年江苏下载客运资格证模拟考试题
- 厦门大学实习报告模板
- 节流孔板孔径及壁厚计算
- 国家电网公司成本管理办法
- 硬件设计checklist
- 光伏发电工程项目施工组织机构与人员配置
- 同分异构体(课堂PPT)
- 膜分离氢器操作指南
- 假期去向表(模板)
- 在银政与银企互惠共赢合作发展座谈会上的发言
- 工业厂房及管理楼水电安装施工方案
- GB 14754-2010 食品安全国家标准 食品添加剂 维生素C(抗坏血酸)
评论
0/150
提交评论