Nghiên cứu hệ thống định vị toàn cầu GPS và ứng dụng trong nghành hàng không - Pdf 35

LỜI CẢM ƠN
Sau 5 năm học tập tại trường, được sự tận tình dạy dỗ của các thầy cô.
Em đã hoàn thành đồ án tốt nghiệp với đề tài “ Nghiên cứu hệ thống định vị toàn
cầu GPS và ứng dụng trong nghành hàng không .”
Em xin gởi lời cám ơn đến ban lãnh đạo nhà trường, quý thầy cô trong
bộ môn Công Nghệ Điên Tử đã truyền đạt những kiến thức quý báu làm cơ sở
cho em hoàn thành đồ án tốt nghiệp.
Em đặc biệt cám ơn cô giáo Nguyễn Thị Dung là người trực tiếp
hướng dẫn em thực hiện đề tài này.
Em cũng xin được cám ơn sâu sắc đến gia đình và bạn bè đã tạo cho em
những điều kiện thuận lợi nhất trong suốt quá trình học tập và thực hiện đề tài.

SVTH: Lê Quang Huy

0


LỜI CAM ĐOAN

Để hoàn thành đồ án tốt nghiệp đúng thời gian quy định và đáp ứng được
yêu cầu đề ra, em đã cố gắng tìm hiểu, học hỏi, tích lũy kiến thức đã học. Em có
tham khảo một số tài liệu đã nêu trong phần “Tài liệu tham khảo” nhưng không
sao chép nội dung từ bất kỳ đồ án nào khác.
Em xin cam đoan đồ án là công trình nghiên cứu của cá nhân nghiên cứu,
xây dựng dưới sự hướng dẫn của cô giáo Nguyễn Thị Dung. Nội dung lý thuyết
trong đồ án có sự tham khảo và sử dụng của một số tài liệu, thông tin được đăng
tải trên các tác phẩm, tạp chí và các trang web theo danh mục tài liệu của đồ án.
Em xin cam đoan những lời khai trên là đúng, mọi thông tin sai lệch em
xin hoàn toàn chịu trách nhiệm trước hội đồng.

Thái Nguyên, ngày tháng 6 năm 2012

1.9. Phương pháp đa truy nhập phân phối trước và đa truy nhập phân phối theo
yêu cầu. ......................................................................................................... 28
1.9.1. Đa truy nhập phân phối trước ............................................................ 28
1.9.2. Đa truy nhập phân phối theo yêu cầu ................................................. 28
1.11. Suy hao trong thông tin vệ tinh ............................................................. 30
1.11.1. Suy hao trong không gian tự do ....................................................... 30

2


1.11.2. Suy hao do tầng đối lưu ................................................................... 31
1.11.3. Suy hao do tầng điện ly .................................................................... 32
1.11.4. Suy hao do thời tiết ......................................................................... 32
1.11.5. Suy hao do đặt anten chưa đúng ...................................................... 33
1.11.6. Suy hao trong thiết bị phát và thu .................................................... 33
1.11.7. Suy hao do phân cực không đối xứng .............................................. 33
1.12. Tạp âm trong thông tin vệ tinh .............................................................. 34
1.12.1. Nhiệt tạp âm hệ thống ..................................................................... 34
1.13. Hiệu ứng Doppler ................................................................................. 41
1.14. Trễ truyền dẫn....................................................................................... 41
CHƯƠNG II: HỆ THỐNG ĐỊNH VỊ TOÀN CẦU GPS ................................... 42
2.1. Sự ra đời của hệ thống GPS ................................................................... 42
2.2. Cấu tạo của hệ thống GPS....................................................................... 48
2.2.1. Chùm vệ tinh ..................................................................................... 48
2.2.2. Hệ thống điều khiển mặt đất .............................................................. 48
2.2.3. Bộ phận người sử dụng ..................................................................... 49
2.3. Nguyên lý hoạt động của hệ thống .......................................................... 49
2.4. Điều chế và giải điều chế GPS ................................................................ 51
2.4.1. Điều chế tín hiệu GPS ....................................................................... 51
2.4.2. Giải điều chế GPS ............................................................................. 52

KẾT LUẬN ....................................................................................................... 71
TÀI LIỆU THAM KHẢO ................................................................................. 72

4


