Matlab語言 關於拉格朗日插值法實現求值

n+1個節點的基函式是n次代數多項式。

基函式在節點取值為0或者1，即第i個節點xi基函式為1，其他都為0

基函式和每個節點有關，與被插值函式無關

**

網上有許多大神寫的**，都大同小異。我這稍微改了一下，寫的囉嗦一點，多多大家指教

function yn=lan(x,y,xi)

%untitled summary of this function goes here
% detailed explanation goes here
m=length(x);
n=length(y);
t=ones(1,xi);%生成基函式
s=0;
if m~=n
error('輸入錯誤，x,y長度不一致');%檢查輸入是否完整
endfor i=1:n;%利用2層for迴圈將公式用**實現
for j=1:n
if j~=i%當j不等於i時，進行連乘
t=t.*(xi-x(j))/(x(i)-x(j));
endend
s=s+t*y(i);%進行連續相加
endyn=s;%得出結果
end

拉格朗日插值與拉格朗日反演

模板 拉格朗日插值 拉格朗日插值法 f x sum limits 我們先把右邊那部分提出來看 ell x prod x x cdots x x cdots x 舉個例子吧 有二次函式上的三點 f 4 10,f 5 5.25,f 6 1 求 f 18 求出三個基本式 ell x ell x ell x...

拉格朗日插值

拉格朗日插值基函式 li x x x 0 x xi 1 x xi 1 x x n x i x0 xi xi 1 xi xi 1 xi xn 拉格朗日差值函式 ln x i 0 nyil i x 其中，x為缺失值對應的下表序號，ln x 為缺失值的插值結果，xi 為缺失值yi 的下表序號。對全部缺失值...

拉格朗日插值

function s larg1 x,y,xi m length x 求出插值節點向量長度 n length y if m n error 向量x與y的長度必須一致 這裡肯定一致，只是為了消除直接選取x,y資料的時候出錯而設定的 end s 0 for i 1 n z ones 1,length x...