数学实验_三_陈雨.docx

实验三 插值与数值积分陈雨 2010012199 基应02【实验目的】1掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析。2掌握用MATLAB及梯形公式、辛普森公式计算数值积分。3通过实例学习用插值和数值积分解决实际问题。【实验内容】第10题10.1问题分析题目给定的数据点是有限的，要想确定更多的数据，就要运用插值方法。而加工断面的面积，则应通过数值积分分别求出上下轮廓线对x轴围城的面积，然后做差求得。10.2模型建立与求解记机翼形状曲线的横轴坐标为x，与x对应的机翼的上轮廓线纵轴值为y1，下轮廓线纵轴值为y2。为了画出机翼形状曲线，对x轴每隔0.1进行插值，于是可用不同插值方法画出与之对应的机翼形状。通过数值积分分别求出上下轮廓线对x轴围城的面积，然后做差求得剖面面积。分别运用拉格朗日、分段线性、三次样条插值方法对原始数据插值。由于拉格朗日方法得到的结果与实际偏差较大，其收敛性难以得到保证，故下面讨论后两种插值方法。完成加工所需数据如下表：xxyy11yy12yy21yy220.00.00000.00000.00000.00000.10.06000.04000.10890.04990.20.12000.08000.21340.09900.30.18000.12000.31370.14740.40.24000.16000.40970.19510.50.30000.20000.50180.24210.60.36000.24000.58980.28840.70.42000.28000.67400.33400.80.48000.32000.75450.37880.90.54000.36000.83140.42301.00.60000.40000.90470.46651.10.66000.44000.97470.50941.20.72000.48001.04130.55151.30.78000.52001.10470.59301.40.84000.56001.16510.63381.50.90000.60001.22250.67391.60.96000.64001.27700.71341.71.02000.68001.32870.75231.81.08000.72001.37780.79041.91.14000.76001.42440.82802.01.20000.80001.46850.86492.11.26000.84001.51040.90122.21.32000.88001.54990.93682.31.38000.92001.58740.97192.41.44000.96001.62291.00632.51.50001.00001.65651.04012.61.56001.04001.68841.07322.71.62001.08001.71851.10582.81.68001.12001.74711.13782.91.74001.16001.77421.16923.01.80001.20001.80001.20003.11.82001.22501.82451.23023.21.84001.25001.84801.25993.31.86001.27501.87041.28893.41.88001.30001.89181.31743.51.90001.32501.91251.34543.61.92001.35001.93251.37273.71.94001.37501.95191.39953.81.96001.40001.97081.42583.91.98001.42501.98941.45154.02.00001.45002.00761.47674.12.02001.47502.02581.50144.22.04001.50002.04391.52554.32.06001.52502.06201.54914.42.08001.55002.08031.57224.52.10001.57502.09891.59474.62.12001.60002.11791.61684.72.14001.62502.13741.63834.82.16001.65002.15751.65944.92.18001.67502.17841.67995.02.20001.70002.20001.70005.12.22501.71502.22251.71965.22.25001.73002.24591.73875.32.27501.74502.27001.75735.42.30001.76002.29481.77545.52.32501.77502.32011.79305.62.35001.79