- 1、本文档共17页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
信计101 郭英贵 学号:201000901031
信计、数应专业10级数值方法计算实习题
要求:1、用Matlab语言或你熟悉的其他语言编写程序,使之尽可能具有通用性;2、根据上机计算实践,对所使用的数值方法的特点、性质、有效性、误差和收敛性等方面进行必要的讨论和分析;3、完成计算后写出实验报告,内容包括:课题名称、解决的问题、采用的数值方法、算法程序、数值结果、对实验结果的讨论和分析等;4、特别说明:严禁抄袭,否则一经发现,所有雷同实验报告最多评为及格。
在区间[-1,1]上分别取用两组等距节点对龙格函数作多项式插值及三次样条插值,对每个值,分别画出插值函数即的图形。
解:(1)多项式插值
①先建立一个多项式插值的M-file;
输入如下的命令(如牛顿插值公式):
function [C,D]=newpoly(X,Y)
n=length(X);
D=zeros(n,n)
D(:,1)=Y
for j=2:n
for k=j:n
D(k,j)=(D(k,j-1)- D(k-1,j-1))/(X(k)-X(k-j+1));
end
end
C=D(n,n);
for k=(n-1):-1:1
C=conv(C,poly(X(k)))
m=length(C);
C(m)= C(m)+D(k,k);
end
②当n=10时,我们在命令窗口中输入以下的命令:
clear,clf,hold on;
X=-1:0.2:1;
Y=1./(1+25*X.^2);
[C,D]=newpoly(X,Y);
x=-1:0.01:1;
y=polyval(C,x);
plot(x,y,X,Y,.);
grid on;
xp=-1:0.2:1;
z=1./(1+25*xp.^2);
plot(xp,z,r)
得到插值函数和f(x)图形:
③当n=20时,我们在命令窗口中输入以下的命令:
clear,clf,hold on;
X=-1:0.1:1;
Y=1./(1+25*X.^2);
[C,D]=newpoly(X,Y);
x=-1:0.01:1;
y=polyval(C,x);
plot(x,y,X,Y,.);
grid on;
xp=-1:0.1:1;
z=1./(1+25*xp.^2);
plot(xp,z,r)
得到插值函数和f(x)图形:
(2)三次样条插值
①先建立一个多项式插值的M-file;
输入如下的命令:
function S=csfit(X,Y,dx0,dxn)
N=length(X)-1;
H=diff(X);
D=diff(Y)./H;
A=H(2:N-1);
B=2*(H(1:N-1)+H(2:N));
C=H(2:N);
U=6*diff(D);
B(1)=B(1)-H(1)/2;
U(1)=U(1)-3*(D(1));
B(N-1)=B(N-1)-H(N)/2;
U(N-1)=U(N-1)-3*(-D(N));
for k=2:N-1
temp=A(k-1)/B(k-1);
B(k)=B(k)-temp*C(k-1);
U(k)=U(k)-temp*U(k-1);
end
M(N)=U(N-1)/B(N-1);
for k=N-2:-1:1
M(k+1)=(U(k)-C(k)*M(k+2))/B(k);
end
M(1)=3*(D(1)-dx0)/H(1)-M(2)/2;
M(N+1)=3*(dxn-D(N))/H(N)-M(N)/2;
for k=0:N-1
S(k+1,1)=(M(k+2)-M(k+1))/(6*H(k+1));
S(k+1,2)=M(k+1)/2;
S(k+1,3)=D(k+1)-H(k+1)*(2*M(k+1)+M(k+2))/6;
S(k+1,4)=Y(k+1);
end
②当n=10时,我们在命令窗口中输入以下的命令:
clear,clc
X=-1:0.2:1;
Y=1./(25*X.^2+1);
dx0= 0.0739644970414201;dxn= -0.0739644970414201;
S=csfit(X,Y,dx0,dxn)
x1=-1:0.01:-0.5;y1=polyval(S(1,:),x1-X(1));
x2=-0.5:0.01:0;y2=polyval(S(2,:),x2-X(2));
x3=0:0.01:0.5; y3=polyval(S(3,:
您可能关注的文档
最近下载
- 备战2023年高考语文一轮复习考点微专题(新高考地区专用)考向28 诗歌鉴赏之语言(含详解).docx VIP
- 飞机交易平台及飞机拆解项目可行性研究报告.doc
- 视听语言PPT全套教学课件.pptx
- 健康评估-河南大学-中国大学MOOC慕课答案.pdf
- 初中音乐人音版《七年级上册青年友谊圆舞曲》课件_1.ppt
- 基于Java的小区物业管理系统的设计与实现.docx VIP
- 普通话课件(完整版)教学文案.ppt
- 【清风语文精品课件】2021高中语文《静女》优质课一等奖.pptx
- 某区南1#矿石泊位升级10万吨级散货泊位工程环境影响报告书.pdf
- 2024高中语文教师课程标准考试模拟试卷及参考答案.docx VIP
文档评论(0)