基于移动最小二乘法的曲线曲面拟合-曾清红

2002102110102020973G19990328051978-(xi,yi),i=1,2,,nxiy=f(x)f()()()f(x)()[1],[2]()()(MLS)()LancasterSalkauskasMLS[3]Belytschko[4]1a(x)p(x)a(x)x2CompactSupportxyxxxw(x)()JOURNALOFENGINEERINGGRAPHICSNo.1230027(MovingLeast-SquaresMLS)(LS)MLSTP391A1003-0158(2004)01-0084-0611.1f(x)∑===miiixxpxpxxf1T)()()()()(aa1T21)](,),(),([)(xaxaxaxmΛ=axT21)](,),(),([)(xpxpxpxpmΛ=km3,],,1[)(T==myxxp2a6,],,,,,1[)(T22==myxyxyxxp2bL2(],,[21nxxxΛ=xL212221niix/⎛⎞⎜⎟⎜⎟=⎝⎠=∑x)∑∑==--=--=nIIIInIIIyxxpxxwyxfxxwJ12T12])()()[(])()[(a3n()fxIyI=xx()IIyy=x()Iw-xxIx()xa33a0)()()(=-=∂∂yxBxxAJaa4yxBxAx)()()(1-=a5∑=-=nIIIIxpxpxxwxA1T)()()()(6)]()(,),()(),()([)(2211nnxpxxwxpxxwxpxxwxB---=Λ7],,,[21TnyyyyΛ=8(5)(1)MLS∑===nIkIkIyxÖyxxf1)()()(9)(xÖkk)()()(],,,[)(1T21xBxAxpxÖknkkk-==Λ100k=(){1}px=Shepard∑==-=nJJIIxxwxxwxÖ1Shepard)()()(11p(x)9f(x)rCp∈swC∈min()rsfC,∈1.2w(xxI)x(x)(1)smaxxw(xxI)2Ixx-w(xxI)C11Cs=xxImaxsss=122⎪⎪⎪⎩⎪⎪⎪⎨⎧-+-+-=)1(0)121(344434)21(4432)(3232sssssssssw126A(x)13maxsx1229()1112x(a)x(b)x(c))(xÖk(d)x34()33.1(MLS)x=[0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0]y=[0,4,5,14,15,14.5,14,12,10,5,4]p(x)=[1,x]T33(a)~(c)MLS13a23b33c3MLSX