Proceedings of the 61st Annual Conference of the Institute of Systems, Control and Information Engineers (ISCIE), Kyoto, May 23-25, 2017 The Visual Se

Similar documents
& Vol.5 No (Oct. 2015) TV 1,2,a) , Augmented TV TV AR Augmented Reality 3DCG TV Estimation of TV Screen Position and Ro

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

{ 8. { CHAPTER 8. Å (sampling time) x[k] =x(kå) u(ú) t t + Å (u[k]) x[k + 1] =A d x[k] +B d u[k] (8:) (diãerence equation) A d =e AÅ ; B d = Z Å 0 e A

28 Horizontal angle correction using straight line detection in an equirectangular image


analog-control-mod : 2007/2/4(8:44) 2 E8 P M () r e K P ( ) T I u K M T M K D E8.: DC PID K D E8. (E8.) P M () E8.2 K P D () ( T ) (E8.2) K M T M K, T

pdf

人文学部研究年報12号.indb

公益社団法人日本都市計画学会都市計画論文集 Vol.53 No 年 10 月 Journal of the City Planning Institute of Japan, Vol.53 No.3, October, 2018 A queueing model for goods d

ばらつき抑制のための確率最適制御

258 5) GPS 1 GPS 6) GPS DP 7) 8) 10) GPS GPS ) GPS Global Positioning System

Untitled

,, 2. Matlab Simulink 2018 PC Matlab Scilab 2


soturon.dvi

IPSJ SIG Technical Report Vol.2016-CE-137 No /12/ e β /α α β β / α A judgment method of difficulty of task for a learner using simple

5D1 SY0004/14/ SICE 1, 2 Dynamically Consistent Motion Design of Humanoid Robots even at the Limit of Kinematics Kenya TANAKA 1 and Tomo

<4D F736F F D B B83578B6594BB2D834A836F815B82D082C88C60202E646F63>

II No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2

1 Fig. 1 Extraction of motion,.,,, 4,,, 3., 1, 2. 2.,. CHLAC,. 2.1,. (256 ).,., CHLAC. CHLAC, HLAC. 2.3 (HLAC ) r,.,. HLAC. N. 2 HLAC Fig. 2

i

28 TCG SURF Card recognition using SURF in TCG play video

THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS TECHNICAL REPORT OF IEICE k

201711grade1ouyou.pdf

Public Pension and Immigration The Effects of Immigration on Welfare Inequality The immigration of unskilled workers has been analyzed by a considerab

通信容量制約を考慮したフィードバック制御 - 電子情報通信学会 情報理論研究会(IT) 若手研究者のための講演会

xia2.dvi

7,, i

I ( ) 1 de Broglie 1 (de Broglie) p λ k h Planck ( Js) p = h λ = k (1) h 2π : Dirac k B Boltzmann ( J/K) T U = 3 2 k BT

??

ver.1 / c /(13)