Mục lục
Hình 1.1. Ba dạng quỹ đạo cơ bản của vệ tinh .....................................................9
Hình 1.2. Vệ tinh quỹ đạo địa tĩnh .................................................................... 10
Hình 1.3. Liên lạc giữa hai trạm mặt đất qua vệ tinh .......................................... 12
Hình 1.4. vị trí các vật thể được xác định qua 4 phép đo ................................... 14
Hình 1.5. định vị điểm bằng vệ tinh ................................................................... 15
Hình 1.6. Miêu tả 3 phép đo bằng 3 vệ tinh ...................................................... 16
Hình 1.7. Miêu tả xác định ví trí ....................................................................... 16
Hình 1.8: Mã giả ngẫu nhiên PRC .................................................................... 17
Hình 1.9. Phương pháp đo giả cự ly................................................................... 19
Hình 1.11. Các sai số do lỗi ............................................................................... 20
hình 1.12. sai số do tầng điện ly......................................................................... 22
Hình 1.13. Nhiễu đa đường................................................................................ 23
Hình 1.14. Truy nhập theo tần số ....................................................................... 24
Hình 1.15. Đa truy nhập phân chia theo thời gian .............................................. 25
Hình 1.16. Đồ thị biểu diễn suy hao do mưa (nhỏ) và do tầng điện ly theo tần số .. 29

Hình 1.17. Sai lệch do đặt anten chưa đúng ....................................................... 33
Hình 1.18. Suy hao trong thiết bị phát và thu ..................................................... 33
Hình 1.19. Các nguồn tạp âm ảnh hưởng đến thông tin vệ tinh .......................... 34
Hình 1.20. Can nhiễu giữa viba và trạm mặt đất và vệ tinh ............................... 38
Hình 1.21. Can nhiễu giữa các hệ thống thông tin vệ tinh .................................. 39
Hình 1.22. Đặc tính vào ra của TWT ................................................................. 40
Hình 1.24. Mức lùi đầu vào và lùi đầu ra ........................................................... 40

lượng và chất lượng, đòi hỏi các dịch vụ của ngành Viễn Thông càng mở rộng.
Trong những năm gần đây thông tin vệ tinh trên thế giới đã có những bước tiến
vượt bậc đáp ứng nhu cầu đời sống, đưa con người nhanh chóng tiếp cận với các
tiến bộ khoa học kỹ thuật.
Sự ra đời của nhiều loại phương tiện tiên tiến như máy bay, tàu vũ trụ đòi
hỏi một kỹ thuật mà các hệ thống cũ không thể đáp ứng được đó là định vị trong
không gian ba chiều, đứng trước sự đòi hỏi đó chính phủ Mỹ đã tài trợ một
chương trình nghiên cứu hệ thống định vị trong vũ trụ.
Được sự hướng dẫn của Cô giáo Nguyễn Thị Dung em đã chọn đề tài
“Nghiên cứu hệ thống định vị toàn cầu GPS và ứng dụng trong hàng không” cho
đồ án tốt nghiệp của mình.
Nội dung của đồ án gồm 3 chương như sau:
Chương 1: Tổng quan về hệ thống thông tin vệ tinh
Chương 2: Hệ thống định vị toàn cầu GPS
Chương 3: Ứng dụng hệ thống GPS trong ngành hàng không
Chương 1: Cung cấp các khái niệm cơ bản về TTVT, về định vị, cơ sở lý
thuyết của phép định vị bằng vệ tinh, phương trình đo mã, thiết lập mô hình toán
học của phép đo cũng như các nguồn gây sai số của phép đo và tìm hiểu về các
phương pháp và những ảnh hưởng của tạp âm và suy hao trong thông tin vệ tinh.
Chương 2: Tìm hiểu về ưu điểm, chức năng, cấu tạo, cấu trúc-đặc tính tín
hiệu, định dạng dữ liệu điện văn GPS và cũng như thông tin dẫn đường của GPS.
Chương 3: Sẽ trình bày ứng dụng của hệ thống GPS vào hàng không,
đồng thời xét hai hệ thống tăng cường GBAS và SBAS.

Sinh viên thực hiện
Lê Quang Huy

6



Quỹ đạo thấp

D/C

Down Converter

Bộ hạ tần

U/C

Up Converter

Bộ nâng tần

HPA

High Power Amplifier

Bộ khuếch đại công suất cao

LNA

Low Noise Amplifier

