本帖最后由 wll 于 2011-11-18 10:26 编辑
& G4 I5 i5 X w4 I
- I' F- r* r: U! C: o& Ff(r)=x1*L1(r)+x2*L2(r)+...+x8*L8(r)
4 _! {" [- [3 F' K; g3 B其中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)
* Y# R0 L' k# c% s4 t! m: d; u3 \如果一项一项的按顺序写下去,可以实现这个式子,但是很繁琐。" s- ^) E9 r2 ]& b! M. z" K+ m5 G
大家有什么办法,用什么函数或者语句可以实现拉格朗日插值多项式?
8 r2 J( z) I9 \* o+ a/ L各位谁有思路给一下提示,先谢过啦^_^* @( _( Y0 T6 h/ X2 n/ U) z7 j
; q: I+ r8 J7 q9 j9 z7 b另: , o, ?( M7 I- r. Q" S
就比如说我想实现这样一个循环
. g! ^/ ^& B$ I% E. mDO/trm:,i,1,5
+ \6 J+ i; I! J! N& Fx(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)))
2 X1 X* _2 N, N: E, n2 |trm:5 V: b) w# V4 Y- A' y
: R" J8 x) h- S2 o: R& @不论i等于几时,分母都有0项,改怎么办呢?0 J, E- `& h" Y; I6 S/ a
|