9 1. (Ti:Al 2 O 3 ) (DCM) (Cr:Al 2 O 3 ) (Cr:BeAl 2 O 4 ) Ĥ0 ψ n (r) ω n Schrödinger Ĥ 0 ψ n (r) = ω n ψ n (r), (1) ω i ψ (r, t) = [Ĥ0 + Ĥint (

2.2 (a) = 1, M = 9, p i 1 = p i = p i+1 = 0 (b) = 1, M = 9, p i 1 = 0, p i = 1, p i+1 = 1 1: M 2 M 2 w i [j] w i [j] = 1 j= w i w i = (w i [ ],, w i [

1 No.1 5 C 1 I III F 1 F 2 F 1 F 2 2 Φ 2 (t) = Φ 1 (t) Φ 1 (t t). = Φ 1(t) t = ( 1.5e 0.5t 2.4e 4t 2e 10t ) τ < 0 t > τ Φ 2 (t) < 0 lim t Φ 2 (t) = 0

SICE東北支部研究集会資料(2012年)

grad φ(p ) φ P grad φ(p ) p P p φ P p l t φ l t = 0 g (0) g (0) (31) grad φ(p ) p grad φ φ (P, φ(p )) xy (x, y) = (ξ(t), η(t)) ( )

ohp_06nov_tohoku.dvi

1 Kinect for Windows M = [X Y Z] T M = [X Y Z ] T f (u,v) w 3.2 [11] [7] u = f X +u Z 0 δ u (X,Y,Z ) (5) v = f Y Z +v 0 δ v (X,Y,Z ) (6) w = Z +

Input image Initialize variables Loop for period of oscillation Update height map Make shade image Change property of image Output image Change time L

第 55 回自動制御連合講演会 2012 年 11 月 17 日,18 日京都大学 1K403 ( ) Interpolation for the Gas Source Detection using the Parameter Estimation in a Sensor Network S. T

ac b 0 r = r a 0 b 0 y 0 cy 0 ac b 0 f(, y) = a + by + cy ac b = 0 1 ac b = 0 z = f(, y) f(, y) 1 a, b, c 0 a 0 f(, y) = a ( ( + b ) ) a y ac b + a y

Virtual Window System Virtual Window System Virtual Window System Virtual Window System Virtual Window System Virtual Window System Social Networking

画像工学特論

LAN LAN LAN LAN LAN LAN,, i

[1.1] r 1 =10e j(ωt+π/4), r 2 =5e j(ωt+π/3), r 3 =3e j(ωt+π/6) ~r = ~r 1 + ~r 2 + ~r 3 = re j(ωt+φ) =(10e π 4 j +5e π 3 j +3e π 6 j )e jωt

I, II 1, A = A 4 : 6 = max{ A, } A A 10 10%


D = [a, b] [c, d] D ij P ij (ξ ij, η ij ) f S(f,, {P ij }) S(f,, {P ij }) = = k m i=1 j=1 m n f(ξ ij, η ij )(x i x i 1 )(y j y j 1 ) = i=1 j

08-Note2-web

(3.6 ) (4.6 ) 2. [3], [6], [12] [7] [2], [5], [11] [14] [9] [8] [10] (1) Voodoo 3 : 3 Voodoo[1] 3 ( 3D ) (2) : Voodoo 3D (3) : 3D (Welc

Dirac 38 5 Dirac 4 4 γ µ p µ p µ + m 2 = ( p µ γ µ + m)(p ν γ ν + m) (5.1) γ = p µ p ν γ µ γ ν p µ γ µ m + mp ν γ ν + m 2 = 1 2 p µp ν {γ µ, γ ν } + m

3D UbiCode (Ubiquitous+Code) RFID ResBe (Remote entertainment space Behavior evaluation) 2 UbiCode Fig. 2 UbiCode 2. UbiCode 2. 1 UbiCode UbiCode 2. 2

(1.2) T D = 0 T = D = 30 kn 1.2 (1.4) 2F W = 0 F = W/2 = 300 kn/2 = 150 kn 1.3 (1.9) R = W 1 + W 2 = = 1100 N. (1.9) W 2 b W 1 a = 0


0401489‐工芸‐医用画像22‐1/12[論文]柳田

c 2006 Yoneda norimasa All rights reserved

Sobel Canny i


構造と連続体の力学基礎

鉄鋼協会プレゼン

( ) x y f(x, y) = ax

24 Depth scaling of binocular stereopsis by observer s own movements

²ÄÀÑʬΥ»¶ÈóÀþ·¿¥·¥å¥ì¡¼¥Ç¥£¥ó¥¬¡¼ÊýÄø¼°¤ÎÁ²¶á²òÀÏ Asymptotic analysis for the integrable discrete nonlinear Schrödinger equation

A Feasibility Study of Direct-Mapping-Type Parallel Processing Method to Solve Linear Equations in Load Flow Calculations Hiroaki Inayoshi, Non-member

IHI Robust Path Planning against Position Error for UGVs in Rough Terrain Yuki DOI, Yonghoon JI, Yusuke TAMURA(University of Tokyo), Yuki IKEDA, Atsus

MmUm+FopX m Mm+Mop F-Mm(Fop-Mopum)M m+mop MSuS+FX S M S+MOb Fs-Ms(Mobus-Fex)M s+mob Fig. 1 Particle model of single degree of freedom master/ slave sy

fiš„v8.dvi



知能と情報, Vol.30, No.5, pp

2007-Kanai-paper.dvi

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H

K E N Z OU

,. Black-Scholes u t t, x c u 0 t, x x u t t, x c u t, x x u t t, x + σ x u t, x + rx ut, x rux, t 0 x x,,.,. Step 3, 7,,, Step 6., Step 4,. Step 5,,.

B ver B

B HNS 7)8) HNS ( ( ) 7)8) (SOA) HNS HNS 4) HNS ( ) ( ) 1 TV power, channel, volume power true( ON) false( OFF) boolean channel volume int

) a + b = i + 6 b c = 6i j ) a = 0 b = c = 0 ) â = i + j 0 ˆb = 4) a b = b c = j + ) cos α = cos β = 6) a ˆb = b ĉ = 0 7) a b = 6i j b c = i + 6j + 8)

