fν
(r, s)t
pν,t(r, s) = ρt(r, s) + c[dtν,t(r)dtν,t(s)] + ǫν,t(r, s)
φν,t(r, s) = ρt(r, s) + c[δtν,t(r)δtν,t(s)]
+λν[ϕν,0(r)ϕν,0(s)] + λνNν(r, s) + εν,t(r, s)
ρt(r, s)
s tτ
r t
λν
Nν(r, s)
(ν, t)r s
dtν,t(r)
dtν,t(s)δtν,t(r)δtν,t(s)
c ϕν,0(r)ϕν,0(s)
r s
ǫν,t(i, j)
εν,t(i, j)
φt(r, s) = ρt(r, s) + c[δtt(r)δtt(s)]
+λ[ϕ0(r)ϕ0(s)] + λN(r, s) + εt(r, s)
a:= b a
b r1
r2s1, s2,...,sn
ϑj:= ϑ(r2, sj)ϑ(r1, sj)
ϑ
φj
t=ρj
t+c[δtt(r2)δtt(r1)]
+λ[ϕ0(r2)ϕ0(r1)] + λaj+εj
t
aj:= Nj
aj
sk
ϑj
k:= ϑjϑk(i=j)
ϑ
j=k
φj
t;k=ρj
t;k+λaj
k+εj
t;k(aj
kZ)
aj
k
ϑj
k
ϑkϑj
r
0
ϑj
k
ϑj
ϑk
ϑj
0:= ϑjϑ0
ϑ0ϑj
ϑ0:= 1
n
n
j=1
ϑj
r
ϑ0
ϑj
0
ϑj
r
0
ϑj
0
ϑj
ϑ0
s0
k
n
j=1
|ϑj
0|2
j=k
|ϑj
k|2
ϑj
0
φj
t;0 =ρj
t;0 +λaj
0+εj
t;0 (aj
0Q)
aj
0
ϑd(ri, sj) :=0i= 1j=k;
ϑj
k
ϑ(ri, sj)
ϑr(ri, sj) :=0i= 1;
ϑj
0
ϑ(ri, sj)
ϑ(1)
c(ri, sj) = ϑ(ri, sj)1
2
2
i=1
ϑ(ri, sj)
= (1)i1
2[ϑ(r2, sj)ϑ(r1, sj)]
= (1)iϑj
2
ϑ(2)
c(ri, sj) = ϑ(1)
c(ri, sj)1
n
n
j=1
ϑ(1)
c(ri, sj)
= (1)i1
2ϑj1
n
n
j=1
ϑj
= (1)iϑj
0
2
ϑc(ri, sj) := (1)iϑj
0
2
ϑ(r, s)
n n
ϑ
E:= R2nϑ
E Eψ
ϑEψ =p
ψ=φ
VψVψ
EVψ
ψ E+
ψ
ϑ|ϑE+
ψ:= (ϑ·V1
ψϑ)E
E+E+
ψ
E0ϑ(ri, sj)
ϕ(sj)ϕ(ri)
En + 1
E0E Ec
E0E+E+
c
EcE+
c
2n(n+ 1)=n1Ec
ϑc
E ϑc
E0
2
i=1
n
j=1
[ϕ(sj)ϕ(ri)](1)iϑj
0=
n
j=1
ϑj
0ϕ(sj)
2
i=1
(1)i
2
i=1
(1)iϕ(ri)
n
j=1
ϑj
0
2
i=1(1)i= 0n
j=1 ϑj
0= 0
ϑc
n1
Pc
E Ec
ϑcϑ Ec
ϑc:= Pcϑ Pc
b:= {ej}n
j=1 Rn
ϑ:=n
j=1 ϑjejϑj
F:= Rn
E+F
:= ϑ()j:= ϑj
S
F E
(Sϑ)(ri, sj) := (1)iϑj
ϑF ϑ:=Sϑ/2
=ϑS
Vψ
ψ:=Sψ Vψ
FVψ
Vψ=SVψS
S
S=VψS
SS=Vψ
SF E+
ϑE+ϑF
(·ϑ)F=ϑ|SϑE+
(·ϑ)F=
n
j=12
i=1
(1)iϑ(ri, sj)ϑj
(·ϑ)F=
2
i=1
n
j=1
ϑ(ri, sj)(Sϑ)(ri, sj)
(·ϑ)F=(ϑ·Sϑ)E=ϑ·[V1
ψVψ]SϑE
(·ϑ)F=ϑ|VψSϑE+
S=VψS
F0ϑF
ϑj
F0F
Fr:= ϑF:n
j=1 ϑj= 0
F0Frn1
Q0QrF F0
Fr
(Q0ϑ)j=ϑ0(Qrϑ)j=ϑjϑ0
ϑ0ϑj
ϑr:= Qrϑ
ϑrbn
ϑj
0Fr
ϑcϑr
ϑc=Sϑr/2
ϑrϑEr
n
j=1
|ϑj
0|2=
n
j=1
|ϑjϑ0|2= inf
ϑoR
n
j=1
|ϑjϑo|2
ZnFrn1
Lrb
LrarLr
n aj
0
ϑr
E Er
FrE
k
F
Fd:= {ϑF:ϑk= 0}
FdRn1Qd
F FdF0
Qd
(Qdϑ)j=ϑjϑk
ϑd:= Qdϑ
bd:={ej}j=kFd
ϑdbdn1
ϑj
kFd
ZnFdn1
Ldbd
adLd
n1aj
kj=k
Ld=QdLdLdQdZnLd
ZnQdZnLd
Ld=QdZnLr:= QrZnQr=QrQd
Lr=QrLd
ϑd
E Ed
FdE
E+Fr
Sr:= QrS
Sr
S Sr
ker SrE0
E0ker Srdim E0=n+ 1
dim(ker Sr)=dim Edim Fr=2n(n1) = n+ 1
E+Fd
Sd:=QdS
SrSdker Sd=E0
FdFr
QrFdFdFr
Rϑd:= Qrϑd
FdFrLdLr
FrFd
Dϑr:= Qdϑr
DR
ϑkϑ0ϑk
(DRϑd)j= (ϑj
kϑ0)(ϑk
kϑ0) = ϑj
k
erj:= Rejj=k Fr
Lrbrd:= Rbd
ϑrFr
ϑj
kϑd=Dϑr
ϑr=Rϑd=R
j=k
ϑj
kej=
j=k
ϑj
kRej=
j=k
ϑj
kerj
arLrn1
aj
kad=Dar
F Lr
arn aj
0
TF Fd
FrϑFrϑ
Fd(ϑ·ϑ)F= (ϑ·Rϑ)F= (Tϑ·ϑ)F
TRF R=T
(Rϑr)j=ϑj
0(j=k); (Rϑr)k= 0
DR DR
(Dϑd)j=ϑj
k(j=k); (Dϑd)k=
j=k
ϑj
k
Vψd
bdψdVψr
R
D
❅■
0
ϑϑ1
,...,ϑn
ϑdϑ1
k,...,ϑn
k
ϑrϑ1
0,...,ϑn
0
F:= Rn
Fd
ekFr
F0
rrr
s
s
s
ad
ar
sLd
rLr
❅■T
ekRn
ϑk
k= 0
Pn
j=1 ϑj
0= 0R
F FrFdR
D
F FdF0FrD
R
F FdFrR=T
D D
LdLr
b
ψrV
ψdFd
VψdV
ψrFr
VψrQrQr
b Vψr=QrVψQT
r=QrVψQr
V
ψrFr
V
ψrϑ=QrV
ψϑ(ϑFr)
FdFr
ϑ
d|ϑdFψ;d+ := (ϑ
d·ϑd+)Fϑd+:= V1
ψdϑd
ϑ
r|ϑrFψ;r+ := (ϑ
r·ϑr+)Fϑr+:= V1
ψrϑr
E+E+
ψ
ψ Fd+Fψ;d+DD
ψ Fr+Fψ;r+RD
ψ
V
ψd=DV
ψrD
V1
ψd=RV1
ψrR
ϑd+=Rϑr+
ϑr+=Dϑd+
ϑ
d|ϑdFd+ = (ϑ
d·Rϑr+)F
ϑ
d|ϑdFd+ =(Rϑ
d·ϑr+)F
ϑ
d|ϑdFd+ =ϑ
r|ϑrFr+
ϑ
r=Rϑ
dϑr=Rϑd
ϑd2
Fd+ =ϑr2
Fr+ ϑr=Rϑd
ϑ E
ϑ+
cϑ
E+
c
ϑ+
c:= P+
cϑ
P+
c
ϑr:=Srϑ ϑ
Srϑ=ϑrE0
ϑ+
cE+
ϑrϑ+
c
SrS+
r
ϑ+
c=S+
rϑr
ϑd=Dϑrϑ+
c=S+
dϑd
ϑ+
cϑrϑd
E+S+
rS+
d
ϑc+
ϑc+:= Sϑr+
ϑ+
c=Vψϑc+
ϑ+
c2
E+= (ϑ+
c·ϑc+)E=ϑr2
Fr+
Sr
S+
r=S
r(SrS
r)1
ϑFrSr=QrS
S
rϑ= (QrS)ϑ=SQ
rϑ=SQrϑ=Sϑ
S=V
ψS
SrS
rϑ=QrSSϑ=QrV
ψϑ=V
ψrϑ
S+
r=S
rV1
ψr
ϑ+
c=S
rV1
ψrϑr=SV1
ψrϑr=V
ψSV1
ψrϑr
ϑ+
c=Vψϑc+
ϑ+
c2
E+=ϑ+
c|ϑ+
cE+
=ϑ+
c|Vψϑc+E+= (ϑ+
c·ϑc+)E
ϑcϑ+
cEc
(ϑ+
c·ϑc+)E= (ϑc·ϑc+)E
= (ϑc·Sϑr+)E
= (c·ϑr+)F
c=ϑr
ϑ+
c2
E+= (ϑr·ϑr+)F=ϑr2
Fr+
ψ
Vψ=η(ri, sj)σ2
ψ
σ2
ψη(r, s)
ψ:=
Vψ= (ηjσ2
ψ)ηj:=η(r1, sj) + η(r2, sj)
ϑj
r+ϑj
r
ϑr+ϑr
ϑj
r+=1
ηjσ2
ψ
(ϑj
rδϑ) (ϑj
rϑj
0)
δϑ:=
n
j=1
µjϑj
rµj:=
1
ηj
n
j=1 1
ηj
ϑ+
c=ησ2
ψϑc+=ησ2
ψSϑr+
ϑr+:= V1
ψrϑr
V
ψrFr
V
ψrϑ=ϑrFrV
ψrϑV
ψϑ
F0
ϑϑr
ηjσ2
ψϑj=ϑj
rδ
δR
ϑj=1
ηjσ2
ψ
(ϑj
rδ)
ϑFrn
j=1 ϑj= 0
δδϑ
η(ri, sj)ηj= 2
µj= 1/nj δϑ= 0
ϑrFr+
ϑr2
Fr+ =
n
j=1
1
ηjσ2
ψ
(ϑj
rδϑ)2
ϑr2
Fr+ = (ϑ+
c·ϑc+)E
=
2
i=1
n
j=1
η(ri, sj)σ2
ψϑ2
c+(ri, sj)
=
n
j=1
2
i=1
η(ri, sj)1
η2
jσ2
ψ
(ϑj
rδϑ)2
=
n
j=1
[η(r1, sj) + η(r2, sj)] 1
η2
jσ2
ψ
(ϑj
rδϑ)2
η(r1, sj) + η(r2, sj) = ηj
ϑr[ψν,t]
[ψν,t]Fr+
[ψν,t]E+
c
Fr+Fd+Er+Ed+Ec+
E+
c
η(ri, sj)
ϑ+
c=ϑc=1
2Sϑr
E+
c
Ec
0
ϑr
ϑ+
c
ϑ
ϑd
ϑc
ϑc+
E
Ed
Er
E+
c
Ec
E0
Sr
Sd
0
ϑ
ϑd
ϑr
F
Fd
Fr
F0
ϑd+
ϑr+
Srϑc+
E0
E
SrSd
ϑcϑ
E0E Ecϑc=Pcϑ
ϑcϑc
ϑr
#rϑd
#d
SrFrE+
c
E
E+ϑ+
c=S+
r#r=P+
cϑ
ϑ+
c=V
ψϑc+ϑc+:= S#r+#r+:= V1
r#r
Srϑc+=c+=SS#r+= 2#r+
S+
dFdE+
cϑ+
c=S+
d#d
ϑ+
cϑc
E+
cEc
ϑ+
c2
E+=ϑr2
Fr+ =ϑd2
Fd+ =
n
j=1
1
2σ2
ψ
|ϑj
r|2
E+E+
c
SdC
η(ri, sj)
k= 1
V1
ψd
κj,j=1
2σ2
ψ
×1
n
n1j=j
1j=jj, j∈{2,...,n}
ϑFd
ϑ2
Fd+ = (ϑ·V1
ψdϑ)F=
n
j=2
ϑj
k(V1
ψdϑ)j
j= 2,...,n
(V1
ψdϑ)j=1
2σ2
ψϑj
k1
n
n
j=2
ϑj
k
=1
2σ2
ψ(ϑjϑk)1
n
n
j=1
(ϑjϑk)
=1
2σ2
ψϑj1
n
n
j=1
ϑj
=1
2σ2
ψ
ϑj
0
n
j=2
ϑj
k(V1
ψdϑ)j=1
2σ2
ψ
n
j=2
(ϑjϑk)ϑj
0
=1
2σ2
ψ
n
j=1
(ϑjϑk)ϑj
0
=1
2σ2
ψ
n
j=1
(ϑj
0ϑk
0)ϑj
0
n
j=1 ϑj
0= 0
ϑ2
Fd+ =
n
j=1
1
2σ2
ψ
|ϑj
0|2(ϑj
0ϑj
r)
t ξt
ρj
t
ξ2;tr2ξ2;t=
ξ2;t+ξt
ρj
t=ρt(r2, sj)ρt(r1, sj)
ρj
t
ρj
t=ρj
t+ (dj
t·ξt)R3
dj
t
sjr2tJt
jth
dj
tJt
ρt=
ρt+Jtξt
ρt;r =
ρt;r +Jt;rξt(Jt;r := QrJt)
t
xt
xt:= (α, ξt)T
αarFr
t1, t2,...tn
X:= (α, ξ1, ξ2,... ,ξn)T
ξnξtnα
t ytt
yt:= pt;r
ρt;r
φt;r
ρt;r
yt=Atxt+
At:= 0Jt;r
λIαJt;r
Y= (y1, y2,... ,yn)T
ynytn
Y=AX +
A
A:=
·J1;r ·· · ··
λIαJ1;r ··· ··
·· J2;r · · ··
λIα·J2;r · ···
·
·
··
·
··
·
·····
·
·
·
·
··
·
··
·
·····
·
·
····· ·Jn;r
λIα··· · ·Jn;r
x(α,
ξ)T
tn
xn|n=xn|n1+Knvn
vn=ynAnxn|n1
xn|n:= (αn,
ξn)xn|n1:= (αn1,0)
Kntnvn
α
Vbα
ˇα Lr
α
Frf(ϑ) := (ϑ·V1
bαϑ)F
ˇα= argmin
αLr
ααV1
bα
Lr
brd
k
α
αd:= bd
αbrd
αd:= Dαbd
α
brd DVbαDVDbα=Vbαd
Lr
Ld
ˇαd= argmin
αdLd
αdαdV1
bαd
ˇα=Rˇαd
k= 1
α
Vbα
αd=DαVbαd=
DVbαDn1
αd
αj
k=αjαk(j= 2,...,n)
Vbα
(n1)2Vbαd
ˇαd
ˇα=Rˇαd
ˇαj= ˇαj
kˇα0
kˇα0
k:=1
n
j=k
ˇαj
k
ˇα
ξn
wn:= ynAnxn|n=Hnvn
Hn:=IAnKn
In
wpwφ
w
w2:= wp2
Fp;r+ +wφ2
Fφ;r+
ψ=p φ
wψ2
Fψ;r+ =
n
jψ=1
cjψ
cjψ:= 1
ηjσ2
ψ
(wj
ψδwψ)2δwψ:=
n
j=1
µjwj
ψ
w2
β=
jpp
βjpejp,
jφφ
βjφejφ
pφ
{1,...,n}
ρj+c[dt(r2)dt(r1)] + ǫj=
pjβjpjp
pj
pφ
βjpβjφ
w
w=Hδv=H δy
y
ejpejφ
zjpzjφ
yset
=yzjψ
zjp:= (erjp,0)
zjφ:= (0,erjφ)
aset
=a+b a
a+b
erj
ej
rj=1/n j=j
ej
rj=11/n
wejpejφ
fjpfjφ
wset
=wHzjψ
fjp:= Hzjp
fjφ:= Hzjφ
w β
:=
jpp
βjpfjp+
jφφ
βjφfjφ
Mβ wM
fjpfjφ
Ω := Ωpφ
θ
θ0
TLOM :=w2/mm= 2(n1)3
tLOM := Fθ(m,,0) θ
Fm,
TLOM< tLOM
r= 1
Ω = Π =
Π
jψ/cjψw2
cmax := max
jψ/cjψ
κ1
Πr:= {jψ/Ω: cjψκcmax}
r rr r
3p5p3φ5φ
c5φ
cjψκ=0.5
n= 7Ω =
3p5p3φ5φ
5φ
jψΠr
jψ/Π
fjψ:=H·
(erjp,0) ψ=p
(0,erjφ)ψ=φ
gjψ:= fjψΠset
={jψ}Π=
Π∪{jψ}
Πjψ
fjψ
r= 1
{g
q}q<r
q<r
ςq,jψ:=g
q|gjψ
:=
ψ=p,φ
g
q;ψ|gjψ;ψFψ;r+
ψg
q;ψ
g
qgjψ;ψςq,jψ
gjψ
set
=gjψςq,jψg
q
ςq,jψ=g
q|fjψ
gjψg
q
q<r
w
gjψhjψ|whjψ
hjψ:=gjψ/gjψ
|hjψ|w|
γjψ:=gjψ|wjψ̺jψ:=gjψ
gjψ|w:=
ψ=p,φ
gjψ;ψ|wψFψ;r+
gjψ2:=
ψ=p,φ
gjψ;ψ2
Fψ;r+
¯¯
ψ
|γjψ|
¯¯
ψ:=arg max
jψΠr
|γjψ|
χ0θ0/2
χ0:=Nθ0/2(0,1)
|γ¯¯
ψ|> χ0m>0
ωr:= ¯¯
ψset
={ωr}r= 1
∪ {ωr}r>1
γ
r:= γωrg
r:= gωrωr
Ω = {ωq}r
q=1
{g
q}r
q=1
Mr
q=1 γ
qg
q
w
β
r:= βωrf
r:= fωr
|γ¯¯
ψ|< χ0TLOM >5tLOM
|γ¯¯
ψ|0TLOM<5tLOM r>1
|γ¯¯
ψ|0m= 0
g
rf
q
uq,r
g
r=
r
q=1
uq,rf
q
uq,r=
1
̺ωr
qq<r
uq,qςq
rq<r
1
̺ωr
q=r
1qruq,r
rth U
β
q
r
q=1 γ
qg
qf
q
r
q=1
γ
qg
q=
r
q=1
β
qf
q
[γ]γ
q
q= 1r[β]
[β] = U[γ]
β
q
set
=β
q+uq,rγ
rq<r
ur,rγ
rq=r
1qr
ww2
wset
=wγ
rg
rw2set
=w2− |γ
r|2
m
mset
=m1
tLOM
set
=Fθ(m,,0)
TLOM
set
=w2/m
TLOM> tLOM rset
=r+1
KK
ω1,...,ωr
xset
=xK[β]
[β]= U[γ]
x
Vbx
set
=Vbx+ [KU][KU]T
{g
q}r
q=1
[γ]
w
w
w
̟0
ζ0
θ0̟0
θ0= 0.001 ̟0= 0.80
ζ0
Fθ(m,,0) θ
ηj= 2j
w2=1
2σ2
p
n
j=1
|wj
p|2+1
2σ2
φ
n
j=1
|wj
φ|2
jpjφ
|wj
p|
p|wj
φ|
φj= 1,...,n
1
303 121238
ˇaj
fν;kj= 1, . . ., nn=
9k=1
j;fνf1f2
10 0
234 868257 496
3625 263196 104
42 502 8961 419 324
512 155 3239 705 967
62 593 1671 303 294
75 973 7734 346 092
89 056 7407 252 801
99 386 8387 332 507
4745
j= 1,...,9
af1,ˇaf2)
k= 1
j= 1,...,n
k= 1
ˇαset
= ˇaf1ˇαset
= ˇaf2
303.39
120.92238.49
2f1
j= 2
• −1f2
j= 4
TLOM 9.13
tLOM (1.26)
TLOM
(f2; 1p) 5.40
(f1; 4φ) 2.63
(f1; 1p) 1.63
(f1; 8p) 1.06
TLOM TLOM
1.06 tLOM
1.30
TLOM 335.09
tLOM 1.26
TLOM
(f1; 2φ) 108.81
(f2; 4φ) 9.42
(f2; 1p)5.41
(f1; 4φ) 2.82
(f1; 1p)1.74
(f1; 8p)1.13
TLOM1.13
tLOM (1.34)
κ= 1κ= 0
κ= 0
κ0.5
βfν;jψ
n= 9
f1
jψ
p4.806 3.304
φ0.043
f2
jψ
p8.755
βf1;2φ2λ1
βf2;4φ≃ −λ2
f1
jψ
p4.806 3.304
φ0.381 0.043
f2
jψ
p8.755
φ0.248