牛顿插值法的matlab实现

shizhe-annibaby

function a=new(c,k) , U9 k4 Q p% f[x,y]=size(c); b=zeros(x,2+k); for i=1:x for j=1:2 & }( B+ Z. s$ d8 A# w b(i,j)=c(i,j); end end ' L0 r0 j& a4 U9 V for j=3:2+k ! d7 k+ K/ {1 @) h1 _: j7 E for i=j-1:x 8 d1 P4 `! `: F7 J4 g b(i,j)=(b(i,j-1)-b(i-1,j-1))/(b(i,1)-b(i-j+2,1)); w `' m+ v9 j) a/ J. J end + h3 a7 {' C! n$ b, W! t" t* }# ? end 5 R) F' @: n4 t' h' u/ M a=b

4 {5 C8 k' H' hfunction fx=newton(c,x,k) ( w, F. N8 O1 j8 M+ J6 c b=new(c,k); + q) O6 J) Z* o L7 [: R [p,q]=size(b); * X- ~( w! w5 I$ | fx=0; # H. R& r9 O# d+ L, s( o4 e for i=1:p 1 l; }7 o3 _" F' F t=b(i,i+1) for j=1:i-1 t=t*(x-b(i-1,1)); end fx=fx+t; ) l/ L3 n; ?3 R# O* X" B end fx

zysh666666

是什么程序呀?C语言吗?看不懂.[em06]

coolistone

顶顶顶顶顶顶顶顶顶顶顶顶得到！

540215460

什么东东呀

lyztt1234

看不太懂                        
- v$ h$ T. m0 v

lyztt1234

有没有手算的算法                  
" N% f# Y( _! k6 @) I

lyztt1234

没有tilia              

lyztt1234

看不太懂啊，怎么                 
  L% B" z  g2 d: w