IPSJ SIG Technical Report Vol.2012-MUS-96 No /8/10 MIDI Modeling Performance Indeterminacies for Polyphonic Midi Score Following and

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

1 filename=mathformula tex 1 ax 2 + bx + c = 0, x = b ± b 2 4ac, (1.1) 2a x 1 + x 2 = b a, x 1x 2 = c a, (1.2) ax 2 + 2b x + c = 0, x = b ± b 2

2011de.dvi

December 28, 2018

kiyo5_1-masuzawa.indd

Z: Q: R: C: sin 6 5 ζ a, b

1 variation 1.1 imension unit L m M kg T s Q C QT 1 A = C s 1 MKSA F = ma N N = kg m s 1.1 J E = 1 mv W = F x J = kg m s 1 = N m 1.

paper.dvi

Title 最適年金の理論 Author(s) 藤井, 隆雄 ; 林, 史明 ; 入谷, 純 ; 小黒, 一正 Citation Issue Date Type Technical Report Text Version publisher URL

2. 2 P M A 2 F = mmg AP AP 2 AP (G > : ) AP/ AP A P P j M j F = n j=1 mm j G AP j AP j 2 AP j 3 P ψ(p) j ψ(p j ) j (P j j ) A F = n j=1 mgψ(p j ) j AP

25 D Effects of viewpoints of head mounted wearable 3D display on human task performance

.,,, [12].,, [13].,,.,, meal[10]., [11], SNS.,., [14].,,.,,.,,,.,,., Cami-log, , [15], A/D (Powerlab ; ), F- (F-150M, ), ( PC ).,, Chart5(ADIns

MPC MPC R p N p Z p p N (m, σ 2 ) m σ 2 floor( ), rem(v 1 v 2 ) v 1 v 2 r p e u[k] x[k] Σ x[k] Σ 2 L 0 Σ x[k + 1] = x[k] + u[k floor(l/h)] d[k]. Σ k x

,,,,., C Java,,.,,.,., ,,.,, i

IPSJ SIG Technical Report Vol.2012-CG-148 No /8/29 3DCG 1,a) On rigid body animation taking into account the 3D computer graphics came

三石貴志.indd

P2P P2P peer peer P2P peer P2P peer P2P i

( ) [1] [4] ( ) 2. [5] [6] Piano Tutor[7] [1], [2], [8], [9] Radiobaton[10] Two Finger Piano[11] Coloring-in Piano[12] ism[13] MIDI MIDI 1 Fig. 1 Syst

Note.tex 2008/09/19( )

II A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )

2 1 κ c(t) = (x(t), y(t)) ( ) det(c (t), c x (t)) = det (t) x (t) y (t) y = x (t)y (t) x (t)y (t), (t) c (t) = (x (t)) 2 + (y (t)) 2. c (t) =

1 Fourier Fourier Fourier Fourier Fourier Fourier Fourier Fourier Fourier analog digital Fourier Fourier Fourier Fourier Fourier Fourier Green Fourier

