信计08级数值方法计算实习题.docVIP

信计08级数值方法计算实习题 要求：1、用Matlab语言或你熟悉的其他算法语言编写程序，使之尽可能具有通用性；2、根据上机计算实践，对所使用的数值方法的特点、性质、有效性、误差和收敛性等方面进行必要的讨论和分析；3、完成计算后写出实验报告，内容包括：课题名称、解决的问题、采用的数值方法、算法程序、数值结果、对实验结果的讨论和分析等；4、特别说明：严禁抄袭，否则一经发现，所有雷同实验报告最多评为及格。 一、下表给出了飞行中鸭子的上部形状的节点数据，试用三次样条插值函数（自然边界条件）和20次Lagrange插值多项式对数据进行插值。用图示出给定的数据，以及和。 0.9 1.3 1.9 2.1 2.6 3.0 3.9 4.4 4.7 5.0 6.0 1.3 1.5 1.85 2.1 2.6 2.7 2.4 2.15 2.05 2.1 2.25 7.0 8.0 9.2 10.5 11.3 11.6 12 12.6 13.0 13.3 2.3 2.25 1.95 1.4 0.9 0.7 0.6 0.5 0.4 0.25 解： x=[0.9 1.3 1.9 2.1 2.6 3.0 3.9 4.4 4.7 5.0 6.0 7.0 8.0 9.2 10.5 11.3 11.6 12 12.6 13.0 13.3]; y=[1.3 1.5 1.85 2.1 2.6 2.7 2.4 2.15 2.05 2.1 2.25 2.3 2.25 1.95 1.4 0.9 0.7 0.6 0.5 0.4 0.25]; %（1）三次样条插值法 xi=0.9:0.01:13.3;yi=interp1(x,y,xi,spline); xi=0.9:0.01:13.3;yi=interp1(x,y,xi,spline); title(试验一--三次样条插值图示) pp=spline(x,y) pp = form: pp breaks: [1x21 double] coefs: [20x4 double] pieces: 20 order: 4 dim: 1 pp.coefs ans = 0.7735 -0.9995 0.7760 1.3000 0.7735 -0.0714 0.3477 1.5000 -2.7894 1.3209 1.0974 1.8500 -0.4585 -0.3528 1.2910 2.1000 0.4489 -1.0405 0.5944 2.6000 0.1738 -0.5018 -0.0225 2.7000 0.0783 -0.0325 -0.5033 2.4000 1.3141 0.0850 -0.4771 2.1500 -1.5812 1.2676 -0.0713 2.0500 0.0431 -0.1555 0.2623 2.1000 -0.0047 -0.0261 0.0808 2.2500 -0.0245 -0.0401 0.0146 2.3000 0.0175 -0.1135 -0.1390 2.2500 -0.0128 -0.0505 -0.3358 1.9500 -0.0201 -0.1003 -0.5319 1.4000 1.2094 -0.1485 -0.7310 0.9000 -0.8279 0.9400 -0.4935 0.7000 0.0122 -0.0535 -0.1389 0.6000 -0.2960 -0.0316 -0.1900 0.5000 -0.2960 -0.3867 -0.3573 0.4000 所以所得方程为 %（2）用拉格朗日法插值 %定义Lagrange程序 function f=Language(x,y,x0) syms t; if(length(x)==length(y)) n=length(x); else disp(xoíyμ???êy2??àμè￡?); return; end f=0.0; for(i=1:n) l=y(i); for(j=1:i-1) l=l*(t-