Tài liệu Lý thuyết Anten pdf

PHẦN LÝ THUYẾT ANTEN CHƯƠNG I : KHÁI NIỆM TRƯỜNG ĐIỆN TỪ -oOo I Đặc điểm sóng điện từ Toàn lý thuyết anten xây dựng sở sóng điện từ Điện trường từ trường đồng thời tồn không gian thống tạo thành trường điện từ Trường điện từ dạng vật chất bản, chuyển động với vận tốc c hệ quy chiếu quán tính chân không; thể tồn vận động qua tương với dạng vật chất khác hạt môi trường chất mang điện Trường điện từ có mang lượng nhiệt độ: W ( B E ( HD) J m3 ) Trường điện từ đặc trưng đại lượng sau : E : vector cường độ điện trường ( D : Vector cảm ứng điện ( C V ) m ) m2 Wb B : Vector cảm ứng từ ( H : Vector cường độ từ trường ( m2 ) A ) m Trường điện từ biểu diễn qua hệ thống phương trình Maxwell rot E = div D B (1.1) t (1.3) rotH J div B Ngoaøi có phương trình liên hệ : D (1.2) t (1.4) D B H J E ( E rH A m2 ( rE ( Vs m2 T) ) C m2 ) (1.5) (1.6) (1.7) J :thông lượng vector mật độ dòng dẫn ( : phân bố điện tích khối ( C m3 : độ dẫn điện môi trường dẫn ( A m2 ) ) m s ) m : hệ số điện thẩm tuyệt đối môi trường ( A ) m 10 F ( ) 36 m Đối với môi trường chân không r 10 :các điện môi thông thường r : không khí r 103 104 : số muối senhet : hệ số từ thẩm môi trường ( H ) m Đối với môi trường chân không r 10 H ( ) m : môi trường thông thường Nguồn tạo trường điện từ dòng điện từ trường Ý nghóa hệ phương trình Maxwell - Phương trình Maxwell (1.1)và (1.2) nêu rõ từ trường điện trường biến thiên gắn bó với có tính chất xoáy - Phương trình Maxwell (1.3)và (1.4)mô tả dạng hình học hai mặt thể điện trường từ trường II Sóng điện từ : Phương trình sóng điện từ có dạng E = Eo cos( t - Tương tự : B = Bo cos( t - x ) V x ) V Eo , Bo , phụ thuộc điều kiện đầu Hàm E ( x,t ) B ( x,t ) hàm sóng, điện trường từ trường lan truyền không gian dạng sóng Vậy sóng điện từ trường điện từ biến thiên truyền không gian Sự lan truyền sóng điện từ thể qua lan truyền lượng điện từ, cường độ trường (sóng E, sóng H )và (sóng A, sóng ) Theo dạng mặt phẳng đồng pha sóng điện từ mà ta có sóng điện từ phẳng, sóng trụ sóng cầu Sóng điện từ phẳng sóng điện từ có mặt đồng pha mặt phẳng, phương truyền sóng phẳng nơi vuông góc với mặt phẳng xác định Sóng điện từ gọi đơn sắc hay hòa vector cường độ điện trường, từ trường biến đổi hình sin theo thời gian với tần số wxác định - Sóng phẳng gọi sóng phẳng đồng vector E,H sóng phụ thuộc tọa độ không gian - Sóng trụ tròn dạng sóng mà trường lượng lan truyền theo chiều bán kính r tỏa từ trục không gian xung quanh hướng vào trục Tính chất sóng điện từ tồn môi trường chất môi trường chân không a Sóng điện từ có chức năng: - Tạo chùm tia công suất theo hướng định trước - Lái chùm tia để khu vực bao phủ sóng - Cho phép đo đạc thông tin góc để xác định hướng b Sóng điện từ sóng không gian Tại điểm không gian, phương E , B vuông góc với phương truyền sóng Ta nói sóng điện từ phẳng thuộc loại sóng điện từ ngang TEM c E , H pha có trị số tiû lệ với d Biên độ điện trưồng từ trường sóng giữ không đổi trình lan truyền e Vận tốc dịch chuyển mặt đẳng pha gọi vận tốc pha v p (hay vận tốc truyền lượng trường (cả điện từ kèm theo nhau) vp (1.8) gọi hệ số pha (rad/m) f Năng lượng thông - Năng lượng sóng điện từ lượng trường điện từ Năng lượng tồn vùng không gian có sóng điện từ - Năng thông lượng truyền đơn vị diện tích đơn vị thời gian Định nghóa :Vectơ Poynting Là vector mật độ dòng công suất điện từ, vector công suất điện từ chảy qua đơn vị diện tích đặt vuông góc với dòng chảy ExH III Sóng điện từ môi trường bán dẫn điện đồng chất: Ta có điện dẫn suất môi trường phương trình Maxwell là: rot H = divB E+ D t rot E = divD hệ thống B t Nếu nguồn trường biến thiên theo qui luật điều hòa jwt E Re( e E (1.9) ) Re phần thực đại lượng phức Từ phương trình liên hệ, lấy đạo hàm (1.9), kết hợp với phương trình liên hệ so sánh với phương trình Maxwell ta có: p ( i F ) ( ) m Từ tương ứng sau khảo sát người ta rút kết luận : Sóng loại sóng ngang Vectơ cường độ điện trường từ trường vuông góc với hai vuông góc với phương truyền sóng Sóng lan truyền bị môi trường hấp thụ làm suy giảm biên độ Mức độ w suy giảm theo hàm mủ đặc trưng thừa số với = c (1.10) hệ số hấp thụ Điện trường từ trường truyền với vận tốc giống v= c/n (m/s) Ở điểm không gian từ trường lệch pha với điện trường góc = arctg n (độ) Biên độ cường độ điện trường từ trường có quan hệ với biểu thức : Hm = n2 120 Em ( A ) m Qua biểu thức hệ số hấp thu (1.10)ï ta thấy điện dẫn tăng bước sóng ngắn hấp thụ nhiều Truyền sóng qua môi trường bán dẫn điện (mặt đất, mặt biển) để giảm bớt tổn hao nên dùng sóng có bước sóng lớn IV Sóng vô tuyến: Đặc điểm sóng vô tuyến: Sóng vô tuyến có chung đặc tính với dạng chuyển động khác Có thể nói cách gần đúng, chuyển động sóng bao gồm liên tục đỉnh sóng đáy sóng với khoảng cách chuyển động theo tốc độ cố định Ví dụ ta nhìn mảnh gỗ mặt nước bị nâng lên hạ xuống sóng qua, gió dòng nước không dịch chuyển hướng Điều sóng tạo biến động xa, tịnh tiến qua môi trường (trong trường hợp biển) với tốc độ cố định thân môi trường không dịch chuyển Khoảng cách hai sóng liên tục gọi bước sóng (kí hiệu ) Một dao động hoàn chỉnh từ đỉnh sóng qua đáy sóng đến đỉnh sóng gọi chu kì Số chu kì sóng qua điểm cố định khoảng thời gian cho sẵn gọi tần số (kí hiệu f), giải thích tần số số chu kì giây gọi Hez (Hz) Rõ ràng số chu kì giây phụ thuộc vào bước sóng tốc độ mà sóng truyền lan (tốc độ kí hiệu c) Sóng có bước sóng 2m chuyển động với tốc độ 10m giây phải dao động với tần số chu kì giây.Mối quan hệ yếu tố sau: Tốc độ tần số x bước sóng Hoặc C = f x (1.11) Bằng cách suy luận từ sóng biển chứng minh thật xa Một đá ném xuống hồ cho thấy sóng truyền lan với tốc độ cố định hướng không bị cản trở Sóng vô tuyến có đầy đủ đặc tính vừa mô tả khác chất chổ sóng vô tuyến sóng điện từ tạo nên trường điện trường từ, mà chấn động học Sóng điện từ có đặc tính với sóng ánh sáng, bước sóng chúng dài yếu tố ảnh hưởng lớn đến đặc tính chúng Khi nói đến chúng người ta thường dùng khái niệm điện khái niệm từ, cường độ chúng đựơc đo vol/ m , công suất đo oat/ m Sự tồn lượng sóng điện từ chứng minh thực cảm ứng dòng điện vật chất mà tiếp xúc Sóng vô tuyến có trục điện trục từ vuông góc Sự định hướng trục không gian gọi phân cực chúng đựơc biểu diễn theo hướng trục điện Bước sóng sóng vô tuyến thường nằm khoảng 20.000m đến 4mm Tốc độ chúng không gian tự cố định khoảng 300 triệu mét 161800 hải lý giây Giá trị sử dụng rộng rãi xét truyền lan sóng khí Vì tốc độ cố định, tần số tăng bước sóng giảm Ví dụ rõ công thức đưa sử dụng để đổi bước sóng tần số: Hãy tìm tần số đài phát sóng vô tuyến phát bước sóng 1500m f c , tần số = 300.000.000 1.500 200.000Hz Tần số sóng vô tuyến thường đo bội số Hz như: kilôhec (Khz) = 103Hz, 1mêgahec (Mhz) = 106 Hz, 1gigahec (Ghz)=109 Hz Đối với khoảng cách thời gian ngắn, thời gian thường tính micrô giây( tức phần triệu giây) Sóng vô tuyến có bước sóng ngắn gọi sóng siêu ngắn (viba ) Từ so sánh đơn giản ta hiểu chất chung sóng vô tuyến Thuộc tính chúng phụ thuộc nhiều vào phương pháp xạ chúng vào không gian Nhưng điều cần nói sóng vô tuyến tập trung lại thành chùm tia theo hướng định bị phản xạ, khúc xạ, tán xạ hay nhiễu xạ giống sóng ánh sáng tùy thuộc vào chất môi trường mà chúng qua mục tiêu mà chúng tiếp xúc Sự định hướng Ta xét thời điểm mà lúc sóng bắt đầu truyền không gian Kích thước cần thiết phản xạ để tập trung sóng vô tuyến vào chùm tia với độ rộng cho trước, phụ thuộc vào bước sóng sử dụng, bước sóng dài phản xạ rộng Vì để có phản xạ có kích thước thích hợp, để nhận chùm tia hẹp phải sử dụng sóng có bước sóng ngắn Với bước sóng 3cm phản xạ rộng fút cho chùm tia rộng khoảng 1.5 độ với phản xạ rộng 10 fút cho chùm tia rộng 0.75 độ Độ xác việc hướng cần thiết mặt phẳng ngang tức phương vị Chúng ta dễ nhận thấy rằng, phương vị chùm tia rộng cường độ yếu Bộ phản xạ phát lượng theo chùm tia hẹp với góc độ định, lượng phát từ tiêu điểm mặt phản xạ phản xạ tập trung tất lượng từ nguồn bên đến phản xạ tiêu điểm theo góc độ lúc phát Điều nói lên anten có tính định hướng cho thu phát Nó có lợi cho độ xác việc hướng mà làm tăng cường độ sóng thu Sự suy giảm cường độ tín hiệu theo khoảng cách: cường độ tín hiệu thu điểm biến đổi thay đổi khoảng cách điểm đến máy phát sau: + Cường độ trường( đo vol/ m ) tỉ lệ nghịch với khoảng cách + Công suất( đo oat/ m ) tỉ lệ nghịch với bình phương khoảng cách ï V Công suất truyền sóng lý tưởng Giả sử nguồn xạ đẳng hướng đặt không gian tự Nghóa môi trường đồng nhất, đồng hướng không hấp thụ có hệ số điện thẩm tương đối Ta tính mật độ thông lượng lượng trường xạ khoảng cách kể từ nguồn giả thuyết lượng xạ phân bố đồng mặt cầu bán kính r Biểu thị công suất xạ W Đơn vị chiều dài m , ta có biểu thức thông lượng lượng qua đơn vị điện tích mặt cầu bán kính r đơn vị thời gian laø: S= r (W / m 2) (1.12) Giá trị biểu thức (1.12) giá trị trung bình vector Poynting mặt cầu ấy: S = En Hn (1.13) Eh (V/m); Hn (A/m) gọi trị hiệu dụng vectơ cường độ điện trường từ trường Trong đơn vị đo lường hợp pháp H = E/120 ( A/m ) (1.13) viết lại S= E2 120 (W/ m ) Eh = 30 p (V/m) r Trong thực tế người ta dùng hệ thống xạ có tính phương hướng Mức độ định hướng đánh giá hệ số phương trình D, hệ số D hệ số đặc trưng cho mật độ tập trung lượng xạ anten theo hướng Có thể hiểu cách đơn giản sau: anten có hướng công suất xạ P có hệ số tính phương hướng hướng D tạo điểm thu hướng cường độ trường có trị giống anten hướng có công suất PD tạo Như việc sử dụng anten có hướng tương đương với việc tăng công suất xạ lên so với anten vô hướng đó: Eh = 90 PD r Trị số biên độ cường độ trường Em = 60 pD (V/m) r Trị số tức thời cường độ trường : E= 60 pD n cos ( t- ) r c 60 pD cos( wt-Kr ) r w : tần số góc sóng K= w c w hệ số sóng : bước sóng không gian tự Trường hợp nguồn xạ đặt không gian tự mà đặt mặt đất dẫn điện lý tưởng, lượng phân bố theo hình cầu, trị số D tăng gấp đôi cường độ trường tăng lên lần VI Phân loại sóng theo vô tuyến điện theo băng sóng theo phương thức lan truyền Các sóng vô tuyến điện chia thành băng sóng Sóng cực dài: sóng có bước sóng lớn 10.000 m (tần số thấp hớn 30 Khz ) Sóng dài: sóng có bước sóng từ 10.000 m đến 1.000 m Sóng trung: sóng có bước sóng từ 1.000 m đến 300Khz đến Mhz ) 100m ( tần số Sóng cực ngắn : sóng có bước sóng 10m đến 1mm đến 300.000 Mhz ) (tần số 30 Mhz ¬ Những phương thức lan truyền sóng vô tuyến điện Mặt cắt diện antena 0.8 0.6 0.4 0.2 x (m) -0.2 -0.4 -0.6 -0.8 -1 -0.8 -0.6 -0.4 -0.2 Y(m) Sự xuất chảo phát xạ phụ cấu hình anten cassegrain giúp tăng hiệu suất anten lên khoảng 10% so với anten parabol chảo phát xạ Hiệu suất anten anten cassegrain thường nằm khoãng 60% đến 70% tùy thuộc vào yêu cầu mức búp sóng phụ Đối với anten Cassegrain hiệu suất anten tính tổng quát từ công thức sau: = 1: B E x hiệu suất xạ anten s :Hiệu p s p suất tràn xạ (spillover ) :hiệu suất sai pha (phase error) x:Hiệu suất B:Hiệu trái phân cực (cross polarization ) suất che chắn E :hiệu suất sai số mặt phát xạ (surface error) Việc tính toán xác tham số đòi hỏi phải giải hàm tích phân phức tạp Hơn , để giải tích phân phức tạp cần phải tìm hàm toán học búp sóng phát xạ loa phát xạ Vì vậy, anten Cassegrain việc tính toán xác hiệu suất anten vấn đề khó khăn Tuy nhiên việc sử dụng phương pháp gần giúp cho việc tìm hiệu suất anten cách nhanh Thep phương pháp hiệu suất xạ biểu diển dạng sau [2] : = sd dl làø hiệu suất tán xạ chảo phát xạ phụ, sd Là hiệu suất xạ chưa tính đến ảnh hưởng tán xạ từ chảo tán xạ phụ Theo [2] ta có biểu thức sau: dl 10log10 s 1.38 sd 0.6 Ds dB o.6 Nếu giả sử hàm phân bố xạ hiệu suất xạ có tính đến ảnh hưởng che chắn chảo phát xạ phụ gây dl , biểu diển công thức sau [2]: db dl B Ds Dm 2 Kết hợp (14), (15) ( 16) ta coù s B= ilb s sd Hiệu suất sai số mặt phát xạ (surface error) thường tính công thức sau E= exp(-4 )2 Trong : độ xác mặt phát xạ, : bước sóng tầng số công tác Cả tham số tính mm Nếu thiết kế cẩn thận thông thường tham số lại thể coi Vì vậy, kết cuối [1]: s B E = ilb s sd E x, p coù 38 Ds 06 10 10 Ds Dm 2 exp 4.2 Độkhuếch đại anten: Đối với anten Cassegrain độ khuếch đại anten tính theo công thức sau [6]: G D2 m Ds Trong hiệu suất anten tính từ công thức (19), tham số lại giải thích ba Ví dụ :tính toán tham số điện cho anten có Dm=1.5 m ví dụ phần Giả sử anten chế tạo xácbề mặt phát xạ =0.7 mm tần số 20GHz áp dụng công thức (19) 38 0.1804 0.15 10 06 2 10 0.1804 exp 1.5 0.7 3.14 15 =93.10%97.13%70.92% = 64.13% Hệ số khuếch đại tính từ công thức (20) sau: G D2 m 2 Ds 64.13% 3.142 1.52 0.18042 0.0152 47.95dB Một điểm cần ý trình tính hiệu suất anten ảnh hưởng hiệu suất sai số mặt phát xạ đến hiệu suất chung anten Thông thường hiệu suất sai số mặt phát xạ thường yếu tố định lớn đến suy giảm hiệu suất anten Tuy nhiên, hiệu suất sai số bề mặt lại phụ thuộc chủ yếu vào độ xác bề mặt phát xạ Vì vậy, việc đạt hiệu suất anten cao đồng nghóa với việc chế tạo bề mặt chảo phát xạ có độ xác cao Hình miêu tả quan hệ giửa hiệu suất anten độ xác bề mặt anten ví dụ Rõ ràng để đạt hiệu suất 70% yêu cầu cần phải co ù =0.6 mm , để đạt hiệu suất cao ví dụ 80% = 0.35 mm Kết luận: Phương pháp tổng hợp dùng phân tích tính toán anten Cassegrain trình Phương pháp đơn giản dể sử dụng chí loại máy tính cầm tay Phương pháp kết tính toán ví dụ áp dụng thực tế vào việc thiết kế anten bám cho vệ tinh q đạo thấp băng tần Ka CHƯƠNG II:GIỚI THIỆU MATLAB -oOo - I Lời giới thiệu MATLAB (Matrix Laboratory) phần mềm ưa chuộng cho lập trình tính toán kỹ thuật Nó phổ biến rộng khắp trường đại học, trung học chuyên nghiệp nhiều nước Và công cụ trợ giúp hữu hiệu cho nhà chuyên môn, sinh viên, kỹ sư, cán kỹ thuật Ngày MATLAB chứng tỏ phần mềm có giao diện cực mạnh nhiều lợi kỹ thuật lập trình Với MATLAB, công việc tính toán trở nên đơn giản nhẹ nhàng nhiều so với nhiều ngôn ngữ lập trình khác nhờ có thiết kế sẵn Toolbox font tiếng việt giúp cho người sử dụng dùng MATLAB công cụ hữu hiệu cho lónh vực chuyên ngành - Control System Toolbox: tảng cho họ Toolbox thiết kế điều khiển MATLAB Nó chứa hàm cho việc mô phỏng, phân tích thiết kế hệ thống điều khiển tự động -Frequency–Domain System Identification Toolbox: gồm M-files cho việc mô hệ thống tính toán sở phép đo đáp ứng tần số hệ thống -Fuzzy Logic Toolbox: cung cấp tập hợp đầy đủ công cụ cho việc thiết kế mô phân tích hệ thống logic mờ -Higher– Oder Spectral Analysis Toolbox: cung cấp công cụ cho việc xử lý tín hiệu dùng phổ bậc cao Các phương pháp đặc biệt hữu dụng cho phân tích tín hiệu có nguồn gốc từ trình phi tuyến hay bị nhiễu phi gaussian (non-gaussian) xâm nhập -Image Processing Toolbox: chứa công cụ cho xử lý ảnh Nó bao gồm công cụ cho việc thiết kế lọc lưu trữ ảnh, nâng cấp ảnh, phân tích thiết kế -Model Predictive Control Toolbox: đặc biệt hữu ích cho ứng dụng điều khiển với nhiều biến input output mà phần lớn có giới hạn, kỹ thuật hóa chất -Mu–Analysis and Synthesis Toolbox: chứa công cụ chuyên môn hóa cho điều khiển tối ưu hóa phân tích tổng hợp , đặc biệt kỹ thuật rôbot cao cấp mà hệ thống đa biến tuyến tính -Signal Processing Toolbox: chứa công cụ xử lý tín hiệu Các ứng dụng bao gồm: audio (đóa Compact, băng digital), video( digital HDTV, xử lý nén ảnh), viễn thông(fax, telephone), y học, địa lý -Non–Linear Control Design Toolbox: cho phép thiết kế hệ thống điều khiển tuyến tính phi tuyến sử dụng kỹ thuật tối ưu hóa sở miền thời gian -Optimization Toolbox: lệnh dùng cho tối ưu hóa hàm tuyến tính phi tuyến tổng quát -Symbolic Math Toolbox: gồm công cụ cho việc tính toán biểu thức -System Identification Toolbox: tập hợp công cụ cho ước lượng nhận dạng(tìm mô hình toán học cho hệ thống vật lý) -Robust Control Toolbox: công cụ phân tích tổng hợp hệ thống điều khiển rôbô Ngoài có Toolbox khác như: NAG Foundation Toolbox, Simulink, Simulink Realtime Workshop, Spline Toolbox, Satics Toolbox Một khả cần phải nhắc đến MATLAB biểu diễn data đồ thị hai chiều lẫn ba chiều lệnh(hàm) đơn giản II Hàm Script File Script file Dùng cửa sổ soạn thảo Matlab chương trình soạn thảo khác Word 97, soạn thảo cửa NC dạng ký tự ASCII; để tạo chương trình Thủ tục theo trình tự dòng lệnh theo cú pháp ngôn ngữ lập trình Matlab lưu với tên –file, phần mở rộng có đuôi m thường gọi Script file.Ví dụ: Hamain.m Thông thường muốn chạy Script file, cửa sổ lệnh Matlab, dấu nhắc ‘>>’ ta gõ tên Script file mà không cần gõ thêm phần mở rộng ‘.m’ sau nhấn ‘Enter’ Script file thực hiện.Ví dụ: >>Hamain nhấn Trong Matlab thông thường Script lưu thư mục /bin Matlab đến 5.2, Matlab 5.3 lưu thư mục Matlab Work Nếu bạn muốn tìm đường dẫn đến Scirpt file dùng lệnh Which.Ví dụ: >>Which Hamain Bạn muốn tìm hiểu Script file để sử dụng dùng lệnh help Ví dụ: >>help Hamain Sẽ Chương trình tạo cửa sổ giới thiệu hường dẫn thầy Lê Cảnh Trung viết Nguyễn Huỳnh Hà Bạn muốn tìm hiểu toàn Script file dùng lệnh type.Ví dụ: >>type Hamain Script file thuận tiện cho việc nhập mảng liệu lớn lưu trữ đóa mềm cần sử dụng lại Hàm M-file M-file trường hợp đặc biệt Script file, thông thường M-file hàm Nghóa M-file có truyền biến xuất kết Hàm M-file khai báo dạng : function < outputs > =fun-name (< input > , < input >, ) Trong input1, input2, … biến (ma trận hay vô hướng ) truyền cho hàm outputs trị trả (ma trận hay vô hướng) Hàm có đặc điểm sau: - Tên hàm tên file lưu trữ đóa phải giống (với phần mở rộng file m) - Các hàng ghi (tính đến hàng lệnh đầu tiên) dùng thông tin trợ giúp từ dấu nhắc Matlab nhập vào lệnh help< fun-name> - Các biến hàm có tác dụng hàm thực thi: Muốn biến hàm biến toàn cục cần có khai báo global trước tên biến - Số biến nhập vào xuất kiểm tra thông qua lệnh nargin nargout Điều có tác dụng to lớn số biến nhập vào nhỏ số biến cần thiết ta cho giá trị mặc định Ví Dụ: Hàm linspace function y = linspace (d1,d2,n) % LINSPACE linearly spaced vector % equally spaced points between x1 and x2 % LINSPACE (x1,x2,N) generates N points between x1 and x2 % % See also LOGSPACE,: % Copyright (c) 1984-94 by the Math Works, Inc if nargin == n = 100; end y =[d1+(0:n-2)*(d2-d1)/(n-1)d2]; Khi từ dấu nhắc Matlab gõ vào help linspace ta thấy: % LINSPACE linearly spaced vector % equally spaced points between x1 and x2 % LINSPACE (x1,x2,N) generates N points between x1 and x2 % % See also LOGSPACE,: % Copyright (c) 1984-94 by the Math Works, Inc Khi sử dụng lệnh linspace (-2,2) hoàn toàn giống linspace (2,2,100) Ta nhận thấy M-file từ ví dụ có - Có tên file linspace - Có phần mở rộng m - Dòng funtion - y đối số (suất kết quả) - dấu ‘=’ - (d1,d2,n) đối số vào - % dòng ghi không thực Khi chạy M-file ý dùng help để đọc thích - Chương trình II Các điều khiển rẽ nhánh vòng lặp Các toán tử Matlab MATLAB có nhiều toán tử bình thường phần mềm tính toán khác Tuy nhiên có khác toán ma trận mảng hay vô hướng Ví dụ phép nhân hay ma trận dấu( *) nhân phần tử hai ma trận dấu(.*) Cách biểu diễn ma trận: dùng dấu ; để ngăn cách hàng, dấu : để mảng 1:10 (các số từ đến 10); 1:2:10(các số 1,3,5, ,9); A(:;1) nghóa øcột tất hàng hay vectơ cột thứ ma trận A ; dấu ‘là toán tử chuyển vị ma trận ; dấu % bắt đầu ghi chú; dấu = phép gán; dấu ~ phủ định; … Điều khiển rẽ nhánh IF: Các dạng lệnh IF sau: a if < biểu thức logic >, < lệnh > end; b if , else end; c if; < lệnh biểu thức logic1 > elseif; … elseif end; Các vòng lặp: Vòng lặp for: - Cấu trúc: for x = , lệnh end; Các lệnh for end thự lần cho cột array Ví dụ: » for n =1:10 x(n) = sin(n*pi/10); end »x x= Columns through 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns through 10 0.5878 0.3090 0.0000 Vòng lặp for không bị dừng cách tăng biến lặp (n) vòng lặp: » for n = 1:10 x(n) = sin(n*pi/10); n = 10; % thay doi bien lap end »x x= Columns through 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns through 10 0.5878 0.3090 0.0000 Ta thấy kết không thay đổi dù có thêm dòng lệnh n = 10; vòng lặp for cần tránh dùng xử lý ma trận để thời gian thi hành giảm - Vòng lặp while: Cấu trúc while: while < biểu thức logic > < lệnh > end; Các lệnh while end thực lặp lại experession Lệnh break: dùng để thoát khỏi vòng lặp (for, while) bất chấp d8iều kiện kết thúc có thỏa hay chưa (thường chung với if) Ví dụ: » ESP =1; » for num = 1:1000 ESP = ESP/2; if (1+ESP)