Doctoral thesis
for the degree of doktor ingeniør Trondheim, May 2006
Norwegian University of Science and Technology
Department of Electrical Engineering, Information Technology, and Cybernetics
Geir Werner Nilsen
TOPICS IN OPEN AND CLOSED LOOP SUBSPACE SYSTEM
IDENTIFICATION: FINITE DATA-
BASED METHODS.
!
" #$ #%%$
Æ
!"#
$%" !
!"# &
!"# ' '
&
( !"#
)
" )
*
* *
( ( *
* +
,
!"# ,
!"# ,
' , *
!"#
* ,
-. !"#
' !"# /
(
*
*
!"# -.
'
!"# 0
) !"#
!"#
-.
!"#
*
$ 1 !#! 234456 7 "/ #
# -
$ 1 ! # ! 2344% 6 8
9
! "#! $ % & '!( %
$ 1 ! # ! 2344%6 8
) *+ ! ,
! % # - '!( %
$ 1 ! # ! 2344%6 :
.,+/ ! %
$ 1 !# ! 23446 ; 7 "
; 0 3$ 5
$ !#234456 *
) * ) !"#
! # 234456 *
)
" ! !
!# ! . .
; * 8 ;
* * * ,&<&
$* = # *
" ; * $ 8 " 2$$86
1 $
- 344
3+ > &
3++ &
3+3 /&
33 "
35 ? @
5+ .' + ++
53 .' 3 +3
55 .' 5 +5
5% .' % +% 5 .' + %+ !"# 34
%++ !"# 34
%+3 .'" " 34
%+5 # 3+
%+% & !"# 3+
%+ ' ' 35
%+ !"# 3%
%3 $%" ! 39
%5 . 3A %% -& 5+
% : * 5
% ! 5@
%+ 8 5A
%3 8 -#>" %
%5 8 %@
%% 8 -#>" %
%9 !"# ) A
%9+ '
) ) 5
%93 : / 5
!" + < * 94
3 < * ) 9%
5 < * , 9@
5+ " * @+ % < * , !"# @3
%+ " " ' @5
%3 ' A4 %5 ; A% ! # " + : A@ 3 !"# +44
3+ !"# +43
5 " ' +45
5+ .' + +4%
53 .' 3 ++@ 55 .' 5 +33
5% .' % +3@
5 .' +5%
% : +5@
: +%5
$ %
& "
' ( !
( !
()
* )
'+ $ , - !) .
/ $ , - !) .
2" !6 &
<
"
0 ! 2+AA6 2+AA96
*
; 0 2+AA96
" ! 2. 06
' /
' 1 2344+6
0 2 06 " !
' . 0 ; 0 2+AA96
" 0 2"
0 6 , 23446 , 2344 6 " !
) &
" ! *
* 0 !
2+AA62+AA96 < 7&2+AA96
< *
1
2344+6
" 0 . 0
B
*
" ! . 0
'
'
, 23446) '
)
*
'
; 3
" ! ;
% ' ; 5
" " ! ; %
" ! !"# :* = 2+AA96 !
# 2+AA62344%6 & !"# '
' &
" ! "%% : "%
!"#
) " !
!"#
" % !
( ' >
: ;
* *
* !# 23445 6
*
" !
( ( *
: !"# ,
"%
: " ! !"# ! # 2344%6
; !"#
'!"#
C 234456 " !
&
:
:#D " +
:#D &
Æ
"+ " !
C 234456 " !
" !
" !
-. 2-.6
' :
7& 2+AAA6
* " !
" !
-.
" )!"#
" % :
"
!"# "
% !"# ,
" % " 5+
"5 !"# -. )
!"# "% :
"
" &
& & /& &
( ( >
- ' *
& - ' Ê
E
2
6
& Ê
Ê
E
2
6
23+6
&
>
E E
2
6
2336
; ./ 2336
E
2
6
2356
'
</
E 23%6
E4 236
2
6
E
236
E
2396
- ./ 23%6 236 !# 2+AA96
/ & Ê
Ê
Ê
E 2
62
6
23@6
" / & 2+AA6
E4 23A6
E 23+46
:
E
F
F
23++6
E
F
F
23+36
Ê
Ê
Ê
E G
G *
'
" )
Æ
,
H
E H
F
F
23+56
E H
F
F
23+%6
Ê
' 2
6 E H
Ê '
' ' ' '
' '
'
H
EH
F
: '
H
E H
F
F
23+6
E H
F
F
23+6
Ê
2
6 E
'
E
,
: '
E
E
26
F 26
23+96
26
26E26
F 23+@6
26
26E2 6
F 23+A6
26 26
I
' 26 , 2 6
) 2
266 2 6
'
2+6
' 2+6
: / /
"
/ < ' Ê
/
E 2
6E
23346
< ' EF&
#26
E 26
E 233+6
E
+
: /
./ 23336 /
"/ .
. 2"..6
E
#22
66#2266
F
22
66 2266
23336
/ 23336 26 Ê
' 2
6 Ê
'
; ;
:/ )
./ 23356 / )
) "/ J .2"J.6
E
#22
6# 22 66
F
22
66 22 66
23356
/ 23356 2 6
) 2
6
) ;
;
: '
*
/ /
: !# 234456
E
4 +
49 +
E
43
43
G
E 4
4
E
+ 4
E 4
4
25+6
-
E3
E '
-
E
!
F
2536
E
F
!
2556
!
!
E"
25%6
2G
G
6 E 44+
2
6 E 44+ $E+444
E4+# +
½
¾
256
3
¿
256
: !
+ > -
E 4
E
B 2B 6 K
L
: B $
: > *
E
M E
F2$ 6 2596
M
E
25@6
E
25A6
E 4
E +44
E +4
I E 3
E+ $ E +4 E 3 $
*
!
Æ
E
25+46
Æ
E
25++6
Æ
E
25+36
Æ$
E $
$ 25+56
) '
Æ
E Æ
F Æ
FÆ$
F
25+%6
Æ
E Æ
FÆ
F
25+6
: ) ) '
E
4@94 4
4444 4@94
E
4449
444+4
E
4434
4
E
4 +
25+6
' ) ' )
) ' .
Æ% E 444+ : 4 + $E+444
E4+# +
) 25+%6 25+6 - 2536
2556
"
E +
2
6E4444+
2
6 E44444+ )
) ) ,
)
:? 7&234456
F
E
F
F
25+96
*
E
F"
25+@6
E 4A E + E 4A E 4 '
"
>
1
./
23++6 23+36 E4A E+ E+@ E+ E4 E+
$E3444 34
;
: ; - 2344%6
< 5+
F(z) +
u y K(z)
H(z)
G(z) e 1
e 2
+ +
r +
<5+I " .' %
26E
&26E F
26E+ 25+A6
E5 E3 E4AAA
826 &26
.' 26 &26
./ 23++623+36
E
+
4 4
E
4
E
E
+ 4
25346
) ' '
" 26 E
./ 23++6 23+36
'
E4AA
E+
E4
E45A9@
E43
E4 253+6
"
:
E
+ +4 4+
49 44 4+
4 4 4@
E
4 4
4 +
+ 4
G
E
44 4+
4+ 4
4 43
E
5 4 4
4 + +
E
4 4
4 4
25336
+ + -
E 443
E3 33 -
E443
E3 - ./
2536 2556 2G
G
6 E
444+ 4
4 444+
2
6 E
444+ 4
4 4444
$E+444 E4+# +
> &
/
* ) 01% 2 .0 10 23 4+
'
Ê %E4+3%
%
F+ 2%+6
(" +, '
( +, + '
5 / +, %
5 2 )
)
E
'
'
'
'
'
'
'
'
'
'
'
'
Ê 2%36
5 62 ,+ ++ /
+, %
5
%
2 4+ '
+ ,
6 2
+, (" , )
¼
+, ( , + )
¼
= * '
= * '
0 /
'
(E+ ;
Ê
*
Ê
1 *
* ' ./ 2%56
' ./ 2%%6
E
Ê
2%56
E
Ê
2%%6
' !
' *
*
E
Ê
2%6
! ' *
+
E
Ê
2%6
' ' *
* '
' ' +
* '
: ./ 23++6 23+36
' ' 2,
6 26
,
E
Ê
2%96
' '
2 6
E K
L Ê
2%@6
*
' '
26
E K
L Ê
2%A6
*
* ) ' 2
6 /
2 6
E
4 4 4
4 4
4
Ê
2%+46
* +
*
* ) ' 2
6 /
26
E
4 4 4
4 4
4
Ê
2%++6
* +
*
. ' ' ,
(
* '
,
& ' * ,
'
/ !"#
, !"#
!# 23445 6 !# 2344%6
!"# - ( . /
- + !"# ' .
- ) '.
) ' - ( - .
* ' ' .
)
'
+ ( .0 0 - /
)
'
' .
' ;
/ ,
/ / E . / E.F+
!
./ 23++623+36 >
.' ' ./ I
*
E,
1
F
+
F
2%+36
1
E K
L Ê
2%+56
' >
'F+./ 2%+36 /
./ 23+5623+%6 ' 1
./
2%+36 .' " " 2.""6 .""
* E G
* F G
+ F G
2%+%6
G
E ,
2,
,
6
,
E,
,
Ê
2%+6
G
E K,
L
G
K
4
L Ê
2%+6
G
E K,
L
G
K
4
L Ê
2%+96
$ *
" # !
> .' ' ./ 2%+36 '
' / 1
1
1
E
1
F
+
F
2%+@6
1
1
.'
' ./ 2%+36 ' ' 4 >
*
EK
2
G
L
+
+
*
FK
2
L
2%+A6
G
E ,
,
2%346
2
E ,
2
,
6 2%3+6
2
E ,
2
,
6 2%336
' &!"#
.' " " ./ 2%+%6
./
2%+A6
+
'
'
Æ
.""
: 2
6 Æ
Æ
0 2+AA%6 >
!"#
+
+
*
2%356
-& ./ 2%+A6 ./ 2%356
*
+
+
*
E K
2
G
L
+
+
*
F 2
+
+
*
F3 2%3%6
3
E
+
+
*
2%36
./ 2%36 ) #
./ 2%3%6
*
+
+
*
E K2
G
2
L
+
*
*
+
+
*
F
+
F3 2%36
> &
+
E
+
2+
+
6
+
2%396
( +
2*
+
+
*
6+
E K2
G
2
L
+
*
*
+
+
*
+
F 3
2%3@6
3
E 3 +
2%3A6
./ 2%3A6 ) & *
./ 2%36 +
$ %&
: ; % ' '
,
( * '
,
& ' * ,
'
!# 2+AA6 ; ./ 2%3@6 ./ 2%346
2%336 ./ 2%3@6
' ' ./ 2%3@6
& ' ! ./ 2%3@6
4
E 2*
+
+
*
6+
2%546
4
' '
4
E,
1
2%5+6
1
E K
,
,
,
L
*
+
+
+
*
+
2%536
' 4
' ',
(
1
/ * 4
/ *
1
26
: ' 1
1
EK
L
1
+
+
*
+
+
+
*
+
2%556
' 1
' '
2 6 ' '
26 "(21
6E(
2K L6
1
Æ 2
6 +
*
' Æ
Æ 2
6 +
*
!"# ?#
'
'
+ $ 3
4
E2*
+
+
*
6+
2%5%6
4
E2*
+
+
*
6+
2%56
' /
4
E
G
4
2%56
. /
3 / - (
; " 0 !2"0!6
4
E+)
2%596
+ Ê
) Ê
Ê
+ EK+ +
L ) E
) 4
4 )
EK
L 2%5@6
) Ê (N)N
4
/ (
) ))
/ N)N +
)
(
5 4 + ,
0*-'2
2'6 ' *
+ +
,
E+2+I.0+I(6E+ 2%5A6
% 4 '
' :
E,
4
)
4
E+
4
)
2%%46
4 *
' 0 ( '
' ,
*
E+2+I0+I(6 2%%+6
4 G
G
E,
2,
,
6
,
E4
)
+
2%%36
9 4
!
4
E *
+
+
*
2%%56
' 4
4
4
E
G
4
F
G
+
2%%%6
! ./ 2%%%6
E G
2%%6
E 4
G
4
2%%6
E +
2%%96
" /
G
2 6
2%%@6
@ 5 $
, # /
& !# 23445 6 !#
2344%6
$ Æ )
' 4
N N 0
*
$%" ! 0 !
2+AA6 !# 2+AA62344%6
$ !"#
) ( ) ;
!"# $%" ! E .F+ $%" !
- * !"# =
$%" ! 3 * !"#
. +
+ 6 $ 3
E *
·½·½ +
*
2%%A6
4
E *
+
*
+
2%46
4
E *
+
*
+
2%+6
3 ,* 7 8$ 6 3
/ & &
+
E+)
2%36
& +
$%" !
(
5 *
) "0! +
% ,
,
,
E + )
½
¾
2%56
,
E ,
2+I.0I6 2%%6
+ 6 ' *
,
4
*
E
,
4
+
F
2%6
E
,
,
K 4L ,
2%6
& ,
,
' *
9 +
E ,
4
*
,
4
2 6+
2%96
@ + 5- ) 6
5 )
)
6
E
2
6 2%@6
2 6
E
2 6 2%A6
A 5 $
: # / ,
: $%" ! 0
?# &* !"#
( $%" ! !"#
' ' /
& & $%" ! !"#
& $%" ! ,
# / !"# &
# / , !#
23445 6
$%" ! )
0 ! 2+AA6 ( $%" !
$%" !
' " %% ( &
$%" ! ) 0
! 2+AA6
!" " " #
"
*
2. 06 = *
. 0
. 0
&
" 0 2" 06
1 2344+6 '
' " 0 ."-
; 0 2+AA96
K+
*
+
L
' > ; 0 2+AA96 1
2344+6K+
*
L
1 2344+6
+
'
B
. 0
; 0 2+AA96
" 0 '
' !# 2+AA62344%6 &
4
E*
+
*
+
+
*
7
2%46
7
' *
" 0 * '7
NN '
' K+
*
L
+
*
+
*
½
¾
+
*
2%+6
NN (
; 0 2+AA96 1 2344+6
" 0
+ $ 3
4
E*
6
½
¾
+
*
+
6
½
¾
+
*
2%36
6
E +
*
+
*
2%56
3 ,* 4
4
E+)
2%%6
5 * "
% ,
,
E+ 2%6
4 '
' :
E,
2+I.0I6
,
20F+I2.F+60I6 2%6
,
$ % &
/ &!"#$%" ! " 0
- . 2-.6
" 5+
; : &
'
' ' ( E3
' . ) ' * -.
( E+ )
)"/ . .
: 2- # > " 6 -#>" #
8E 8E N /N -#>"
N /N -#>"
E34 E544
N /N -#>"
< %+
)
!"# - E4 . E / E
$%" !
) .E - E 4 " 0 )
.E@/ EA - E4
< %3
& $%" ! !"#
/
& $%" ! & !"# " 0
* $%" ! !"#
0 100 200 300 400 500 600 700 800 900 1000
−2
−1 0 1 2
u k 1
0 100 200 300 400 500 600 700 800 900 1000
−5 0 5
y k
Discrete time
<%+I
.' + "5+
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) PEM
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) DSR
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) N4SID
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) SIV
< %3I . ;
/ -#>"
.' +" 5+
< %5
"
0 100 200 300 400 500 600 700 800 900 1000
−2
−1 0 1 2
r k 1
0 100 200 300 400 500 600 700 800 900 1000
−0.5 0 0.5
u k
0 100 200 300 400 500 600 700 800 900 1000
−2
−1 0 1 2
y k
Discrete time
<%5I # "
.'
+" 5+
!"# - E4 .E% / E "
$%" !
) .E - E 4 " 0 )
. E / E - E 4
<%%
. " 0
$%" ! !"# "
*
'
/
< % !"# - E 4
.E% / E $%" ! ) . E
- E4 " 0 ) .E / E - E4
/ <
% < %% $ " 0
:
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) PEM
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) DSR
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) N4SID
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) SIV
< %%I . ;
/ -#>" "
.' +
" 5+
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) PEM
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) DSR
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) N4SID
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) SIV
< %I . ;
/
.' + "5+
<%
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) PEM
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) DSR
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) N4SID
0.6 0.8 1 1.2
−0.5 0 0.5
λ(A) SIV
< %I . ;
/
4 .' +" 5+
< % !"#
- E4.E / E9 $%" ! ) .E9
- E 4 " 0 ) . E + / E +
- E 4 :
$ $%" ! !"# " 0
. 0
" 0 * $%" ! !"#
" 0
; 0 2+AA96 1 2344+6
#
/ '
> +
NN &
/
! # 23445 6 $%" !
$%" !
" / & $%" !
!"#
' *
$%" ! !"#
'
' ( )
' "55 .' 5
? 7&234456 -:#" .
* ? 7& 234456
* *
!"#
: -#>" 2- # > " 6
< %9
0 200 400 600 800 1000 1200 1400 1600 1800 2000
−5 0 5
r k 2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
−5 0 5
u k
0 200 400 600 800 1000 1200 1400 1600 1800 2000
−5 0 5
y k
Discrete time
<%9I "
)
.' 5" 55
< %9 (
8 1 %
; 34
<%@ ' ? 7& 234456
= !"#
? 7& 234456
0 0.2 0.4 0.6 0.8 1
−0.5
−0.4
−0.3
−0.2
−0.1 0 0.1 0.2 0.3 0.4
0.5 λ(A) PEM
0 0.2 0.4 0.6 0.8 1
−0.5
−0.4
−0.3
−0.2
−0.1 0 0.1 0.2 0.3 0.4
0.5 λ(A) DSR
< %@I !"# ( E + - E 4 . E +
/ E .F+ 1 %
.' 5 "55
234456 ) 9 !"# /
/ E34 <%A
0 0.2 0.4 0.6 0.8 1
−0.5
−0.4
−0.3
−0.2
−0.1 0 0.1 0.2 0.3 0.4 0.5
λ(A) PEM
0 0.2 0.4 0.6 0.8 1
−0.5
−0.4
−0.3
−0.2
−0.1 0 0.1 0.2 0.3 0.4 0.5
λ(A) DSR
< %AI !"# ( E + - E 4. E +
/ E34 1 % .'
5" 55
$ < %A
!#23445 6
' !"#
/ :
!"# )"/ .
. / E .F+ !"#
"
E 2:6 :E 4%94 ( E+ - E 4
. E +4 / E ++ ) "/ . .
<%+4
0 0.2 0.4 0.6 0.8 1
−0.5
−0.4
−0.3
−0.2
−0.1 0 0.1 0.2 0.3 0.4
0.5 λ(A) PEM
0 0.2 0.4 0.6 0.8 1
−0.5
−0.4
−0.3
−0.2
−0.1 0 0.1 0.2 0.3 0.4
0.5 λ(A) DSR
< %+4I ;
"
E2:6 :E4%94 -. (
E+
!"# ( E+ - E4 .E+4 / E++
.' 5" 55
-. / * *
!"#
-.
" % !"#
'
; 34
? 7& 234456 $ +44
; *
* -:#" .
.
*
* )+ + ! )
' '
!# 23445 6
"53
!"# -. 2-
. 6 ' !
' '
'
-. -. '
; $ ! # 2344% 6
$ !# 2344%6 '
!"#
-.
-#>" 2- # >
" 6
<%++
0 100 200 300 400 500 600 700 800 900 1000
0.8 0.9 1 1.1
r k
0 100 200 300 400 500 600 700 800 900 1000
10 20 30 40
u k
0 100 200 300 400 500 600 700 800 900 1000
0.8 0.9 1 1.1
y k
Discrete time
<%++I "
)
.' 3 "53
' (
!
$+
: / : E 493 4+
I
"
E+F4+2:6 2%96
+44
) < %+3
; )
I (E3 - E4.E++ / E+3 -.
(
E+ <%+3
!"# -.
> !"#
-.
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) PEM
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) DSR
< %+3I ;
"
E+F4+2:6 :E493
!"# ( E 3 - E4 .E ++ / E+3 -.
(
E+ .' 3" 53
%+ 2"6
;
!"# %3 -. %+
%3 < %+3
!"# -.
%+I ;
"
E
+F4+2:6 : E 493 !"# ( E3 - E4. E++
/ E+3 .' 3" 53
!"# -+ -3
K6 0L K6 0L
K4@@3A44++@ L K4@@5 44++@L
" K44++A44+%@L K44+3 44+%@L
4 *&
9 6 $ $
$
' . !"#
%3I ;
"
E
+F4+2:6 : E 493 -.
(
E+ .' 3 "53
-. - + - 3
K6 0L K6 0L
K4A455444+5L K4@%4 444+5L
" K44+A@4445L K4435 4445L
'
/ / E. / E.F+
# /
!"# /:
) .
/ / E.F+
"
53 $ ; /
; 4343 $//
/ /
/
"/ . .
< %+5 : . .
; : .
<<%+5 :
<%+5
.
$ *& (
6 $ $ $
< %+% -. !"# :
!"# . E++
/ E .F+ !"# -.
/
0 0.5
1 1.5
2 2.5
3
2 4 6 8 10 12 14 16
0 0.5 1 1.5 2 2.5 3 3.5
ω L
V 1
<%+5I "/ . . :
.
"
E +F4+2:6 !"# ( E 3 - E 4 / E .F+
.' 3 "53
'+ 6
E
4 +
E
E
+ 4
2%@6
. ) '
$ EK
L 2%A6
$EK49 +944 444+4 4444L 2%946
< %+ !"#
-. ;
./ 2%96 : E493 !"#
(E3 - E4 .E++ / E+3 !"#
0 0.5 1 1.5 2 2.5 3 0
0.5 1 1.5 2 2.5 3
V 1
ω PEM DSR
< %+%I "/ . . :
"
E +F4+2:6 .' 3 " 53 -.
(
E + !"# ( E 3 - E 4 . E ++
/ E+3 -. !"# *
-. )
!"#
%5 %% !"# -.
0 20 40 60 80 100
−0.8
−0.75 a 21
Simulation number DSR PEM True
0 20 40 60 80 100
1.68 1.7 1.72 1.74 1.76 1.78 1.8 1.82
a 22
Simulation number DSR PEM True
0 20 40 60 80 100
−1.8
−1.2
−0.6 x 10 −3
b 11
Simulation number DSR PEM True
0 20 40 60 80 100
−1
−0.8
−0.6
−0.4
−0.2 0 x 10 −3
b 21
Simulation number DSR PEM True
<%+I ./ 2%A6
; .' 3"53
"
E +F4+ 2:6
: E 493 !"# ( E 3 - E 4 . E ++ / E +3 -.
(
E+ -.
!"# *
%5I ;
.' 3 " 53
"
E+F4+2:6 :E493 !"#
( E3 - E4 .E++ / E+3
!"#
49@4 +9%+3 444+5 44449
" 4443 444 4444+ 4444+
%%I ;
.' 3 " 53
"
E+F4+2:6 : E493 -.
-.
49@A +9%5@ 444+4 4444
" 44+9 44+@ 4444+ 4444+
' ( )* !
$+
: -#>" 2- # > " 6 0:-;;-<=
>2
' E4+
2 *
6 : /
/ *
/
6
Ê
!
+#
6
E 2#+6F4+2#
9"'
K4 L6 2%9+6
+44
) < %+
; )
4@9 !"#
I ( E 3 - E 4 . E ++ / E +3 -.
(
E+ <%+ -#>"
!"# -.
% !"#
% -.
4 *&
9 6 (& $ $
/ !"#
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) PEM
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) DSR
< %+I ; .' 3 " 53
-#>" 6
E
2#+6F4+2#
9"'
K4 L6 E43 !"# ( E3
- E 4 . E ++ / E +3 -. (
E +
%I ;
.' 3 " 53 -#>"
6
E2#+6F4+2#
9"'
K4 L6
E43!"#( E3- E4.E++ / E+3
!"# -+ -3
K6 0L K6 0L
K4@@9@4455@L K4@9+4455@L
" K44+4544335L K4434+44335L
).
/ / E.F+
"53 7
;
-#>" ; AA ( / :
AA ( P 44+ 4AA "/ . .
/
/
< %+9 . .
;
.
< <%+9
. N N
) . !"#
%I ;
.' 3 " 53 -#>"
6
E 2#+6F4+2#
9"'
K4 L6
E43 -. (
E+
-. - + - 3
K6 0L K6 0L
K4@A4544+9+L K4@A4 44+9+L
" K44+5%443+4L K4435 443+4L
$ *& (
6 (& $$
<%+@ -. !"# -#>"
!"# . E++ / E.F+ -.
(
E+ / !"#
-. (
/
'+ 6
< ./ 2%A6
./ 2%946 <%+A
!"# -.
;
./ 2%9+6 E 43 !"# ( E3 - E 4
. E ++ / E +3
%+ %+ !"# -.
%9I ;
.' 3 " 53 -#>"
6
E 2#+6F4+2#
9"'
K4 L6 E 43
!"# (E3- E4.E+3 / E+5
!"#
499@% +93A 444+4 4444
" 44+9 44++ 4444+ 4444+
0.2 0 0.6 0.4
1 0.8
0 5 10 15 0
0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
L
B V 1
< %+9I "/ . .
. -#>"
6
E 2#+6F4+2#
9"'
K4 L6 .' 3 " 53
!"# ( E3 - E4 / E.F+
%@I ;
.' 3 " 53 -#>"
6
E 2#+6F4+2#
9"'
K4 L6 E 43
-. (
E+
-.
493 +9%A3 444+4 4444
" 44+%3 44+% 4444+ 4444+
$ !"#
-. (
!"#
'
'" (
!"#
'
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0
1 2 3 4 5 6 x 10 −3
V 1
B DSR
PEM
< %+@I "/ . .
-#>" 6
E
2#+6F4+ 2#
9"'
K4 L6 .' 3 " 53
-. (
E+ !"# ( E3
- E 4 . E ++ / E +3 -. !"#
*
0 20 40 60 80 100
−0.8
−0.75 a 21
Simulation number DSR PEM True
0 20 40 60 80 100
1.68 1.7 1.72 1.74 1.76 1.78 1.8 1.82
a 22
Simulation number DSR PEM True
0 20 40 60 80 100
−1.8
−1.2
−0.6 x 10 −3
b 11
Simulation number DSR PEM True
0 20 40 60 80 100
−1
−0.8
−0.6
−0.4
−0.2 0 x 10 −3
b 21
Simulation number DSR PEM True
< %+AI ./ 2%A6
; .' 3 " 53
-#>" 6
E 2#+6 F4+
2#
9"'
K4 L6 E43 !"# (E3- E4.E++
/ E+3 -. (
E+ -.
!"# *
"53
$
"
E+ 2%936
*
> * B
' /
!; !;
./ 23+@6 E+
2+6E44+4
$+
: / :E493 3
3+
E32:6 2%956
+44
)
< %34 ;
!"# I ( E 3 - E 4 . E ++ / E +3
-. (
E + < %34
!"#
-.
> !"#
-. "%+
!"#
-. 8
(
!"# # -.
%A 2"6
;
!"# %+4 -. %A
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) PEM
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) DSR
<%34I ; .' 3"53
"
E+
3+
E32:6 :E 493 !"# ( E3 - E4 . E++
/ E+3 -. (
E+
!"# -.
%AI ; .'
3"53 "
E+
3+
E 32:6 : E493
!"# (E3 - E4 .E++ / E+3
!"# -+ -3
K6 0L K6 0L
K4@@3944+3+ L K4@@ 44+3+L
" K44++A44++L K44+3 44++L
4 *&
9 6 $ $
$
7*"%+/
!"# / :
) .
/ / E.F+ ;
43 43 $/ /
"/ . . /
< %3+ : . .
; : .
%+4I ;
.' 3" 53 "
E+
3+
E32:6: E493
-. (
E+
-. - + - 3
K6 0L K6 0L
K4@@@@44+99L K4@5A4 44+99L
" K443+4443%@L K4+334 443%@L
0 <%3+ <%+5
: <
%3+ <%+5
.
$ *& (
6 $ $ $
<%33-. !"# :
!"# .E++ / E.F+
!"# -. /
&* <%+%
'+ 6
./ 2%A6 ./ 2%946
<%35 !"#
-. ;
./ 2%936
./ 2%956 : E 493 !"#
(E3- E4.E++ / E+3
7* < %+
!"# (
-. )
!"#
0 0.5
1 1.5
2 2.5
3
0 5
10 15
0 0.5 1 1.5 2 2.5 3 3.5
ω L
V 1
<%3+I "/ . . :
. "
E+
3+
E32:6 .' 3" 53
!"# ( E3 - E4 / E.F+
%++ %+3!"# -.
%++I ;
.' 3 "53 "
E+
3+
E32:6:E493
!"# ( E3 - E4 .E++ / E+3
!"#
49@4 +9%+3 444++ 4444
" 4443 444 4444+ 44444
' ( )*
!"# )
-#>"
)
0 0.5 1 1.5 2 2.5 3 0
0.5 1 1.5 2 2.5 3
V 1
ω DSR
PEM
<%33I "/ . . :
"
E +
3+
E 32:6 .' 3 " 53
-. (
E+ !"# ( E3
- E 4 . E ++ / E +3 -. !"#
*
0 20 40 60 80 100
−0.8
−0.75 a 21
Simulation number DSR PEM True
0 20 40 60 80 100
1.68 1.7 1.72 1.74 1.76 1.78 1.8 1.82
a 22
Simulation number DSR PEM True
0 20 40 60 80 100
−1.8
−1.2
−0.6 x 10 −3
b 11
Simulation number DSR PEM True
0 20 40 60 80 100
−1
−0.8
−0.6
−0.4
−0.2 0 x 10 −3
b 21
Simulation number DSR PEM True
<%35I ./ 2%A6
; .' 3"53
"
E+ 3+
E
3 2:6 : E 493 !"# ( E 3 - E 4 . E ++
/ E+3 -. (
E+ -.
!"# *
%+3I ;
.' 3" 53 "
E+
3+
E32:6: E493
-.
-.
49%93 +939@ 444++ 4444
" 4+499 4+3@ 4444+ 4444+
$+
./ 2%936
3+
Ê
!
+ #
3+
E32#
9"'
K4 L6 2%9%6
+44
)
< %3% ;
!"# I ( E3 - E4. E++ / E+3 -.
(
E+ < %+
-#>" !"#
-.
!"#
<%34
%+5!"#
%+% -.
4 *&
9 6 (& $ $
7* " %+ /
!"# / -#>"
).
/
/ E.F+ ; AA(/
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) PEM
0.7 0.8 0.9 1 1.1
−0.1
−0.08
−0.06
−0.04
−0.02 0 0.02 0.04 0.06 0.08
0.1 λ(A) DSR
<%3%I ; .' 3"53
"
E + -#>"
3+
E 32#
9"'
K4 L6 E 43 !"#
( E 3 - E 4 . E ++ / E +3 -.
(
E+
%+5I ;
.' 3 " 53 "
E +
-#>" 3+
E 3
2#
9"'
K4 L6 E43 !"# (E3- E4.E++
/ E+3
!"# -+ -3
K6 0L K6 0L
K4@@445%AL K4@9%@445%AL
" K44+++4435AL K44+944435AL
/
< %3 : . .
; : .
<%3
.N N & -#>"
< %+9
).!"#
)
$ *& (
6 (& $ $
<%3 -. !"# -#>"
%+%I ;
.' 3 " 53 "
E +
-#>" 3+
E 3
2#
9"'
K4 L6 E43 -.
(
E+
-. - + - 3
K6 0L K6 0L
K4@A4944+9@L K4@@A 44+9@L
" K44+@443+4L K4433@ 443+4L
/ !"# -.
-#>" (
<%+5
< %3 >
/
'+ 6
./ 2%A6 ./ 2%946
<%39 !"#
-. ;
./ 2%936
./ 2%9%6 E43 !"# ( E3
- E4.E++ / E+3
%+ %+!"# -.
:
<%35
, -) - % ).
"%+
/
).
8 '
0.2 0 0.4
0.8 0.6 1
0 5
10 15
0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
L B
V 1
<%3I "/ . .
. "
E+ -#>"
3+
E32#
9"'
K4 L6 .'
3 "53 !"# ( E3 - E4 / E.F+
* '
.
; * " %+
) !"# ) 7)
) )
< !"# 37C ) .
"
/
"
' '
) .
/
' ( )
%+I ; .'
3"53 "
E+
-#>" 3+
E32#
9"'
K4 L6
E43 !"# (E3 - E4 .E+3 / E+5
!"#
49993 +9+5 444+4 4444
" 44+5% 44+5 4444+ 44444
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1
2 3 4 5 6 7 x 10 −3
V 1
B
< %3I "/ . .
"
E+ -#>"
3+
E32#
9"'
K4 L6.'
3"53 -. (
E+
!"# ( E3 - E4 .E++ / E+3 -.
!"# *
0 20 40 60 80 100
−0.8
−0.75 a 21
Simulation number DSR PEM True
0 20 40 60 80 100
1.68 1.7 1.72 1.74 1.76 1.78 1.8 1.82
a 22
Simulation number DSR PEM True
0 20 40 60 80 100
−1.8
−1.2
−0.6 x 10 −3
b 11
Simulation number DSR PEM True
0 20 40 60 80 100
−1
−0.8
−0.6
−0.4
−0.2 0 x 10 −3
b 21
Simulation number DSR PEM True
<%39I ./ 2%A6
; .' 3"53
"
E + -#>"
3+
E 3 2#
9"'
K4 L6 E 43 !"# ( E 3
- E 4 . E ++ / E +3 -. (
E +
-. !"#
*
%+I ; .'
3"53 "
E+
-#>" 3+
E32#
9"'
K4 L6
E43 -. (
E+
-.
49 +9%A 444+4 4444
" 44++9 44++A 4444+ 44444
' !"#
+ !
! ! , - ! -
)
!"# 37C )
) "53
) .
/
'
./
2%96 : E 493 ) !"#
) . / < %3@
"/ J .
).
/
E.
F+
./ 2%96: E493!"# 37C
(E3- E4.E++ /E+3
).
)
'
; ).
.
"
)
+ .% &
)
./
2%A6 ./ 2%946
2 4 6 8 10 12 14 0
0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045
V 2
L z
<%3@I "/ J.
).
"
E+F4+2:6 :E493 .' 3"
53 !"# 37C ( E3 - E4.E++ /E+3
< %3A !"#
!"# 37C +44 ;
./ 2%96 : E 493 !"#
( E 3 - E 4 . E ++ / E +3 !"# 37C
( E3 - E4.
"
E++ /
"
E+3.
E /
E9
8!"# 37C
)
) !"#
!"# 37C
%5 %+9!"# !"# 37C
!"# 37C
> '
' <%54
.
< %5+
"/ J .
). / E. F+ -#>"
0 20 40 60 80 100
−0.8
−0.75 a 21
Simulation number DSR DSR2LJ PEM True
0 20 40 60 80 100
1.68 1.7 1.72 1.74 1.76 1.78 1.8 1.82
a 22
Simulation number DSR DSR2LJ PEM True
0 20 40 60 80 100
−1.8
−1.2
−0.6 x 10 −3
b 11
Simulation number DSR DSR2LJ PEM True
0 20 40 60 80 100
−1
−0.8
−0.6
−0.4
−0.2 0 x 10 −3
b 21
Simulation number DSR DSR2LJ PEM True
<%3AI ./ 2%A6
; .' 3"53
"
E+F4+2:6
: E493 !"# ( E3 - E 4 .E++ / E+3
!"# 37C ( E3- E4.
"
E++/
"
E+3 .
E
/
E9 -. (
E+ -.
!"# * !"# 37C
0 20 40 60 80 100
−0.8
−0.75 a 21
Simulation number DSR DSR2LJ PEM True
0 20 40 60 80 100
1.68 1.7 1.72 1.74 1.76 1.78 1.8 1.82
a 22
Simulation number DSR DSR2LJ PEM True
0 20 40 60 80 100
−1.8
−1.2
−0.6 x 10 −3
b 11
Simulation number DSR DSR2LJ PEM True
0 20 40 60 80 100
−1
−0.8
−0.6
−0.4
−0.2 0 x 10 −3
b 21
Simulation number DSR DSR2LJ PEM True
<%54I ./ 2%A6
; .' 3"53
-#>" 6
E 2#+6F4+
2#
9"'
K4 L6 E 43 !"# ( E3
- E4 .E++ / E +3 !"# 37C ( E3 - E4
.
"
E++/
"
E+3 .
E /
E9 -.
(
E+ -. !"#
* !"# 37C
%+9I ;
.' 3 " 53
"
E+F4+2:6 :E493 !"# 37C
( E 3 - E 4 .
"
E ++ /
"
E +3 .
E /
E 9
!"#
49@4 +9%+3 444+4 4444
" 4443 444 4444+ 4444+
./ 2%9+6 E 43 !"# 37C ( E 3
- E4 .E++ /E+3
2 4 6 8 10 12 14
0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045
V 2
L z
<%5+I "/ J .
) .
-#>"
6
E 2#+6F4+2#
9"'
K4 L6
E 43 .' 3 " 53 !"# 37C ( E 3 - E 4
.E++ /E+3
; <%5+ <%3@
.
'
' )
Æ
!
;
2" !6
<
&
; %
) .
) *
) :
' ) 8
:
= *
* *
! # 23445 6 *
" !
( ( *
' " %%
" 5+ !"# -.
$%" ! ( E 3 ' . )
' * -. $%" !
( E + . !"# )
"/ . . / E.F+ - E4 ;
+44 <
I '
26
2+6 ,
26 2+6
/"( % ) %
:
/ ./ 23+%6
,
E
2+6
* - NN
./ 2+6
(
N/N-#>" <%5 55
!"# < %% 5%
-. !"# $%" ! *
./ 2+6 <+ -. $%" !
( E + !"# - E 4 . E %
/ E
< +
!"# -.
$%" ! -.
!"#
4 < 3
-. $%" ! ( E + !"#
- E4.E+5 / E+%
$
*
' < B
'
< 5 N /N -#>" "
*
0.6 0.8 1 1.2
−0.5 0 0.5
PEM
λ (A)
0.6 0.8 1 1.2
−0.5 0 0.5
DSR
0.6 0.8 1 1.2
−0.5 0 0.5
N4SID
1.5 2 2.5 3 3.5
−1 0 1
H d (1)
1.5 2 2.5 3 3.5
−1 0 1
1.5 2 2.5 3 3.5
−1 0 1
0 0.5 1
−0.5 0 0.5
λ (A−KD)
0 0.5 1
−0.5 0 0.5
0 0.5 1
−0.5 0 0.5
0.5 1 1.5 2
−1 0 1
H s (1)
0.5 1 1.5 2
−1 0 1
0.5 1 1.5 2
−1 0 1
<+I . / -#>""
55
* .' +" 5+ -.
$%" ! ( E + !"#
- E4 .E% / E
0.6 0.8 1 1.2
−0.5 0 0.5
PEM
λ (A)
0.6 0.8 1 1.2
−0.5 0 0.5
DSR
0.6 0.8 1 1.2
−0.5 0 0.5
N4SID
1.5 2 2.5 3 3.5
−1 0 1
H d (1)
1.5 2 2.5 3 3.5
−1 0 1
1.5 2 2.5 3 3.5
−1 0 1
0 0.5 1
−0.5 0 0.5
λ (A−KD)
0 0.5 1
−0.5 0 0.5
0 0.5 1
−0.5 0 0.5
0.5 1 1.5 2
−1 0 1
H s (1)
0.5 1 1.5 2
−1 0 1
0.5 1 1.5 2
−1 0 1
<3I . / -#>""
55
* .' +" 5+
4 -. $%" !
(E+ !"# - E4 .E+5 / E+%