matlab數值分析??插值法

?

1 拉格朗日插值法

function yh=lagrange(x,y,xh) n=length(x); m=length(xh); yh=zeros(1,m); for j=1:m;for i=1:nxp=x([1:i-1 i+1:n]);yh(j)=yh(j)+y(i)*prod((xh(j)-xp)./(x(i)-xp)); %注意區分yh和yend end

調用程序

x=[11,12,13]; y=[2.3979,2.4849,2.5649]; xh=11.75; yh=lagrange(x,y,xh)

2 牛頓插值法

?

function yh=newtonPol(x,y,xh) n=length(x); p(:,1)=x; p(:,2)=y; for j=3:n+1p(1:n+2-j,j)=diff(p(1:n+3-j,j-1))./(x(j-1:n)-x(1:n+2-j))'; %求差商表 （注意這里有一個 ’ 符號，與差商表不一樣的地方） end q=p(1,2:n+1)'; %求牛頓法的系數--取第一行 yh=0; m=1; yh=q(1); for i=2:nm=q(i);for j=2:im=m*(xh-x(j-1)); %求牛頓法中各多項式值(xh-x0)…(xh-xn-1)endyh=yh+m;%求和 end

調用程序?

x=[11,12,13]; y=[2.3979,2.4849,2.5649]; xh=11.75; yh= newtonPol(x,y,xh)

?

?

《新程序員》：云原生和全面數字化實踐50位技術專家共同創作，文字、視頻、音頻交互閱讀

總結

以上是生活随笔為你收集整理的2021-01-07 matlab数值分析  插值法 拉格朗日插值法 牛顿插值法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。