Project: IEEE P802.15 Working Group for Wireess Persona Area Networks N (WPANs( WPANs) Tite: [MATLAB Simuation Program for TSV-channe mode] Date Submitted: [September 18, 2006] Source: [Hiroshi Harada, Ryuhei Funada, Hirokazu Sawada, Chang-soon Choi, Yozo Shoji, Shuzo Kato] Company [NICT] Address[3-4 Hikari-no-oka, Yokosuka-shi, Kanagawa 239-0847, Japan] Voice:[+81-46-847-5074] FAX:[+81-46-847-5440] E-Mai:[harada@nict.go.jp, funada@nict.go.jp, sawahiro@nict.go.jp, shoji@nict.go.jp, cschoi@nict.go.jp,shu.kato@nict.go.jp] Re: [] Abstract: [Proposing a MATLAB Simuation Program for TSV-channe mode] Purpose: [To be considered in 15.3c transmission performance by computer simuation] Notice: Reease: The contributors acknowedge and accept that this contribution becomes the property of IEEE and may be made pubicy avaiabe by P802.15. Side 1
MATLAB Simuation Program for TSV-channe mode Hiroshi Harada, Ryuhei Funada, Hirokazu Sawada,Chang-Soon Choi, Yozo Shoji, Shuzo Kato (NICT) Side 2
Summary of this document Finished to prepare MATLAB simuation program for TSV-channe mode Expain the fowchart of the MATLAB mode Show comparison of experimenta and simuated resuts Summarize avaiabe LOS / NLOS channe modes by the MATLAB-based TSV channe mode Side 3
β = α September 2006 2 μ PL D CIR: Definition of fina TSV mode L 1 M 1 () t = β δ () t + α, m δ ( t T τ, m ) δ ( ϕ Ψ ψ m ) h, = 0 m= 0 (Compex impuse response) τ γ k 2π 2h1h Γ0 exp j λ f D 2 D 2 Gt1Gr 1 + Gt 2Gr2 [ 1 δ ( m) ] G (, Ψ + ψ ), α Uniform [ 0,2π ) T Γ, m =Ω 0, 0e e m r, m, m Two-path response Path number of G ti and G ri (1: direct, 2 : refrect) Arriva rate: Poisson process p p ( T T ) = Λ exp[ Λ( T T )], 1 1 [ ( )], m > 0 ( τ τ, ( m 1) ) = λ exp λ τ τ, ( m 1) > 0 PL: Path oss of the first impuse response t: time[ns] δ( ): Deta function = custer number, m= ray number in -th custer, L = tota number of custers; M = tota number of rays in the -th custer; T = arriva time of the first ray of the -th custer; τ,m = deay of the m-th ray within the -th custer reative to the firs path arriva time, T ; Ω 0 = Average power of the first ray of the first custer Ψ Uniform[0,2π); arriva ange of the first ray within the -th custer ψ,m = arriva ange of the m-th ray within the -th custer reative to the first path arriva ange, Ψ 1 2 Two-path parameters (4) S-V parameters (7) D h h Uniform : Distance between Tx and Rx Uniform : Height of Uniform : Height of μ Average of distance between Tx and Rx Γ D 0 : Refection coefficient Γ Γ 0 0 Tx Rx 1: LOS Desktop environment (incident ange π 2) 0: Other LOS environment Γ : custer decay factor 1/ Λ : custer γ : ray decay factor 1/ λ : ray σ : custer ognorma standard deviation 1 σ : ray ognorma standard deviation 2 σ φ arriva rate arriva rate : Ange spread of ray within custer (Lapace distribution) K Antenna parameters (2) Gt Gr = L 1 M 1 ( θ, φ) : Antenna gain of Tx ( θ, ι) : Antenna gain of Rx = 0 m= 0 Rician factor (2) k :Sma Rician effect in each custer α 2, m δ ( t T τ ) δ ( ϕ Ψ ψ ) G ( 0, Ψ + ψ ), m 2 β, m r, m Side 4
Impuse response LOS penetration wave component Reative Ampitude β Ω 0 Rician factor (ΔK) Sma Rician factor (Δk) S-V mode response Γ,Λ,γ,λ Time of Arriva This response can be aso obtained in TSV mode by setting Γ0 =0 Side 5
Exampes of parameters for TSV mode TSV Sma S-V mode oriented parameters Number Mode Rician of custer factor Parameter Ω 0 (D) [db] k (Δk) Γ [ns] 1/Λ [ns] γ [ns] 1/λ [ns] σ 1 custer σ 2 ray σ φ [deg] N Tx:60 Rx:60 3.46 D- 98.4 3.97 22.3 21.1 17.2 2.68 7.27 4.42 38.1 3 Side 6
Function cas Main sv_params_tg3c sv_mode_ct_tg3c sv_cnvrt_ct_tg3c Side 7
Fowchart (overview) Main A input cm_num, number of channe reaizations (N), and minimum time resoution ca function sv_params_tg3c and get parameter regarding antenna beam-width, S-V parameters, and direct-component-reated parameters if necessary Cacuate and pot out performance metrics store off channe reaizations ca function sv_mode_ct_tg3c, and generate N impuse responses, using antenna gain determined by AOA information which are generated by the Monte Caro method done ca function sv_cnvrt_ct_tg3c and generate discrete time output from continuous time input decimate impuse response if necessary A Side 8
Fowchart of sv_mode_ct_tg3c yes function sv_mode_ct_tg3c generate N impuse responses yes LOS? no direct LOS component generated by SV mode no cacuate first custer arriva time cacuate direct LOS component and set the component at time 0 first custer arrives at time 0 no Tr: arriva time of ray in the k-th custer L: Number of custers to be generated k 0 k k+1, Tr 0, and set k-th custer s AOA N channes generated? no k=<l yes yes done no Tr<10*gam yes set ray AOA and then obtain antenna gain from the AOA cacuate ray arriva time Tr determine ray ampitude and phase using antenna gain and AOA Tr Tr+ Tr Side 9
Comparison of experimenta and simuated resuts Reative power [db] -70-80 -90-100 -110 0 20 40 60 80 100 Time of arriva [ns] LOS component (a) Experimenta resut Antenna height Tx: 170 mm Rx: 150 mm Beam width: 60 deg Distance: 3m S-V custer Reative power [db] -70-75 -80-85 -90-95 -100-105 Beam width: 60 deg Assumed distance: 3m S-V custers -110 0 20 40 60 80 100 Time of arriva [ns] LOS component Simuation data is a snap-shot. (b) Simuation resut Average RMS deay spread Experimenta resuts 10.6[ns] Side 10 Simuated resuts 9.2 [ns]
Summary of avaiabe LOS / NLOS channe modes by MATLAB based TSV-channe mode Office Residentia Desktop Library LOS Avaiabe (NICT) Avaiabe (NICT) Avaiabe (NICT) Avaiabe (IMST/Inte) NLOS Avaiabe (NICTA) N/A N/A N/A Measurement and anaysis to get TSV parameters are finished by NICT. MATLAB program is now avaiabe by using anayzed parameters. Measurement is finished by NICT. Anaysis to get TSV parameters is aso possibe within a coupe of week. MATLAB program wi be prepared by using anayzed parameters. Side 11
Summary Finished to prepare MATLAB simuation program for TSV-channe mode Expained the fowchart of the program Showed comparison of experimenta and simuated resuts Performance is amost simiar to the experimenta one Summarized avaiabe LOS / NLOS channe modes by MATLAB based TSV-channe mode NLOS as we as LOS can be covered by the proposed MATLAB program. Side 12