phạm toàn thắng Bộ 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 Ngành: điện tử viễn thông điện tử viễn thông tối ưu hóa quỹ công suất đường truyền thông tin vệ tinh phạm toàn thắng 2004 - 2006 Hµ Néi 2006 Hµ néi 2006 Mơc lục Danh mục chữ viết tắt 18T 18T Danh môc bảng 18T 18T Danh mục hình vẽ 18T 18 T Mở đầu 18T 18T Ch¬ng 1: Tỉng quan 1.1 Giíi thiƯu 1.2 Các nghiệp vụ thông tin vệ tinh 10 1.3 Lịch sử phát triĨn cđa th«ng tin vƯ tinh 11 1.4 Tình hình sử dụng vị trí quỹ đạo vệ tinh địa tĩnh 14 1.5 Mét sè hƯ thèng vƯ tinh trªn thÕ giíi 17 1.6 Các quy định ITU 18 18T 18 T 18T 18T 18T 18T 18T 18 T 18T 18T 18T 18T 18T 18T Chương 2: đường truyền vệ tinh băng tần sử dụng 21 2.1 Đặc điểm đường truyền vÖ tinh 21 2.1.1 Tuyến thông tin vệ tinh 21 2.1.2 Các đặc tính anten 23 2.1.2.1 Độ khuếch đại anten 23 2.1.2.2 §é më hiƯu dụng độ khuếch đại anten 24 2.1.3 Công suất phát công suất thu anten 26 2.1.4 Công suất tạp ©m 29 2.1.4.1 Nhiệt độ tạp âm 29 2.1.4.2 Nhiệt độ tạp ©m cđa m¸y thu 30 2.1.5 Phân hệ không gian 32 2.1.6 Phân hệ mặt đất 35 2.2 Phân bổ tần số thông tin vệ tinh 37 2.2.1 Nguyên tắc phân bổ tần số cho thông tin vệ tinh ITU 37 2.2.2 Phân bổ băng tần cho c¸c nghiƯp vơ 38 2.3 Thông số đặc trưng mạng vệ tinh băng tần 40 2.3.1 Băng tần C 40 2.3.2 Băng tần X 41 2.3.3 Băng tần Ku 42 2.3.4 Băng tần Ka 43 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T Ch¬ng 3: phèi hợp tần số quỹ đạo vệ tinh địa tĩnh 46 3.1 Quy trình phối hợp 46 3.1.1 Giai đoạn xuất b¶n tríc 46 3.1.2 Giai đoạn phối hợp 47 3.1.3 Giai đoạn thông báo, ghi ấn định tần số vào Bảng tần số chñ 51 18T 18 T 18T 8T 18T 18T 18T 18T 18T Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 18T 3.2 Các phương pháp xác định cần thiết phải thực phối hợp 53 3.2.1 Phương pháp giới hạn cung quỹ đạo phối hợp 55 3.2.2 Phương pháp giới h¹n tû sè ∆T/T 56 3.2.3 Phương pháp giới hạn tỷ số C/I 57 18T 18 T 18T 18 T 18T 18 T 18T 18T ch¬ng 4: Tính toán quỹ công suất đường truyền 59 4.1 Khái niệm quỹ công suất đường truyền 59 4.2 Quỹ công suất đường lên (C/N ) u 61 4.3 Quỹ công suất đường xuống (C/N ) d 63 4.4 Quỹ công suất đường truyền cho phát đáp truyền thống 63 4.4.1 Kết hợp đường lên đường xuống (C/N ) ud 64 4.4.2 ảnh hưởng tạp âm khác tới quỹ công suất đường truyền 65 4.4.3 Quỹ công suất toàn tuyến 66 4.4.4 Suy hao khÝ quyÓn 67 4.4.5 Mét sè chó ý 68 4.5 Qịy c«ng st ®êng trun cho bé ph¸t ®¸p cã xư lý 69 18T 18 T 18T 18 T 18T R R R 18T R 18T R R R R 18T 18T 18 T 18T R R R R 18T 18T 18 T 18T 18T 18T 18T 18T 18T 18T 18T Chương 5: tối ưu hóa quỹ công suất ®êng trun 71 5.1 Tỉng quan 71 5.2 TÝnh to¸n suy hao ma 72 5.2.1 C¸c bước xác định suy hao mưa 72 5.2.2 Xác định tỷ lệ suy hao theo phân cực tần số 77 5.2.3 Mô hình xác định suy hao đặc trng ma 78 5.2.4 X¸c ®Þnh ®é cao ma 81 5.3 TÝnh to¸n suy hao khÝ quyÓn 81 5.3.1 Suy hao đặc trưng 82 5.3.2 Suy hao ®êng trun 83 5.4 Tèi u hãa 86 5.4.1 Các thông số đầu vào 86 5.4.2 Bài toán thuận 87 5.4.3 Bài toán nghịch 89 5.4.3.1 C¸c bíc thùc hiƯn 89 5.4.3.2 Lưu đồ thuật toán 91 5.5 Chương trình phÇn mỊm tèi u hãa 91 5.5.1 Bài toán áp dụng băng tÇn C 94 5.5.2 Bài toán áp dụng băng tần Ku 98 5.5.3 KÕt luËn 101 18T 8T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T Lêi kÕt 103 18T 18T Tµi liƯu tham kh¶o 104 18T 18T phụ lục Mà chương trình Visual Basic Net 106 18T Ph¹m Toàn Thắng Luận văn cao học ĐTVT 2004-2006 18T Danh mục chữ viết tắt BER Bit Error Rate BR Bureau Radiocommunication BSS CATV CS CV Radiocomunication Bureau International Frequency Information Circular Broadcasting Satellite Service Cable Television Constitution Convention DTH Direct To Home BRIFIC EIRP FS FSS HEO HPA ISS ITU LEO LHCP LNA LNB LNC MR MEO MS MSS Equivelent Isotropic Radiation Power Fixed Service Fixed Satellite Service High Earth Orbit High Power Amplifier Inter-Satellite Service International Telecommunication Union Low Earth Orbit Left Hand Circular Polarization Low Noise Amplifier Low Noise Block Low Noise Converter Master Register Medium Earth Orbit Mobile Service Mobile Satellite Service Tû sè lỗi bít Văn phòng thông tin vô tuyến ITU Đĩa CD chứa thông tin cập nhật phát triển mạng vệ tinh giới Nghiệp vụ quảng bá qua vệ tinh Truyền hình cáp Hiến chương ITU Công ước ITU Dịch vụ truyền hình vệ tinh thu trực tiếp hộ gia đình Công suất xạ đẳng hướng tương đương Nghiệp vụ cố định Nghiệp vụ cố định qua vệ tinh Quỹ đạo tầm cao Bộ khuếch đại công suất cao Nghiệp vụ liên vệ tinh Liên minh viễn thông quốc tế Quỹ đạo tầm thấp Phân cực tròn trái Bộ khuyếch đại tập âm thấp Khối tạp âm thấp Bộ chuyển đổi tạp âm thấp Bảng tần số chủ Quỹ đạo tầm trung Nghiệp vụ di động Nghiệp vụ di động qua vệ tinh Phạm Toàn Thắng Luận văn cao häc §TVT 2004-2006 OBP PFD PP RoP RR RRB On-Board Processing Power Flux Density Plenipotentiary Right Hand Circular Polarization Rules of Procedures Radio Regulations Radio Regulation Board TVRO TeleVision Receive Only RHCP TWTA UPC VSAT WRC Travelling Wave Tube Amplifier Up-Link Power Control Very Small Aperture Terminal World Radiocommunication Conference Xử lý trạm Mật độ thông lượng công suất Hội nghị toàn quyền ITU Phân cực tròn phải Quy định thủ tục thực Thể lệ vô tuyến điện Uỷ ban Thể lệ vô tuyến điện Trạm mặt đất thu tín hiệu truyền hình từ vệ tinh Bộ khuếch đại kiểu ống sóng chạy Điều khiển công suất phát lên Trạm có góc mở nhỏ Hội nghị thông tin vô tuyến giới Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 Danh mục bảng Bảng 1.1 Số lượng antenna vệ tinh dùng gia đình 10 Bảng 2.1 Các băng tần thông dụng cho thông tin vƯ tinh 39 B¶ng 2.2 Bảng phân bổ tần số cho thông tin vệ tinh 39 Bảng 2.3 Các tham số vệ tinh băng tần C 40 Bảng 2.4 Các tham số trạm mặt đất băng tần C 41 Bảng 2.5 Các tham số vệ tinh băng tần Ku điển hình 42 Bảng 2.6 Các tham số trạm mặt đất băng tần Ku điển hình 42 Bảng 2.7 Các tham số vệ tinh băng tần Ka điển hình 43 Bảng 2.8 Các tham số trạm mặt đất băng tần Ka điển hình 45 Bảng 3.1 Cung quỹ đạo phối hợp cho hệ thống vệ tinh địa tĩnh 56 Bảng 3.2 Tỷ số bảo vệ C/I đường vào sóng mang FSS 58 Bảng 4.1: Quan hệ đại lỵng (C/N), (C/N ), (C/T), (E b /N ) 61 Bảng 5.1 Bảng hệ số k phụ thuộc tần số 79 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 18T R R R 18T R R R 18T Danh mục hình vẽ Hình 1.1 Vị trí vệ tinh địa tĩnh băng tần C 15 Hình 1.2 Vị trí vệ tinh địa tĩnh băng tần Ku 16 H×nh 2.1 Sơ đồ tuyến thông tin vệ tinh 21 Hình 2.2 Công suất xạ anten phát 23 Hình 2.3 Độ mở hiệu dụng anten thu 25 H×nh 2.4 Suy hao kh«ng gian tù 28 Hình 2.5 Tạp âm nhiệt hƯ thèng thu 31 H×nh 3.1 Mối quan hệ phương pháp phối hợp kỹ thuật cho nghiệp vụ cố định qua vệ tinh 54 Hình 5.1 Sơ đồ thể đường truyền vƯ tinh víi c¸c 73 tham sè đưa vào trình dự đoán suy hao mưa 73 Hình 5.2 Bản đồ lượng mưa (mm/h) bị vượt thời gian 0.01% trung bình năm 76 Hình 5.3 Suy hao đặc trưng khÝ quyÓn 84 Hình 5.4 Bản đồ nhiệt độ trung bình năm bề mặt trái đất (o C) 85 Hình 5.5 Lưu đồ thuật toán tối ưu hóa 92 H×nh 5.6 Giao diƯn chương trình 93 Hình 5.7 Modul tối ưu hóa quỹ công suất ®êng trun - Input 93 H×nh 5.8 Modul tối ưu hóa quỹ công suất đường truyền Quick Result 94 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T 18T P P 18T 18 T 18 T 18T 18T Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 18 T 18 T Mở đầu Cùng với tiến khoa học công nghệ, thông tin vệ tinh đà góp phần đặc biệt quan trọng vào phát triển chung ngành viễn thông toàn giới ®· chøng tá u thÕ vỵt tréi cđa nã nhiều dịch vụ Trong năm gần đây, hệ thống cáp quang đà trở thành phương tiện truyền dẫn ứng dụng băng thông rộng Tuy nhiên, với tính chất đặc thù mình, thông tin vệ tinh tiếp tục khẳng định vai trò thay thế, đặc biệt dịch vụ thông tin quảng bá, di động, an toàn cứu nạn, hàng hải viễn thông cho vùng sâu, vùng xa, hải đảo Chính lý đó, hầu có công nghiệp dịch vụ viễn thông phát triển xây dựng hệ thống vệ tinh viễn thông riêng Ngoài ra, khả phủ sóng cung cấp dịch vụ toàn cầu, hệ thống vệ tinh viễn thông quốc tế khu vực INTELSAT, INMARSAT, INTERSPUTNIC, EUTELSAT đà thiết lập ngày phát triển mạnh mẽ Nhiều nước phát triển đà trọng phát triển hệ thống vệ tinh viễn thông để cung cấp dịch vụ nước tham gia vào thị trường thông tin vệ tinh quốc tế Không bỏ lỡ hội Việt nam nước vậy, với dự án vệ tinh viễn thông mang tên VINASAT Sự phát triển thông tin vệ tinh đà dẫn đến tình trạng khai thác tối đa nguồn tài nguyên tần số quỹ đạo vệ tinh Quỹ đạo vệ tinh địa tĩnh đà trở nên tắc nghẽn, đặc biệt băng tần không quy hoạch Các nước bắt đầu phát triển thông tin vệ tinh khó có khả tìm vị trí quỹ đạo không gây can nhiễu để phóng vệ tinh Các quy trình phối hợp tần số cho vị trí quỹ đạo địa tĩnh ngày chặt chẽ yêu cầu tính toán can nhiễu hệ thống ngày khắt khe phức tạp Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 Do tham gia vào lĩnh vực thông tin vệ tinh muộn nước khác, nên Việt nam phải chịu thách thức tình hình nói Trong luận văn này, tác giả đề cập đến số kết nghiên cứu trình thực nhiệm vụ phối hợp tần số quỹ đạo vệ tinh Việt nam Ngoài vấn đề chung dịch vụ thông tin vệ tinh, quy định Liên minh viễn thông quốc tế (ITU) phối hợp tần số quỹ đạo vệ tinh địa tĩnh, nhiễu hệ thống vệ tinh viễn thông, ứng dụng công nghệ để giảm yêu cầu chống nhiễu hệ thống vệ tinh, tác giả sâu vào việc tính toán tối ưu hóa quỹ công suất đường truyền thông tin vệ tinh để nâng cao khả hoàn thành phối hợp tần số để khai thác hiệu dự án VINASAT Tác giả xin trân trọng cảm ơn hướng dẫn khoa học tận tình PGS TS Trần Văn Cúc; cảm ơn giúp đỡ thầy giáo, cô giáo khoa Điện tử viễn thông, Trường Đại học Bách khoa Hà nội; cảm ơn hỗ trợ đồng nghiệp Phòng Quy hoạch tần số Phối hợp quốc tế, Cục tần số vô tuyến điện trình thực luận văn Hà nội, tháng 10 năm 2006 Phạm Toàn Thắng Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 Chương 1: Tổng quan 1.1 Giới thiệu Thông tin vệ tinh đà trở thành phần thiếu mạng viễn thông toàn cầu Một điều dễ dàng nhìn thấy phổ biến thông tin vệ tinh ngày nhiều gia đình lắp đặt chảo anten parabol để thu tÝn hiƯu trun h×nh trùc tiÕp tõ vƯ tinh, sè liệu cụ thể xem Bảng 1.1 Ngoài ra, lượng lớn lưu lượng thoại liệu truyền tải qua thông tin vệ tinh Trong nhiều trường hợp, thông tin vệ tinh giải pháp hiệu nhÊt: vÝ dơ nh viƯc ¸p dơng hƯ thèng VSAT IP để triển khai mạng điện thoại đến hầu hết xà vùng sâu vùng xa nước ta; việc sử dụng thông tin vệ tinh làm hệ thống dự phòng cho hệ thống cáp quang biển có cố; hay tình tìm kiếm cứu nạn có thiên tai lũ lụt, tai nạn máy bay, tàu biển Xu hướng phát triển thông tin vệ tinh dựa vào đặc tính riêng biệt mà phương tiện truyền thông khác không có: vùng phục vụ rộng lớn, khả đa truy nhập tính linh họat việc thiết lập đường truyền Thông tin vệ tinh cung cấp dịch vụ vùng rộng lớn bề mặt trái đất, kết nối lúc nhiều khách hàng, cụm dân cư vùng địa lý khác không kể vùng núi cao, rừng sâu hay hải đảo, điều mà thực phương tiện truyền thông khác Khả đa truy cập thể việc quảng bá phân tán file liệu từ điểm tới nhiều điểm giới thời điểm Khả triển khai nhanh linh hoạt việc thiết lập đường truyền việc lắp đặt trạm mặt đất thu phát Ngoài ra, khách hàng tự di chuyển trạm thu phát Phạm Toàn Thắng Luận văn cao häc §TVT 2004-2006 108 txtWMonth.Text = Math.Round(MonthAvail, 2) 'hien thi tai phan Quick Result, khong lien quan gi den tinh toan ' Calculate rainfade, co lien quan den thong so Availability RainFadeUp = Rainfade(txtUpRainFall.Text, txtUpFreq.Text, txtESUpLat.Text, txtESUpLong.Text, txtUpAltitude.Text, txtSatLong.Text, txtAvail.Text) RainFadeDown = Rainfade(txtDownRainFall.Text, txtDownFreq.Text, txtESDownLat.Text, txtESDownLong.Text, txtDownAltitude.Text, txtSatLong.Text, txtAvail.Text) ' Calculate G/T of receiving E/S under clear sky TLNA = (Math.Exp(txtLNBF.Text / 10 * Math.Log(10)) - 1) * 290 l = Math.Exp(txtDownCLoss.Text / 10 * Math.Log(10)) Ttotal = txtAntTemp.Text / l + 290 * (1 - / l) + TLNA ' noise temperature at the input of LNA GTES = GesDown - txtDownCLoss.Text - 10 * Math.Log(Ttotal) / Math.Log(10) ' including couplink loss ' Calculate reduction of G/T of receiving E/S under rain condition Train = 290 * (1 - / Math.Exp(RainFadeDown / 10 * Math.Log(10))) TtotalR = Ttotal + Train GTreduction = 10 * Math.Log(TtotalR) / Math.Log(10) - 10 * Math.Log(Ttotal) / Math.Log(10) ' Power control here If chkPowerControl.Checked = True Then If txtPowerUpControl.Text >= RainFadeUp Then Pwcompensation = RainFadeUp Else Pwcompensation = txtPowerUpControl.Text End If Else Pwcompensation = ' khong co pw control End If ' Xac dinh bang tan If (txtUpFreq.Text > 4) And (txtUpFreq.Text < 7.5) Then Bands = " C" End If If (txtUpFreq.Text > 12.9) And (txtUpFreq.Text < 17.8) Then Bands = " Ku" End If If (txtUpFreq.Text > 24) And (txtUpFreq.Text < 35) Then Bands = " Ka" End If If (txtUpFreq.Text > 7.7) And (txtUpFreq.Text < 9) Then Bands = " X" End If 'Calculate values of C/No at saturation point and the clear sky condition UpEIRPSatu = Math.Round(txtSatSFD.Text + 10 * Math.Log(4 * Math.PI * (Dup * 1000) ^ 2) / Math.Log(10) + TotalAtLossUp, 2) UpCNoSatu = Math.Round(UpEIRPSatu - UpPathLoss + txtSatGT.Text + 228.6 - txtUpMis.Text - TotalAtLossUp, 2) DownCNoSatu = Math.Round(txtSatEIRP.Text - DownPathLoss + GTES + 228.6 - txtDownMis.Text TotalAtLossDown, 2) TotalCNoSatu = Math.Round(10 * Math.Log10(1 / (1 / Math.Exp(UpCNoSatu / 10 * Math.Log(10)) + / Math.Exp(DownCNoSatu / 10 * Math.Log(10)))), 2) ' Calculate bandwidth require If optBPSK.Checked = True Then M=2 ModulationMode = "BPSK" Else If optMPSK.Checked = True Then M = txtM.Text ModulationMode = txtM.Text + "-PSK" Else If optMQAM.Checked = True Then M = txtM.Text ModulationMode = txtM.Text + "-QAM" End If Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 109 End If End If ' determine the index of modulation scheme Inforate = txtinforate.Text * (1 + txtOverhead.Text) ' Mbps transrate = Inforate / txtFEC.Text Symbolrate = Math.Round(transrate / (M / 2), 6) OccBW = Symbolrate * txtSpreadfactor.Text ' MHz NoiseBW = 10 * Math.Log10(OccBW * 1000000) ' dBHz MinimumBW = Symbolrate * txtCarrierfactor.Text ' MHz 'minimumbw = OccBW * txtCarrierfactor ' MHz dung hay la bieu thuc tren dung n = MinimumBW / txtBWstep.Text n = Math.Round(n, 0) If n * txtBWstep.Text >= MinimumBW Then n=n Else n=n+1 End If AlloBW = n * txtBWstep.Text 'MHz NcarrierBW = txtSatBW.Text / AlloBW ' number of carriers per transponder by BW usage ' Optimizer ' Uplink calculation ' set up IPBO per carrier deltaP = 0.01 ' dB IPBOC = 100 ' later it should plus transponder IPBO to become real IPBOC OPBOC = 100 ' later it should plus transponder OPBO to become real IPBOC ' >> vong lap tinh truoc kiem tra sau can duoc chen vao day 'Dim CNoUpCS As Double, CNoUpRUp As Double, CNoUpRDown As Double 'Dim CIoCoUpCS As Double, CIoCoUpRUp As Double, CIoCoUPRDown As Double 'Dim CIoAdjUpCS As Double, CIoAdjUpRUP As Double, CIoAdjUPRDown As Double 'Dim CNoIoUpCS As Double, CNoIoUpRUp As Double, CNoIoUpRDown As Double 'Dim CNIUpCS As Double, CNIUpRUp As Double, CNIUpRDown As Double 'Dim EbNoIoUpCS As Double, EbNoIoUpRUp As Double, EbNoIoUpRDown As Double Do ' begin of the loop UpEIRP = UpEIRPSatu - txtSatIPBO.Text - IPBOC ' C/No (Thermal) CNoUpCS = UpEIRP - txtUpMis.Text - UpPathLoss - TotalAtLossUp + txtSatGT.Text + 228.6 'Clear Sky CNoUpRUp = UpEIRP - txtUpMis.Text - UpPathLoss - TotalAtLossUp + txtSatGT.Text + 228.6 - RainFadeUp + Pwcompensation 'Rain at Uplink CNoUpRDown = UpEIRP - txtUpMis.Text - UpPathLoss - TotalAtLossUp + txtSatGT.Text + 228.6 'Rain at Downlink 'C/Io Co-channel Interference CIoCoUpCS = txtUpCoCI.Text + NoiseBW CIoCoUpRUp = txtUpCoCI.Text + NoiseBW - RainFadeUp + Pwcompensation CIoCoUpRDown = txtUpCoCI.Text + NoiseBW ' C/Io Adjacent satellite interference CIoAdjUpCS = txtUpAdjSCI.Text + NoiseBW CIoAdjUpRUp = txtUpAdjSCI.Text + NoiseBW - RainFadeUp + Pwcompensation CIoAdjUpRDown = txtUpAdjSCI.Text + NoiseBW ' C/(No+Io) : for all kinds of noise CNoIoUpCS = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoUpCS / 10 * Math.Log(10)) + / Math.Exp(CIoCoUpCS / 10 * Math.Log(10)) + / Math.Exp(CIoAdjUpCS / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoIoUpRUp = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoUpRUp / 10 * Math.Log(10)) + / Math.Exp(CIoCoUpRUp / 10 * Math.Log(10)) + / Math.Exp(CIoAdjUpRUp / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoIoUpRDown = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoUpRDown / 10 * Math.Log(10)) + / Math.Exp(CIoCoUpRDown / 10 * Math.Log(10)) + / Math.Exp(CIoAdjUpRDown / 10 * Math.Log(10)))) / Math.Log(10), 2) ' C/(N+I) for all kinds of noise CNIUpCS = CNoIoUpCS - NoiseBW CNIUpRUp = CNoIoUpRUp - NoiseBW CNIUpRDown = CNoIoUpRDown - NoiseBW ' Eb/(No+Io) for all kinds of noise EbNoIoUpCS = CNoIoUpCS - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) EbNoIoUpRUp = CNoIoUpRUp - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) Phạm Toàn Thắng Luận văn cao học ĐTVT 2004-2006 110 EbNoIoUpRDown = CNoIoUpRDown - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) ' Downlink Calculation 'Dim CNoDnCS As Double, CNoDnRUp As Double, CNoDnRDown As Double 'Dim CIoCoDnCS As Double, CIoCoDnRUp As Double, CIoCoDnRDown As Double 'Dim CIoAdjDnCS As Double, CIoAdjDnRUP As Double, CIoAdjDnRDown As Double 'Dim CNoIoDnCS As Double, CNoIoDnRUp As Double, CNoIoDnRDown As Double 'Dim CNIDnCS As Double, CNIDnRUp As Double, CNIDnRDown As Double 'Dim EbNoIoDnCS As Double, EbNoIoDnRUp As Double, EbNoIoDnRDown As Double ' Calculate C/No (thermal) CNoDnCS = txtSatEIRP.Text - txtSatOPBO.Text - OPBOC - txtDownMis.Text - DownPathLoss - TotalAtLossDown + GTES + 228.6 CNoDnRUp = txtSatEIRP.Text - txtSatOPBO.Text - OPBOC - RainFadeUp - txtDownMis.Text - DownPathLoss TotalAtLossDown + GTES + 228.6 + Pwcompensation CNoDnRDown = txtSatEIRP.Text - txtSatOPBO.Text - OPBOC - txtDownMis.Text - DownPathLoss TotalAtLossDown - RainFadeDown - GTreduction + GTES + 228.6 ' Calculate C/Io Co-channel Interference CIoCoDnCS = txtDownCoCI.Text + NoiseBW CIoCoDnRUp = txtDownCoCI.Text + NoiseBW - RainFadeUp + Pwcompensation CIoCoDnRDown = txtDownCoCI.Text + NoiseBW ' Calculate C/Io Adjacent satellite interference CIoAdjDnCS = txtDownAdjSCI.Text + NoiseBW CIoAdjDnRUP = txtDownAdjSCI.Text + NoiseBW - RainFadeUp + Pwcompensation CIoAdjDnRDown = txtDownAdjSCI.Text + NoiseBW ' Calculate C/(No+Io) CNoIoDnCS = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoDnCS / 10 * Math.Log(10)) + / Math.Exp(CIoCoDnCS / 10 * Math.Log(10)) + / Math.Exp(CIoAdjDnCS / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoIoDnRUp = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoDnRUp / 10 * Math.Log(10)) + / Math.Exp(CIoCoDnRUp / 10 * Math.Log(10)) + / Math.Exp(CIoAdjDnRUP / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoIoDnRDown = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoDnRDown / 10 * Math.Log(10)) + / Math.Exp(CIoCoDnRDown / 10 * Math.Log(10)) + / Math.Exp(CIoAdjDnRDown / 10 * Math.Log(10)))) / Math.Log(10), 2) 'Calculate C/(N+I) CNIDnCS = CNoIoDnCS - NoiseBW CNIDnRUp = CNoIoDnRUp - NoiseBW CNIDnRDown = CNoIoDnRDown - NoiseBW ' Calculate Eb/(No+Io) EbNoIoDnCS = CNoIoDnCS - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) EbNoIoDnRUp = CNoIoDnRUp - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) EbNoIoDnRDown = CNoIoDnRDown - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) ' calculate total link to optimize ' Calculate C/No total CNoTCS = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoUpCS / 10 * Math.Log(10)) + / Math.Exp(CNoDnCS / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoTRUp = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoUpRUp / 10 * Math.Log(10)) + / Math.Exp(CNoDnRUp / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoTRDown = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoUpRDown / 10 * Math.Log(10)) + / Math.Exp(CNoDnRDown / 10 * Math.Log(10)))) / Math.Log(10), 2) ' Calculate C/Io Transponder Intermod Interference CIoIntCS = txtSatIM_CI.Text + NoiseBW CIoIntRUp = txtSatIM_CI.Text + NoiseBW - RainFadeUp + Pwcompensation CIoIntRDown = txtSatIM_CI.Text + NoiseBW ' Calculate C/(No+Io) CNoIoTCS = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoIoUpCS / 10 * Math.Log(10)) + / Math.Exp(CNoIoDnCS / 10 * Math.Log(10)) + / Math.Exp(CIoIntCS / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoIoTRUp = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoIoUpRUp / 10 * Math.Log(10)) + / Math.Exp(CNoIoDnRUp / 10 * Math.Log(10)) + / Math.Exp(CIoIntRUp / 10 * Math.Log(10)))) / Math.Log(10), 2) CNoIoTRDown = Math.Round(10 * Math.Log(1 / (1 / Math.Exp(CNoIoUpRDown / 10 * Math.Log(10)) + / Math.Exp(CNoIoDnRDown / 10 * Math.Log(10)) + / Math.Exp(CIoIntRDown / 10 * Math.Log(10)))) / Math.Log(10), 2) 'Calculate C/(N+I) Phạm Toàn Thắng Luận văn cao häc §TVT 2004-2006 111 CNITCS = CNoIoTCS - NoiseBW CNITRUp = CNoIoTRUp - NoiseBW CNITRDown = CNoIoTRDown - NoiseBW ' Calculate Eb/(No+Io) EbNoIoTCS = CNoIoTCS - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) EbNoIoTRUp = CNoIoTRUp - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) EbNoIoTRDown = CNoIoTRDown - 10 * Math.Log(Inforate * 1000000) / Math.Log(10) ' Calculate Net Eb/(No+Io) NEbNoIoTCS = EbNoIoTCS - txtMargin.Text NEbNoIoTRUp = EbNoIoTRUp - txtMargin.Text NEbNoIoTRDown = EbNoIoTRDown - txtMargin.Text 'Cac bien kiem tra vong lap EmarginCS = NEbNoIoTCS - txtEbNo.Text 'kiem tra Eb/No tinh toan so voi yeu cau dich vu EmarginRup = NEbNoIoTRUp - txtEbNo.Text 'kiem tra Eb/No tinh toan so voi yeu cau dich vu EmarginRDown = NEbNoIoTRDown - txtEbNo.Text 'kiem tra Eb/No tinh toan so voi yeu cau dich vu If (IPBOC