NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI

NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI

BỘ GIAO THÔNG VẬN TẢI BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC HÀNG HẢI VIỆT NAM LÊ THỊ TOÁN ĐỒ ÁN TỐT NGHIỆP ĐẠI HỌC NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI HẢI PHÒNG - 2015 BỘ GIAO THÔNG VẬN TẢI BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC HÀNG HẢI VIỆT NAM LÊ THỊ TOÁN ĐỒ ÁN TỐT NGHIỆP ĐẠI HỌC NGHIÊN CỨU ĐIỀU KHIỂN NGHẼN MẠNG TRONG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI NGÀNH: KỸ THUẬT ĐIỆN TỬ TRUYỀN THÔNG; MÃ SỐ: D52027 CHUYÊN NGÀNH: ĐIỆN TỬ - VIỄN THÔNG Người hướng dẫn: ThS Trương Thanh Bình HẢI PHÒNG – 2015 LỜI CẢM ƠN Để hoàn thành đồ án em nhận giúp đỡ nhiệt tình Ths Trương Thanh Bình bạn lớp DTV52 – DH1 Đầu tiên , em xin gửi lời biết ơn sâu sắc đến Ths Trương Thanh Bình , người dành nhiều thời gian hướng dẫn giúp em hoàn thành đồ án Em xin cảm ơn Ban giám hiệu trường Đại học Hàng Hải thầy cô khoa Điện – Điện Tử tạo điều kiện cho em học tập hoàn thành khóa học Em xin cảm ơn bạn lớp DTV52 – DH1 giúp đỡ em nhiều trình học tập Do kiến thức em hạn chế nên đồ án em nhiều thiếu sót Em mong nhận ý kiến góp ý thầy cô bạn để đồ án em hoàn thiện LỜI CAM ĐOAN Em xin cam đoan đồ án tốt nghiệp em thực Nếu có vấn đề xảy em xin chịu hoàn toàn trách nhiệm Hải Phòng , ngày 18 tháng 11 năm 2015 Sinh viên thực Lê Thị Toán MỤC LỤC DANH MỤC CÁC CHỮ VIẾT TẮT VÀ KÝ HIỆU ATM AON BA CPG DWDM DR FDL HDTV IBT IP JET JIT NACK OCS OBS OADM O/E/O OLT OPS OXC OSC QOS RFD RAM SONET SCU TAG TAW WADM Asynchronous Transfer Mode - Chế độ truyền không đồng All Optical Network - Mạng toàn quang Burst Assembler - Sắp xếp chùm Control Packet Generalor - Máy phát gói điều khiển Dense WDM - WDM mật độ cao Delay Reservation - Đăng ký trễ Fiber Delay Line - Đường dây trễ quang High Difinition Television - Truyền hình độ phân giải cao In-Band-Terminater - Phương thức IBT Internet Protocol - Giao thức Internet Just Enough Time - Giao thức JET Just In Time - Giao thức JIT Negative Acknowledget - Gói xác nhận Optical Circuit Switch - Chuyển mạch kênh quang Optical Burst Switch - Chuyển mạch chùm quang Optical Add-Drop Multiplexer - Bộ ghép kênh xen – rẽ quang Optical/Electricity/Optical - Chuyển đổi quang/điện/quang Optical Line Termination - Thiết bị đầu cuối đường dây quang Optical Packet Switch - Chuyển mạch gói quang Optical Cross-Connect -Kết nối chéo quang Optical Supervisory Channels - Kênh giám sát quang Quality of service - Chất lượng dịch vụ Reserve-a-Fixed-Delay - Đăng ký trễ cố định Random Access Memory Bộ nhớ truy cập tạm thời Synchronous Optical Network - Mạng quang đồng Switch Control Unit - Khối điều khiển chuyển mạch Tell And Go - Phương thức TAG Tell And Wait - Phương thức TAW Wavelenghth Add-Drop Multiplexer - Ghép kênh thêm bớt theo bước sóng DANH MỤC CÁC HÌNH VẼ Số hình 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 Tên hình Trang Hình 1.1 Sự phát triển mạng thông tin quang Hình 1.2 Sơ đồ mạng quang định tuyến bước sóng Hình 1.3 Sơ đồ khối chuyển mạch gói quang Hình 2.1.Sơ đồ nguyên lý chuyển mạch chùm quang Hình 2.2 Sơ đồ kiến trúc mạng chuyển mạch chùm quang 10 OBS Hình 2.3 Bộ định tuyến lõi 11 Hình 2.4 Bộ định tuyến biên 12 Hình 2.5 Kiến trúc OBS dạng mắc lưới 13 Hình 2.6 Kiến trúc OBS dạng vòng nút 14 Hình 2.7 Kiến trúc nút OBS 15 Hình 2.8 Sơ đồ thiết bị đầu cuối đường dây OLT 16 Hình 2.9 Vai trò xen / rẽ quang 17 Hình 2.10 Các chức OXC 18 Hình 2.11 sơ đồ khuếch đại quang 19 Hình 2.12 Phương thức TAG 20 Hình 2.13 Phương thức TAW 21 Hình 2.14 Quá trình hoạt động giao thức JIT 22 Hình 2.15 Quá trinhg hoạt động giao thức JET 23 Hình 2.16 Giải tranh chấp việc sử dụng đệm 24 FDL Hình 2.17 Giải tranh chấp việc biến đổi bước 25 sóng Hình 2.18 Giải tranh chấp việc làm lệch hướng 27 Hình 3.1 Mạng OBS với kĩ thuật làm lệch hướng Hình 3.2 Lược đồ thuật toán Hình 3.3 (a) Trường ưu tiên gói điều khiển; (b) Lớp ưu tiên đầu vào Hình 3.4 Ảnh hưởng định tuyến lệch hướng Hình 3.5.Phương pháp định tuyến lệch hướng Hình 3.6 cấu trúc kiểm tra Wh = Hình 3.7 Lưu lượng đầu vào nút nguồn Hình 3.8 Sơ đồ mạng ảo Hình 3.9 Kết so sánh có không sử dụng FDL Hình 3.10 Kết sử dụng FDL khác Hình 3.11 Kết so sánh thay đổi khoảng thời gian tạo chùm Hình 3.12 Kết sử dụng FDL = 60 ms thay đổi thời 29 30 31 33 34 37 37 39 39 40 41 42 gian tạo chùm MỞ ĐẦU Hiện nay, kinh tế toàn cầu phát triển mạnh mẽ kéo theo nhu cầu sử dụng mạng thông tin ngày cao người dùng Do đó, đòi hỏi phải có hệ thống truyền lượng liệu lớn với tốc độ truyền cao băng tần rộng Một giải pháp đưa việc sử dụng sợi quang kết hợp với kỹ thuật ghép kênh phân chia theo bước sóng WDA để cung cấp tốc độ truyền liệu nhanh nhằm đáp ứng tốc độ phát triển vượt bậc lưu lượng mạng Mục tiêu tương lai xây dựng mạng mà liệu cần gửi gửi hoàn toàn miền quang gói điều khiển gửi miền điện Vì thế, công nghệ chuyển mạch quang đề xuất công nghệ chuyển mạch gói quang OPS, chuyển mạch kênh quang WRN công nghệ chuyển mạch chùm quang OBS Trong đó, chuyển mạch chùm quang OBS cân mạng chuyển mạch gói quang OPS mạng chuyển mạch kênh WRN Chuyển mạch chùm quang đáp ứng yêu cầu truyền lượng liệu lớn với tốc độ truyền cao, công nghệ hứa hẹn tương lai Trong đồ án này, em tập trung tìm hiểu mạng chuyển mạch quang OBS, sâu vào tìm hiểu mô phương pháp làm lệch hướng để giải nghẽn mạng OBS Mục đích tìm cách giải tốt xảy tranh chấp liệu cho xác suất chùm suy hao nhỏ Nội dung đồ án gồm chương : CHƯƠNG I : GIỚI THIỆU VỀ MẠNG THÔNG TIN QUANG CHƯƠNG II : MẠNG CHUYỂN MẠCH CHÙM QUANG OBS CHƯƠNG III : GIẢI PHÁP ĐIỀU KHIỂN NGHẼN TRONG MẠNG OBS BẰNG PHƯƠNG PHÁP LÀM LỆCH HƯỚNG ĐI Do thời gian làm đồ án có hạn kiến thức hạn chế nên đồ án em nhiều thiếu sót cần bổ sung Do em mong thầy cô thông cảm hướng dẫn cho em, bạn đọc quan tâm đến vấn đề đóng góp ý kiến để đồ án em hoàn thiện Em xin gửi lời cảm ơn chân thành đến thầy giáo hướng dẫn - Ths Trương Thanh Bình, người tận tình bảo, bổ sung kiến thức cho em, giúp em hoàn thành tốt đồ án theo Em xin chân thành cảm ơn thầy cô môn Điện Tử - Viễn Thông truyền đạt kiến thức cho em năm qua giúp đỡ em thời gian làm đồ án Hải Phòng, ngày 18 tháng 11 năm 2015 Sinh viên thực Lê Thị Toán 10 i=7:12; Ti=B.*C./(lamda.*L.*(i-n0).*(i-n0+1)); Ti1=B1.*C./(lamda.*L.*(i-n0).*(i-n0+1)); Ti2=B2.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a=3.8462; b=1.3514; n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); Pl=p.*exp(-i.*a.*Ti) Pl1=p.*exp(-i.*a.*Ti1) Pl2=p.*exp(-i.*a.*Ti2) %VE TONG HOP DO THI w=plot(i,Pl,'b*-',i,Pl1,'ko ',i,Pl2,'r+:'); grid on ; title('so sanh ket qua'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI DO THI TONG HOP TREN HE TOA DO % - Executes on button press in NUT3 function NUT3_Callback(hObject, eventdata, handles) % hObject handle to NUT3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT3 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',1); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); 62 set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO C = 6; L = 1; lamda = 1000; %KHI SU DUNG FDL = 60 ms B2=60; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti2=B2.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a=3.8462; b=1.3514; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); Pl2=p.*exp(-i.*a.*Ti2); %VE HAM Pl THEO i m=plot(i,Pl2,'r+:'); grid on ; title('ket qua FDL=60ms'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); 63 %HIEN THI TREN HE TOA DO % - Executes on button press in NUT5 function NUT5_Callback(hObject, eventdata, handles) % hObject handle to NUT5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT5 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',1); set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO C = 6; L = 1; lamda = 1000; %KHI KHONG SU DUNG FDL B5=0; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti5=B5.*C./(lamda.*L.*(i-n0).*(i-n0+1)); 64 a5=3.5714; b5=1.3889; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b5./(a5+b5)).^i).*((b5./ (a5+b5)).^(n-i)); Pl5=p.*exp(-i.*a5.*Ti5) %VE HAM Pl THEO i l=plot(i,Pl5,'k*:'); grid on ; title('khong FDL voi 1/a=320 va 1/b=680'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT7 function NUT6_Callback(hObject, eventdata, handles) % hObject handle to NUT7 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT7 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); set(handles.NUT6,'Value',1); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); 65 set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO C = 6; L = 1; lamda = 1000; %KHI KHONG SU DUNG FDL B6=0; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti6=B6.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a=3.8462; b=1.3514; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); Pl6=p.*exp(-i.*a.*Ti6) %VE HAM Pl THEO i q=plot(i,Pl6,'r+:'); grid on ; title('khong FDL voi 1/a=300 va 1/b=700'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT7 function NUT7_Callback(hObject, eventdata, handles) % hObject handle to NUT7 (see GCBO) 66 % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT7 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',1); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO C = 6; L = 1; lamda = 1000; %KHI KHONG SU DUNG FDL B7=0; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti7=B7.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a7=4.1667; b7=1.3158; %CHON n THOA DIEU KIEN n>i n=25; 67 p=factorial(n)./(factorial(i).*factorial(n-i)).*((b7./(a7+b7)).^i).*((b7./ (a7+b7)).^(n-i)); Pl7=p.*exp(-i.*a7.*Ti7); %VE HAM Pl THEO i x=plot(i,Pl7,'bo:'); grid on ; title('khong FDL voi 1/a=280 va 1/b=720'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT8 function NUT8_Callback(hObject, eventdata, handles) % hObject handle to NUT8 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT8 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',1); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO 68 C = 6; L = 1; lamda = 1000; B6=0; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti6=B6.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a5=3.1250; b5=1.4706; a=3.333; b=1.4268; a7=3.5714; b7=1.3889; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); p5=factorial(n)./(factorial(i).*factorial(n-i)).*((b5./(a5+b5)).^i).*((b5./ (a5+b5)).^(n-i)); p7=factorial(n)./(factorial(i).*factorial(n-i)).*((b7./(a7+b7)).^i).*((b7./ (a7+b7)).^(n-i)); Pl6=p.*exp(-i.*a.*Ti6); Pl5=p5.*exp(-i.*a5.*Ti6); Pl7=p7.*exp(-i.*a7.*Ti6); %VE HAM Pl THEO i l=plot(i,Pl6,'r+:',i,Pl5,'k*:',i,Pl7,'bo:'); grid on ; title('so sanh ket qua'); xlabel('luu luong tai'); 69 ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT9 function NUT9_Callback(hObject, eventdata, handles) % hObject handle to NUT9 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT9 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',1); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO C = 6; L = 1; lamda = 1000; %KHI SU DUNG FDL = 60 ms B6=60; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; 70 Ti6=B6.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a5=3.1250; b5=1.4706; a=3.333; b=1.4268; a7=3.5714; b7=1.3889; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); Pl6=p.*exp(-i.*a.*Ti6); %VE HAM Pl THEO i q=plot(i,Pl6,'k*:'); grid on ; title('FDL = 60 voi 1/a=280 va 1/b=720'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT10 function NUT10_Callback(hObject, eventdata, handles) % hObject handle to NUT10 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT10 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); 71 set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',1); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',0); %THONG SO C = 6; L = 1; lamda = 1000; %KHI SU DUNG FDL = 60 ms B7=60; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti7=B7.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a7=4.1667; b7=1.3158; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b7./(a7+b7)).^i).*((b7./ (a7+b7)).^(n-i)); Pl7=p.*exp(-i.*a7.*Ti7); %VE HAM Pl THEO i x=plot(i,Pl7,'r+:'); grid on ; title('FDL=60 voi 1/a=260 va 1/b=740'); xlabel('luu luong tai'); 72 ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT11 function NUT11_Callback(hObject, eventdata, handles) % hObject handle to NUT11 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT11 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',1); set(handles.NUT12,'Value',0); C = 6; L = 1; lamda = 1000; %KHI SU DUNG FDL = 60 ms B6=60; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti6=B6.*C./(lamda.*L.*(i-n0).*(i-n0+1)); 73 a5=3.1250; b5=1.4706; a=3.333; b=1.4268; a7=3.5714; b7=1.3889; %CHON n THOA DIEU KIEN n>i n=25; p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); p5=factorial(n)./(factorial(i).*factorial(n-i)).*((b5./(a5+b5)).^i).*((b5./ (a5+b5)).^(n-i)); p7=factorial(n)./(factorial(i).*factorial(n-i)).*((b7./(a7+b7)).^i).*((b7./ (a7+b7)).^(n-i)); Pl6=p.*exp(-i.*a.*Ti6); Pl5=p5.*exp(-i.*a5.*Ti6); Pl7=p7.*exp(-i.*a7.*Ti6); %VE HAM Pl THEO i x=plot(i,Pl7,'bo:'); grid on ; title('so sanh ket qua'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO % - Executes on button press in NUT12 function NUT12_Callback(hObject, eventdata, handles) % hObject handle to NUT12 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hint: get(hObject,'Value') returns toggle state of NUT12 74 set(handles.NUT1,'Value',0); set(handles.NUT2,'Value',0); set(handles.NUT3,'Value',0); set(handles.NUT4,'Value',0); set(handles.NUT5,'Value',0); set(handles.NUT6,'Value',0); set(handles.NUT7,'Value',0); set(handles.NUT8,'Value',0); set(handles.NUT9,'Value',0); set(handles.NUT10,'Value',0); set(handles.NUT11,'Value',0); set(handles.NUT12,'Value',1); C = 6; L = 1; lamda = 1000; %KHI SU DUNG FDL = 60 ms B6=60; n0 = 10^3.*C./(lamda.*L); %CHON i THOA DIEU KIEN i>n0+1 i=7:12; Ti6=B6.*C./(lamda.*L.*(i-n0).*(i-n0+1)); a5=3.1250; b5=1.4706; a=3.333; b=1.4268; a7=3.5714; b7=1.3889; %CHON n THOA DIEU KIEN n>i n=25; 75 p=factorial(n)./(factorial(i).*factorial(n-i)).*((b./(a+b)).^i).*((b./(a+b)).^(n-i)); p5=factorial(n)./(factorial(i).*factorial(n-i)).*((b5./(a5+b5)).^i).*((b5./ (a5+b5)).^(n-i)); p7=factorial(n)./(factorial(i).*factorial(n-i)).*((b7./(a7+b7)).^i).*((b7./ (a7+b7)).^(n-i)); Pl6=p.*exp(-i.*a.*Ti6); Pl5=p5.*exp(-i.*a5.*Ti6); Pl7=p7.*exp(-i.*a7.*Ti6); %VE HAM Pl THEO i q=plot(i,Pl6,'r+:',i,Pl5,'k*:',i,Pl7,'bo:'); grid on ; title('so sanh ket qua'); xlabel('luu luong tai'); ylabel('xac suat chum suy hao'); %HIEN THI TREN HE TOA DO 76