Transcription:

The Visual Servo Control of Drone in Consideration of Dead Time,, Junpei Shirai and Takashi Yamaguchi and Kiyotsugu Takaba Ritsumeikan University Abstract Recently, the use of drones has been expected in many applications such as exploration in dangerous areas, security surveillance, and logistics Accurate position control of the drones is required in these applications This work is concerned with the remote position tracking control of a quadrotor-type drone via visual feedback control The transient response or stability of the drone may be significantly deteriorated by the adverse effect of the time delays due to wireless communication and image processing In this work, we construct the LQI optimal servo control system with state augmentation to cope with the time delays The effectiveness of the proposed method is verified by an experiment 1, () 4,,,,, 1],,,,, PC,, LQI,, 8], 2 21, x, y, z, X, Y, Z 1 1: 2:,, θ, ϕ, ψ 2,,, R(Θ) = R(ψ)R(ϕ)R(θ) R( ) 17 第 61 回システム制御情報学会研究発表講演会 (217 年 5 月 23 日 25 日, 京都 )

22 2] m ξ + mge z = F ξ (1) η = J 1 (τ C(η, η) η) (2) F ξ Z F = E z u R(Θ) F ξ = R(Θ)F (3) 1 1: m ξ = (x, y, z) g g = 98m/s 2 ] E z = ] η = (ψ, ϕ, θ),, J = diag(i x, I y, I z ) τ=(τ ψ, τ ϕ, τ θ ) C u Z () 23 (η =, η =, u = mg) mẍ = mgϕ (4) I x ϕ = τϕ mÿ = mgθ (5) I y θ = τθ m z = u mg (6) I z ψ = τψ (7) I x, I y, I z (4) (7),, (y ) (5) 24 u τ (Parrot AR Drone2) PC η d = (ψ d, ϕ d, θ d ) u τ (8) PD I y θ(t) = b (θ(t) θ d (t L 1 )) a θ(t) mÿ(t) = mgθ K ẏ(t) (8) (5) 2 (8) 2 K WiFi L 1 WiFi PC L 2 25 (8) t (8) ẋ(t) = A c x(t) + B c u(t L 1 ) y(t) = C c x(t) θ(t) x(t) = θ(t) y(t), u(t) = θ d(t), y(t) = y(t), ẏ(t) 1 A c = b a 1, B c = b, g K v ] C c = 1 (9a) (9b) a = a /I y, b = b /I y, K v = K /m 1 1 L 1, L 2 t l 1, l 2 L 1 = l 1 t, L 2 = l 2 t 1 sine sweep 2

ts] u(k t + τ) = u(k t), τ, t), (k =, 1, 2, ) xk + 1] = Axk] + Buk l 1 ] yk] = Cxk] (1a) (1b) t A = e Ac t, B = e Ac t e Acτ dτb c, C c = C f(t) t = k t, k =, 1, 2, fk] := f(k t) 2 (i) (A, B) (C, ] A) A I B (ii) rank C 3,, 3, G(z) 31 3: 4 (y ) 4 y rk] (i) lim rk] yk] = k 4: (ii) (i) (ii) LQ LQI 32 l 1 l 2 3] yk l 2 ] yk l 2 + 1] x od k] := (11) yk 1] ] x od k], Xk] := xk] (1), Xk + 1] = A ob Xk] + B ob uk l 1 ] yk l 2 ] = C ob Xk] (12a) (12b) I p 1 A ob = I p, B ob =, C ob = C A B T

