本帖最后由 wll 于 2011-11-18 10:26 编辑
! @0 L) \6 i* c8 v( H
5 ]% v+ z! [2 L/ |7 L4 sf(r)=x1*L1(r)+x2*L2(r)+...+x8*L8(r)6 o' B. V3 q' k3 I& M
其中Li(r)=(r-r1)...(r-r(i-1))(r-r(i+1))...(r-rn)/(ri-r1)...(ri-r(i-1))(ri-r(i+1))...(ri-rn)
b2 e' u! _0 Z7 F* i如果一项一项的按顺序写下去,可以实现这个式子,但是很繁琐。
) ?$ j" W7 _7 _! {* o大家有什么办法,用什么函数或者语句可以实现拉格朗日插值多项式?$ j5 F9 o4 p8 T- n3 a6 d, B
各位谁有思路给一下提示,先谢过啦^_^
. w, c" H' u! L1 D. Z5 {# K2 @
8 ~& p0 c9 Y& A! ]另:
' A% Q; C6 d2 B1 U5 \, T就比如说我想实现这样一个循环
% ~* O6 h4 h; ADO/trm:,i,1,5
) C' c, U/ o3 Y- n) _! Zx(i)=n(1)*(e(i)-m(1))*(e(i)-m(2))*(e(i)-m(3))/((m(i)-m(1))*(m(i)-m(2))*(m(i)-m(3)))' m! N6 W. j5 Z! W: |
trm:7 m0 k# i- n9 L
2 ^; F5 l7 o5 P. |+ Q: W不论i等于几时,分母都有0项,改怎么办呢?8 V" a% R1 ]6 y9 E- m6 N0 |: [
|