giáo dục đào tạo Trường đại học bách khoa hà nội - luận văn thạc sỹ khoa học xác định ảnh hưởng nhiễu từ vệ tinh n-sat-heo tới hệ thống 3G việt nam băng tần 2.6 GHz ngành: Điện tử-Viễn thông lê văn tuấn Người hướng dÉn Khoa häc: PGS.Ts Hå Anh Tuý Hµ néi - 2005 lời cảm ơn Để thực luận này, đà nhận hướng dẫn, giúp đỡ tận tâm thày cô đồng nghiệp Trước hết, xin trân thành cảm ơn PGS TS Hồ Anh Tuý - người thầy đà dìu dắt từ năm Đại học, đà hướng dẫn tận tình cho Tôi trình thực đề tài Các gợi ý góp ý mặt khoa học PGS đà giúp định hướng việc lựa chọn bước thực nhằm hoàn thành tốt ®Ị tµi ®iỊu kiƯn eo hĐp vỊ thêi gian Tôi bày tỏ lòng biết ơn giúp đỡ thày cô giáo Khoa Đào tạo Bồi dưỡng sau Đại học Trường Đại học Bách khoa Hà nội, đặc biệt thày, cô giáo tham gia giảng dạy lớp cao học Điện tử- Viễn thông khoá 2003-2005 đà cung cấp kiến thức, tạo tiền đề cho thực luận văn Tôi xin cảm ơn đồng nghiệp đà tạo điều kiện thuận lợi, hỗ trợ trình Tôi thực luận văn Cuối cùng, xin chân thành cảm ơn gia đình bạn bè đà động viên, quan tâm giúp đỡ hoàn thành luận văn Hà nội, tháng 10 năm 2005 Lê Văn Tuấn Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn Mục lục Trang lời cảm ơn Môc lôc .2 Danh môc hình vẽ Danh mục Bảng biểu .5 PhÇn I : Tỉng quan PhÇn I : Tæng quan Ch¬ng 1: Néi dung, mục đích đề tài .6 1.1 T×nh h×nh sư dơng băng tần 2500-2690MHz giới Việt nam 1.2 Nội dung, mục đích luận văn Ch¬ng 2: Tỉng quan vỊ hƯ thèng N-SAT-HEO vµ hƯ thèng 3G .10 2.1 HÖ thèng vÖ tinh phi ®Þa tÜnh N-SAT-HEO cđa NhËt .10 2.1.1 CÊu h×nh hƯ thèng 10 2.1.2 C¸c tham sè hÖ thèng 12 2.1.2.1 Băng tần 12 2.1.2.2 Ph©n cùc 12 2.1.2.3 §iỊu chÕ sãng mang 12 2.1.2.4 HƯ sè c¾t (roll-off factor) 12 2.1.2.5 Tham số quỹ đạo .12 2.1.2.6 Anten vÖ tinh .14 2.2 Sơ lược cấu trúc hệ thống di động thÕ hÖ ba (3G hay IMT-2000 ) 15 2.2.1 Các dịch vụ 3G 15 2.2.2 C¸c giao diƯn v« tun cđa 3G (IMT-2000) .16 2.2.2.1 Tiêu chuẩn CDMA băng rộng (Wide band CDMA-WCDMA) 16 2.2.2.2 Tiêu chuẩn CDMA đa sóng mang (CDMA-2000 ) .17 2.2.2.3 Tiªu chuÈn CDMA song c«ng theo thêi gian ( CDMA TDD ) 18 2.2.2.4 Tiêu chuẩn TDMA sóng mang đơn (TDMA Single Carrier) 18 2.2.2.5 Tiªu chuÈn FDMA/TDMA 19 2.2.2.6 Lé tr×nh nâng cấp lên 3G Việt nam 19 2.2.3 Các loại cell .19 2.2.4 Băng tần 20 Phần II: Phương pháp xác định ảnh hưởng nhiễu từ hệ thống vệ tinh tíi hƯ thèng 3G 22 Chương 3: ảnh hưởng cđa nhiƠu ®Õn hƯ thèng 3G .22 3.1 NhiƠu hƯ thèng CDMA 22 3.1.1 Kh¸i niƯm vỊ nhiƠu 22 3.1.1.1 Khái niệm phân lo¹i nhiƠu 22 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn 3.1.1.2 Nhiễu cộng nhiễu nhân -tû sè tÝn hiƯu/nhiƠu 23 3.1.2 Kỹ thuật trải phổ CDMA 25 3.1.2.1 Giíi thiƯu chung 25 3.1.2.2 Băng thông tín hiệu trải phổ .27 3.1.2.3 M· ngÉu nhiªn 27 3.1.2.4 HƯ thèng m¸y ph¸t, m¸y thu tr¶i phỉ trùc tiÕp (DSSS) .29 3.1.2.4.1 M¸y ph¸t DSSS-BPSK .29 3.1.2.4.2 M¸y thu DSSS-BPSK 30 3.1.2.5 Mật độ phổ công suất 32 3.1.2.6 Điều khiển công suÊt hÖ thèng CDMA 34 3.1.3 NhiƠu hƯ thèng CDMA .34 3.1.3.1 Víi nhiƠu t¹p ©m Gausse tr¾ng céng (AWGN), No 36 3.1.3.2 Nhiễu băng hẹp 36 3.1.3.3 Nhiễu băng réng 37 3.1.3.4 Nhiễu đồng kênh (từ người sư dơng kh¸c) 37 3.2 ¶nh hëng cđa nhiƠu tíi m¹ng 3G 39 3.2.1 Tæng quan .39 3.2.2 Phương pháp tính suy giảm vùng phủ sãng cđa hƯ thèng 3G nhiƠu 40 3.2.2.1 Tính toán hệ số tăng tạp âm nhiệt đường lên CDMA 41 3.2.2.2 Tính toán hệ sè suy gi¶m vïng phđ cđa sector 44 3.2.2.3 ¶nh hëng cđa nhiƠu tíi mét site 3G cã nhiỊu sector 45 Ch¬ng 4: Phương pháp tính toán mức nhiễu .49 4.1 Phương pháp tính mức nhiễu từ hệ thống N-SAT-HEO vào trạm gốc 3G .49 4.1.1 ThuËt to¸n .49 4.1.2 HƯ sè ph©n cách phân cực .54 4.1.3 C¸c hƯ sè suy hao kh¸c 54 4.1.4 Tính toán độ tăng ích búp sóng phụ cđa anten tr¹m gèc 3G 54 4.2 Xác định vị trí tức thời vệ tinh phi địa tĩnh không gian 55 4.2.1 Ba ®Þnh luËt Kepler .55 4.2.2 Các tham số quỹ đạo vÖ tinh 58 4.2.3 Xác định vị trí vệ tinh kh«ng gian 60 4.2.3.1 Xác định toạ độ vệ tinh .60 4.2.3.2 Xác định Ly giác nút lên u(t) theo thời gian 64 4.2.4 Xác định góc ngẩng góc phương vị 69 Phần III: Xác định ảnh hưởng nhiễu với ®iỊu kiƯn thùc tÕ ViƯt nam .71 Chương 5: Thực xác định ảnh hưởng nhiƠu víi ®iỊu kiƯn thùc tÕ cđa ViƯt nam 71 5.1 Thuật toán chương tr×nh thùc hiƯn phÐp tÝnh 71 5.2 C¸c tham sè sư dơng .74 5.2.1 HÖ thèng vÖ tinh N-SAT-HEO 74 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn 5.2.2 Hệ thống 3G 75 5.3 Chương trình dùng bảng tính excel .77 5.3.1 PhÇn số liệu đầu vào 77 5.3.2 Phần tính toán 79 5.3.3 Phần hiển thị kết 82 5.4 Ch¬ng tr×nh Visual Basic 84 5.4.1 Giao diện chương trình 85 5.4.2 Phần mà chương trình .88 5.4.3 Hiển thị kết chi tiÕt .88 Chương 6: Kết quả, nhận xét kiến nghị 90 6.1 Kết tính cho số địa điểm cụ thể t¹i ViƯt nam 90 6.2 NhËn xÐt 91 6.3 KiÕn nghÞ 93 tài liệu tham khảo .95 Phô lôc 1: Mà chương trình VB 97 Phụ lục 2: Chương trình Excel 112 Danh mục hình vẽ Hình 1: Bảng phân bổ tần số băng 2.6 GHz ITU .7 H×nh 2:HƯ thèng vƯ tinh N-SAT-HEO 11 Hình 3: Vết di chuyển mặt đất cđa vƯ tinh N-SAT-HEO 14 Hình 4: Lộ trình nâng cấp lên 3G mạng di ®éng ViƯt nam 19 H×nh 5: loại cell 3G (IMT-2000) 20 Hình 6: Các phương án quy hoạch tần số cho IMT-2000 .21 Hình 7:Chuỗi PN với chu kỳ N=15 28 Hình 8:Sơ đồ khối máy phát DSSS BPSK dạng tín hiệu .29 Hình 9: sơ đồ khối máy thu DSSS-BPSK dạng tÝn hiƯu .31 H×nh 10: phỉ c«ng st hƯ thèng DSSS-BPSK 33 Hình 11: Các loại nhiễu hệ thèng CDMA .35 Hình 12: Phổ tín hiệu băng hẹp gây nhiễu cho mạng S-Fone Tp.HCM .37 Hình 13: Nhiễu đồng kênh (từ người sử dụng khác) .38 Hình 14: Nhiễu làm giảm khoảng c¸ch phđ sãng 40 Hình 15:Biến thiên hệ số tăng tạp âm đường lên data throughput 43 Hình 16: ¶nh hëng cđa nhiƠu tíi cell cã sector 46 Hình 17: ảnh hưởng nhiễu tới trạm gốc ni=0.5, 1, dB 47 Hình 18: Các tham số quỹ đạo vệ tinh 58 Hình 19:Tam giác cầu ABC 61 Hình 20: Xác định toạ độ vệ tinh 62 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn Hình 21:Lưu đồ thuật toán 73 Hình 22: Modun tham số quỹ đạo 78 H×nh 23: Modun tham số phát xạ 79 H×nh 24: Modun xác định góc nút u(t) 80 Hình 25: Modun tính toạ độ, khoảng cách vệ tinh 80 H×nh 26: Modun tÝnh gãc ngÈng, gãc phương vị 81 H×nh 27: Modun tÝnh Off-axis Gain cđa anten Sector 81 H×nh 28: Modun tÝnh Isat/Nth modun tính độ suy giảm vùng phủ sóng 82 Hình 29: Bảng kết .82 Hình 30: Vết quỹ đạo vệ tinh mặt ®Êt .83 Hình 31: Đồ thị suy giảm vùng phủ sóng với pfd=-113dBW/m2.MHz 84 Hình 32: Giao diện chương tr×nh chung 86 Hình 33: Giao diện chương trình tính Isat/Nth .86 H×nh 34: Bảng kết tổng hợp 11 tỉnh khảo sát 89 Hình 35: Suy giảm vùng phủ sóng với mức pfd khác 90 Hình 36: Suy giảm vùng phủ 1BTS Tp.HCM vơi mức pfd khác 91 Danh mục Bảng biểu Bảng 1:Các tham số đường lên hệ thống 3G 42 Bảng 2: Số lượng người sử dụng gây hệ số tăng tạp âm .43 B¶ng 3: Sù suy giảm vùng phủ sóng theo giá trị Isat/Nth 48 B¶ng 5: Tham sè hƯ thèng N-SAT-HEO 74 B¶ng 6: Tham sè hƯ thèng 3G 76 Bảng 7: Toạ độ địa điểm đặt BTS 76 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn Phần I : Tổng quan Chương 1: Nội dung, mục đích đề tài 1.1 Tình hình sử dụng băng tần 2500-2690MHz giới Việt nam Do đặc tính truyền lan không phụ thuộc biên giới hành sóng vô tuyến điện, nên việc sử dụng hệ thống thông tin vô tuyến nước nói chung Việt nam nói riêng phải tuân theo quy định quốc tế, cụ thể Thể lệ Vô tuyến điện (Radio Regulations-RR) Liên minh ViƠn th«ng Qc tÕ (the International Telecommunications Union-ITU ) Thể lệ Vô tuyến điện Hội nghị Thông tin v« tun ThÕ giíi (the Worl Radiocommunications ConferenceWRC) th«ng qua sửa đổi Các WRC thường tổ chức theo chu kỳ từ đến năm lần Thông tin di động năm gần đà có sù bïng nỉ vỊ sè lỵng ngêi sư dơng víi tốc độ phát triển lên tới hàng chục % năm số thuê bao điện thoại di động đà vượt qua số thuê bao điện thoại cố định Bên cạnh đó, xu chuyển đổi từ hệ thống thông tin động hệ (2G) lên hệ thống thông tin di động hệ 3, hay IMT-2000 theo cách gọi ITU diễn mạnh mẽ Hai mạng Vinaphone MobiFone đà triển khai dịch vụ GPRS công nghệ GSM, mạng S-Fone triển khai công nghệ CDMA-2000 1x tới mạng di động Công ty viễn thông điện lực triển khai công nghệ CDMA2000 1x EV-DO Thấy trước nhu cầu phát triển mạng di động 3G, từ năm 2000, ITU đà thông qua việc sử dụng thêm băng tần 2500-2690MHz (2.6GHz) cho hệ thống toàn cầu Trong Quy hoạch phổ tần quốc gia Việt nam, băng tần đà xác định giành cho việc triển khai hệ thống 3G tương lai gần Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn Hội nghị thông tin vô tuyến giới năm 2000 (WRC-2000) đà thông qua định cho phép số nước khu vực Châu Thái bình dương ( Khu vực theo phân chia ITU) triển khai hệ thống Phát số (sound digital broadcasting) qua vệ tinh phi địa tĩnh (non-GSO satellites) băng tần 2630-2655MHz Hội nghị Thông tin vô tuyến giới năm 2003 (WRC-03) tiếp tục mở rộng, cho phép nước sử dụng thêm băng tần 26052630MHz loại hình vệ tinh địa tĩnh (GSO satellite) để cung cấp loại hình dịch vụ phát số qua vệ tinh nói Hình 1: Bảng phân bổ tần số băng 2.6 GHz ITU Ngoài ra, yếu tố lịch sử, băng tần 2500-2690MHz sử dụng cho loại hình thông tin vô tuyến cố định (viba, hệ thống truyền hình MMDS, ), hệ thống thông tin di động mặt đất (hệ thống vô tuyến điều hành taxi, ) thông tin vệ tinh cố định, dịch vụ truyền hình qua vệ tinh Hình trình bày Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn bảng loại hình thông tin vô tuyến phép sử dụng băng tần 25202655MHz Thể lệ Vô tuyến điện 2004 ITU Các quy định liên quan ®Õn c¸c hƯ thèng ph¸t sè qua vƯ tinh nêu điều khoản 5.418 Quy định sử dụnng IMT-2000 nêu điều khoản 5.384A Điều khoản 5.418 ghi Hàn quốc, ấn độ, Nhật bản, Pakistan Thái lan băng tần 2535-2655MHz phân bổ cho nghiệp vụ Thông tin quảng bá (phát thanh) qua vệ tinh hệ thống phụ trợ mặt đất nghiệp vụ Việc sử dụng hạn chế cho hệ thống phát số Điều khoản 5.384A ghi Các băng tần, hay số đoạn băng tần băng, 1710-1885 MHz 2500-2690MHz, xác định để nước sử dụng cho hÖ thèng IMT-2000 ” Nh vËy, khu vùc có số nước sử dụng băng tần 2605-2655MHz cho hƯ thèng ph¸t sè qua vƯ tinh, mét sè níc sư dơng cho hƯ thèng th«ng tin di động hệ Trên sở quy định ITU, Nhật có kế hoạch triển khai hƯ thèng ph¸t sè qua vƯ tinh phi địa tĩnh (non-GSO BSS) với tên gọi NSAT-HEO băng tần 2630-2655MHz dự kiến cung cấp dịch vụ vào năm 2008 Hàn quốc có dự định xây dựng hệ thống khác làm việc băng tần 2605-2630 MHz Với Việt nam, băng tần 2500-2690MHz Đài truyền hình Việt nam sử dụng cho hệ thống truyền hình MMDS Tuy nhiên, thời gian tới, theo Quy hoạch phổ tần quốc gia băng tần dự kiến sử dụng cho hệ thống thông tin di động 3G (IMT-2000) Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn 1.2 Nội dung, mục đích luận văn Dù hệ thèng vƯ tinh N-SAT-HEO chØ cung cÊp dÞch vơ ë Nhật bản, hệ thống vệ tinh có công suất lớn, sử dụng băng tần trùng với băng tần hệ thống 3G nên khả búp sóng phơ cđa hƯ thèng vƯ tinh nµy vÉn cã møc công suất lớn, gây ảnh hưởng tới hệ thống 3G nước khu vực, có Việt nam Hiện nay, Nhóm nghiên cứu (Study Group) ITU-R nghiên cứu ảnh hưởng hệ thống tới hệ thống thông tin vô tuyến mặt đất băng tần 2.6 GHz để xác lập điều kiện mức công suất phát xạ phụ mà hệ thống phải tuân theo ITU đà thống đưa số khuyến nghị (Recommendation) phương pháp tính (methodology) tham số đầu vào làm sở cho việc tính toán chi tiết Mục tiêu luận văn tìm hiểu phương pháp tính toán, xác định mức độ ảnh hưởng hệ thống vệ tinh quỹ đạo cao N-SAT-HEO Nhật tới hệ thống 3G, mà cụ thể mức độ suy giảm vùng phủ sóng trạm BTS (Base Transceiver Station), từ áp dụng vào điều kiện Việt nam Kết tính sử dụng để xác định mức công suất phát xạ tối đa chấp nhận từ hệ thống N-SAT-HEO tới lÃnh thổ Việt nam để làm sở đàm phán với Nhật xây dựng quan điểm Việt nam vấn đề ITU Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn Dim ItoN(0 To 10000, To 11) As Double, ItoN1(0 To 10000, To 11) As Double, ItoN2(0 To 10000, To 11) As Double, ItoN3(0 To 10000, To 11) As Double Dim ItoNmax(1 To 11) As Double, ItoNmin(1 To 11) As Double Dim DeltaA(1 To 10000, To 11) As Double, DeltaA1(1 To 10000, To 11) As Double, DeltaA2(1 To 10000, To 11) As Double, DeltaA3(1 To 10000, To 11) As Double Dim DeltaAmax(1 To 11) As Double, DeltaAmin(1 To 11) As Double Dim CoverageLossMax As Double, CoverageLossMin As Double Dim IoverNmax As Double, IoverNmin As Double Dim Go As Double, G1 As Double, G2 As Double, G_phi As Double Dim G_sector2 As Double, G_sector3 As Double, Phi3 As Double Dim Nth As Double Dim rad As Double, deg As Double Dim ConvertX As Double, ConvertY As Double Dim PointX As Double, PointY As Double Dim PointXo As Double, PointYo As Double Dim PointXo1 As Double, PointYo1 As Double, PointYo2 As Double Dim PointX1 As Double, PointY1 As Double, PointY2 As Double ' Khai bao cac Hang so Const re = 6378 Const Muy = 398600 Te = 360 / 86164 ' van toc goc trung binh cua trai dat rad = Pi / 180 deg = 180 / Pi '********Ve ban va he truc toa do** ‘ Load ban so Image1.Refresh 98 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn Image1.Picture = LoadPicture("C:\SatCalculation\Picture\WorldColor.bmp") Image1.Left = 400 Image1.Top = Label1.Move Image1.Left - 250, Image1.Height / - 50 Label2.Move Image1.Left - 250, Image1.Height / - 50 - Image1.Height / Label3.Move Image1.Left - 250, Image1.Height / - 50 - * Image1.Height / Label4.Move Image1.Left - 250, Image1.Height / - 50 - * Image1.Height / Label5.Move Image1.Left - 250, Image1.Height / - 50 - * Image1.Height / Label6.Move Image1.Left - 300, Image1.Height / - 50 + Image1.Height / Label7.Move Image1.Left - 300, Image1.Height / - 50 + * Image1.Height / Label8.Move Image1.Left - 300, Image1.Height / - 50 + * Image1.Height / Label9.Move Image1.Left - 300, Image1.Height / - 50 + * Image1.Height / Label15.Move 200 + * Image1.Width / 18, Image1.Height + 50 Label16.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label17.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label18.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label19.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label20.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label21.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label22.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label23.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label24.Move 300 + * Image1.Width / 18, Image1.Height + 50 Label25.Move 300 + 10 * Image1.Width / 18, Image1.Height + 50 Label26.Move 300 + 11 * Image1.Width / 18, Image1.Height + 50 Label27.Move 300 + 12 * Image1.Width / 18, Image1.Height + 50 Label28.Move 300 + 13 * Image1.Width / 18, Image1.Height + 50 Label29.Move 300 + 14 * Image1.Width / 18, Image1.Height + 50 Label30.Move 300 + 15 * Image1.Width / 18, Image1.Height + 50 99 LuËn văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn TuÊn Label31.Move 300 + 16 * Image1.Width / 18, Image1.Height + 50 Label37.Move 300 + 17 * Image1.Width / 18, Image1.Height + 50 Label45.Move 300 + 18 * Image1.Width / 18, Image1.Height + 50 '*************NHAP SO LIEU DAU VAO ************************** Apogee = re + Text1.Text Perigee = re + Text2.Text Inc = (Text3.Text) * rad ' tinh bang don vi rad LongSat0 = Text4.Text ' Kinh cua Ascending Node, degree Omega0 = (Text5.Text) * rad 'Quy dao ellipse, argument cua perigee, cho truoc va tinh theo don vi rad ActiveArcStart = Text6.Text * 3600 ' thoi diem ve tinh bat dau hoat dong, s ActiveArcEnd = Text7.Text * 3600 ' thoi diem ve tinh ket thuc hoat dong, s ' NHAP Dia danh cac Dia phuong ma ta khao sat Diadanh(1) = "Lang Son" BTSlong(1) = 106.92 * rad BTSlat(1) = 21.83 * rad Diadanh(2) = "Dien Bien" BTSlong(2) = 103.03 * rad BTSlat(2) = 21.38 * rad Diadanh(3) = "Hai Phong" BTSlong(3) = 106.67 * rad BTSlat(3) = 20.8 * rad Diadanh(4) = "Ha Noi" BTSlong(4) = 105.88 * rad BTSlat(4) = 21.02 * rad 100 Luận văn tèt nghiƯp cao häc §TVT 2003- 2005 Diadanh(5) = "Hue" BTSlong(5) = 107.58 * rad BTSlat(5) = 16.47 * rad Diadanh(6) = "Da Nang" BTSlong(6) = 108.23 * rad BTSlat(6) = 16.07 * rad Diadanh(7) = "Nha Trang" BTSlong(7) = 109.17 * rad BTSlat(7) = 12.25 * rad Diadanh(8) = "Tp.HCM" BTSlong(8) = 106.72 * rad BTSlat(8) = 10.77 * rad Diadanh(9) = "Vung Tau" BTSlong(9) = 107.07 * rad BTSlat(9) = 10.35 * rad Diadanh(10) = "Can Tho" BTSlong(10) = 105.77 * rad BTSlat(10) = 10.05 * rad Diadanh(11) = "Bac Lieu" BTSlong(11) = 105.75 * rad BTSlat(11) = 9.27 * rad Go = Text8.Text PFD = Text9.Text 101 Lê Văn Tuấn Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Nth = Text10.Text ' TINH CAC THAM SO THU CAP ‘ Canh bao cao Vien diem < cao Can diem If Apogee < Perigee Then MsgBox ("Apogee must be greater than or equal Perigee") Else E = (Apogee - Perigee) / (Apogee + Perigee) End If ' Xac dinh goc Doi so cua can diem If Omega0 * deg >= 180 Then UkApogee = Omega0 - Pi ' rad Else UkApogee = Omega0 + Pi End If 'Ban kinh lon cua ellipse a = (Apogee + Perigee) / 'km ' Chu ky quy dao Period = * Pi * Sqr(a * a * a / Muy) 'seconds ' Ro Ro = * Apogee * Perigee / (Apogee + Perigee) ' He so anpha Anpha = (1 - E * E) ^ 1.5 / (2 * Pi) ' Buoc lay mau thoi gian DeltaT = Period / 10000 'seconds ' Xac dinh dieu kien dau- thoi diem ve tinh di qua Ascending Node ts = Uk(0) = 102 Lê Văn Tuấn Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn PUk(0) = Anpha / (1 + E * Cos(-Omega0)) ^ '***** XAC DINH TOA DO VE TINH, GOC NGANG, GOC PHUONG VI ****** For k = To 10000 'Lay chu ky quy dao, chia nho 10000 diem ' Tinh Ly giac goc nut u(t) ts = ts + DeltaT ' quet theo thoi gian voi buoc nhay DeltaT ' tinh goc quet Uk tai thoi diem tk Uk(k) = Uk(k - 1) + / (10000 * PUk(k - 1)) 'Ham Xac xuat xuat hien ve tinh khoang Uk + dUk PUk(k) = Anpha / (1 + E * Cos(Uk(k) - Omega0)) ^ 'Tinh toa tai thoi diem tk cua Ve tinh 'Vi cua ve tinh, rad LatSat(k) = asin(Sin(Inc) * Sin(Uk(k))) ' Goc giua Ascending Node va vi tri tuc thoi cua Ve tinh If Cos(Inc) * Sin(Uk(k)) > Then So = acos(Cos(Uk(k)) / Cos(LatSat(k))) ' rad Else So = * Pi - acos(Cos(Uk(k)) / Cos(LatSat(k))) 'rad End If 'Kinh tai thoi diem tk cua Ve tinh tinh theo deg LongSat(k) = So * deg + LongSat0 - 7.292115856 * 10 ^ -5 * ts * deg ' Xac dinh vi tri cua Vien diem (Apogee) va thoi gian ve tinh di qua If Uk(k - 1) UkApogee Then kApogee = k - LongSatApogee = LongSat(k - 1) LatSatApogee = LatSat(k - 1) 103 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn End If ' Xac dinh khoang cach Tam trai dat - Ve tinh: Ru tai thoi diem tk Ru(k) = Ro / (1 + E * Cos(Uk(k) - Omega0)) ' XAC DINH GOC NGANG, GOC PHUONG VI MA TRAM BTS QUAN SAT VE TINH For j = To 11 X = Ru(k) * Cos(LatSat(k)) * Cos(LongSat(k) * rad - BTSlong(j)) Y = Ru(k) * Cos(LatSat(k)) * Sin(LongSat(k) * rad - BTSlong(j)) Z = Ru(k) * Sin(LatSat(k)) x1 = X * Cos(BTSlat(j)) + Z * Sin(BTSlat(j)) y1 = Y Z1 = -X * Sin(BTSlat(j)) + Z * Cos(BTSlat(j)) ' Goc ngang BTS quan sat ve tinh Elevation(k, j) = Atn((x1 - re) / Sqr(y1 ^ + Z1 ^ 2)) * deg ' Goc phuong vi BTS quan sat ve tinh If Z1 > Then Azimuth(k, j) = Atn(y1 / Z1) * deg Else Azimuth(k, j) = Atn(y1 / Z1) * deg + 180 End If Next Next ' Ket thuc vong lap k, KET THUC TINH TOA DO VE TINH, GOC NGANG, GOC PHUONG VI ' ******* VE VET QUY DAO - Ground Track ****** 104 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn ConvertY = Image1.Height / 180 ConvertX = Image1.Width / 360 PointXo = Round(LongSat(1) * ConvertX, 0) PointYo = Round(Image1.Height / - LatSat(1) * ConvertY * 180 / Pi, 0) Image1.Height / ' ' Bat dau qua trinh ve For k = To 10000 PointX = Round(LongSat(k) * ConvertX, 0) PointY = Round(Image1.Height / - LatSat(k) * ConvertY * 180 / Pi, 0) If Abs(PointX - PointXo) < Image1.Width / Then ' khong ve duong noi giua 360 va ' Trong doan active arc ve mau DO - RED If k >= kApogee - Round(16200 / DeltaT, 0) And k ItoN1(k, j) Then ItoNmin(j) = ItoN1(k, j) End If ' Xac dinh he so suy giam vung phu song cua tung sector va cua ca BTS DeltaA1(k, j) = (1 + 10 ^ ((ItoN1(k, j) - Noiserise) / 10)) ^ (-20 / 35.2) DeltaA2(k, j) = (1 + 10 ^ ((ItoN2(k, j) - Noiserise) / 10)) ^ (-20 / 35.2) DeltaA3(k, j) = (1 + 10 ^ ((ItoN3(k, j) - Noiserise) / 10)) ^ (-20 / 35.2) DeltaA(k, j) = (DeltaA1(k, j) + DeltaA2(k, j) + DeltaA3(k, j)) / ' Xac dinh he so suy giam vung phu song thap nhat va cao nhat If DeltaAmax(j) < DeltaA(k, j) Then DeltaAmax(j) = DeltaA(k, j) End If If DeltaAmin(j) > DeltaA(k, j) Then DeltaAmin(j) = DeltaA(k, j) End If Next 'j Next 'k CoverageLossMax = DeltaAmax(1) CoverageLossMin = DeltaAmax(1) IoverNmax = ItoNmax(1) IoverNmin = ItoNmax(1) ‘ Xac dinh ty so Isat/Nth max tren ca nuoc For j = To 11 'Xac dinh ty so Isat/Nth lon nhat 11 dia diem khao sat If IoverNmax < ItoNmax(j) Then 107 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn IoverNmax = ItoNmax(j) End If ' Xac dinh muc suy giam vung phu song lon nhat If CoverageLossMax > DeltaAmax(j) Then CoverageLossMax = DeltaAmax(j) End If ' Xac dinh muc suy giam vung phu song be nhat If CoverageLossMin < DeltaAmax(j) Then CoverageLossMin = DeltaAmax(j) End If Next ' Hien thi ket qua tren giao dien Text11.Text = IoverNmax Text12.Text = CoverageLossMax Text13.Text = CoverageLossMin Text14.Text = 100 - CoverageLossMax * 100 '************ Xuat du lieu Excel*************** Dim cot As String Set XL = New Excel.Application XL.Visible = False XL.Workbooks.Add XL.Sheets.Add.Name = "Ket qua tong hop" XL.Cells(1, 3) = "KET QUA CHUONG TRINH TINH ANH HUONG CUA NHIEU" XL.Cells(3, 2) = "Coverage Reduction" XL.Cells(3, 4) = "Increased BTS" XL.Cells(3, 6) = "Isat/Nth" XL.Cells(3, 7) = "Isat/Nth" 108 Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 XL.Cells(4, 2) = "max,%" XL.Cells(4, 3) = "min,%" XL.Cells(4, 4) = "max,%" XL.Cells(4, 5) = "min,%" XL.Cells(4, 6) = "max,dB" XL.Cells(4, 7) = "min,dB" XL.Cells(4, 8) = "degrees" XL.Cells(4, 9) = "degrees" XL.Cells(4, 10) = "Long" XL.Cells(4, 11) = "Lat" XL.Cells(3, 8) = "Elevation" XL.Cells(3, 9) = "Azimuth" XL.Cells(3, 10) = "Satellite position" XL.Cells(17, 1) = "MAX" XL.Cells(18, 1) = "MIN" XL.Cells(17, 2) = CoverageLossMax XL.Cells(18, 3) = CoverageLossMin XL.Cells(17, 4) = 100 - 100 * CoverageLossMax XL.Cells(18, 5) = 100 - 100 * CoverageLossMin XL.Cells(17, 6) = IoverNmax XL.Cells(18, 7) = IoverNmin ‘ Xuat so lieu cac cell cua Excel For j = To 11 XL.Cells(4 + j, 1) = Diadanh(j) XL.Cells(4 + j, 2) = DeltaAmax(j) XL.Cells(4 + j, 3) = DeltaAmin(j) XL.Cells(4 + j, 4) = 100 - 100 * DeltaAmin(j) XL.Cells(4 + j, 5) = 100 - 100 * DeltaAmax(j) XL.Cells(4 + j, 6) = ItoNmax(j) 109 Lê Văn Tuấn Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Lê Văn Tuấn XL.Cells(4 + j, 7) = ItoNmin(j) XL.Cells(4 + j, 8) = Elevation(l(j), j) XL.Cells(4 + j, 9) = Azimuth(l(j), j) XL.Cells(4 + j, 10) = LongSat(l(j)) XL.Cells(4 + j, 11) = LatSat(l(j)) * deg Next XL.Sheets.Add.Name = "Ket qua chi tiet" j=1 XL.Cells(1, 1) = "HA NOI" XL.Cells(2, 1) = "tk" XL.Cells(2, 2) = "Reduced" XL.Cells(3, 2) = "Coverage" XL.Cells(1, 4) = "DA NANG" XL.Cells(2, 4) = "tk" XL.Cells(2, 5) = "Reduced" XL.Cells(3, 5) = "Coverage" XL.Cells(1, 7) = "Tp HCM" XL.Cells(2, 7) = "tk" XL.Cells(2, 8) = "Reduced" XL.Cells(3, 8) = "Coverage" ‘ Hien thi ket qua tinh ve tinh Active Arc For k = kApogee - Round(16200 / DeltaT, 0) To kApogee + Round(12600 / DeltaT, 0) Step 10 XL.Cells(3 + j, 1) = k / 10000 XL.Cells(3 + j, 2) = DeltaA(k, 4) XL.Cells(3 + j, 4) = k / 10000 110 Luận văn tốt nghiệp cao học §TVT 2003- 2005 XL.Cells(3 + j, 5) = DeltaA(k, 6) XL.Cells(3 + j, 7) = k / 10000 XL.Cells(3 + j, 8) = DeltaA(k, 8) j=j+1 Next End Sub ' KET THUC TINH TOAN Private Sub Command2_Click() Form12.Hide Unload Form12 End Sub Private Sub Command3_Click() XL.Visible = True End Sub 111 Lê Văn Tuấn Luận văn tốt nghiệp cao học ĐTVT 2003- 2005 Phụ lục 2: Chương trình Excel 112 Lê Văn Tuấn ... phương pháp nhằm đánh giá mức nhiễu từ hệ thống vệ tinh N- SAT- HEO hệ thống IMT -20 00 băng t? ?n 26 30 - 26 55MHz sử dụng để xác đính mức độ ảnh hưởng cđa hƯ thèng vƯ tinh tíi hƯ thèng 3G trường hợp hệ thống. .. gây ảnh hưởng tới hệ thống 3G n? ?ớc khu vực, có Việt nam Hi? ?n nay, Nhóm nghi? ?n cứu (Study Group) ITU-R nghi? ?n cứu ảnh hưởng hệ thống tới hệ thống thông tin vô tuy? ?n mặt đất băng t? ?n 2. 6 GHz để xác. .. b? ?n kÝnh tíi 3km 2. 1 .2 C¸c tham sè hệ thống 2. 1 .2. 1 Băng t? ?n Như đà trình bày tr? ?n, băng t? ?n hệ thống 26 30 - 26 55MHz, t? ?n số phát tới máy thu Băng thông 25 MHz 2. 1 .2. 2 Ph? ?n cực Ph? ?n cực sử dụng