![常州大学数值分析第三章_第1页](http://file4.renrendoc.com/view/2c02d97167ed286ac44aaf44a0ce2106/2c02d97167ed286ac44aaf44a0ce21061.gif)
![常州大学数值分析第三章_第2页](http://file4.renrendoc.com/view/2c02d97167ed286ac44aaf44a0ce2106/2c02d97167ed286ac44aaf44a0ce21062.gif)
![常州大学数值分析第三章_第3页](http://file4.renrendoc.com/view/2c02d97167ed286ac44aaf44a0ce2106/2c02d97167ed286ac44aaf44a0ce21063.gif)
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第三章作业1.设节点xo=0,x1=n/8,x2=n/4,x3=3n/8,x4=n/2,试适当选取上述节点,用拉格朗日插值法分别构造cosx在区间[0,n/2]上的一次、二次、四次差值多项式P1(x),P2(x)和P4(x),并分别计算P1(n/3),P2(n/3)和P4(n/3).解:xOx1x2x3x484x4P(兀/3)=0.500147.解:选取x=0,x=1,x=2,x=3为节点0123>>T0=[0.00.5];x=[123]';y=[1.252.753.5]';x0=2.8;T=aitken(x,y,x0,T0)0.00000.50000.00000.00000.00001.00001.25002.60000.00000.00002.00002.75000.00000.50000.00000.00000.00001.00001.25002.60000.00000.00002.00002.75003.65004.49000.00003.00003.50003.30003.23003.4820161)拉格朗日差值.选取函数点八、、P(x)^―0.636620x+11P(兀/3)«0.3333331(2)选择x0=0,x2=n/4,x4=n/2的节点y0=cosx0=1,y2=cosx2=0.707106,y4=cosx4=0,可得、(x-x)(x-x)(x-x)(x-x)P(x)=y12+y0210(x-x)(x-x)0(x-x)1(x-x)(x-0(x-x)1(x-x)+(x-x)(x-x)2(x-x)(x-x)即P(x)«-1.482067x2+5.511124x+1'%(兀/3)«1.145968(3)选择x0=0,,x1=n/8,x2=n/4,x3=3n/8,x4=n/2的节点y0=cosx0=1,y1=cosx1=0.923879,y2=cosx2=0.707106,y3=cosx3=0.382683,y4=cosx4=0可得P4(x)=£y(H)i=0j=0,声%Xj,得P3(x)=1+0.0031x-0.5154x0110+0.0242x3+0.02y=cos(x)+sin(x),xg[一兀,兀]x0二-pi:0.5*pi:pi;y0=cos(x0);x=-pi:0.05*pi:pi;iflength(x0)~=length(y0)error('Thelengthofx0mustbeequaltoitofy0');endw=length(x0);n=w-1;L=zeros(w,w);fork=1:n+10110V=1;forj=1:n+1ifk~二jifabs(x0(k)-x0(j))<epserror('DividedbyZero,therearetwonodesarethesame');endV=conv(V,poly(x0(j)))/(x0(k)-x0(j));endendL(k,:)=V;
endC=y0*L;Y=polyval(C,x);y=cos(x)+sin(x);r=y-Y;plot(x,Y,'r—,,x,y,,b-,,x,r,,k-.','LineWidth',2);legend(,Lagrangepolynomial,,,Theoriginalf(x),,,Error,,0)2)牛顿差值12)牛顿差值.选取函数y=1+担,]日一5,5]>>x0=-5:1:5;〉〉y0=1./(1+x0.*x0);>>x=-5:0.1:5;>>iflength(x0)~=length(y0)>>error('Thelengthofx0mustbeequaltoitofy0');end>>n=length(x0);>>D=zeros(n,n);>>D(:,1)=y0';>>forj=2:nfork=j:nifabs(x0(k)-x0(k-j+1))<epserror('DividedbyZero,therearetwonodesarethesame');end>>D(k,j)=(D(k,j-1)-D(k-1,j-1))/(x0(k)-x0(k-j+1));endendC=D(n,n);fork=(n-1):-1:1C=conv(C,poly(x0(k)));m=length(C);C(m)=C(m)+D(k,k);endY=polyval(C,x);y=1./(1+x.*x+20*x);r=y-Y;>plot(x,Y,'r--',x,y,'b-',x,r,'k-.','LineWidth',2);>legend('Newtonpolynomial','Th
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 消化内科护理工作上半年总结
- 撞锤机器项目成效分析报告
- GPS导航设备细分市场深度研究报告
- 丝线绕线机产品市场需求分析报告
- 蜡笔项目成效分析报告
- 烧烤用手套产品商业计划书
- 橡胶成型机项目成效分析报告
- 斯诺克台球架杆产品相关项目实施方案
- 乘客登乘用可移动非金属梯产品相关项目建议书
- 库房整改方案8篇
- 小学一下语文无纸化测试题例
- 2023年的18项医疗核心制度
- 2022年嵊州市国企招聘考试真题及答案
- 工业区位因素及其变化教学设计
- 王志魁《化工原理》课后思考题参考答案
- 第5章不完全竞争市场-1-垄断-学生课件
- 少儿美术网红教案课件- 3岁-6岁 《夏日荷塘》
- 办公场所风险分级管控清单
- 2023年浙江绍兴市越城区国有企业专场招聘笔试参考题库附带答案详解
- 2023年云南省执业药师继续教育全科答案
- 2023ESC急性肺栓塞诊断和管理指南中文完整版
评论
0/150
提交评论