Bộ khuếch đại tạp âm thấp

PRN

Pseudo Random Code


Distance Measuring Equipment

Đài đo cự ly

Instrument Loading Sysem

Là hệ thống thiết bị mặt đất

S

P

ILS
VOR

SBAS

GBAS
GMS

Very High Frequency Omi
Range

Đài phát sóng dài phát vô
hướng

Đài vô tuyến phương vị mặt đất

Satellite Based Augmentation



Hệ thống mặt đất

AS

Aircraft System

Hệ thống trên máy bay

CNMP

Multipath error confidence
bounds

Gới hạn nhiễu đa đường

TROP

Troposphere delay

UDRE

User Differential Range Error

GIVE

Grid Ionosphere Vertical Error

Lỗi lưới dọc ion



Hệ thống trắc địa học toàn cầu
1984

UTC

Universal Time Coordinated

Là một hệ thời gian

TLM

Telemetry

Là từ đầu tiên của khung dữ liệu

HOW

Hand-over

Là từ thứ 3 của khung dữ liệu

International Civil Aviation

Tổ chức hàng không dân dụng

Organization

quốc tế


TDMA

Time Division Multiplex Access

8

Đa truy nhập phân chia theo
thời gian


CHƯƠNG I: TỔNG QUAN VỀ HỆ THỐNG THÔNG TIN VỆ TINH
1.1. Khái quát về định vị
Thật vậy, từ xa xưa con người đã biết vận dụng nhiều phương pháp như:
- Định vị cổ điển như quan sát theo dõi các ngọn núi, tòa nhà cao...
- Định vị quan sát các chòm sao và hành tinh trên vũ trụ: như sao Bắc đẩu
để xác định vị trí của mình tuy là không chính xác lắm.
Ngày nay, nhờ vào sự tiến bộ của khoa học kỹ thuật và các phương pháp
định vị hiện đại như:
- Định vị quán tính
- Định vị vô tuyến mặt đất
- Định vị vô tuyến không gian người ta đã tính toán và đo được các thông
số của quá trình chuyển động (vị trí, vận tốc, thời gian,…) chính xác hơn.
1.2. Nguyên lý của hệ thống thông tin vệ tinh
Sau khi được phóng vào vũ trụ, vệ tinh trở thành trạm thông tin ngoài trái
đất. Nó có nhiệm vụ thu tín hiệu dưới dạng sóng vô tuyến từ một trạm ở trái đất,
khuếch đại rồi phát trở về trái đất cho một trạm khác.
Có hai quy luật chi phối quỹ đạo của các vệ tinh bay xung quanh quả đất là:
- Mặt phẳng quỹ đạo bay của vệ tinh phải cắt ngang tâm Trái đất.
- Qủa đất phải là trung tâm của bất kỳ quỹ đạo nào của vệ tinh.


Hình 1.2. Vệ tinh quỹ đạo địa tĩnh

10


1.2.3.1. Quỹ đạo địa tĩnh GEO (Geosychronous Earth Orbit)
Vệ tinh địa tĩnh là vệ tinh được phóng lên quỹ đạo tròn ở độ cao khoảng
36.000km so với đường xích đạo, vệ tinh loại này bay xung quanh quả đất một
vòng mất 24h. Do chu kỳ bay của vệ tinh bằng chu kỳ quay của trái đất xung quanh
trục của nó theo hướng đông cùng với hướng quay của trái đất, bởi vậy vệ tinh
dường như đứng yên khi quan sát từ mặt đất, do đó nó được gọi là vệ tinh địa tĩnh.
Bởi vì một vệ tinh địa tĩnh có thể đảm bảo thông tin ổn định liên tục nên
có nhiều ưu điểm hơn vệ tinh quỹ đạo thấp dùng làm vệ tinh thông tin.
Nếu ba vệ tinh địa tĩnh được đặt ở cách đều nhau bên trên xích đạo thì có thể
thiết lập thông tin liên kết giữa các vùng trên trái đất bằng cách chuyển tiếp qua một
hoặc hai vệ tinh. Điều này cho phép xây dựng một mạng thông tin trên toàn thế giới.
1.2.3.2. Quỹ đạo thấp LEO (Low Earth Orbit)
Độ cao điển hình của dạng quỹ đạo này là 160 đến 480 km, nó có chu kỳ
90 phút. Thời gian quan sát thấy vệ tinh khoảng dưới 30 phút. Việc bố trí các vệ
tinh LEO gần nhau có thuận lợi là thời gian để dữ liệu phát đi đến vệ tinh và đi
về là rất ngắn. Do khả năng thực hiện nhanh của nó, tác dụng tiếp sức tương hỗ
toàn cầu giữa các mạng và loại hình hội thoại vô tuyến truyền hình sẽ có hiệu quả
và hấp dẫn hơn. Nhưng hệ thống LEO đòi hỏi phải có khoảng 60 vệ tinh loại này
mới bao trùm hết bề mặt địa cầu
1.2.3.3. Quỹ đạo trung bình MEO (Medium Earth Orbit)
Vệ tinh MEO ở độ cao từ 10.000km đến 20.000 km, chu kỳ của quỹ đạo là 5
đến 12 giờ, thời gian quan sát vệ tinh từ 2 đến 4 giờ. Ứng dụng cho thông tin di động
hay thông tin radio. Hệ thống MEO cần khoảng 12 vệ tinh để phủ sóng toàn cầu.
1.3. Đặc điểm của thông tin vệ tinh
Trong thời đại hiện nay, thông tin vệ tinh được phát triển và phổ biến

