- 1、本文档共42页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值计算B大作业精要
课 程 设 计
课程名称: 数值计算B 设计题目: 数值计算B大作业 学 号: 姓 名:
完成时间:
题目一:多项式插值
某气象观测站在8:00(AM)开始每隔10分钟对天气作如下观测,用三次多项式插值函数(Newton)逼近如下曲线,插值节点数据如上表,并求出9点30分该地区的温度(x=10)。
x 1 2 3 4 5 6 7 8 y 22.5 23.3 24.4 21.70 25.2 28.5 24.8 25.4 二、数学原理
假设有n+1个不同的节点及函数在节点上的值(x,y),……(x,y),插值多项式有如下形式:
(1)
其中系数(i=0,1,2……n)为特定系数,可由插值样条(i=0,1,2……n)确定。
根据均差的定义,把x看成[a,b]上的一点,可得
f(x)= f()+f[]()
f[x, ]= f[]+f[x,] ()
……
f[x, ,…x]= f[x, ,…x]+ f[x, ,…x](x-x)
综合以上式子,把后一式代入前一式,可得到:
f(x)= f[]+f[]()+ f[]()()+
…+ f[x, ,…x]()…(x-x)+ f[x, ,…x,]= N(x)+
其中
N(x)= f[]+f[]()+ f[]()()+
…+ f[x, ,…x]()…(x-x) (2)
= f(x)- N(x)= f[x, ,…x,] (3)
=()…(x-x)
Newton插值的系数(i=0,1,2……n)可以用差商表示。一般有
[] (k=0,1,2,……,n ) (4)
把(4)代入(1)得到满足插值条件N(i=0,1,2,……n)的n次Newton插值多项式
N(x)=f()+f[]()+f[]()()+……+f[]()()…().
其中插值余项为:
介于之间。
三、程序设计
function [y,A,C,L]=newdscg(X,Y,x,M)
% y为对应x的值,A为差商表,C为多项式系数,L为多项式
% X为给定节点,Y为节点值,x为待求节点
n=length(X); m=length(x); % n为X的长度
for t=1:m
z=x(t); A=zeros(n,n);A(:,1)=Y;
s=0.0; p=1.0; q1=1.0; c1=1.0;
for j=2:n
for i=j:n
A(i,j)=(A(i,j-1)- A(i-1,j-1))/(X(i)-X(i-j+1));
end
q1=abs(q1*(z-X(j-1)));c1=c1*j;
end
C=A(n,n);q1=abs(q1*(z-X(n)));
for k=(n-1):-1:1
C=conv(C,poly(X(k)));
d=length(C);C(d)=C(d)+A(k,k);
end
y(k)= polyval(C, z); %输出y值
end
L(k,:)=poly2sym(C); %输出多项式
syms M,X=[1,3,5,7];Y=[22.5,24.4,25.2,24.8];x=10;
[y,A,C,L]=newdscg(X,Y,x,M)
y =
21.7313
A =
22.5000 0 0 0
24.4000 0.9500 0 0
25.2000 0.4000 -0.1375 0
24.8000 -0.2000 -0.1500 -0.0021
C =
-0.0021 -0.1187 1.4521 21.1688
L =
- x^3/480 -
您可能关注的文档
- 人教版二年级语文下册《三个儿子》课件PPT.ppt
- 教师职场礼仪培训.ppt
- 人教版二年级语文下册《玩具柜台前的孩子》课件PPT.ppt
- 人教版二年级语文上册《称赞》课件上课1.pptx
- 人教版二年级语文上册25古诗两首课件.ppt
- 公路工程交竣工验收-课件(PPT演示).ppt
- 人教版二年级语文下册各单元知识归纳.docx
- 教育心理学基本知识及重点(自己整理).doc
- 人教版二年级语文下册语文园地三第二课时.ppt
- 教育心理学思考题.docx
- 2024年部编版一年级数学(下册)期中试卷及答案(完整).pdf
- 2024年步步高高考语文一轮复习(部编新高考版)板块九 整本书阅读《红楼.pdf
- 2024年黑龙江省专升本考试大学艺术概论测试题含解析 .pdf
- 2024年有关幼儿园科研工作计划(三篇).pdf
- 2024年税务师之涉税服务实务过关检测试卷A卷附答案.pdf
- 2024浙江高考历史模拟试题含答案 .pdf
- 2024年监理工程师之交通工程监理案例分析通关题库(附答案) .pdf
- 2024年第一次广东省普通高中学业水平合格性考试化学仿真模拟卷01(考试.pdf
- 2024年高职单招环境科学专业技能考试题库及参考答案.pdf
- 2024年海上风电项目经营分析报告 .pdf
最近下载
- 《大观念下初中跨学科大单元课程开发的实践研究》课题研究方案.doc
- 妇女权益保障法讲座讲稿四篇.docx
- 企业主要负责人安全述职报告PPT.pptx
- 水浒传回目(全120回).pdf
- 逆变器eg-芯片EG8010串口通信使用说明.pdf
- GBZT213-2008血源性病原体职业接触防护导则-出版.pdf
- 学习贯彻党的创新理论情况,看学了多少、学得怎样,有什么收获和体会四个检视对照检查材料2篇文2024年.docx VIP
- 集中带量药品采购与使用的精细化管理系统及方法.pdf VIP
- 《TSG ZF001-2006 《安全阀安全技术监察规程》》.pdf
- Schneider Electric施耐德HVX12kV (U) 中压真空断路器中文操作手册安装和用户指南(中文).pdf
文档评论(0)