- 1、本文档共6页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数学实验“等距节点插值,Hermite插值,分段插值(线性,二次,三次)”实验报告(内含matlab程序)
西京学院数学软件实验任务书
课程名称 数学软件实验 班级 数0901 学号 0912020107 姓名 李亚强 实验课题 等距节点插值,Hermite插值,分段插值(线性,二次,三次) 实验目的 熟悉等距节点插值,Hermite插值,分段插值(线性,二次,三次) 实验要求 运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成 实验内容 等距节点插值,Hermite插值,分段插值(线性,二次,三次) 成绩 教师 实验十六实验报告
实验名称:等距节点插值,Hermite插值,分段插值(线性,二次,三次)。
实验目的:进一步熟悉等距节点插值,Hermite插值,分段插值(线性,二次,三次)。
实验要求:运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成程序设计。
实验原理:
等距节点插值:
差分分为前向差分、后向差分和中心差分三种,它们的记法及定义如下所示:
阶前向差分公式
阶后向差分公式
阶中心差分公式
其中: -前向差分; -后向差分; -中心差分。
假设,为了方便计算,构造差分表()。
这里只说明前向牛顿插值,其多项式可表示为如下形式:
其中为步长,,且的取值范围为。
埃尔米特插值:
埃尔米特插值法满足在节点上等于给定函数值,而且在节点上的导数值也等于给定的导数值,对于有高阶导数的情况,埃尔米特插值多项式比较复杂,在实际应用中,常常遇到的是函数值与一阶导数值给定的情况,在这种情况下,个节点的埃尔米特插值多项式的表达形式如下所示:
其中
分段插值:
给定插值节点、节点函数值及对应的导数值,则满足下面条件
的分段埃尔米特插值函数的表达式如下所示:
实验内容:
%等距节点插值
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;
end
f = y(1);
y1 = 0;
xx =linspace(x(1),x(n),(x(2)-x(1)));
if(xx ~= x)
disp(节点之间不是等距的!);
return;
end
for(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;
end
f0=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;
end
else
disp(x和y的维数不相等!);
return;
end
for 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));
end
f0=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的导数的维数不相等!);
retu
您可能关注的文档
- 数学八年级下资源与评价答案.doc
- 数学分析.docx
- 数学分析3.doc
- 数学分析习题.docx
- 数学分析教学大纲.doc
- 数学分析研究性学习.ppt
- 数学分析绪论.ppt
- 数学分析论文.docx
- 数学分析课程教学大纲.doc
- 数学分析论文.doc
- 山东省烟台市牟平区(五四制)2023-2024学年七年级下学期期中考试道德与法治试卷(含答案).pdf
- 山东省临沂市临沭县2024年四年级英语第二学期期中考试模拟试题含答案.pdf
- 实验报告家兔动脉血压的神经体液调节影响尿生成的因素.pdf
- 国家开放大学一网一平台电大《教育学》形考任务1-4网考题库及答案.pdf
- 小学数学四年级下册易错题期末综合复习-人教版(含答案).pdf
- 四川省成都市成都市第七中学2024届高三第二次联考英语试题理试题含解析.pdf
- 小学四年级数学下册期末复习试卷及答案.pdf
- 声乐社团活动策划书.pdf
- 大学生创新创业训练计划项目.pdf
- 安徽芜湖市初中英语七年级上册Unit 5知识点复习(含答案解析).pdf
文档评论(0)