牛顿插值法的matlab实现

shizhe-annibaby

function a=new(c,k) [x,y]=size(c); b=zeros(x,2+k); ) n! U4 d: [! Y+ q3 [' ` for i=1:x for j=1:2 b(i,j)=c(i,j); ' R6 h& C3 I# B v) _ end 9 o8 L% P8 h5 }1 |- c7 l% m% j end 9 f4 m/ G1 l9 X( i* R for j=3:2+k ; @7 V: q/ f( ^* l ` for i=j-1:x 4 m( P& E2 D, p8 G) }8 C3 K b(i,j)=(b(i,j-1)-b(i-1,j-1))/(b(i,1)-b(i-j+2,1)); end # |( n& H/ P6 A7 q: E% f( g% g end ' N+ b- i. q1 g! v+ a0 T a=b

, |8 S- S6 J" B

9 S- N1 q% f$ i: Lfunction fx=newton(c,x,k) 8 D# C9 p8 Q$ x* w- d/ E b=new(c,k); $ U9 ~5 F; H4 T* u* t B [p,q]=size(b); 2 N2 P) ^. X9 O$ P7 j* |+ X* Z1 z . C; \3 M7 p6 g3 c& c* u fx=0; ; q2 K4 d% `1 y; n0 A for i=1:p ' d. U1 \6 o3 E! j$ `) T" L! ^ t=b(i,i+1) for j=1:i-1 t=t*(x-b(i-1,1)); 8 o: |6 I) H" o2 k+ S- A end ! t+ x* _2 o2 s4 O5 m/ \) S2 Y fx=fx+t; / e/ v) \5 j' O9 S- @2 p! i end fx

zysh666666

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

coolistone

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

540215460

什么东东呀

lyztt1234

看不太懂                        

lyztt1234

有没有手算的算法                  
" @# `. ~' Q' q* J

lyztt1234

没有tilia              

lyztt1234

看不太懂啊，怎么                 
4 z& P& W2 s$ U