Nghiên cứu quy hoạch phủ sóng LTE cho khu vực Thái Nguyên

QU TR ỜNG I HỌ N NG NGHỆ NGUYỄN THỊ HINH NGHIÊN ỨU QUY HO H PHỦ SÓNG LTE HO KHU VỰ THÁI NGUYÊN LU N V N TH NG NGHỆ IỆN T H N I – 2014 S - VIỄN TH NG QU TR ỜNG I HỌ N NG NGHỆ ************** NGUYỄN THỊ HINH NGHIÊN ỨU QUY HO H PHỦ SĨNG LTE HO KHU VỰ THÁI NGUN Ngành: ơng nghệ iện Tử - Viễn Thông huyên ngành: Kỹ thuật iện tử Mã số: 60 52 02 03 LU N V N TH S NG NGHỆ IỆN T - VIỄN TH NG Người hướng dẫn khoa học: PGS.TS Nguyễn Quốc Tuấn H N I – 2014 LỜI Tôi xin m o n nh ng g vi t không s o h p, nh ng k t qu ông ố t luận v n M O N t o y hồn tồn h nh thống mơ ph ng tài liệu ngu n gố i luận v n h i m i h nh th tr h thông tin sử t ng ng n r ràng Tơi xin hồn tồn hịu tr h nhiệm n u u hiệu s o h p k t qu t tài liệu kh H N n y th n năm TÁ GIẢ NGUYỄN THỊ HINH i LỜI ẢM Trong th i gi n h s gi p tập, nghiên r t tận t nh hu Viễn thông Tr ng ih ng i ã tận t nh h ng sỹ suốt th i gi n v m ng xin tử - Viễn thông, nh ng ng h n ng l nhận ih quố gi àN i qu qu th y ô, nh hị nh ng g p i ã ng h n t i kho kịp th i iện h, gi p u luận v n ng xin ày t l ng i t n s u sắ m ãn l iện tử – m hoàn thành luận v n tốt nghiệp th ông nghệ ã ng tr nh t o Th M th y ô gi o kho m n s u sắ t i P S TS Nguyễn Quố Tu n, m suốt qu tr nh nghiên Ngoài r , m o n gi p m n ih u hoàn thiện luận v n, m ã nhận ông nghệ, u tiên m xin gửi l i N m suốt qu tr nh h sỹ t i kho n gi n , tập hoàn thành iện tử – Viễn thơng, ố gắng hồn thiện luận v n nh, ih ng t t ông nghệ nhiệt t nh m nh, nhiên tr nh kh i nh ng thi u s t, r t mong nh ng ng g p qu m xin h n thành u qu th y ô n m n H N n y 10 t n năm HỌ VIÊN NGUYỄN THỊ HINH ii NH M K hiệu K HIỆU V TỪ VI T T T T Vi t Tắt One Generation Cellular 1G Second Generation Cellular 2G Third Generation Cellular 3G Four Generation Cellular 4G 3GPP ACK BCCH BCH BW CDMA CP DL Third Generation Patnership Project knowl g m nt T n hiệu x nhận Broadcast Control Channel Broadcast Channel Band Width Code Division Multiple Access Cycle Prefix Downlink Enhance Data rates for GSM EDGE Evolution EUTRAN Evolved UMTS Terrestrial Radio Access eNodeB Enhance NodeB EPC Evolved Packet Core FDMA Frequency Division Multiple Access FDD FrequencyDivision Duplexing FEC Forward Error Correction GSM Global System for Mobile GERAN GSM/EDGE Radio Access Network GPRS General Packet Radio Service GI Guard Interval HDTV High Definition Television HSOPA High Speed OFDM Packet Access HSPA High Speed Packet Access HSS Home Subscriber Server PAP Peak-to-Average Power Ratio P2P Point to Point PDSCH Physical Downlink Shared Channel iii Nghĩ ệ thống thông tin i ng th hệ th nh t ệ thống thông tin i ng th hệ th h i ệ thống thông tin i ng th hệ th ệ thống thông tin i ng th hệ th t D n h p t th hệ T n hiệu x nhận Kênh iều khiển qu ng Kênh qu ng B ng thông truy ập ph n hi th o mã Tiền tố l p ng xuống Tố liệu t ng ng ho m ng SM i ti n M ng truy nhập vô n i ti n No B i ti n M ng l i g i truy ập ph n hi th o t n số h p kênh ph n hi th o t n số Sử l i h i ti p ệ thống i ng toàn u M ng truy nhập vô n GSM/EDGE Dị h v g i vô n thông ng Kho ng o vệ Tivi ph n gi i o Truy ập g i OFDM tố o Truy nhập g i tố o Qu n l thuê o Tỷ số ông su t ỉnh ông su t trung nh iểm n iểm Kênh vật l hi sẻ ng xuống PUCCH PBCH PCCH PCH QoS RLC Physical Uplink Control Channel Physical Broadcast Channel Paging Control Channel Paging Channel Quality of Services Radio Link Control RRC RB RE RSRP RSRQ Radio Resource Control Resource Block Resource Element Reference Signal Receive Power Receive Power Reference Signal Receive Quality RS Reference Signal SDR Software - Defined Radio SNR Signal to Noise Ratio SCFDMA Single Carrier Frequency Division multiple Access SMS Short Message Service SUMIMO Single User Multi Input Multi Output TDMA Time Division Multiple Access TDMA Time Division Multiple Access TTI Time Transmit Interval TDD Time Division Duplexing TPC Transmit Power Command UMB Ultra Mobile Broadband UL Uplink UTRAN UTMS Terrestrial Radio Access Networks UTMS Universal Telecommunication Mobile System UE User Equipment VHE Virtual Home Environment Voice IP VoIP Voice IP WCDMA Wideband Code Division Multiple Access WAP Wireless Applicaion protocol iv Kênh vật l iều khiển ng lên Kênh vật l qu ng Kênh iều khiển tin nhắn Kênh tin nhắn h t l ng ị h v iều khiển k t nối vô n iều khiển tài nguyên vô n Khối tài nguyên Thành ph n tài nguyên ông su t thu t n hiệu th m kh o h t l ng thu t n hiệu th m kh o T n hiệu th m kh o Ph n mềm nhận ng vô n Tỷ số t n hiệu nhiễu truy ập ph n hi th o t n số tr gi o n s ng m ng Tin nhắn ngắn n us r- ng vào ng r truy ập ph n hi th o th i gi n truy ập ph n hi th o th i gi n Kho ng th i gi n phát h p kênh ph n hi th o th i gi n Lệnh ông su t ph t Di ng ng r ng mở r ng ng lên M ng truy nhập vô n m t t ệ thống thông tin i ng Thi t ị ng i ng (Di ng) Môi tr ng nhà o Tho i sử ng P truy ập ph n hi th o mã ng r ng i o th ng ng không y M L L M O N L M N ii D N M M L K V T TẮT iii iv D N M B N vii D N M L UV T N M N viii U Ư N TỔN QU N Á T 1 Lị h sử ph t triển hệ thống thông tin i 1 Th hệ (First 1 Th hệ (Thir 12 N ng n r tion) 1 Th hệ (S on 114 N T ÔN T N D n r tion) n r tion) ông nghệ i i thiệu ông nghệ LT M t số tr ng 1 Tố 1212 m ng LT liệu ỉnh trễ Thông l ng ng i ng Kh n ng i ng S linh ho t ng thông ph t n số V ng o ph 10 2 T ng qu n ng ownlink LT 10 2 Mô h nh ng xuống LT 10 2 Ph ng th truy nhập OFDM 12 2 Ki n trúc khung gói tin 13 2 Kênh vật l 123 Ư l ownlink LT 15 ng kênh truyền LT 18 Mô h nh t n hiệu 18 1232 tr thử nghiệm l v ng kênh 18 1233 n ng 20 Mô h nh kênh truyền Ư N QUY O ịnh 211 M N LT 23 m ng không ịnh nt n M MO 21 y 23 truy ập m ng LT 24 1 ịnh u vào 2112 ur ịnh 1 Qu tr nh ịnh LT 24 LT 25 LT 25 2 Quy ho h hi ti t m ng LT 26 Quy ho h v ng v ng ph s ng 29 Quy ho h ung l 241 ịnh 242 Ư l ng 33 thông l ng i ng l u l ng truy ập y u tố v 2.4.3 T nh to n sit ng trung ung l nh 33 t qu 35 ng 36 Ư N ỨN D N QUY O P Ủ SÓN LT CHO VÙNG THÁI NGUYÊN 37 31 i i thiệu 37 Quy ho h ung l 321 số liệu 3.2.2 L u l ng ng 39 thống kê khu v ll 42 3 Quy ho h v ng ph x 331 kh o s t 40 ịnh vị tr sit 45 i i thuật M-PAM (Modified – PAM) 48 3 h n v ng 49 333X ịnh số no B 51 3 Nhận x t gi i thuật kh n ng p K T LUẬN V T P L UT L ƯỚN P ÁT TR ỂN Ề T MK ng th t 54 55 O 56 57 vi NH M B ng 1 Tỷ lệ BẢNG ng thông RBs 15 B ng 2.1 : Các giá trị K sử d ng cho tính tốn vùng ph sóng 32 B ng 2.2: Thông l ng ll ng xuống trung bình c a LTE 34 B ng 3.1: Thống kê m t số dịch v d liệu 41 vii NH M Hình 1.1: Mơ hình hệ thống H NH ẢNH ng xuống LTE v i MIMO 2x2 11 nh 2: OFDM OFDM 13 nh 1.3: C u trúc frame LTE 13 nh 4: Fr m lo i 14 nh 5: Kênh ng xuống LTE 15 nh 6: K hiệu tham chi u downlink LTE sử d ng tiền tố l p normal 17 nh 7: Q tr nh truyền ng downlink LTE 17 Hình 1.8: Tùy ch n cân b ng cho LTE, miền th i gian miền t n số 20 nh 9: Mô h nh truyền nh 10: nt n M MO 21 mô h nh M MO ng downlink LTE 22 Hình 2.1: Quá trình quy ho ch m ng khơng dây nói chung 23 nh 2: ịnh c m ng LTE 26 Hình 2.3: Quy ho h l u l ng ph s ng l u l Hình 2.4 : Thuật tốn tối u gi a dung l ng 28 ng diện tích vùng ph sóng 29 Hình 2.5: Các tham số c a mơ hình Walfisch-Ikegami 31 nh 1: S vùng Thái Nguyên lân cận theo Google map 38 nh 2: S vùng Thái Nguyên lân cận theo Google satellite 39 nh 3: Thông l ng tho i mô ph ng theo gi i thuật lập lịch 44 Hình 3.4: Công c Universal map download 46 Hình 3.5: File b n số vùng thành phố Thái Nguyên 47 Hình 3.6: B n quy ho ch ph sóng vùng Thái Nguyên 47 Hình 3.7: B n rút g n quy ho ch ph sóng vùng Thái Nguyên 48 Hình 3.8: Số hóa b n Hình 3.9: B n Hình 3.10: B n Hình 3.11: Khu v khu v c xét 50 số thành phố Thái Nguyên ph cận 50 quy ho ch m ng thông tin i ng thành phố Thái Nguyên 53 t tr m 53 viii ) S u ã lo v i h n t ti n hành k t h p h ng l i thành nh BMP n ng M ps om in r Univ rs l M ps Downlo r Hình 3.4: Công c Universal map download ể t o fil n t mở M pinfo s u số * t t t xã, ph h h n R gist r nhập t iểm ng k googl ng ị khu v vào Fil -> Op n th m typ R st r m g ng nh àn nh n thành phố Th i Nguyên t l u nh h n Fil s of Ti p t t iểm (nên h n rth ) uối fil * t n Nguyên ã h n 46 ng t ng k vị tr nh số khu v u nh qu t ễ nh ểl y trung tâm thành phố Th i Hình 3.5: File đồ số vùng thành phố Th i Ngun [9] Hình 3.6: Bản đồ quy hoạch phủ sóng vùng Th i Nguyên [9] 47 Hình 3.7: Bản đồ rút gọn quy hoạch phủ sóng vùng Thái Nguyên [9] 3.3.1 G ả t uật M-PAM (Modified – PAM) i i thuật M – P M gi i thuật M –P M ã sử t ởng ng gi i thuật ịnh th m số k (Trong P M ph i x gi i thuật gi i thuật P M v i iểm m i quy ho h m ng viễn thông ể x ịnh th m số k tr ng nh m x h y gi i thuật) ịnh vị tr sit ng v ng ph s ng  Giải thuật dùng để x c định tham số k [6] Sử l ng gi i thuật quy ho h v ng ph s ng LT ng v ng ph n thi t th o quy ho h v ng ph số v ng ph quy ho h ung l n thi t th o ng  k1= Số l v )/( iện t h ể t nh to n số ng v ng ph th o quy ho h v ng ph = (t ng iện t h khu m i  k2 = Số l ll) ng v ng ph th o quy ho h ung l bao khu v )/(số thuê o mà m i ll ph ng = (t ng số thuê v )  k = Max(k1, k2) h t nh iện t h khu v kh o h số thuê o mà m i ng trên)  Giải thuật M – PAM [6] Bước 1: X ịnh số l ng l ng 48 ll mong muốn (k) ll ph v (th m Bước 2: Ch n ng u nhiên k iểm Bước 3: Ph n phối hỉ ph thu non – m oi vào m oi uy nh t, Bước 4: T nh to n gi trị T m oi n non – m oi TCcurrent  ng ể làm m oi l n ho m oi , m t non – medoid m oi g n n nh t u tiên (t ng ài liên k t a nó) k   dis (medoid , non  medoid i medoid n j i j ) Qu trình l p (nh m h n r nh ng iểm th h h p nh t ể làm medoid):  Sử ng v ng l p for l ng vào nh u ể h n tr  Kiểm tr t ng p m oi , non – m oi ; th y ng h p kiểm tr i t nh h t h ng:  medoid  non – medoid  non – medoid  medoid  Ti n hành l i Tm 3, T nh l i gi trị TCi , j ho tr ng h p p i,j s o ho gi trị TCi , j nh nh t  Ti n hành so s nh: N u TCi , j < TCcurrent th t h p nhận s th y L p l i qu tr nh l p trên, iều kiện ể k t th i v ng l p là: TCi , j  TCcurrent S u ã x th o ịnh tr ng h p tối u, t ph n phối m oi , hiển thị k t qu Vị tr m oi non – medoid h nh vị tr th m kh o dành cho EnodeB 3.3.2 C ọn vùn D n l thông tin i h n n t vẽ h nh vuông iểm mong muốn ng thành phố 49 m th o trở thành on ng BS m ng Hình 3.8: Số hóa đồ khu vực xét [9] Hình 3.9: Bản đồ số thành phố Th i Nguyên ph cận Tập liệu non-m io s thuật n t t h r thành tập og m:T ên iểm nh [km] [ n /km2] ) nh iểm n n th o l h nt y liệu n số n số h hỉ r h nh 50 m io s – liệu thông tin (kho ng hiều vào iện t h này ài quãng n ,t ng, mật thể hi n ng ho gi i ht i n r v ng 3.3.3 X c địn số EnodeB  T nh to n số no B v ng ph thông số t nh to n : sử 95% p ng ho khu v ng s tor v ng ph v i yêu ô thị, u v ng ph t n k nh s tor v ng ph R = 30km, iện t h sector v ng ph : S = k R2 S = 1,95.(30)2 = 1,755 km2 Trong : K th m kh o u h nh tr m ng gi trị s u [6]: – Sector Ommi K 2,6  Số v ng ph 1,3 n thi t ể k1 = (t ng iện t h khu v th khu v Số thuê o (gi sử nhà ung i thi t sử u v ng ph : yêu 146,7  83,5  k1  90 1,755 u ung l n kh o s t th o ng ng : 393,294 Số thuê bao o n ho m i nhà m ng : 255,65 = 85,22 p) ng m u t i sử ng 4/12 ( nghĩ t n số sử ng hi thành 12 nh m t n số n ịnh vị tr tr m gố ), kênh l u l ng oS 1%  Tr thể m i v ng ph 30T v i p 20,337 rl ng gi sử m i thuê thể ph ng rl ng th m i ll o hi m 0,033 rl ng th m i ll 20,337 = 616 thuê bao 0,033 v  Số thuê o m t s tor v ng ph  Số v ng ph n thi t ể k = (t ng số thuê = o yêu 2,6 hi m kho ng 65%  65%*393,294 = 255,65 (thuê bao) t ung m – Sector 1,95 )/ iện t h m i v ng ph =  T nh to n số no B D n số – Sector o m thể ph v : 616*3=1848 thuê o ung l o n)/(Số thuê 85,22 =0,046 1848  k2 = k = max (k1 , k ) = k = 51 o ng : o tối m i s tor v ng ph )  Số no B yêu u – ông su t ph t trung  Diện t h trung V ik=9 Th o tr m i s tor S = 60,55 = 60,55km2 S 60,55 = = 30,56 km2 R = 10,88km 1,95 1,95 Mà S = 1,95.R2  R2 =  B n k nh trung nh nh ho EnodeB nh m i s tor 10,88km ng xuống : PinMS = PoutEnodeB – LdupEnodeB – L f  EnodeB + Ga EnodeB – Lslan_ EnodeB – L V i : PinMS ông su t thu LduplBTS : suy h o l MS = – 92dBm song ông = 0,5 B L fBTS = 3dB GaBTS : sử ng nt n ịnh h ng = 17 Bi LslantBTS = 3,5dBi L : suy h o ng truyền ng v i n k nh v ng ph th o ông th s u: L = A + Blgf – 13,82lg H b + a H m + (44,9 – 6,55lg H b )logR V i iều kiện : f = 900M z, H b = 30m o nt n ph t a H m = (1,1lgf – 0,7) H m – (1,56lgf – 0,8) (m) H m : hiều o nt n = 1,5m R = 10,88km bán kính cell A = 69,55 B = 26,16 L = 104,75dB Th y gi trị vào ông th t t nh PoutEnodeB = PinMS + LduplBTS + L f  EnodeB – Ga EnodeB + Lslant EnodeB + L PoutEnodeB = –92 + 0,5 + – 17 + 3,5 + 104,75 = 2,75dB Vậy ông su t ph t trung  K t qu l h nh ho h n vị tr ng tr nh t nh th o gi i thuật P M v i c EnodeB PoutEnodeB = 2,75dB no B gi i thuật M-PAM l th m số t nh nh 52 h n vị tr no B y ho k t qu nh hỉ hình 4.4 ng i y Vị tr t tr m no B th o t nh to n u h m m u x nh ậm, vuông Nh vậy, thành phố Th i Nguyên g m ph m i khu v thể ng thành phố Th i Nguyên t tr m no B (vị tr t nh to n tối u số th o nh tiêu h v ng ph ận quy ho h hi thành s tor (khu v ), h m vuông, x nh) vị tr th m kh o ể x y y vị tr ng tr m no B gi i thuật Hình 3.10: Bản đồ quy hoạch mạng thông tin di động thành phố Thái Nguyên Hình 3.11: Khu vực đ t trạm 53 3.3.4 N ận xét ả t uật v k ả năn Quy ho h y u tố ràng u p dụn t ực tế i h i ph i xử l m t khối l nên việ p ng p ng ph h p v i th gi i thuật t l thuy t ng  Nên r ph n hi ti t gi o iểm u mật thể p o nh lị h, tiễn t on on n th tiễn n s th y i n: ng v liệu mà gi i thuật xử l nên ng quố l T nên th m kh o ng, quận nh m ể ễ àng x T nh ng kinh nghiệm t th u t m hiểu gi i thuật, t nhận th y r ng : ể gi i thuật tốt h n th h n n tiễn S u m t th i gi n nghiên p th m số gi i thuật ể ph n t h xử l số liệu v n ề t t y u Tuy nhiên qu tr nh quy ho h tiễn nên việ ng l n ng ịnh t i tr ng l u l liệu ho nhà l n, ối t , tr n ghi ng t i m i iểm ng khu ng h , ệnh viện, n v i khu u trung t m mu sắm  Về kh n ng p ng th tiễn th gi i thuật r t h u ng m ng viễn thông ho m t khu v ng việ x y r ng l n nh m t tỉnh, m t v ng h y m t miền  x y iv i m ng ng thêm ph i t nh n sẵn sở h t ng, hệ thống h ố ịnh 54 liệu tr m BS th muốn K T LU N V H ỚNG PHÁT TRIỂN Ề T I ông nghệ LT m t ông nghệ m i, ã ng ti p t u triển kh i toàn th gi i, v i kh n ng truyền t i tố m ng n gi n, sử thống tr ( ng ng t n hiệu qu hoàn tồn t SM& W DM trở thành hệ thống thơng tin i ông nghệ LT ) u t nghĩ th Trong luận v n ã ề ập m t tr ng t m g m o ki n tr ng th h v i hệ m t m ng toàn P LT ng toàn n thi t nghiên thể ng l i V việ t m hiểu t h t ng qu n ông nghệ LT , ph n : - T m hiểu qu tr nh ph t triển hệ thống thông tin i ng ho n nay, gi i thiệu ông nghệ LT - T m hiểu ki n tr sử m ng LT , thành ph n h nh gi o th ng - T nh h nh triển kh i LT th gi i t i V T N M - So s nh hiệu su t hệ thống M MO-OFDM (tuy n xuống) S -FDMA (tuy n lên) v i hệ thống MCMC-CDMA- ph ng th truyền ng ng nh kh n ng ng i th r ối v i LT ph t triển - Tuy nhiên, v th i gi n h n h mà luận v n hỉ m i h n ng l i m t m hiểu LT th sánh SC-FDM v i M M -CDMA, minh h t i khu v ng Th i Nguyên h h r n so qu tr nh quy ho h m ng LT nh gi hiệu su t h y kh n ng sử m ng - Nghiên u th v i M M - DM , th mô ph ng ể so s nh hệ thống M MO-OFDM v i M MO 2x2, t mở r ng h n v i M MO 4x4, 8x8 LT m t ông nghệ ph t triển s u so v i W M X, nh ng v i t nh tuyệt v i mà n m l i, nên n y ã th gi i ã ng h l h n ể triển kh i r t nhiều nhà h t o thi t ị ng ã ti n hành t h h p ông nghệ LT vào s n phẩm th nhà m ng nhà m ng l n u uối m nh T i việt n m ng ã ti n hành thử ngiệm ông nghệ LT ng ã nh ng k t qu kh qu n Do vậy, việ nắm ông nghệ LT h t s n thi t 55 t T I LIỆU TH M KHẢO [1] Harri Holma, Antti Toskala,LTE for UMTS – OFDMA and SC-FDMA Based Radio Access, John Wiley & Sons Ltd, 2009 [2] Agilent Technologies, 3GPP Long Term Evolution: System Overview, Product Development,and Test Challenges, 2009 [3] kuno, M Wruli h, M Rupp, “System level simulation of LTE networks”, in Pro 2010 71st V hi ul r T hnology onf r n , Taipei, Taiwan, May 2010 Available at: http://publik.tuwien.ac.at/files/PubDat_184908.pdf [4] Elias Yaacoub, Hussein AI- s i, n Z h r D wy, “Algorithms for the LTE Uplink”, Research paper, 2009 [5] Dr J y sh Kot h , J son Wong, “LTE: MIMO Techniques in 3GPPLTE, Freescale Semiconductors”, Pr s nt tion, Nov 5, 2008 [6] Nguyễn Quố Tu n, “Bài giảng quy hoạch mạng viễn thông”, Tr h [7] ông nghệ, Nguyễn h ih Quố gi ih Quố gi àN i [8] Thống kê số liệu t i tập oàn VNPT tỉnh Th i Nguyên, 2010 [9] Google maps [10] i ng i àN i o Quyền, “Bài giảng thông tin di động nâng cao”, Tr ông nghệ, ng w sit th m kh o : www.Thongtincongnghe.com www.Vntelecom.org www.Tapchibcvt.gov.vn www.Tudiencongnghe.net www.Xahoithongtin.com 56 PH L Giải thuật k-medoids thuật to n phân nhóm thuật to n có liên quan đến k-medoids chọn điểm liệu c c trung tâm function [IDX, Cluster, Err] = kmedoid2(data, NC, maxIter, varargin) % Thực phân nhóm s d ng thuật to n PAM % C m ban đầu phải lựa chọn từ c c điểm liệu có sẵn % Kích thước liệu dsize = size(data); % Số điểm liệu L = dsize(1); % Số c c thuộc tính NF = dsize(2); % Kích thước c m csize = [NC, NF]; if (L < NC) error('Qu số liệu để làm k clusters'); end if(nargin > 5) error('Usage : [IDX,C,error] = kmedoid(data, NC, maxIter, [init_cluster], [init_idx])'); elseif(nargin == 5) vsize1 = size(varargin{1}); vsize2 = length(varargin{2}); if(isequal(vsize1,csize)) Cluster = varargin{1}; else error('Kích thước c m khởi tạo không "); end 57 if(vsize2 == NC) IDX = varargin{2}; else error('Kích thước x c cho số c m ban u "); end elseif(nargin == 4) rror(''B n ph i ung p h i ối số t y h n: init_ lust r, init_i x'); elseif(nargin == 3) % khơng có c m ban đầu cung cấp IDX = randint(NC,1,L)+1; Cluster = data(IDX,:); % Khởi tạo c c c m ban đầu ngẫu nhiên else display('Usage : [IDX,C] = kmedoid(data, NC, maxIter, [init_cluster], [init_idx]'); error('Chức cần đối số'); end % Mảng để lưu trữ chi phí cho lần l p Err = zeros(maxIter,1); Z = zeros(NC, NF, L); for i = 1:L Z(:,:,i) = repmat(data(i,:),NC,1); end FIRST = 1; total_cost = 0; for Iteration = 1:maxIter %disp(Iteration); if(FIRST) % Đầu tiên, tính to n c c chi phí liên quan đến c c c m ban đầu C = repmat(Cluster,[1,1,L]); %Kho ng B = sqrt(sum((abs(Z-C).^2),2)); 58 h u li B1 = squeeze(B); [Bmin,Bidx] = min(B1,[],1); cost = zeros(1,NC); for k = 1:NC cost(k) = sum(Bmin(Bidx==k)); end total_cost = sum(cost); % Đ t lại cờ FIRST FIRST = 0; % Nếu lần đầu else % thay đổi trọng tâm c m % có thay đổi, chấp nhận Nếu khơng từ chối while(1) % tìm số ngẫu nhiên khoảng từ Tidx = randint(1,1,L)+1; nL % tr nh c ch điểm trung tâm if(isempty(find(IDX==Tidx))) chọn break; end end % thay đổi ngẫu nhiên c m pos = randint(1,1,NC) + 1; % i l i OLD_IDX = IDX; nh s h hỉ số IDX(pos) = Tidx; %g n c c trung tâm c m PCluster = Cluster; Cluster = data(IDX,:); % Đối với điểm liệu, tìm c c c m gần với % tính to n c c chi phí thu từ liên kết 59 C = repmat(Cluster,[1,1,L]); %đường Quang trung %B = sum(abs(Z-C),2); %khoảng c ch Euclide B = sqrt(sum((abs(Z-C).^2),2)); % Row - cluster (NC) % Column - each data point 1, 2, L B1 = squeeze(B); % Đối với điểm liệu, tìm c c c m gần khu vực % size(Bmin) = 1xL [Bmin,Bidx] = min(B1,[],1); cost = zeros(1,NC); for k = 1:NC cost(k) = sum(Bmin(Bidx==k)); end previous_cost = total_cost; total_cost = sum(cost); if(total_cost > previous_cost) % cost increases % lựa chọn hơn, ph c hồi lại gi trị ban đầu IDX = OLD_IDX; total_cost = previous_cost; Cluster = PCluster; end end Err(Iteration) = total_cost; end if(nargout == 0) disp(IDX); elseif(nargout > 3) rror(' h n ng ho ph p hỉ k t qu end return 60 u r ");