6GHz(14GHz)

Điều
chế

Nâng
tần

Khuếch đại
tạp âm tháp

Khuếch đại
công suất

Hạ
tần

Giải điều
chế

Hình 1.3. Liên lạc giữa hai trạm mặt đất qua vệ
Trong đó vệ tinh
tinh đóng vai trò lặp lại tín hiệu truyền giữa các trạm mặt đất,
thực chất kỹ thuật thông tin vệ tinh là kỹ thuật truyền dẫn mà trong đó môi
trường truyền dẫn là không gian vũ trụ với khoảng cách đường truyền khá dài.

12


Tại đây ta cũng gặp lại một số vấn đề đối với một bài toán truyền dẫn, đó là các


3,400  7,075

8,82  4,41

X

7,025  8,425

4,41  3,56

Ku

10,90  18,10

2,75  1,66

Ka

17,70  36,00

1,95  0,83

Bảng 1.1. Tần số sử dụng trong thông tin vệ tinh

13


Hiện nay, băng C và băng Ku được sử dụng phổ biến nhất, băng C (4/6
GHz) nằm ở khoảng giữa cửa sổ tần số, suy hao ít do mưa, trước đây được dùng

chính xác vị trí của vệ tinh rj(t) tại một thời điểm nào đó là rất khó khăn. Nhiệm
vụ dự đoán quỹ đạo hay lịch thiên văn (ephemeris) của vệ tinh đòi hỏi phải có
kiến thức đặc biệt về động lực học vệ tinh mà người vận hành hệ thống cần phải
quan tâm.
Vệ tinh thứ j
(đo được) eijij
rj (đã biết)
Anten thứ i P

Mặt đất

(cần tìm) Ri
C (Gốc tọa độ)
(Tâm trái đất)
Hình 1.5. định vị điểm bằng vệ tinh
Giả sử ta bỏ qua sai số đồng hồ máy thu trên vật thể i và đo được cự ly từ
vật thể i đến vệ tinh 1 là i1, nghĩa là vật thể i đang nằm trên một mặt cầu (S1) có
tâm là vệ tinh 1 (C1) và bán kính là i1. Tiếp theo ta thực hiện phép đo cự ly từ
vật thể i đến vệ tinh 2 và nhận được kết quả là i2, điều này cho chúng ta biết
rằng vật thể i không chỉ nằm trên mặt cầu (S1) mà còn nằm trên mặt cầu (S2) cách
vệ tinh 2 (C2) một khoảng cách là i2. Nói cách khác, vật thể i sẽ nằm trên đường
tròn (O) do 2 mặt cầu (S1), (S2) cắt nhau tạo ra. Nếu chúng ta tiếp tục đo được cự
ly từ vật thể i đến vệ tinh 3 là i3 thì vị trí chính xác của nó là một trong hai giao
điểm P1, P2 của mặt cầu (S3) với đường tròn (O), như ở hình 1.5.
Như vậy, bằng các phép đo cự ly từ vật thể i đến 3 vệ tinh, ta có thể xác
định được 2 vị trí có thể có của nó trong không gian. Để xác định vị trí nào là vị
trí thật ta có thể thực hiện 1 phép đo bổ sung, tuy nhiên 1 trong 2 vị trí tính được
từ phép đo này sẽ cho một kết quả không phù hợp (hoặc là ở rất xa trái đất, hoặc
là chuyển động với vận tốc vô cùng lớn) và do đó có thể bỏ qua mà không cần


