实验一 三次样插值的三弯矩法.docxVIP

实验一三次样条插值的三弯矩法一、实验目的用三次样条插值的三弯矩法,编制第一与第二种边界条件的程序.已知数据,,及边界条件,编程计算的三次样条插值函数.具体要求为:输出用追赶法解出的弯矩向量及的值;画出的图形,图形中描出插值点及分别用‘o’和‘*’标记.二、实验原理（1）由,按公式, ,求出.（2）由边界条件确定和。（3）用追赶法解线性方程组求出（4）由确定,得到S(x).三、实验结果所用数据:x=[-5,-3.8884,-3.0142,-2.3216,-1.7433,-1.3771,-1.1318,-0.86821,-0.58598,-0.4314,-0.0293,0.4314,0.58598,0.86821,1.1318,1.3771,1.7433,2.1242,2.3216,3.0142,3.8884,5];y=[0.03846,0.06204,0.09915,0.1565,0.2476,0.3453,0.4384,0.5702,0.7444,0.8431,0.9991,0.8431,0.7444,0.5702,0.4384,0.3453,0.2476,0.1814,0.1565,0.09915,0.06204,0.03846];d2s1=0.0084;d2sn=0.0084;%第二种边界条件t=[-4.166,-3.233,-2.495,-1.888,-1.469,-1.193,-0.9341,-0.6565,-0.47,-0.1298,0.3162,0.5473,0.7977,1.066,1.316,1.652,2.029,2.272,2.841,3.67];用追赶法求得的弯矩量为要计算的20个节点出的值为M0=0.0084M1=0.019204M2=0.047836M3=0.10891M4=0.2393M5=0.38381M6=0.48683M7=0.50559M8=0.081512M9=-0.55789M10=-2.2251M11=-0.64145M12=0.097766M13=0.50104M14=0.48825M15=0.3828M16=0.24171M17=0.14866M18=0.11118M19=0.047302M20=0.019321M21=0.0084S(-4.166)=0.054449S(-3.233)=0.087288S(-2.495)=0.13838S(-1.888)=0.219S(-1.469)=0.3167S(-1.193)=0.41267S(-0.9341)=0.53401S(-0.6565)=0.69895S(-0.47)=0.81911S(-0.1298)=0.9833S(0.3162)=0.90801S(0.5473)=0.76957S(0.7977)=0.61123S(1.066)=0.46809S(1.316)=0.36609S(1.652)=0.26819S(2.029)=0.1954S(2.272)=0.16229S(2.841)=0.11017S(3.67)=0.069093要计算的20个节点出的一阶导数值为:要计算的20个节点出的二阶导数值为:dS(-4.166)=0.024928dS(-3.233)=0.049505dS(-2.495)=0.095911dS(-1.888)=0.18189dS(-1.469)=0.29464dS(-1.193)=0.40603dS(-0.9341)=0.53267dS(-0.6565)=0.63919dS(-0.47)=0.63031dS(-0.1298)=0.25502dS(0.3162)=-0.51127dS(0.5473)=-0.6502dS(0.7977)=-0.59727dS(1.066)=-0.46734dS(1.316)=-0.35247dS(1.652)=-0.2376dS(2.029)=-0.15484dS(2.272)=-0.11968dS(2.841)=-0.068625dS(3.67)=-0.034911d2S(-4.166)=0.016506d2S(-3.233)=0.04067d2S(-2.495)=0.093623d2S(-1.888)=0.20668d2S(-1.469)=0.34754d2S(-1.193)=0.46113d2S(-0.9341)=0.5009d2S(-0.6565)=0.18747d2S(-0.47)=-0.39823d2S(-0.1298)=-1.8084d2S(0.3162)=-1.0374d2S(0.5473)=-0.087206d2S(0.7977)=0.40029d2S(1.066)=0.49144d2S(1.316)=0.40907d2S(1.652)=0.27689d