%runge函式插值計算
%參考教材：《數值分析》李乃成，梅立泉，科學出版社
%《計算方法教程》第二版 凌永祥，陳明逵
%runge函式f(x)=(1+25x^2)^(-1),(-1clear;clc;close all;
format short;
fprintf('輸入n：(5,10,15,20)\n');
n=input('');
ni=0:n;
%format short;
fprintf('插值節點xi');xx=-1+2/n*i
fprintf('插值節點xi處的函式值f(xi)');fx=(1+25*xx.^2).^(-1)
yy=fx;
%----------求newton插值多項式的各階差商-------
%注意：陣列下表為1...n+1，區別於課本中的0...n
for k=1:n
for i=n+1:-1:k+1
yy(i)=(yy(i)-yy(i-1))/(xx(i)-xx(i-k));
endend
fprintf('newton插值多項式的各階差商');yy
syms x y
%--------用秦九韶演算法計算多項式的表示式
y=yy(n+1);
for i=n:-1:1
y=y*(x-xx(i))+yy(i);
endfprintf('newton插值多項式的表示式:');simplify(y)
%----------（自然）三次樣條插值多項式-
matlab 插值函式
matlab中的插值函式為interp1，其呼叫格式為 yi interp1 x,y,xi,method 其中x，y為插值點，yi為在被插值點xi處的插值結果 x,y為向量，method 表示採用的插值方法，matlab提供的插值方法有幾種 method 是最鄰近插值，linear 線性插值 spl...
Matlab插值函式
x 0 2 pi y sin x xx 0 0.5 2 pi interp1對sin函式進行分段線性插值，呼叫interp1的時候，預設的是分段線性插值 y1 interp1 x,y,xx figure plot x,y,o xx,y1,r title 分段線性插值 臨近插值 y2 interp1 ...
MATLAB 牛頓插值法程式
function p newton x,xi,yi,n ton n,n 是乙個二位陣列，用來儲存newton插值多項式的表。ton第一列儲存的是yi,函式值。x 11.5 n 5 x1 10 for i 1 1 n xi 1,i x1 ton i,1 yi 1,i ton i,1 log x1 x1...