của vệ tinh khác, vì vậy các vệ tinh có thể sử dụng cùng tần số mà không làm
nhiễu lẫn nhau. Không những vậy, việc sử dụng mã PRC này còn giúp cho quá
trình xử lý và khuếch đại tín hiệu dựa trên lý thuyết thông tin được thực hiện dễ
dàng hơn, giúp tối ưu hóa anten thu và tiết kiệm chi phí.

Hình 1.8: Mã giả ngẫu nhiên PRC
1.6.1.2. Giả cự ly
Là cự ly đo được giữa vệ tinh và máy thu khi kể đến các sai số đồng hồ
(đồng hồ máy thu và vệ tinh) cũng như các nguồn sai số khác (sai số do lịch thiên
văn, do tầng điện ly, do tầng đối lưu, …). Các cự ly đo được trên hình 1.5 ở trên
chính là các số đo giả cự ly cần thiết để xác định vị trí của vật thể i được tính
theo công thức sau:
prij = ij + c.T
Trong đó :
prij - giả cự ly giữa vật thể i và vệ tinh thứ j ;
ij - cự ly thật giữa vật thể i và vệ tinh thứ j ;
c - vận tốc ánh sáng (3x108 m/s);
T - các nguồn sai số.

1.6.1.3. Đo cự ly bằng sóng xung và sóng liên tục

17

(1.2)


Các hệ thống đo cự ly thường dùng các tín hiệu xung hoặc các tín hiệu
sóng liên tục. Mỗi phương pháp đều có những ưu khuyết điểm riêng và đều có
thể sử dụng trong phép đo một chiều hoặc hai chiều. Hệ thống định vị vô tuyến
toàn cầu là hệ thống đo cự ly một chiều có khả năng sử dụng cả hai loại: sóng

thu tạo ra
Thời trễ
(Giả cự ly)

t
Thời trễ được xác định bằng cách sử dụng nguyên lý tương quan
tín hiệu ngẫu nhiên trong máy thu tương quan (correlator).
Hình 1.9. Phương pháp đo giả cự ly
1.6.2. Các nguồn gây sai số trong phép đo

Như chúng ta đã biết để xác định vị trí của một vật thể, ta cần phải tính
toán được khoảng cách từ nó đến 4 vệ tinh dựa vào phép đo khoảng thời gian
truyền tín hiệu sóng điện từ từ các vệ tinh đến vật thể này. Do đó độ chính xác
của đồng hồ vệ tinh và đồng hồ máy thu là các thông số rất quan trọng ảnh hưởng
đến các phép đo cự ly cần thiết. Ngoài ra, các yếu tố khác như tầng điện ly, tầng
đối lưu, quỹ đạo vệ tinh, độ ồn của máy thu, nhiễu đa đường (multipath)... cũng
góp phần gây ra các sai số không nhỏ trong các phép đo cự ly này.

Tín hiệu truyền từ vệ tinh đến
máy thu bị sai số và bị trễ

Quỹ đạo vệ tinh
Đồng hồ vệ tinh

28.000km

Đồng hồ máy thu
Nhiễu đa đường
Nhiễu máy thu
200km


20


