建模实验拉格朗日插值.docx

建模实验 插值方法 【实验目的】1．掌握用MATLAB计算拉格朗日、分段线性、三次样条三种插值的方法，改变节点的数目，对三种插值结果进行初步分析；2．掌握用MATLAB及梯形公式、辛普森公式计算数值积分；3．通过实例学习用插值和数值积分解决实际问题。【实验内容】题1（3）选择一些函数，在n个节点上（n不要太大，如5~11）用拉格朗日、分段线性、三次样条三种插值的方法，计算m个插值点的函数值（m要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较，适当增加n，再作比较，由此作初步分析。下列函数供选择参考：（3）y=cos10x,-2≤x≤2解：当n=8，m=50时，在matlab中编程如下：x=-2:0.5:2;y=(cos(x)).^10;a=-2:4/50:2;yy=(cos(a)).^10;y1=lagr(x,y,a);y2=interp1(x,y,a);y3=interp1(x,y,a,spline);[a;yy;y1;y2;y3]数值输出如下： -2 0 0 0 0-1.92 2.1987e-005 -0.99289 0 0.0019267 -1.84 1.771e-006 -1.2394 0 0.0027391 -1.76 5.538e-008 -1.078 8.0997e-005 0.0027867 -1.68 2.3645e-010 -0.74009 5.6075e-005 0.0022632 -1.6 4.5057e-016 -0.37345 3.1153e-005 0.0013621 -1.52 1.1388e-013 -0.063296 6.2305e-006 0-1.44 1.4242e-009 0.1509 0 -0-1.36 1.6085e-007 0.26129 0 -0.0016693 -1.28 3.7541e-006 0.28038 0 -0.0021436 -1.2 3.9028e-005 0.23209 0.0012721 -0.0020271 -1.12 0 0.14495 0.0016113 -0.001126 -1.04 0.0011051 0.047167 0.0019505 0.0007532 -0.96 0.0038502 -0.03686 0.023626 0.0039155 -0.88 0.011026 -0.088593 0.066638 0.011223 -0.8 0.026946 -0.096401 0.10965 0.028094 -0.72 0.057684 -0.055717 0.15266 0.06006 -0.64 0.11022 0.031477 0.19567 0.11265 -0.56 0.19061 0.15766 0.23869 0.19139 -0.48 0.30145 0.3111 0.30011 0.30176 -0.4 0.43942 0.4774 0.41676 0.44195 -0.32 0.59394 0.64102 0.5334 0.59603 -0.24 0.74766 0.78685 0.65005 0.74643 -0.16 0.87937 0.90159 0.7667 0.8756 -0.08 0.96847 0.97483 0.88335 0.96597 0 1 1 1 1 0.08 0.96847 0.97483 0.88335 0.96597 0.