拉格朗日等方法测温度.docVIP

实验目的 1、掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析。 2、 掌握用MATLAB及梯形公式、辛普森公式计算数值积分。 3、 通过实例学习用插值和数值积分解决实际问题。 二、实验内容 1、在1-12的11小时内，每隔1小时测量一次温度，测得的温度依次为：5，8，9，15，25，29，31，30，22，25，27，24。试估计每隔1/10小时的温度值。 计算过程： x=1:12; y=[5,8,9,15,25,29,31,30,22,25,27,24]； plot(x,y,k.,markersize,11) axis([1 12 0 40]) grid; hold on t=1:0.1:12; u=spline(x,y,t); plot(t,u,r--) fprintf(u=%.2f ,u) u=5.00 u=5.62 u=6.15 u=6.60 u=6.97 u=7.26 u=7.50 u=7.68 u=7.82 u=7.93 u=8.00 u=8.05 u=8.10 u=8.14 u=8.18 u=8.24 u=8.31 u=8.42 u=8.57 u=8.76 u=9.00 u=9.31 u=9.67 u=10.11 u=10.60 u=11.17 u=11.80 u=12.49 u=13.26 u=14.09 u=15.00 u=15.97 u=17.00 u=18.06 u=19.14 u=20.23 u=21.29 u=22.31 u=23.29 u=24.19 u=25.00 u=25.71 u=26.32 u=26.85 u=27.30 u=27.68 u=28.01 u=28.30 u=28.55 u=28.78 u=29.00 u=29.21 u=29.42 u=29.63 u=29.84 u=30.04 u=30.24 u=30.44 u=30.63 u=30.82 u=31.00 u=31.18 u=31.33 u=31.45 u=31.52 u=31.52 u=31.44 u=31.27 u=30.98 u=30.56 u=30.00 u=29.29 u=28.46 u=27.55 u=26.59 u=25.62 u=24.69 u=23.82 u=23.06 u=22.44 u=22.00 u=21.77 u=21.73 u=21.86 u=22.12 u=22.49 u=22.94 u=23.44 u=23.97 u=24.50 u=25.00 u=25.45 u=25.83 u=26.17 u=26.45 u=26.67 u=26.84 u=26.96 u=27.02 u=27.04 u=27.00 u=26.91 u=26.78 u=26.59 u=26.36 u=26.08 u=25.75 u=25.38 u=24.97 u=24.50 u=24.00 输出的图形为： 结论：利用样条插值法可以把问题解决，还可以用分段线性插值法和拉格朗日插值法。 2、已知 ，利用数值积分求π值，要求精确到小数点后10位。 计算过程及结果： x=0:0.01:1; y=4*1./(1+x.*x); a=trapz(x,y); pi=vpa(a,11) pi = 3.1415759869 结果分析：用梯形法很快就求出了π的近似结果，当然此值的近似结果还可以用辛普生方法（quad）和高阶方法(quadl)求得。过程及结果如下： f=inline(4*1./(1+x.*x)); a=quad(f,0,1); pi=vpa(a,11) pi = 3.1415926829 f=inline(4*1./(1+x.*x)); a=quadl(f,0,1); pi=vpa(a,11) pi = 3.1415926829