Giả sử rằng, tại một thời điểm nào đó, đồng hồ máy thu có sai số 1ms và
do đó gây ra sai số khoảng cách 300.000m. Nếu các khoảng cách đến tất cả các
vệ tinh được đo chính xác vào cùng một thời điểm thì tất cả khoảng cách này đều
bị lệch 300.000m. Vì vậy, ta có thể xem sai số đồng hồ máy thu là một trong các
ẩn số cần tìm và đó cũng là lý do mà tại sao khi xác định vị trí ta cần thực hiện
các phép đo cự ly đến 4 vệ tinh, nghĩa là cần 4 phương trình để giải ra 4 ẩn số (3
ẩn số vị trí x, y, z và 1 ẩn số thời gian là sai số đồng hồ máy thu), và từ đó giúp ta
có thể sử dụng đồng hồ rẻ tiền và gọn nhẹ hơn trong máy thu.
Chú ý rằng việc xem sai số đồng hồ máy thu là một ẩn số chỉ hợp lệ nếu
chúng ta thực hiện các phép đo cự ly đến các vệ tinh chính xác vào cùng một thời
điểm. Nếu các phép đo này không xảy ra đồng thời thì đối với mỗi phép đo ta sẽ
có một sai số đồng hồ khác nhau. Thực hiện các phép đo đồng thời đến 4 vệ tinh,
ta không những tính toán được vị trí 3 chiều của mình mà còn xác định được sai
số của đồng hồ máy thu với độ chính xác rất cao. Một đồng hồ tiêu biểu có độ
trôi khoảng 1000ns mỗi giây nhưng bằng phương pháp trên ta có thể điều chỉnh
thời gian máy thu đạt độ chính xác bằng với đồng hồ GPS và biến đồng hồ máy
thu rẻ tiền này trở thành một đồng hồ nguyên tử có độ chính xác cao. Máy thu
hiệu chỉnh đồng hồ của nó mỗi giây và cung cấp một tín hiệu thời chuẩn cho các
ứng dụng bên ngoài. Nếu chúng ta đặt máy thu tại một vị trí chính xác đã biết thì
ta chỉ cần theo dõi 1 vệ tinh để tính toán sai số đồng hồ máy thu và điều chỉnh
nó.4 vệ tinh là số lượng tối thiểu mà chúng ta cần để tính toán vị trí và thời gian.
Càng sử dụng nhiều vệ tinh thì kết quả đo nhận được càng chính xác hơn.
1.7. Sai số quỹ đạo vệ tinh
Như đã thảo luận ở các phần trên, độ chính xác của vị trí cần tính toán
cũng phụ thuộc vào cách xác định vị trí chính xác của các vệ tinh (được xem là
các điểm tham chiếu). Quỹ đạo của các vệ tinh liên tục được theo dõi từ nhiều


mang khác nhau để đo độ trễ sai lệch giữa 2 tín hiệu này và từ đó loại bỏ được
ảnh hưởng của tầng điện ly. Đó chính là lý do tại sao mà tất cả các vệ tinh GPS
truyền thông tin bằng 2 tần số L1, L2. Máy thu chính xác (máy thu 2 tần số) chủ
yếu phục vụ cho quân sự theo dõi cả 2 tín hiệu L1, L2 và thực hiện các kỹ thuật
phức tạp để trích ra các tín hiệu mã và sóng mang nhằm loại bỏ ảnh hưởng của
tầng điện ly. Máy thu không chính xác (máy thu đơn tần) phục vụ chủ yếu trong
dân sự chỉ theo dõi 1 tín hiệu L1. Đây là 1 trong những đặc điểm phân biệt chính
giữa 2 loại máy thu này.
1.7.2. Sai số do tầng đối lưu
Lớp thấp hơn của tầng khí quyển chứa đựng hơi nước được gọi là tầng đối
lưu, gây tác động làm chậm cả tín hiệu mã lẫn tín hiệu sóng mang. Ta không thể
loại bỏ ảnh hưởng của tầng đối lưu bằng cách sử dụng hệ thống 2 tần số.
Phương pháp duy nhất để loại bỏ ảnh hưởng của tầng đối lưu là tiến hành
phép đo lượng hơi nước, nhiệt độ, áp suất của nó và áp dụng một mô hình toán
học để có thể tính toán độ trễ gây ra bởi tầng này.
1.7.3. Nhiễu đa đường

Hình 1.13. Nhiễu đa đường
Khi đo khoảng cách đến mỗi vệ tinh, ta giả sử rằng tín hiệu vệ tinh được
truyền thẳng từ vệ tinh đến anten của máy thu. Nhưng trong thực tế ngoài tín hiệu
trực tiếp này anten máy thu còn nhận được các tín hiệu phản xạ đến từ mặt đất và
các vật thể gần anten qua nhiều đường gián tiếp khác nhau, xen nhiễu vào tín
hiệu trực tiếp, gây ra sai lệch về thời điểm đến của tín hiệu thực sự.

23



Nhờ tải bản gốc

Tài liệu, ebook tham khảo khác

Music ♫

Copyright: Tài liệu đại học © DMCA.com Protection Status