Chapter 15 - Kỹ thuật thông tin vô tuyến

Kỹ thuật thông tin vô tuyến

Chapter 15 1. City has 10 macro-cells each cell has 100 users ∴ total number of users = 1000 Cells are of size 1 sqkm maximum distance traveled to traverse = √ 2km ∴ time = √ 2 30 = 169.7s In the new setup number of cells = 10 5 microcells total number of users = 1000 × 100 2 users time = √ 2×10 30×10 3 =1.69s ∴ number of users increases by 10000 and handoff time reduces by 1/100 2. See Fig 1 (0,0) (0,1) v u (1,1) (1,0) 60 Figure 1: Problem 2 D 2 = (j20) 2 + (i20) 2 − 2(j20)(i20) cos(2π/3) ⇒ D = 2a  i 2 + j 2 + ij = √ 3R  i 2 + j 2 + ij 3. diamond shaped cells, R= 100m D min = 600m D = 2KR K = D 2R = 600 2×100 = 3 N = K 2 = 9 (a) number of cells per cluster = N = 9 (b) number of channels per cell = total number/N = 450/9 = 50 4. (a) R=1km D=6km N = A cluster A cell = √ 3D 2 /2 3 √ 3R 2 /2 = 1 3 (D/R) 2 = 1 3 6 2 = 12 number of cells per cluster = N = 12 D SHADED CELLS USE SAME FREQUENCY Figure 2: Problem 3 (b) number of channels in each cell = 1200/12 = 100 (c)  i 2 + j 2 + ij = 2 √ 3 ⇒ i = 2, j = 2 5. R=10m D=60m γ I = 2 γ 0 = 4 M = 4 for diamond shaped cells SIR a = R −γ I MD −γ 0 = R −2 4D −4 = 32400 SIR b = R −4 4D −4 = 324 SIR c = R −2 4D −2 = 9 SIR a > SIR b > SIR c 6. γ = 2 BP SK P b = 10 −6 → P b = Q( √ 2γ b ) ⇒ γ b = SIR 0 = 4.7534 B = 50M Hz each user100KHz = B s SIR = 1 M  D R  γ M=6 for hexagonal cells a 1 = 0.167 a 2 = 3 N > 1 a 2  SIR 0 a 1  2/γ ⇒ N ≥ 9.4879 ∴ N = 10 C u = 50 7. G = 100 ξ = 1 λ = 1.5 With no sectorization SIR = 1 ξ 3G (N c − 1)(1 + λ) = 4.7534 N c = 26.2450 = 26 With sectorization, interference is reduced by a factor of 3 N c = 76.7349 = 76 8. SINR = G  N c −1 i=1 X i +N α = p(X i = 1) N ∼ G(0.247N c , 0.078N c ) P o ut = p(SIR < SIR 0 ) (a) P out = p  G  N c −1 i=1 X i +N < SIR 0  = p   N c −1 i=1 X i + N > G SIR 0  (b) X =  N c −1 i=1 X i then X ∼ Bin(α, N c − 1) p(x + N > G/SIR 0 ) =  N c −1 N=0 p(n + N > G/SIR 0 |x = n)p(x = n) p(x = n) =  N c − 1 n  α n (1 − α) N c −1−n p(x + N > G/SIR 0 ) = N c −1  N=0 p(N > G/SIR 0 − n|x = n)p(x = n) = N c −1  N=0 p  N − 0.247N c √ 0.078N c > G SIR 0 − n − 0.247N c √ 0.078N c |x = n  = N c −1  N=0 Q  G SIR 0 − n − 0.247N c √ 0.078N c  p(x = n) (c) N c = 35 α = 0.5 SIR 0 = 5 G = 150 p = 0.0973 MATLAB for i = 1:length(n) pn(i) = (factorial(Nc-1)./(factorial(n(i)).*factorial(Nc-1-n(i)))) . *alpha.^n(i)*(1-alpha).^(Nc-1-n(i)); end sump = 0; for i = 1:length(n) f = ((G/sir0)-n(i)-.247*Nc)/(sqrt(.078*Nc)); sump = sump + .5*erfc((f)/sqrt(2))*pn(i); end (d) If x can be approximated as Gaussian then x ∼ G((N c − 1)α, (N c − 1)α(1 − α)) x + N ∼ G(0.247N c + (N c − 1)α, 0.078N c + (N c − 1)α(1 − α)) p(x + N > G/SIR 0 ) = Q  G SIR 0 − (0.247N c + (N c − 1)α)  0.078N c + (N c − 1)α(1 − α)  (e) p= 0.0969 (very accurate approximation!) 9. define γ k = g k P k n k + ρ  k=j g kj p j k, j ∈ {1, . . . K} where, g k is channel power gain from user k to his base station n k is thermal noise power at user k’s base station ρ is interference reduction factor (ρ ∼ 1/G) g kj is channel power gain from j th interfering transmitter to user k’s base station p k is user k’s Tx power p j is user j’s Tx power define a matrix F such that F kj =  0 k = j γ  k g kj ρ g k k = j k, j ∈ {1, . . . K} u =  γ  1 n 1 g 1 , γ  2 n 2 g 2 , . . . , γ  K n K g K  If Perron Ferbinius eigenvalue of F is less than 1 , then a power control policy exists. The optimal power control policy is given to be P  = (I − F ) −1 u 10. Matlab D = 2:.01:10; R = 1; gamma = 2; Pdes = R^(-gamma); for i = 1:length(D) Pint = 6*(.2*(D(i)-R)^(-gamma)+.2*(D(i)-R/2)^(-gamma)+.2*(D(i))^(-gamma) . +.2*(D(i)+R/2)^(-gamma)+.2*(D(i)+R)^(-gamma)); Pintbest = 6*((D(i)+R)^(-gamma)); Pintworst = 6*((D(i)-R)^(-gamma)); ASE(i) = log(1+Pdes/Pint)/(pi*(.5*D(i))^2); ASEbest(i) = log(1+Pdes/Pintbest)/(pi*(.5*D(i))^2); ASEworst(i) = log(1+Pdes/Pintworst)/(pi*(.5*D(i))^2); end 2 3 4 5 6 7 8 9 10 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 D ASE(D) 5 Case ASE Best Case ASE Worst Case ASE Figure 3: Problem 10 11. P t = 5W B = 100KHz N 0 = 10 −16 W/Hz P r = P t K  d 0 d  3 d 0 = 1, K = 100 (a) D=2R 2 users share the band available Each user gets 50KHz BASESTATION USER R=1Km Figure 4: Problem 11a (b) P b = 1 4γ b ⇒ 10 −3 = 1 4γ b ⇒ γ b = 250 If D(n) = 2nR, number of users that share band = 2(2(n-1)+1) ∴ each user gets 100KHz 2(2n−1) = B u (n) interference is only from first tier SIR(n) = P t K  d 0 R  3 N 0 2 B u (n) + 2  P t K  d 0 √ R 2 +D(n) 2  3  > 25 using Matlab , n = 4, SIR = 261.9253, D = 8R Matlab Pt = 5; R = 1000; sigma_2 = 1e-16; n = 1; D = 2*n*R; Bu = (100/(2*(2*n-1)))*1e3; K = 100; d0 = 1; Pdes = Pt*K*(d0/R)^3; Pint = 2*(Pt*K*(d0/sqrt(R^2+D^2))^3); Npower = sigma_2*Bu; sir = Pdes/(Npower+Pint); while sir < 250 n = n+1; D = 2*n*R; Bu = (100/(2*(2*n-1)))*1e3; K = 100; d0 = 1; Pdes = Pt*K*(d0/R)^3; Pint = 2*(Pt*K*(d0/sqrt(R^2+D^2))^3); Npower = sigma_2*Bu; sir = Pdes/(Npower+Pint); end (c) ASE = (R 1 +R 2 )/B 2km×2km R 1 = R 2 = B u (1) log(1 + SIR(1)) B u (1) = 50KHz, SIR(1) = 5.5899 ASE = 0.6801bps/Hz/km 2 12. B=100KHz N 0 = 10 −9 W/Hz K = 10 P = 10mW per user (a) 0 ≤ α ≤ 1 α is channel gain between cells. See Matlab If α is large, interference can be decoded and subtracted easily so capacity grows with α as high SNR’s (beyond an α value) . For low SNR values (α less than a value) c decreases with increase in α as interference is increased which cannot be easily decoded due to low SNR. MATLAB CODE: B = 100e3; sigma_2 = 1e-9; P = 10e-3; K = 10; ss = .001; alpha = 0:.01:1; theta = 0:ss:1; for i = 1:length(alpha) capvec = log2(1+(K*P*(1+2*alpha(i)*cos(2*pi*theta)).^2)/(sigma_2*B)); C(i) = (1/K)*sum(capvec)*ss; end 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.82 0.84 0.86 0.88 0.9 0.92 0.94 0.96 0.98 1 1.02 α C(α)/B B = 100 KHz Figure 5: Problem 12a (b) C(K)↓ as K↑ because as the number of mobile per cell increases system resources get shared more and so per user capacity C(K) has to fall. MATLAB CODE: BB = 100e3; sigma_2 = 1e-9; P = 10e-3; K = 1:.1:30; ss = .001; alpha = .5; theta = 0:ss:1; for i = 1:length(K) capvec = log2(1+(K(i)*P*(1+2*alpha*cos(2*pi*theta)).^2)/(sigma_2*B)); C(i) = (1/K(i))*sum(capvec)*ss; end 0 5 10 15 20 25 30 0 1 2 3 4 5 6 K C(K)/B B = 100 KHz Figure 6: Problem 12b (c) as transmit power P ↑, capacity C ↑ but gets saturated after a while as the system becomes interference limited. MATLAB CODE: B = 100e3; sigma_2 = 1e-9; P = [0:.1:100]*1e-3; K = 10; ss = .001; alpha = .5; theta = 0:ss:1; for i = 1:length(P) capvec = log2(1+(K*P(i)*(1+2*alpha*cos(2*pi*theta)).^2)/(sigma_2*B)); C(i) = (1/K)*sum(capvec)*ss; end 0 10 20 30 40 50 60 70 80 90 100 0 0.2 0.4 0.6 0.8 1 1.2 1.4 P (in mW) C(P)/B B = 100 KHz Figure 7: Problem 12c . Chapter 15 1. City has 10 macro-cells each cell has 100 users ∴ total number of. SIR 0 − n − 0.247N c √ 0.078N c  p(x = n) (c) N c = 35 α = 0.5 SIR 0 = 5 G = 150 p = 0.0973 MATLAB for i = 1:length(n) pn(i) = (factorial(Nc-1)./(factorial(n(i)).*factorial(Nc-1-n(i)))) .