x id k] uk l 1 ] uk l 1 + 1] x id k] := uk 1] 33 LQI (i),(ii) LQI W k + 1] := W k] + Ek], Ek] := rk l 2 ] yk l 2 ] Xk] X m k] = W k], x id k] X m k + 1] = ÂX mk] + ˆBuk] + r 1] T yk l 2 ] = ĈX mk] A ob B ob C ob 1 1 Â =, ˆB = 1 ] Ĉ = 1 1 (13a) (13b) uk] = K X Xk] + K W W k] + K xid x id k] + u (14) K X, K W, K xid, u (13a), r 2(ii) (i) x, u x, ū ] x = ū A I C ] ] B r (15) W r X ū X m = W, X = r, x id = (16) x id ū x (14) k u = K X X KW W Kxid x d + ū (17) Xk] X m k] = X m k] X m = W k] x id k] Xk] = Xk] X, W k] = W k] W, x id k] = x id k] x id, ũk] = uk] ū X m k + 1] = Â X m k] + ˆBũk] (18) (14), (17) ũk] = K X m k], K = K X K W K xid ] (ii) J J = = (Ek] 2 + Q w W k] 2 + Rũk] 2) k= ( Xk] W k] k= x d k] T Xk] ) Q W k] + Rũk] 2 x d k] Q = diag(, C T C, Q w, ) 3 Q w R (19) K = (R + ˆB T P ˆB) 1 ˆBT P Â (2) P Riccati P = ÂT P Â ÂT P ˆB(R + ˆB T P ˆB) 1 ˆBT P Â + Q (21) 2, 3, W

, J = X T m]p X m ], W = P XW T X] + P W xid x id ] (22) P W W P XW, P W W Riccati P 34 2(i) (C ob, A ob ) yk l 2 ] ] yk l x1 k] 2 + 1] Xk] = = x 2 k] yk 1] xk] (12a) x 1 k + 1] = A 12 x 2 k] (23) x 2 k + 1] = A 22 x 2 k] + B 2 uk l 1 ] (24) I p ] I A11 A 12 p A ob = = A 21 A 22 Ip, C A ] B1 B ob = = B 2 B x 1 k] = yk l 2 ] x 2 k] x 2 k] 5] 43 (23), (24) (A z ) γ γ 1: γ 2 P L γ (, 1) γq L γq T L γ 2 P L A T 22P L A 22 A T 12Q L Q T L A 12 ] > P L Q L L = P 1 L Q L A z = A 22 LA 12 γ 4 4 yz- z = 1m] y 1m] On-Off x, z x = 2m] z = 1m] 15s], 3s] 2 1/3 333s] 3 LQ Q w = 1, R = 3, γ = 86 25% 1s] 2: AR Drone 2 (Parrot ) 515cm] 52cm] 472 g] 15rad] 35,rpm] HD 72p 3fps] zk + 1] = A z zk] + A z Lx 1 k] + B 2 uk l 1 ] ˆx 2 k] = zk] + Lx 1 k] A z := A 22 LA ] 12, Xk] ˆXk] x 1 k] = L, ˆx 2 k] 3: I x = 33, I y = 33, I z = 62kg m 2 ] L 1 = 133, L 2 = 566 s] a = 12, b = 627, K v = 28567

41 8 6 4 2 2 4 6 8 1 12 14 16 18 2 5: (), (l 1 = l 2 = ) 5 42 6 15 1 5-5 5 1 15 2 25 3 2 1-1 (a) (y ) -2 5 1 15 2 25 3 (b) 6: () 6(a) 6(b) 5,, LQI anti-windup 6] 1] T Pobkrut, T Eamsa-ard, T Kerdcharoen: Sensor drone for aerial odor mapping for agriculture and security services, Proc of 13th ECTI-CON, 216 2] LR Garcia Carrillo et al: Quad Rotorcraft Control: Vision-based Hovering and Navigation, Springer-Verlag, 213 3] F Liao, K Takaba, T Katayama: Design of an optimal preview servomechanism for discrete-time systems in a multirate setting, Dynamics of Continuous, Discrete and Impulsive Systems: Series B, vol 1, pp 727-744, 23 4] LV Santana, AS Brandao, M Sarcinelli-Filho: Outdoor waypoint navigation with the ARDrone quadrotor, Proc of Int Conf on Unmanned Aircraft Systems (ICUAS), pp 33-311, 215 5] 1999 6] Anti-Windup, 32, pp139-144, 23 7] :,, 21 8] :, 59, 4p, 215