数值分析课程设计报告书-三次样条插值的三弯矩法.docVIP

数值分析课程设计报告书 院系名称 ： 学生姓名 ： 专业名称 ： 班 级 ： 时间 ：  实验一 三次样条插值的三弯矩法 一、实验目的 已知数据,,及边界条件,求的三次样条插值函数.要求输出用追赶法解出的弯矩向量及的值.画出的图形,图形中描出插值点及分别用‘o’和‘*’标记. 二、实验原理 1.用追赶法求解第二类边界条件的三弯矩方程: 其中. 2.得出样条函数表达式: . 3.计算. 三、实验结果 所用数据: x=[-2.223,-1.987,-1.8465,-1.292,-1.2266,-1.1056,-0.8662,-0.6594,-0.2671,-0.0452,0.5385,1.2564,1.4398,1.5415,1.7646,1.9678,2.236]; y=[0.83995,1.1696,1.3141,1.6992,1.7312,1.7847,1.8708,1.9262,1.9881,1.9997,1.9511,1.7169,1.618,1.5543,1.3871,1.191,0.81662]; d2s1= -4.5000;d2sn= -4.8967; %第二种边界条件 t=[-2.223,-1.9443,-1.6656,-1.3869,-1.1083,-0.82956,-0.55088,-0.27219,0.0065,0.28519,0.56387,0.84256,1.1212,1.3999,1.6786,2.236]; ;(指定计算点) 计算结果: 用追赶法求得的弯矩量为 要计算的16个节点出的值为 M0=-4.5000 M1=-1.9282 M2=-0.5336 M3=-0.8000 M4=-0.4188 M5=-0.4663 M6=-0.4264 M7=-0.3343 M8=-0.3677 M9=-0.3113 M10=-0.3705 M11=-0.4185 M12=-0.6533 M13=-0.7253 M14=-1.0132 M15=-1.2968 M16=-4.8967 S(-2.223)=0.83995 S(-1.9443)=1.2163 S(-1.6656)=1.4617 S(-1.3869)=1.6485 S(-1.1083)=1.7836 S(-0.82956)=1.8818 S(-0.55088)=1.9487 S(-0.27219)=1.9876 S(0.0065)=2 S(0.28519)=1.9865 S(0.56387)=1.9462 S(0.84256)=1.8767 S(1.1212)=1.7767 S(1.3999)=1.6411 S(1.6786)=1.4567 S(2.236)=0.81662 要计算的16个节点出的一阶导数值为: 要计算的16个节点出的二阶导数值为: S1(-2.223)=1.8267 S1(-1.9443)=1.058 S1(-1.6656)=0.76267 S1(-1.3869)=0.57108 S1(-1.1083)=0.41615 S1(-0.82956)=0.29348 S1(-0.55088)=0.18877 S1(-0.27219)=0.089715 S1(0.0065)=-0.0028801 S1(0.28519)=-0.095035 S1(0.56387)=-0.19691 S1(0.84256)=-0.30324 S1(1.1212)=-0.41474 S1(1.3999)=-0.56694 S1(1.6786)=-0.76939 S1(2.236)=-1.8916 S2(-2.223)=-4.5 S2(-1.9443)=-1.5044 S2(-1.6656)=-0.62053 S2(-1.3869)=-0.75438 S2(-1.1083)=-0.46529 S2(-0.82956)=-0.41006 S2(-0.55088)=-0.34356 S2(-0.27219)=-0.36728 S2(0.0065)=-0.31653 S2(0.28519)=-0.34481 S2(0.56387)=-0.37221 S2(0.84256)=-0.39084 S2(1.1212)=-0.40946 S2(1.3999)=-0.60221 S2(1.6786)=-0.90224 S2(2.236)=-4.8967 画出的图形为: 四、实验分析 通过实验结果我们，知道三弯矩法求出满足初始条件的三次样条函数，与其他插值函数的构造相比，三次样条插值法的计算量要小得多。其运行速度比较快结果精确，在解决大量数据时可以采用，在工程，计算领域比较广泛。 实验二 基于正交多项式的最小二乘法拟合 一、实验目