matlab中，已知原函數(shù)和插值點(diǎn)，怎么求三次拉格朗日插值多項(xiàng)式？函數(shù)YY=Lagrange（x1，Y1，XX）%，本程序?yàn)長agrange 1插值，其中x1，Y1%為插值節(jié)點(diǎn)及節(jié)點(diǎn)上的函數(shù)值，輸出為

matlab中，已知原函數(shù)和插值點(diǎn)，怎么求三次拉格朗日插值多項(xiàng)式？

函數(shù)YY=Lagrange（x1，Y1，XX）%，本程序?yàn)長agrange 1插值，其中x1，Y1%為插值節(jié)點(diǎn)及節(jié)點(diǎn)上的函數(shù)值，輸出為插值點(diǎn)XX的函數(shù)值，%XX可為矢量。Syms xn=length（x1）for I=1:NT=x1t（I）=[]l（I）=prod（（x-t）/（x1（I）-t））%l向量用于存儲插值基函數(shù)endu=sum（l.*Y1）P=simplify（U）%P是簡化的拉格朗日插值函數(shù)（字符串）YY=subs（P，x，XX）clfplot（x1，Y1，“ro”，XX，YY，“*”）