Tài liệu Khóa luận Tối ưu thiết kế anten loga – chu kỳ bằng thuật toán di truyền: Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẬI HỌC CÔNG NGHỆ
Lưu Thị Hoa Linh
TỐI ƯU THIẾT KẾ ANTEN LOGA – CHU KỲ
BẰNG THUẬT TOÁN DI TRUYỀN
KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC HỆ CHÍNH QUY
Ngành: Điện tử - Viễn thông
HÀ NỘI - 2005
Lưu Thị Hoa Linh 1
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Lưu Thị Hoa Linh 2
ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẬI HỌC CÔNG NGHỆ
Lưu Thị Hoa Linh
TỐI ƯU THIẾT KẾ ANTEN LOGA – CHU KỲ
BẰNG THUẬT TOÁN DI TRUYỀN
KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC HỆ CHÍNH QUY
Ngành: Viễn thông
Cán bộ hướng dẫn: GS TSKH Phan Anh.
Cán bộ đồng hướng dẫn: ThS Lê Quang Toàn
HÀ NỘI - 2005
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
LỜI CẢM ƠN
Em xin chân thành cảm ơn thầy giáo GS. TSKH Phan Anh và Th.S Lê
Quang Toàn đã tận tụy hướng dẫn, giúp đỡ em trong quá trình thực hiện khóa
luận này.
Em xin gửi lời cảm ơn đến các thầy cô của trường ĐH Công nghệ - ĐHQG
Hà Nội đã trang bị kiến thức v...
53 trang |
Chia sẻ: hunglv | Lượt xem: 1259 | Lượt tải: 1
Bạn đang xem trước 20 trang mẫu tài liệu Khóa luận Tối ưu thiết kế anten loga – chu kỳ bằng thuật toán di truyền, để tải tài liệu gốc về máy bạn click vào nút DOWNLOAD ở trên
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẬI HỌC CÔNG NGHỆ
Lưu Thị Hoa Linh
TỐI ƯU THIẾT KẾ ANTEN LOGA – CHU KỲ
BẰNG THUẬT TOÁN DI TRUYỀN
KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC HỆ CHÍNH QUY
Ngành: Điện tử - Viễn thông
HÀ NỘI - 2005
Lưu Thị Hoa Linh 1
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Lưu Thị Hoa Linh 2
ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẬI HỌC CÔNG NGHỆ
Lưu Thị Hoa Linh
TỐI ƯU THIẾT KẾ ANTEN LOGA – CHU KỲ
BẰNG THUẬT TOÁN DI TRUYỀN
KHÓA LUẬN TỐT NGHIỆP ĐẠI HỌC HỆ CHÍNH QUY
Ngành: Viễn thông
Cán bộ hướng dẫn: GS TSKH Phan Anh.
Cán bộ đồng hướng dẫn: ThS Lê Quang Toàn
HÀ NỘI - 2005
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
LỜI CẢM ƠN
Em xin chân thành cảm ơn thầy giáo GS. TSKH Phan Anh và Th.S Lê
Quang Toàn đã tận tụy hướng dẫn, giúp đỡ em trong quá trình thực hiện khóa
luận này.
Em xin gửi lời cảm ơn đến các thầy cô của trường ĐH Công nghệ - ĐHQG
Hà Nội đã trang bị kiến thức và tạo điều kiện thuận lợi để em có thể hoàn thành
khóa luận một cách tốt nhất.
Con xin gửi lời cảm ơn sâu sắc đến cha mẹ, người đã ủng hộ, động viên,
chăm sóc, quan tâm đến từng bước đi của con.
Cuối cùng, tôi xin cảm ơn bạn bè tôi, những người đã giúp đỡ, động viên
tôi hoàn thành bản khóa luận.
Hà Nội, ngày 6 tháng 6 năm 2005.
Sinh viên
Lưu Thị Hoa Linh
Lưu Thị Hoa Linh 3
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Lời mở đầu
Nói đến một hệ thống vô tuyến là nói tới việc sử dụng sóng điện từ. Và thành
phần quan trọng trong việc thu và phát sóng điện từ không có gì khác chính là
anten. Anten được sử dụng trong nhiều lĩnh vực như trong các hệ thống truyền
hình, phát thanh, điều khiển, vệ tinh, …Tương ứng với những mục đích sử dụng
khác nhau đó thì các anten cũng được thiết kế và tối ưu theo các phương pháp khác
nhau nhằm đạt được kết quả mong muốn cuối cùng. Đã có nhiều phương pháp,
nhiều công trình nghiên cứu khoa học tham gia vào quá trình tối ưu này. Các
phương pháp có thể kể đến như: Phương pháp Gradien, phương pháp cổ điển,
phương pháp di truyền, …
Thuật toán di truyền không chỉ sử dụng đơn thuần trong sinh học như ta
thường nghĩ mà nó có thể giải quyết rất tốt các bài toán về điện từ hay bất cứ bài
toán nhiều tham số nào khác và trong thiết kế anten cũng không phải là trường hợp
ngoại lệ. Trong phạm vi nghiên cứu tôi đã tìm hiểu về phương pháp này để áp dụng
cho bài toán tối ưu thiết kế anten Yagi. Ngoài việc xây dựng nên một thuật toán
tính toán hoàn chỉnh tối ưu hầu hết tất cả các tham số anten, tôi còn thực hiện xây
dựng nên giao diện phần mềm tiện lợi cho người sử dụng và thực hiện thiết kế hoàn
chỉnh anten có thể sử dụng được.
Hà nội 08/05/05
Lưu Thị Hoa Linh 4
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
CHƯƠNG 1: ANTEN LOGA – CHU KỲ
1.1. Dải thông tần và dải tần công tác của anten.
1.1.1. Dải thông tần.
Dải thông tần là một đặc tính quan trọng của anten. Đó là một dải tần số
mà trong giới hạn ấy anten có thể đảm bảo quá trình bức xạ hoặc thu phổ của tín
hiệu không bị méo dạng.
Thông thường trở kháng vào của mỗi anten là một hàm số của tần số. Do
đó, nếu anten làm việc với tín hiệu có phổ rộng (tín hiệu xung, số, tín hiệu vô
tuyến truyền hình...) thì ứng với mỗi tần số khác nhau của phổ, biên độ tương đối
của dòng điện đặt vào anten (anten phát) hoặc sức điện động thu được (anten thu)
sẽ biến đổi, làm thay đổi dạng phổ của tín hiệu. Khi dùng fide tiếp điện cho anten,
sự biến đổi trở kháng vào của anten theo tần số sẽ dẫn đến tình trạng lệch phối
hợp trở kháng và xuất hiện sóng phản xạ trong fide. Khi một tín hiệu có phổ rộng
truyền qua fide thì ứng với mỗi tần số khác nhau sẽ có sự trễ pha khác nhau và
gây ra méo dạng tín hiệu. Vì vậy, tốt nhất là phải đảm bảo được trong suốt dải tần
số làm việc và constVAR = 0VAX = .
Ngoài ra vì đặc tính phương hướng của anten cũng phụ thuộc tần số, nên
khi anten làm việc với tín hiệu có phổ rộng thì biên độ tương đối của cường độ
trường bức xạ (hoặc thu được) đối với các tần số khác nhau của phổ cũng biến đổi
và gây méo dạng tín hiệu. Thường thì ảnh hưởng của yếu tố này không lớn lắm và
trong thực tế, độ rộng dải tần của anten chủ yếu được quyết định bởi đặc tính phụ
thuộc của trở kháng vào anten với tần số.
1.1.2. Dải tần công tác.
Có nhiều trường hợp chúng ta đòi hỏi anten không chỉ làm việc ở một tần
số mà nó có thể làm việc ở một số tần số khác nhau. Ứng với mỗi tần số khác
nhau ấy anten phải đảm bảo được những chỉ tiêu kỹ thuật nhất định về đặc tính
phương hướng, trở kháng vào, dải thông tần . . . Dải tần số mà trong giới hạn đó
anten làm việc với các chỉ tiêu kỹ thuật đã cho gọi là dải tần công tác của anten.
Chỉ tiêu kỹ thuật này có thể khác nhau đối với từng loại anten cụ thể.
Lưu Thị Hoa Linh 5
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Căn cứ theo dải tần số công tác, có thể phân loại anten thành 4 nhóm:
- Anten dải tần hẹp (anten điều chuẩn):
0
10%f
f
∆ <
- Anten dải tần tương đối rộng : max
min
10% 50%f
f
< <
- Anten dải tần rộng : max
min
1,5 4
1 1
f
f
< <
- Anten dải tần siêu rộng : max
min
4
1
f
f
>
Tỷ số của tần số cực đại và cực tiểu của dải tần công tác max
min
f
f
gọi là hệ số
bao trùm dải sóng.
1.2. Phương pháp mở rộng dải tần số của anten chấn tử.
Đối với mọi loại anten chấn tử thì yếu tố quyết định dải tần số công tác của
anten là sự phụ thuộc của trở kháng vào anten với tần số, còn đồ thị phương
hướng của anten thường có hướng bức xạ ít thay đổi trong một dải tần khá rộng.
Để giảm sự phụ thuộc vào tần số của trở kháng vào chấn tử, có thể áp dụng
các biện pháp sau:
- Giảm nhỏ trở kháng sóng của chấn tử,
- Biến đổi từ từ thiết diện của chấn tử,
- Hiệu chỉnh trở kháng vào của chấn tử.
Ta lần lượt khảo sát các biện pháp trên:
1.2.1. Giảm nhỏ trở kháng sóng của chấn tử.
Anten chấn tử có thể được xem như một mạch dao động mà trở kháng đặc
trưng của mạch chính là trở kháng sóng của chấn tử, điện trở tổn hao của mạch là
điện trở bức xạ của chấn tử. Hệ số phẩm chất của anten chấn tử được định nghĩa
bằng công thức:
0
A
AQ A R
ρ
∑
= (1.1)
Lưu Thị Hoa Linh 6
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Tương tự như các mạch dao động thường, để mở rộng dải tần cần giảm hệ
số phẩm chất , theo như công thức (1.1) thì việc này đồng nghĩa với việc giảm
trở kháng sóng
AQ
Aρ . Trở kháng sóng Aρ được tính theo công thức:
120 ln 1A
l
a
ρ ⎛ ⎞= −⎜⎝ ⎠⎟ (1.2)
Nếu anten có chiều dài l cố định thì việc giảm trở kháng sóng tức là tăng
đường kính a của chấn tử bằng một số cách sau: tăng đường kính của dây dẫn
làm chấn tử (a), chấn tử có thể được làm bằng tấm kim loại bản rộng (b), sắp xếp
các dây dẫn nhỏ thành bản lớn hoặc khối lớn (c).
Khi tính toán trở kháng sóng theo công thức (1.2) thì trong trường hợp (b)
cần thay bán kính tương đương bằng ¼ độ rộng của tấm kim loại. (
4td
ha = ) . Còn
trong trường hợp sắp xếp các dây dẫn nhỏ thành lồng, thì bán kính tương đương
của chấn tử được tính bằng công thức:
n
td
nra a
a
= (1.3)
Trong đó: a: bán kính thực của lồng dây dẫn,
n: số dây dẫn hợp thành lồng,
r: bán kính dây dẫn nhỏ.
Việc thực hiện giảm trở kháng sóng theo phương pháp này có thể mở rộng
dải tần số để có hệ số bao trùm dải sóng khoảng 2.5
1
khi đảm bảo hệ số sóng chạy
trong fiđe cung cấp không nhỏ hơn 0.3. Trên thực tế hệ số sóng chạy lớn hơn
( ) thì chỉ có thể mở rộng dải tần số khoảng 0.9k ≥ 1 5 % .±
1.2.2. Biến đổi từ từ thiết diện của chấn tử.
Khi nối hai đường truyền sóng trong đó có hai dạng lan truyền sóng khác
nhau (như từ ống dẫn sóng tròn sang ống dẫn sóng hình chữ nhật) thì ta có thể
dùng thiết bị chuyển tiếp dạng sóng. Còn khi nối hai đường truyền sóng có trở
kháng sóng khác nhau (về kích thước và tiết diện) thì ta có thể dùng bộ chuyển
tiếp biến đổi dần kích thước của hai đường truyền cho phù hợp.
Lưu Thị Hoa Linh 7
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Anten cũng có thể được xem như một thiết bị chuyển tiếp giữa hệ thống
fiđe tiếp điện và môi trường truyền sóng (không gian tự do), anten biến đổi sóng
điện từ ràng buộc trong fiđe thành sóng tự do truyền không gian. Do đó bản thân
anten phải có dạng kết cấu chuyển tiếp để giảm sự phụ thuộc của trở kháng vào
anten với tần số. Tức là kích thước của anten được biến đổi một cách từ từ.
β
Hình 1.2: Chấn tử đối xứng hình chóp.
Trong trường hợp fiđe tiếp điện bằng dây song hành đối xứng thì sự chuyển
tiếp từ đường dây vào không gian bên ngoài có thể được thực hiện dứi dạng hai
chóp đối xứng (hình 1.1), lúc này ta có chấn tử đối xứng hình chóp.
Nếu cặp chấn tử hình chớp này dài vô tận thì hệ thống này sẽ hình thành
một đường truyền sóng đồng nhất không có phản xạ, sóng truyền lan là sóng cầu.
Trở kháng sóng của đường truyền là hằng số, và bằng trở kháng vào của
anten.
276 cotg
2
βρ ⎛= ⎜⎝ ⎠
⎞⎟ (1.4)
Khi chóp có độ dài hữu hạn thì sóng sẽ bị phản xạ lại một phần ở đầu cuối,
phần còn lại bức xạ ra không gian. Trở kháng vào của anten sẽ là một số phức,
phụ thuộc vào tần số. Nhưng nếu chọn vị trí góc β phù hợp thì thì trở kháng vào
của anten sẽ ít phụ thuộc vào tần số trog một dải tần khá rộng. Thực tế anten có
thể làm việc với hệ số bao trùm dải sóng max
min
4 /1f
f
= . Trong dải tần này hệ số sóng
chạy trong fiđe không nhỏ hơn 0,5.
Lưu Thị Hoa Linh 8
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
1.2.3. Hiệu chỉnh trở kháng vào của chấn tử.
Theo phương pháp này thì chấn tử cần được thiết lập thế nào để kết cấu của
nó sẽ bao gồm hai bộ phận mà điện kháng vào của mỗi bộ phận ấy có dấu ngược
nhau và có thể bù nhau trong dải tần số công tác. Ví dụ như nếu một bộ phận có
trở kháng dung tính i
Cω
⎛−⎜⎝ ⎠
⎞⎟ ) thì phần kia sẽ có trở kháng cảm tính (i Lω . Hình 1.3
là sơ đồ ả có hiệu chỉnh trở kháng vào. Hai nhánh của chấn tử là ae và bf, đoạn
cod hình thành một đường dây nhánh mắc song song với chấn tử.
o
Nếu ta chọn kích thước của các chấn tử phù hợp thì thành phần điện kháng
của trở kháng vào sẽ nhỏ, điều này được thể hiện qua sơ đồ tương đương - hình
1.3b: hai nhánh của chấn tử tương đương với một đoạn dây song hành ngắn mạch.
Nếu chọn kích thước của các đoạn dây thích hợp thì đoạn dây nhánh sẽ có
điện kháng vào cảm tính, còn đoạn hở mạch ce – df sẽ có điện kháng vào dung
tính, sao cho hai điện kháng này có thể bù trừ cho nhau trong một dải tần nhất
định. Đồng thời đoạn ac – bd được chọn để trở kháng vào tại cd sau khi biến đổi
về ab sẽ có phần thực lớn và phần kháng nhỏ.
1.3. Mở rộng dải tần công tác theo nguyên lý tương tự. Áp dụng cho anten
lôga - chu kỳ.
Nguyên lý tương tự của điện động học có thể được phát biểu như sau: Nếu
biến đổi đồng thời bước sóng công tác và tất cả các kích thước của anten theo một
tỷ lệ giống nhau thì các đặc tính của anten như đồ thị phương hướng, trở kháng
a)
c a b d f
o
a b
c d
e f
b)
Hình 1.3: Sơ đồ anten có hiệu chỉnh trở kháng vào
Lưu Thị Hoa Linh 9
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
vào ... sẽ không biến đổi. Hệ số tỷ lệ này được gọi là tỷ lệ xích của phép biến đổi
tương tự.
Dựa vào nguyên lý trên có thể thiết lập các anten không phụ thuộc tần số
bằng cách cấu tạo anten từ nhiều khu vực có kích thước hình học khác nhau. Kích
thước hình học của các khu vực ấy tỷ lệ với nhau theo một hệ số nhất định Khi
anten làm việc với một bước sóng nào đó thì sẽ chỉ có một khu vực của anten
tham gia vào quá trình bức xạ. Khu vực này gọi là miền bức xạ của anten. Khi
bước sóng công tác thay đổi thì miền bức xạ của anten sẽ dịch chuyển đến khu
vực mà tỷ lệ kích thước hình học của các phần tử bức xạ với bước sóng giống như
lúc trước.
1.3.1. Nguyên lý cấu tạo của anten lôga - chu kỳ.
Anten được tạo bởi tập hợp các chấn tử có kích thước và khoảng cách
khác nhau và được tiếp điện từ một đường fiđe song hành chung như hình 1.4,
các chấn tử nhận dòng từ fiđe theo cách tiếp điện chéo:
α
1l 2l
1nl −
nl
fiđe
Hình 1.4: Kết cấu của anten loga - chu kỳ
Kích thước của các chấn tử và khoảng cách giữa chúng biến đổi dần theo
một tỉ lệ, tỉ lệ này được gọi là chu kỳ của kết cấu:
11 2
2 3
... n
n
ll l
l l l
τ −= = = = (1.5)
Đặc tính kết cấu của anten lôga - chu kỳ được xác định bởi hai thông số
chính là τ và góc α .
Nếu máy phát làm việc ở tần số of nào đó, tần số này lại là tần số cộng
hưởng của một trong các chấn tử thì trở kháng của chấn tử đó sẽ là điện trở thuần.
Lưu Thị Hoa Linh 10
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Các chấn tử khác vẫn còn thành phần điện kháng, giá trị của điện kháng càng lớn
khi độ dài của chấn tử này khác càng xa với chấn tử cộng hưởng, tức là chấn tử
này càng xa chấn tử cộng hưởng. Chấn tử cộng hưởng được kích thích mạnh nhất.
Các chấn tử không cộng hưởng có dòng điện chạy qua nhỏ nên trường bức
xạ của anten được quyết định chủ yếu bởi bức xạ của của chấn tử cộng hưởng và
một số chấn tử lân cận đó. Những chấn tử này tạo nên miền bức xạ của anten.
Dòng điện trong các chấn tử của miền bức xạ có được do tiếp nhận trực tiếp từ
fiđe và hình thành do cảm ứng điện trường của chấn tử cộng hưởng.
Các chấn tử ở phía trước chấn tử cộng hưởng có chiều dài nhỏ hơn, sẽ có
dung kháng vào, dòng cảm ứng trong chấn tử này chậm pha so với dòng trong các
chấn tử có độ dài hơn nó. Và ngược lại, các chấn tử ở phía sau chấn tử cộng
hưởng có chiều dài lớn hơn, sẽ có cảm kháng vào, dòng cảm ứng trong chấn tử
này sớm pha so với dòng trong các chấn tử có độ dài hơn nó. Các chấn tử nhận
dòng từ fiđe theo cách tiếp điện chéo nên 2 chấn tử kề nhau có dòng điện lệch pha
nhau 180 cộng với góc lệch pha do truyền sóng trên đoạn fiđe mắc giữa 2 chấn tử
đó. Ta nhận được dòng tổng hợp trong các chấn tử của miền bức xạ có góc lệch
pha giảm dần theo chiều giảm kích thước anten.
o
Các chấn tử có quan hệ pha như trên, nên chấn tử đứng trước chấn tử cộng
hưởng sẽ thỏa mãn điều kiện chấn tử dẫn xạ, còn chấn tử tử dứng sau sẽ thỏa mãn
điều kiện của chấn tử phản xạ. Bức xạ của anten chủ yếu do chấn tử cộng hưởng
quyết định và sẽ được định hướng theo trục anten, về phía các chấn tử nhỏ hơn.
Nếu tần số máy phát giảm đi, còn là ofτ (τ nhỏ hơn 1) thì vai trò của chấn
tử cộng hưởng sẽ được dịch chuyển sang chấn tử có độ dài lớn hơn kế đó, và
ngược lại, nếu tần số tăng lên bằng ofτ thì chấn tử cộng hưởng sẽ chuyển sang
chấn tử ngắn hơn kế đó.
Chấn tử cộng hưởng với tần số 1l 1f , ta có 11 2
l λ= . Nếu tần số máy phát
giảm xuống ' ' 11f f
λτ λ τ= → = , chấn tử cộng hưởng mới có độ dài là
'
' 1 1
22 2
ll lλλ τ τ= = = = .
Ở các tần số 1 1nnf fτ −= (1.6)
Lưu Thị Hoa Linh 11
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
thì các chấn tử cộng hưởng có độ dài tương ứng là: 1 1n n
ll τ −= (1.7)
trên anten cũng sẽ xuất hiện miền bức xạ mà chấn tử phản xạ có độ dài chính là ln.
Trong đó: là số thứ tự các chấn tử, n
nf là cộng hưởng của chấn tử thứ , n
nl là độ dài của chấn tử thứ . n
Miền bức xạ của anten sẽ dịch chuyển khi tần số công tác thay đổi, nhưng hướng
bức xạ cực đại vẫn giữ nguyên.
Lấy loga hai vế của biểu thức (1.6) ta có
1ln ( 1) ln lnnf n fτ= − + (1.8)
Ta thấy khi biểu thị tần số trên thang đo logarit thì tần số cộng hưởng của
anten sẽ được lặp lại qua các khoảng giống nhau là lnτ , chính vì thế mà người ta
gọi anten là anten Lôga-chu kỳ.
Khi anten hoạt động ở tần số cộng hưởng thì các thông số điện như đồ thị
phương hướng, trở kháng vào,… sẽ không có sự thay đổi. Nhưng ứng với các tần
số trung tâm giữa các tần số cộng hưởng 1 2 2 3 1, ,..., n nf f f f f f−÷ ÷ ÷ , các tần số của
anten sẽ bị thay đổi nhỏ. Ta cũng có thể cấu tạo anten sao cho trong khoảng giưa
2 tần số kề nhau các thông số biến đổi trong một giới hạn chấp nhận được.
Đồ thị phương hướng của anten được xác định bởi số lượng chấn tử của
miền bức xạ tác dụng, thông thường là khoảng 3 5÷ , và bởi tương quan biên độ và
pha của dòng điện trong các chấn tử ấy. Các đại lượng này lại phụ thuộc vào các
thông số hình học τ và α của kết cấu anten.
Với α xác định, tăng τ thì số chấn tử thuộc miền bức xạ tác dụng cũng
tăng, do đó đồ thị phương hướng hẹp lại. Nhưng nếu tăng τ quá lớn thì đặc tính
phương hướng lại xấu đi vì lúc đó kích thước miền bức xạ tác dụng giảm do các
chấn tử qúa gần nhau. Giữ nguyên τ , giảm α đến một giới hạn nhất định nào đó
sẽ làm hẹp đồ thị vì khi đó khoảng cách giữa các chấn tử lại tăng và do đó tăng
kích thước miền bức xạ tác dụng.
Các giá trị giới hạn của τ và α thường là ax m in0 , 95; 10mτ α≈ ≈ ° .
Lưu Thị Hoa Linh 12
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
1.4. Đặc điểm kết cấu anten loga chu kỳ và phương pháp tính toán
Để đảm bảo đồ thị phương hướng của anten trong mặt phẳng thẳng đứng
không biến đổi khi thay đổi tần số công tác, anten được đặt nghiêng một góc ∆ so
với mặt đất, sao cho độ cao tương đối của mỗi chấn tử so với mặt đất là đại lượng
không đổi:
1 2
1 2
n
n
HH H
λ λ= = =L λ (1.9)
Khi đó độ cao tương ứng của miền bức xạ tác dụng – di chuyển dọc theo
anten khi biến đổi tần số công tác – không thay đổi và do đó hướng bức xạ cực đại
trong mặt phẳng thẳng đứng sẽ không biến đổi.
Nếu coi mỗi đoạn dây truyền sóng giữa hai chấn tử là một mạng 4 cực, còn
mỗi chấn tử chấn tử tương đương với một trở kháng có giá trị bằng tổng trở vào
của chấn tử đó (có kể đến ảnh hưởng tương hỗ với các phần tử hệ thống), ta sẽ có
sơ đồ tương đương anten loga chu kỳ được vẽ như hình 1.5.
Theo hình 1.5 thì quan hệ điện áp ở cửa các tầng có thể được viết như sau:
Tầng I:
1
0
1 2
1
V
R V
V V
V V V
⎧ =⎨ = =⎩ (điện áp của nguồn cấp điện)
Tầng II:
2 1
1
2 3
2
V R
R V
V V V
V V V
⎧ = =⎨ = =⎩
Tầng N:
1
1
N N
V R N
N
R N
V V V
V V
−
−⎧ = =⎨ =⎩
Kí hiệu ,nV RnI I và là các dòng
điện và điện áp cửa vào và cửa ra của
tầng thứ n, còn là điện áp ra của
tầng thứ n.
,nV RV V
n
n
n
nV
Phương trình mạch điện đối với mạng tứ cực thứ n được viết như sau:
n n n n
V VV V VR R
n n n n
R RV V RR R
I y V y V
I y V y V
= +
= + (1.10)
Trong đó:
0nV
n
n V
VV n
V V
Iy
V =
= là dẫn nạp vào của đoạn dây truyền sóng khi đầu kia nối tắt
Lưu Thị Hoa Linh 13
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
0nR
n
n R
VR n
V V
Iy
V =
= là dẫn nạp truyền đạt của đoạn dây truyền sóng khi đầu kia nối tắt
Từ lý thuyết đường dây, ta xác định được các dẫn nạp vào và dẫn nạp truyền đạt:
( )
( )
1 cth
1 cosech
n n
VV RR
o
n n
VR RV
o
y y d
y y d
γρ
γρ
= =
= = −
(1.11)
Với:
oρ là trở kháng của đường dây,
d là độ dài của đoạn dây truyền sóng,
γ là hằng số truyền lan phức.
VN
Vn
NIII
a)
1I
1V 1Z∼ 0
V
N-1 2V 2Z
2I 2NI −
N 1NV − 1NZ − NV NZ
NI1NI −
II
I
b)
nI
nV
n
rI
n
n
vI
1n n
R VV V V
+
n= =
n
VV
n
RV
c)
Hình 1.5: Sơ đồ nguyên lý của anten
Lưu Thị Hoa Linh 14
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Từ (1.7) ta có phương trình mạch điện đối với các tầng:
Tầng 1:
1 1 1 1 1 1 1
0 1
1 1 1 1 1 1 1
0 1
V VV V VR R VV VR
R RV V RR R RV RR
I y V y V y V y V
I y V y V y V y V
= + = +
= + = +
Tầng 2:
2 2 2 2 2 2 2
1 2
2 2 2 2 2 2 2
1 2
V VV V VR R VV VR
R RV V RR r RV RR
I y V y V y V y V
I y V y V y V y V
= + = +
= + = +
……………………………………
Tầng N: 1
1
N N N N N N N
V VV V VR R VV N VR N
N N N N N N N
R RV V RR R RV N RR N
I y V y V y V y V
I y V y V y V y V
−
−
= + = +
= + = +
(1.12)
Tại các nút, ta có phương trình mạch điện được viết như sau:
( )
(
2 1
1
3 2
2
2
........................
V R
V R
N
n
)
I I I
I I I
I I
= − +
= − +
= −
(1.13)
Từ việc xác định dòng điện trong các chấn tử theo các công thức (1.10) và
(1.11), ta có thể đưa anten loga-chu kỳ về mô hình đơn giản gồm các chấn tử có
độ dài thay đổi đặt song song cách nhau những khoảng cách nhất định dọc theo
trục z ở các vị trí có tọa độ (hình 1.3). Mỗi chấn tử được tiếp điện bởi một
nguồn riêng biệt có sức điện động .Các kích thước và tọa độ được xác
định khi cho trước các thông số của kết cấu như chu kỳ
nl
nz
nV nl nz
τ các góc mở α .
Lưu Thị Hoa Linh 15
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
nz
Nz
0θ =
N
x
n
ln 2
1 l2
l1 0
z2z
1z
Hình 1.6: Mô hình đơn giản của anten lôga chukỳ
Hệ phương trình Kirchhoff đối với hệ thống N chấn tử ghép khi có tính đến
ảnh hưởng tương hỗ của các phân tử được viết dưới dạng:
11 12 1 1 1
21 22 2 2 2
1 2
...
...
... ... ... ... ... ...
...
N
N
N N NN N N
Z Z Z I V
Z Z Z I V
Z Z Z I V
⎡ ⎤ ⎡ ⎤ ⎡ ⎤⎢ ⎥ ⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥× =⎢ ⎥ ⎢ ⎥ ⎢ ⎥⎢ ⎥ ⎢ ⎥ ⎢ ⎥⎣ ⎦ ⎣ ⎦ ⎣ ⎦
(1.14)
Dựa vào các phương trình (1.12), (1.13), và (1.14) ta có các bước tính toán
đối với anten lôga – chu kỳ như sau:
Bước 1: Thay (1.12), (1.13) vào (1.14) sẽ nhận được hệ mới gồm N phương
trình. Giải hệ phương trình này được N nghiệm V V . Thay các nghiệm này
vào (1.10), (1.11) sẽ xác định được dòng điện trong các chấn tử:
1 2, ,..., NV
1 2, ,..., NI I I .
Tính trở kháng tương hỗ theo công thức: ij ij ijZ R iX= + trong đó ijR và ijX
xác định theo công thức sau:
Lưu Thị Hoa Linh 16
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
( ) ( )
1 ij ij2
ij 1 2
1 2
2
ij -2
ij ij ij ij ij
2 230 sin .sin {[sin sin
2 2
2sin .cos . ] . .sin cos sin sin
2
s2.sin .cos
2
j
j
l i i
z z
ji
l
zi
x y
i
l lS H S HklklR kr
r r
S HklkR S d S
R
klkR
ρ θ ϕ θ ϕ
−
−
⎛ ⎞ ⎛+ + + −⎜ ⎟ ⎜⎡ ⎤⎛ ⎞⎛ ⎞ ⎝ ⎠ ⎝= − +⎢ ⎥⎜ ⎟⎜ ⎟⎝ ⎠ ⎝ ⎠⎣ ⎦
+⎛ ⎞ ⎡ ⎤− +⎜ ⎟ ⎣ ⎦⎝ ⎠
⎛ ⎞+ −⎜ ⎟⎝ ⎠
∫ kr
⎞⎟⎠
+ +
1 2
ij
1 2
in sin .cos }.sin d
2
jlkr kr k S S
r r
θ⎡ ⎤ ⎛ ⎞− −⎜ ⎟⎢ ⎥ ⎝ ⎠⎣ ⎦
(1.15)
Phần ảo của trở kháng tương hỗ:
( ) ( )
1 ij ij2
ij 1 2
1 2
2
ij -2
ij ij ij ij ij
2 230 sin .sin {[cos cos
2 2
2cos .cos . ] . .sin cos sin sin
2
cos .cos
22.
j
j
l i i
z z
ji
l
zi
x y
i
l lS H S HklklX kr
r r
S HklkR S d S
R
klkR
R
ρ θ ϕ θ ϕ
−
−
⎛ ⎞ ⎛+ + + −⎜ ⎟ ⎜⎡ ⎤⎛ ⎞⎛ ⎞ ⎝ ⎠ ⎝= − +⎢ ⎥⎜ ⎟⎜ ⎟⎝ ⎠ ⎝ ⎠⎣ ⎦
+⎛ ⎞ ⎡ ⎤− +⎜ ⎟ ⎣ ⎦⎝ ⎠
⎛ ⎞⎜ ⎟⎝ ⎠+ −
∫ kr
⎞⎟⎠
+ +
1 2
ij
1 2
cos cos .cos }.sin d
2
jlkr kr k S S
r r
θ
⎡ ⎤⎢ ⎥ ⎛ ⎞⎢ ⎥− −⎜ ⎟⎢ ⎥ ⎝ ⎠⎢ ⎥⎣ ⎦
(1.16)
Ở đây:
( )
( )
( )
( )
ij ij
ij ij
1/ 22
ij
1/ 222
ij
1/ 222
1 ij
1/ 222
2 ij
.sin .cos
sin sin
/ 2
/ 2
x
y
x y
z
z i
z i
x x y y z z
S S
S S
S d S
R S H
r S H l
r S H l
S S i S i S i
θ ϕ
θ ϕ
ρ
ρ
ρ
ρ
=
=
⎡ ⎤= + +⎢ ⎥⎣ ⎦
⎡ ⎤= + +⎢ ⎥⎣ ⎦
⎡ ⎤= + + +⎢ ⎥⎣ ⎦
⎡ ⎤= + + −⎢ ⎥⎣ ⎦
= − +
Bước 2: theo giá trị các dòng đã tính được, ta tìm được hàm phương hướng
của anten trong hai mặt phẳng chính theo công thức:
Lưu Thị Hoa Linh 17
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Mặt phẳng H (mặt phẳng từ trường, yOz):
( ) ( ) Hos
1
1 cos
21
sin
2
n
n
N
n ikz cH
n
n
kl
f I kl e
θθ
−
= −∑ (1.17)
Mặt phẳng E (mặt phẳng điện trường, xOz):
( ) ( ) Eos
E1
cos sin cos
2 21
sin os
2
n
Hn n
N
n ikz cE
n
n
kl kl
f I kl c
e θ
θ
θ
θ
⎛ ⎞ −⎜ ⎟⎝ ⎠= −∑ (1.18)
Các góc ,H Eθ θ là góc hợp bởi hướng khảo sát và trục Oz trong mỗi mặt
phẳng E và H. Do cách mắc chéo nhau nên dòng điện trong 2 chấn tử kề nhau có
dấu ngược nhau, vì thế có số hạng (-1)n trong công thức.
Hệ số định hướng được xác định theo công thức:
( )max 2 2
0 0
4
, sin d dH E H H Em
D
F
π π
π
θ θ θ θ θ
=
∫ ∫
(1.19)
Trong đó ( ) ( ) ( )2 2, 2H E HmF f f Eθ θ θ= + θ là tổng bình phương hàm phương hướng
Bước 3: tính tổng trở vào của các chấn tử:
n
n
n
VZ
I
= (1.20)
Và tổng trở của anten:
1 1
1 1
1 0
VA
VV VR
V VZ
I y V y V
= = + 1 (1.21)
Bài toán xác định các thông số tối ưu của anten lôga – chu kỳ được giải
quyết bằng cách lặp lại nhiều lần các bước 1 và 2, đến chừng nào đạt được các chỉ
tiêu chất lượng tốt nhất (tùy thuộc vào hệ số định hướng).
Lưu Thị Hoa Linh 18
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
CHƯƠNG 2: THUẬT TOÁN DI TRUYỀN
2.1. Giới thiệu về thuật toán di truyền.
Trong thiên nhiên, mỗi cá thể có các tính chất và đặc điểm riêng biệt thể
hiện ra ngoài môi trường được gọi là kiểu hình. Các kiểu hình này được quyết
định bởi kiểu gen của chính cá thể đó. Kiểu gen rất đa dạng giữa các loài với nhau
và ngay trong cùng một loài, sự đa dạng này dẫn đến sự đa dạng về kiểu hình
trong một quần thể sinh học.
Thế giới tự nhiên phát triển và tiến hóa qua các thế hệ được là do tuân theo
quy luật chọn lọc tự nhiên: Cá thể nào có kiểu hình thích nghi với môi trường hơn
thì có khả năng tồn tại cao hơn và truyền kiểu gen mang sự thích nghi môi trường
này cho thế hệ con cháu thông qua quá trình sinh sản, còn cá thể nào yếu hơn thì
sẽ bị loại bỏ.
Trong quá trình sinh sản có cả di truyền và biến dị nên con cái không hoàn
toàn giống cha mẹ. Kiểu gen tiến hóa và hoàn thiện để thế hệ tiếp theo có thể phù
hợp với sự thay đổi của môi trường sống. Quá trình tiến hóa này được lặp lại
nhiều lần.
Các kỹ sư có thể sử dụng khoa học tự nhiên của sự thích nghi để thiết kế các
sản phẩm tốt hơn. Phương pháp gradient và phương pháp tìm kiếm ngẫu nhiên hội
tụ khá nhanh chóng tới cực tiểu khi thuật toán gần tới cực tiểu đó. Vì thế với
những bài toán có các thông số nhỏ biến thiên liên tục, và không gian nghiệm nhỏ
thì các phương pháp này tỏ ra có hiệu quả. Nếu số lượng thông số tăng lên thì
chất lượng của lời giải phụ thuộc vào dự đoán ban đầu. Dự đoán ban đầu rơi vào
vùng không gian chất lượng thấp thì việc tìm kiếm lời giải tốt sẽ rất khó khăn.
Trong khi thiết kế và tổng hợp của anten, mục đích của việc tìm kiếm cấu
trúc bức xạ đáp ứng được các tiêu chuẩn chất lượng như độ tăng ích, mức cực đại
phụ lớn nhất, độ rộng chùm, trở kháng đầu vào, và kích thước vật lý. Ngoại trừ
các cấu trúc anten đơn giản nhất, tất cả anten khác đều có rất nhiều các biến số
thiết kế tác động đến chất lượng anten. Các cấu trúc khác nhau của anten cũng
gây khó khăn cho việc dự đoán ban đầu cho từng thiết kế. Với những vấn đề này
thì thuật toán di truyền được đưa ra nhằm giải quyết vấn đề tìm kiếm toàn cục cho
lời giải thỏa mãn một số tiêu chuẩn chất lượng đặc biệt.
Lưu Thị Hoa Linh 19
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
2.2. Khái quát về thuật toán di truyền và các bước tiến hành.
2.2.1. Khái quát về thuật toán di truyền
Thuật toán di truyền bắt chước sự tiến hóa và tổ hợp gen trong tự nhiên nên
gen là thành phần cơ bản của thuật toán. Thuật toán này mã hóa mỗi thông số trở
thành một chuỗi nhị phân, được gọi là gen, và tập hợp các gen là một nhiễm sắc
thể. Các nhiễm sắc thể lại được đánh giá bằng các hàm mục tiêu, trải qua quá
trình chọn lọc tự nhiên, lai ghép và biến dị, để đạt được lời giải tối ưu cuối cùng.
Các nhiễm sắc thể đó được xếp loại từ phù hợp nhất đến ít phù hợp nhất, tùy
thuộc vào hàm mục tiêu tương ứng. Phù hợp hay không là do người lập trình quy
định. Những nhiễm sắc thể không phù hợp bị loại bỏ, để lại những nhiễm sắc thể
tốt nhất ban đầu.
Các cá thể được duy trì sẽ trở thành cha mẹ, và trao đổi vật chất di truyền để
tạo ra hai con cái mới, đủ để bù lại những nhiễm sắc thể đã bị loại bỏ. Vì vậy sau
mỗi vòng lặp tổng số nhiễm sắc thể còn lại là một hằng số. Sự biến dị dẫn đến sự
thay đổi ngẫu nhiên nhỏ trong một nhiễm sắc thể.
Hàm mục tiêu lại được đánh giá cho thế hệ con cái và các nhiễm sắc thể biến
dị, sau đó cả quá trình được lặp lại. Thuật toán chỉ dừng lại sau một lượng xác
định vòng lặp hoặc khi thu được một lời giải có thể chấp nhận được.
2.2.2. Các bước tiến hành
Hình 2.1 là lưu đồ của thuật toán di truyền, bắt đầu bằng việc định nghĩa
nhiễm sắc thể như một hệ anten có các giá trị của thông số được tối ưu hóa. Nếu
nhiễm sắc thể có Npar thông số (một bài toán tối ưu có N chiều), được cho bởi
ar1 2 3
, , , ,
pN
p p p pL , nhiễm sắc thể sẽ được viết dưới dạng:
ar1 2 3 pN
N ST p p p p⎡ ⎤= ⎣ ⎦L (2.1)
Mỗi nhiễm sắc thể có một hàm mục tiêu, được tìm bởi việc đánh giá một
hàm f tại vị trí . Biểu diễn lại hàm mục tiêu như sau: ( ar1 2 3, , , , pNp p p pL )
( )a r1 2 3c o s t= , , , , pNf p p p pL (2.2)
Lưu Thị Hoa Linh 20
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Các thông số np (với n = 1, 2, …, Npar) có thể rời rạc hoặc liên tục. Ta giới
hạn các thông số liên tục bằng cách mã hóa chúng thành các chuỗi nhị phân:
(2.3) [ ] 1W
1
2
nL
m
n
m
q b m −
=
= ∑ Q
Trong đó:
nq = giá trị lượng tử hóa của np
nL = số mức lượng tử cho nq
Wb = hệ anten bao gồm chuỗi nhị phân biểu diễn lại bằng nq
Q = mức lượng tử lớn nhất hoặc một nửa giá trị lớn nhất có thể của nq
thiết lập mã hóa hoặc giải
mã các thông số
tạo M nhiễm sắc thể ngẫu nhiên
Đánh giá hàm mục tiêu
cho các nhiễm sắc thể
Những nhiễm
sắc thể còn lại
lai ghép
loại bỏ những
nhiễm sắc thể không
phù hợp
biến dị
xếp hạng các nhiễm sắc thể
Dừng
chưa
Hoàn thành?
rồi
Hình 2.1: Lưu đồ của thuật toán di truyền
Lưu Thị Hoa Linh 21
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Các thông số được mã hóa nhị phân ( ) có thể chỉ thể hiện một vài giá trị của nq
np . Ví dụ như np thể hiện 8 giá trị của suất điện trở, có thể biểu diễn như sau: nq
Ω==
Ω==
Ω==
800111
200001
100000
n
n
n
q
q
q
LLL
Các thuật toán di truyền ở đây chỉ làm việc với các thông số đã được mã hóa,
mà không làm việc với chính các thông số đó. Mỗi khi hàm mục tiêu được đánh
giá, thì các nhiễm sắc thể phải được giải mã đầu tiên. Phương trình (2.3) cho ta
thấy một hệ anten có các thông số được lượng tử. Nhiễm sắc thể này có tổng cộng
argbit pbit pN N N= × bit.
Sau khi đặt ra sơ đồ lập mã và giải mã các thông số, một danh sách các
nhiễm sắc thể ngẫu nhiên được tạo ra. Mỗi nhiễm sắc thể có một hàm
mục tiêu kết hợp, được tính toán từ phương trình (2.2).
8NSTN =
Bước tiếp theo trong thuật toán là xếp hạng các nhiễm sắc thể từ tốt nhất
đến tồi nhất theo mục tiêu của nhiễm sắc thể đó. Tại đây, các nhiễm sắc thể không
được chấp nhận bị loại bỏ. Chấp nhận được hay không là do người lập trình định
nghĩa tùy thuộc vào chất lượng yêu cầu. Nếu x nhiễm sắc thể đứng đầu được chọn
(x là số chẵn) thì bị loại bỏ. NSTN − x
Ta loại bỏ 50% thì chỉ còn nhiễm sắc thể cho việc lai ghép. Hai
nhiễm sắc thể bất kỳ đều có thể lai ghép với nhau. Có thể ghép đôi nhiễm sắc thể
đứng đầu và nhiễm sắc thể đứng cuối bảng, hoặc ghép đôi ngẫu nhiên: nhiễm sắc
thể số 1 với , NST số 2 với
/ 2NSTN
/ 2NSTN / 2 1NSTN − . v.v. Với mỗi cặp ghép đôi, hai con
cái mới được tạo ra từ việc trao đổi vật chất di truyền của cha mẹ: chọn một điểm
tương giao chéo ngẫu nhiên, các cặp bit nhị phân ở bên phải của điểm tương giao
chéo được đổi để định dạng lứa con cái.
Lưu Thị Hoa Linh 22
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Nếu điểm tương giao chéo ngẫu nhiên giữa các bit của nhiễm sắc thể được
chọn như trong ví dụ sau, các nhiễm sắc thể mới được định dạng bằng cách:
Cha mẹ thứ nhất (NST số 6)
Cha mẹ thứ hai (NST số 2)
Con thứ nhất
Con thứ hai
}
}
10101 0010101
11100 1100101
101011100101
111000010101
4243
4243
Sau khi cho nhiễm sắc thể ghép cặp và lai ghép, danh sách
cặp cha mẹ và con cái lời giải cho tổng số nhiễm sắc thể (giống
lượng nhiễm sắc thể ban đầu)
/ 2NSTN / 2NSTN
/ 2NSTN NSTN
Tại thời điểm này, biến dị ngẫu nhiên làm cho một tỉ lệ nhỏ các bit trong
danh sách nhiễm sắc thể có sự thay đổi - biến dị là việc thay 1 thành 0 hoặc
ngược lại. Một bit được lựa chọn ngẫu nhiên cho biến dị từ tổng số gbit NSTN N×
bit trong tất cả các nhiễm sắc thể. Việc tăng lượng bit bị biến dị làm tăng sự tự
do của thuật toán ra ngoài phạm vi hiện thời của không gian thông số. Sự tự do
này trở nên quan trọng hơn khi thuật toán sắp hội tụ thành lời giải cụ thể. Sự
biến dị không xảy ra ở vòng lặp cuối cùng, vì lúc này các nhiễm sắc thể đã gần
như giống nhau.
Sau khi biến dị, hàm mục tiêu tương ứng của các thế hệ con và các nhiễm
sắc thể bị biến dị được tính toán, và quá trình lại được lặp lại. Số thế hệ phát triển
phụ thuộc vào lời giải có phù hợp không, hoặc số vòng lặp đã thực hiện. Sau một
thời gian, tất cả nhiễm sắc thể và hàm mục tiêu tương ứng trở nên giống nhau,
ngoại trừ trường hợp biến dị. Khi đó thuật toán cần được dừng lại.
2.2.3. Thuật toán di truyền cho tối ưu và thiết kế anten
Hình 2.2 cho ta thấy biểu đồ các bước thực hiện của thuật toán di truyền.
Hầu hết tất cả các bước là chung cho tất cả các thuật toán di truyền. Những khác
biệt đáng chú ý là phương pháp ước lượng pha anten momen ở đầu vòng lặp.
Thuật toán di truyền bắt đầu tìm kiếm với một quần thể lớn của những cá
thể được tạo ngẫu nhiên. Sự ngẫu nhiên này cho ra một quần thể đa dạng gồm
những lời giải có thể đại diện được cho một nhóm lớn của toàn bộ không gian lời
Lưu Thị Hoa Linh 23
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
giải. Mỗi cá thể trong quần thể đại diện cho một thiết kế anten. Những tính chất
của mỗi thiết kế là một lời giải sử dụng phương pháp của momen và sau đó định
một cá thể phù hợp nhất dựa vào thiết kế trùng lặp bao nhiêu so với mong muốn.
Bắt đầu với lựa chọn quần thể
ngẫu nhiên của anten
Kết thúc
Áp dụng toán tử di truyền
để tạo thế hệ anten mới
Sắp xếp anten theo
chất lượng hàm mục tiêu
Kiểm tra tiêu chuẩn
cuối cùng
Tái tạo các anten bằng
phương pháp momen
Hình 2.2: Sơ đồ khối cho thuật toán di truyền.
Bằng cách sử dụng: “cá thể phù hợp thì tồn tại”, các cá thể phù hợp hơn
có khả năng được lựa chọn lớn hơn. Với ý tưởng là việc trao đổi thông tin di
truyền giữa hai lời giải tốt thì sẽ tạo ra một lời giải tốt hơn. Trong quần thể cũng
có một vài sự thay đổi nhỏ làm cho cá thể có thể bị biến dị. Thay đổi ngẫu nhiên
của một tỉ lệ nhỏ của quần thể giúp duy trì thế giới tự nhiên trong việc nghiên
cứu. Những quần thể mới được tạo ra và được đánh giá cho nhiều thế hệ cho tới
khi quần thể tập trung lại thành một lời giải giống nhau hoặc một vài tiêu chí
dừng khác được gặp.
Lưu Thị Hoa Linh 24
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
2.2.4. Hàm mục tiêu và tỉ lệ phù hợp
Mục đích của quá trình thiết kế là phát triển một loại anten đáp ứng một vài
yêu cầu đặc tính chất lượng. Vài đặc tính định nghĩa chất lượng anten là: cực đại
phụ, độ rộng chùm, tỉ lệ trước – sau, kích cỡ, tăng ích và trở kháng đầu vào. Chất
lượng thiết kế được thể hiện theo toán học bằng hàm mục tiêu.
Giả sử hàm mục tiêu cho anten thiết kế x một hệ số tăng ích G và trừ đi
phần sai lệch giữa phần thực của trở kháng với 75Ω và phần ảo khác 0. Trong
khóa luận này không thực hiện tối ưu phần ảo của trở kháng.
object( ) a ( ) b 75 Re( ( ))x G x Z x= − − (2.4)
Những hằng số a, b, c là ảnh hưởng của mỗi số hạng tới toàn bộ hàm mục tiêu.
2.2.5. Sự chọn lọc
Sau khi đánh giá quần thể, ta chọn các cá thể tạo thế hệ tiếp theo. Trong
các cá thể được lựa chọn, có một số được sao chép đơn giản sang quần thể mới.
Số khác được lai ghép với cá thể nào đó và con cháu của chúng trở thành một
phần của quần thể mới. Trong suốt quá trình lựa chọn, cá thể i được giả định là có
thể có Pi cách lựa chọn, dựa vào tỉ lệ phù hợp của cá thể trên tổng số phù hợp của
N cá thể trong quần thể, Pi được tính như sau:
∑
=
= N
k
k
i
i
f
fP
1
(2.5)
Dạng chuẩn của sự lựa chọn (thường được gọi là bánh xe xổ số) cấp một
phần của một bánh xe cho mỗi cá thể vì thế vào mọi thời điểm vòng quay này
luôn được thực hiện. Nếu bánh xe quay được N lần thì cá thể có phù hợp trung
bình có thể được chon khi mà cá thể đó có độ phù hợp gấp 2 lần giá trị trung bình;
tuy nhiên N mẫu thường là quá ít mẫu để phân loại lựa chọn như mong đợi. Vài cá
thể có thể được chọn thường xuyên hơn trong khi một vài cá thể trung bình thì
không bao giờ được chọn.
Lưu Thị Hoa Linh 25
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Một trong các phối hợp lựa chọn Brindle [3] là chọn cá thể dựa trên cơ sở
sự có mặt được mong đợi của chúng trong quần thể N cá thể.
i ie P N= (2.6)
Mỗi cá thể được đảm bảo sẽ có ít nhất (làm tròn tới dưới số nguyên gần nhất)
mẫu đại diện trong quần thể mới. Phần của sau đó được áp dụng kiểu tương tự
bánh xe xổ số để lấp đầy các phần còn lại trong quần thể.
ie
ie
2.2.6. Mô tả anten như nhiễm sắc thể
Cách thức biểu diễn một cá thể trong thuật toán di truyền như thế nào là rất
quan trọng. Trong thuật toán di truyền, các nhiễm sắc thể là sự thể hiện tóm tắt lới
giải của bài toán. Vì thế mỗi lời giải riêng của quần thể được thể hiện bởi một
nhiễm sắc thể. Mỗi loại anten có các thông số đặc trưng cho kết cấu khác nhau.
Những thông số này được mã hóa nhị phân rồi ghép lại với nhau thành xâu đơn.
Xét ví dụ minh họa:
L3 L1 L2
S1 S2
Hình 2.3 : Một anten có 3 chấn tử
Một anten 3 chấn tử như hình 2.3, các chấn tử có chiều dài là L1, L2, L3,
khoảng cách giữa chúng là S1, S2, S3. Các thành phần này được miêu tả như là
một nhiễm sắc thể có giá trị nhị phân như hình 2.4 rồi sau đó ghép lại với nhau
thành xâu đơn:
e n c o d e d 1 2 LB
L l l l= L (2.7)
Lưu Thị Hoa Linh 26
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
1010010000000111 11000111 1100100011001011
Lmin Lmax Smin Smax
L1 S1
Chuyển thành giá trị
nhị phân
Các giá trị rời rạc được mã hóa 8 bit
Nhiễm sắc thể
Kết hợp tất cả các giá trị để tạo thành nhiễm sắc thể
L1 S1 L2 S2 L3 S4
00000111
0 1 2 3 255
11001011
0 1 2 3 255
10111000
Hình 2.4: Biểu diễn anten như một nhiễm sắc thể
Trong đó BL là số bit sử dụng để biễu diễn chiều dài L. BL được tính bằng
công thức sau:
max min
2
res
log ( 1)L
L LB
L
−= + (2.8)
Với Lmin, Lmax là khoảng mà chiều dài chấn tử có thể thay đổi được và Lres
là độ phân giải của sự rời rạc hóa. Nếu BL không nguyên thì nó sẽ được làm tròn
lên giá trị nguyên kế tiếp. Biễu diễn nhị phân của chiều dài được giải mã như sau:
LN
k-1
min res
k=1
2kL L L l= +∑ (2.9)
Nhiễm sắc thể được định dạng bằng cách ràng buộc tất cả các thông số mã
hóa với nhau thành một chuỗi nhị phân đơn. Kích thước của không gian tìm kiếm
được định nghĩa bằng độ dài của nhiễm sắc thể. Sự hội tụ của thuật toán di truyền
thành lời giải nhanh hơn trong không gian tìm kiếm nhỏ hơn, nên cần có sự cân
bằng giữa khoảng tham số và độ phân giải của các thông số với tốc độ hội tụ.
Lưu Thị Hoa Linh 27
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
2.2.7. Các toán tử di truyền
Thông tin di truyền được nhân giống từ một thế hệ sang thế hệ sau bằng
toán tử di truyền. Có 3 cách cơ bản của toán tử di truyền là: sao chép
(reproduction), tương giao chéo (crossover) và biến dị (mutation). Sự sao chép, là
toán tử đơn giản nhất, sao chép trực tiếp một cá thể từ thế hệ cũ sang thế hệ mới
(giống như quá trình sinh sản vô tính).
Toán tử tương giao chéo chỉ phức tạp hơn sao chép một chút: thông tin di
truyền giữa 2 nhiễm sắc thể được đổi chỗ bằng việc cắt những nhiễm sắc thể tại
những nơi được chọn ngẫu nhiên và đổi đoạn cuối của nhiễm sắc thể để tạo nên
những nhiễm sắc thể mới. Hình 2.5 minh họa quá trình này, nêu ra một cách đơn
giản cho nhiễm sắc thể trao đổi thông tin di truyền trong việc tìm kiếm nhiễm sắc
thể tốt hơn. Tỉ lệ của các cá thể được chọn để tham gia vào tương giao chéo trong
suốt mỗi thế hệ được ghi rõ bởi thông số (xác suất tương giao chéo). pcrossover
Cha mẹ A Cha mẹ B
Sau vài thế hệ, quần thể dường như bao gồm nhiều bản sao đồng nhất của
một cá thể rất phù hợp. Một cách tổng quát, điều này là tốt vì con cái của cá thể
này cũng phù hợp; tuy nhiên, lứa con cái rất giống cha mẹ chúng có thể dẫn đến
hạn chế trong việc tìm kiếm của không gian lời giải. Để tạo sự đa dạng trong quần
thể, toán tử biến dị được đưa vào. Nếu một bit được chọn để gây biến dị, giá trị
của bit đó chỉ đơn giản là thay đổi từ 0 thành 1 hoặc từ 1 thành 0. Áp dụng biến dị
Con B
01011 110001 00110 011001
01011 011001 00110 110001
thế hệ cũ
thế hệ mới
Con A
Hình 2.5 : Tương giao chéo
Lưu Thị Hoa Linh 28
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
vào một tỉ lệ nhỏ của quần thể, tính toàn cục của việc tìm kiếm vẫn được
đảm bảo.
crossoverp
Thuật toán di truyền có xác suất tự nhiên làm cho thu được các lời giải
khác nhau cho cùng một chương trình tính toán, trừ khi bài toán này đủ đơn giản
để tìm thấy lời giải tối ưu. Thuật toán di truyền hướng tới việc tìm kiếm lời giải
thích hợp và tối ưu nó. Ví dụ như một hàm phù hợp của anten bao gồm một số
hạng nhằm cực đại hóa tăng ích và số hạng khác để hướng tới trở kháng bằng
75Ω. Nếu ta lấy anten từ một vài lần chạy chương trình khác nhau, thì có thể một
vài anten sẽ có giá trị tăng ích cao với trở kháng đầu vào không như mong đợi,
trong khi đó, các anten khác sẽ có giá trị tăng ích nhỏ hơn nhưng giá trị trở kháng
đầu vào tốt hơn. Vì thế, ta cho chạy chương trình nhiều lần và lựa chọn anten
thích hợp nhất để ứng dụng.
Lưu Thị Hoa Linh 29
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
CHƯƠNG 3: ỨNG DỤNG THUẬT TOÁN DI TRUYỀN
Trong chương này tôi sẽ nêu chi tiết hơn về việc áp dụng thuật toán di
truyền trong tối ưu hóa anten lôga - chu kỳ. Dựa trên những lý thuyết về anten
lôga - chu kỳ và thuật toán di truyền, tôi đã mô phỏng các tính toán trên bằng
phần mềm MatLab.
3.1. Các tham số tối ưu cho anten lôga - chu kỳ.
Chương trình còn tính toán tối ưu anten lôga - chu kỳ các tham số : chu kỳ
kết cấu τ , góc α , hệ số định hướng và trở kháng vào.
Các tham số khác có thể được thay đổi trong file TTDT.m. Trong file này,
nhằm mục đích nghiên cứu tôi đã đặt các tham số như sau:
- tần số lớn nhất fmax = 600 (MHz),
- tần số nhỏ nhất fmin = 400 (MHz),
- tần số hoạt động F0 = 500 (MHZ),
- số chấn tử cuat anten là 6,
- bán kính chấn tử là 0.4 (mm),
- số anten thử nghiệm (số cá thể trong quần thể) là 8
- chu kỳ kết cấu τ giới hạn giữa 0.85 0.95÷
- góc α giới hạn giữa 10 20ο ο÷ .
3.2. Các bước thực hiện tính toán.
3.2.1. Khởi tạo ngẫu nhiên quần thể gồm một số nhiễm sắc thể.
Ta khởi tạo ngẫu nhiên quần thể bằng lệnh:
Gene = round(rand(row,cell));
Trong đó Gene là một ma trận ngẫu nhiên gồm các bit nhị phân có row
hàng, và cell cột. Mỗi hàng là một nhiễm sắc thể đại diện cho một cá thể trong
quần thể.
Lưu Thị Hoa Linh 30
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.2.2. Thực hiện mã hóa các tham số của từng anten.
Vì bài toán tối ưu anten lôga - chu kỳ có hai tham số là chu kỳ kết cấu τ ,
góc α nên số gen trong mỗi nhiễm sắc thể là 2. Mỗi tham số được mã hóa bằng
một chuỗi nhị phân tương ứng với một đoạn gen trong nhiễm sắc thể, với nhiễm
sắc thể thứ ic thì công việc này được thực hiện qua các lệnh
cost = Gene(ic,:);
for jb = 1:No_bit
ckkc(1,jb) = cost(1,jb);
goc_d(1,jb) = cost(1,jb+No_bit);
end
3.2.3. Giải mã các gen sang dạng thập phân và tính theo giới hạn.
Ta giải mã các gen sang dạng thập phân bằng hàm decode trong file
decode.m
function deco = decode(g_in,No_bit)
deco = 0;
for ic = 1:No_bit
deco = deco+g_in(1,ik)*2^(No_bit-ic);
end
Tham số vào là g_in cần được giải mã, No_bit là số bit mã hóa.
Các tham số sau khi được giải mã sang dạng thập phân thì cần được tính
dựa trên giới hạn ban đầu của chương trình:
chuky = a0+chuky*(b0-a0)/(2^No_bit);
alpha = a1+alpha*(b1-a1)/(2^No_bit);
Lưu Thị Hoa Linh 31
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.2.4. Thực hiện tính toán mỗi thiết kế anten.
Việc tính toán hàm mục tiêu được thực hiện trên file loga.m dựa vào hai
tham số τ và α vừa tìm được để tìm hàm mục tiêu của chúng theo các bước sau:
Bước 1: Tính chiều dài các chấn tử theo chu kỳ, góc và chiều dài của chấn
tử lớn nhất.
L(so_chan_tu) = lmax;
for n = so_chan_tu-1:-1:1
L(n) = L(n+1)*chky;
end
Bước 2: Tính tọa độ của mỗi chấn tử (tọa độ chính là khoảng cách giữa một
chấn tử nào đó với chấn tử đầu tiên). Tính khoảng cách giữa mỗi chấn tử. Khoảng
cách này dùng cho việc tính trở kháng tương hỗ.
x(1) = (L(1)/2)/tan(alfa);
for n = 2:so_chan_tu
x(n) = x(n-1)/chky;
end
for n = 1:so_chan_tu
for m = 1:so_chan_tu
if n ~ = m
k_cach(n,m) = abs(x(n)-x(m));
else k_cach(n,m) = 0;
end
end
end
Lưu Thị Hoa Linh 32
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Bước 3: Tính các dẫn nạp, trở kháng riêng, trở kháng tương hỗ của các
chấn tử theo các hàm DanNap(n, m ,k), trk_rieng(bk, L(n)) và trk_th(L(n), L(m),
k_cach(n,m)). Sau đó thế các giá trị này vào ma trận để tính các hệ số của ma trận.
Từ (1.10) và (1.11) ta thấy:
1 1 2
1
1 2( )
( j = 1,2 ; i,k = 1,2,...N)
k k k
j j k j k
i i
i
I y V y V
I I I
−
+
= +
= − +
Thay vào (1.12) chúng ta có :
1
1
1 2
1
1 1
11 12 1 21 1 22
1
( )
( )
N
ij j i
j
N
j j
ij i
j
N
j j j j
ij j j j j i
Z I V
Z I I V
Z y V y V y V y V V
i
=
+
=
+ +
+ −
=
=
− + =
− + + +
∑
∑
∑ =
jy
Các hệ số của ma trận bậc N mới sẽ có các hệ số là:
1 1
1 12 11 22 1 21
1
1 21 0
( )
. .
( i,j = 1,2,...,N)
j j j
ij ij ij ij ij
i i
a Z y Z y y Z
b Z y V
δ + +− += + + + +
=
Trong đó: 0 1 0i iNZ Z += =
1 0Nijy
+ =
Chương trình tính toán này sẽ được nêu trong phần phụ lục.
Bước 4: Tính trở kháng đầu vào ở dạng phức:
trovao = v(1) / dong(1);
Bước 5: Tính hệ số định hướng theo công thức (1.15), (1.16) và (1.17) qua
hàm bf_hfh (goc1, goc2) được nêu trong phần phụ lục.
Lưu Thị Hoa Linh 33
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.2.5. Loại bỏ những nhiễm sắc thể không thỏa mãn theo tiêu chí đặt ra.
Ta đặt ra tiêu chí là hàm mục tiêu phải lớn nên sẽ loại bỏ những nhiễm sắc
thể có hàm mục tiêu thấp nhất. Xếp hạng các nhiễm sắc thể theo hàm mục tiêu
tăng dần và loại bỏ 50%.
muctieu = - muctieu ;
[muctieu,ind] = sort(muctieu);
Gene = Gene(ind(1:round(row/2)),:);
3.2.6. Sử dụng các toán tử di truyền.
Chọn ngẫu nhiên điểm tương giao chéo và thay đổi các bit để tạo ra thế hệ
con mới (như lý thuyết phần 2.2.7) bằng các lệnh:
cross = ceil((cell-1)*rand(round(row/2),1));
for ic = 1:2:(round(round(row/2)))
Gene(round(row/2)+ic,1:cross) = Gene(ic,1:cross);
Gene(round(row/2)+ic,cross+1:cell) = Gene(ic+1,cross+1:cell);
Gene(round(row/2)+ic+1,1:cross) = Gene(ic+1,1:cross);
Gene(round(row/2)+ic+1,cross+1:cell) = Gene(ic,cross+1:cell);
end
3.2.7. Kiểm tra điều kiện dừng.
Nếu chưa đến vòng lặp cuối sử dụng toán tử biến dị để tạo ra biến dị ngẫu
nhiên trong quần thể. Ngược lại, tại vòng lặp cuối thì không thực hiện biến dị, lúc
này các nhiễm sắc thể đã gần như giống nhau, ta lấy nhiễm sắc thể đầu tiên để
giải mã, sau đó kết thúc quá trình tính toán và hiện kết quả ra màn hình Command
Window của phần mềm Matlab (bao gồm chu kỳ kết cấu, góc, hệ số định hướng,
và hàm mục tiêu).
Lưu Thị Hoa Linh 34
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.3. Một số kết quả.
3.3.1. Trường hợp 1: chỉ tối ưu hệ số định hướng, không tối ưu trở kháng.
3.3.1.1 Kết quả 1.
Ngoài các tham số chung như nêu ở mục 3.1. thì trong lần tính toán này sử
dụng 10 bit mã hóa, 8 vòng chạy. Ta thu được kết quả sau:
Bảng 1: sử dụng 10 bit mã hoá, 8 vòng chạy, không tối ưu phần thực trở kháng
Chu kỳ kết cấu (chuky) 0.8707
Góc (alpha) o12.6647
Hệ số định hướng (hsdh) 3.7520 dBi
Thời gian chạy chương trình khoảng 6 phút.
3.3.1.2 Kết quả 2.
Lần tính toán này sử dụng 10 bit mã hóa, tăng số vòng chạy lên 10.
Bảng 2: sử dụng 10 bit mã hoá, 10 vòng chạy, không tối ưu phần thực trở kháng
Chu kỳ kết cấu (chuky) 0.8911
Góc (alpha) o12.5290
Hệ số định hướng (hsdh) 3.8014 dBi
Nhận xét: So với kết quả đầu tiên thì hệ số định hướng tăng, chu kỳ kết cấu tăng,
nhưng vẫn nhỏ hơn trong kết quả 2. Vì tăng số vòng lặp nên thời gian tính toán
cũng lâu hơn, khoảng 7.5 phút.
Lưu Thị Hoa Linh 35
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.3.1.3 Kết quả 3.
Lần này ta tăng số bit mã hóa lên 12, và cho chạy chương trình 10 vòng.
Bảng 3: sử dụng 12 bit mã hoá, 10 vòng chạy, không tối ưu phần thực trở kháng
Chu kỳ kết cấu (chuky) 0.9013
Góc (alpha) o11.5152
Hệ số định hướng (hsdh) 3.9244 dBi
Nhận xét: Chương trình chạy lâu hơn, khoảng 9 phút tuy nhiên cho ta dết quả tốt
hơn, hệ số định hướng cao hơn 2 kết quả ở trên.
3.3.2. Trường hợp 2: Tối ưu hệ số định hướng với hệ số là 30, đồng thời tối ưu
phần thực của trở kháng.
3.3.2.1 Kết quả 4.
Vẫn giữ nguyên số bit mã hóa và số vòng chạy. Mục đích tối ưu là tối ưu
hệ số định hướng với hệ số là 30, có tối ưu phần thực của trở kháng.
Bảng 4: sử dụng 10 bit mã hoá, 8 vòng chạy, có tối ưu phần thực trở kháng
Chu kỳ kết cấu (chuky) 0.9305
Góc (alpha) o12.6258
Hệ số định hướng (hsdh) 3.6957 dBi
Phần thực trở kháng vào (real(trovao)) 33.1995 Ω
Nhận xét: Khi giữ nguyên số vòng chạy, số bit mã hóa (so với kết quả 1), thêm tối
ưu hóa phần thực của trở kháng vào thì hệ số định hướng có giảm đi một chút.
Phần thực của trở kháng vẫn còn khác khá nhiều so với 75Ω .
Lưu Thị Hoa Linh 36
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.3.2.2 Kết quả 5.
Lần tính toán này sử dụng 10 bit mã hóa, 10 vòng chạy.
Bảng 5: sử dụng 10 bit mã hoá, 10 vòng chạy, có tối ưu phần thực trở kháng
Chu kỳ kết cấu (chuky) 0.9391
Góc (alpha) o12.2016
Hệ số định hướng (hsdh) 3.7357 dBi
Phần thực trở kháng vào (real(trovao)) 48.2741Ω
Nhận xét: Trường hợp này thời gian tính toán dài hơn, phần thực của trở kháng
gần với hơn, hệ số định hướng cũng cao hơn trường hợp trên. 75Ω
3.3.2.3 Kết quả 6.
Trong lần tính toán này sử dụng 12 bit mã hóa, 10 vòng chạy.
Bảng 6: sử dụng 12 bit mã hoá, 10 vòng chạy, có tối ưu phần thực trở kháng
Chu kỳ kết cấu (chuky) 0.9413
Góc (alpha) o11.1235
Hệ số định hướng (hsdh) 3.8652 dBi
Phần thực trở kháng vào (real(trovao)) 72.6374Ω
Nhận xét: Trường hợp này công việc tính toán diễn ra lâu nhất, nhưng lại cho ra
kết quả chính xác nhất chu kỳ kết cấu gần với giới hạn trên 0.95, góc gần với giới
hạn dưới 10 , phần thực của trở kháng gần với o 75Ω .
3.4. Nhận xét tổng quát các kết quả đạt được.
Phương pháp sử dụng thuật toán di truyền mang tính ngẫu nhiên rất lớn do
việc tạo ra một quần thể ngẫu nhiên ngay từ thế hệ đầu và việc chọn điểm biến dị
Lưu Thị Hoa Linh 37
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
trong các thế hệ con tiếp theo cũng dược thực hiện một cách ngẫu nhiên. Do đó
thuật toán di truyền chậm hơn các phương pháp tối ưu cục bộ khác.
Thông thường khi các tham số đầu vào là không thay đổi thì qua những lần
mô phỏng khác nhau, kết quả bài toán thu được là khác nhau. Với mỗi bài toán
đặt ra thường thì ta phải chạy chương trính tính toán nhiều lần thì mới có thể tìm
được kết quả tốt nhất.
Qua các ví dụ trên ta có thể đi đến một số nhận xét như sau:
- Số bit biểu diễn một tham số và số vòng lặp: khi mã hóa với số bit
càng lớn hoặc đặt số vòng lặp lớn thì độ chính xác càng cao nhưng tốc độ hội tụ
chậm và thời gian tính toán sẽ càng cao.
- Có nhiều cách để chọn lọc nững cá thể phù hợp với tiêu chí đạt ra.
Trong khoá luận sự chọn lọc được thực hiện bằng cách sắp xếp giá trị của hàm
mục tiêu từ lớn tới nhỏ và loại bỏ đi 50% các cá thể có hàm mục tiêu thấp nhất.
- Sự hội tụ của thuật toán: Việc xác định khi nào bài toán tối ưu là rất
khó vì không thể điều khiển dừng thuật toán trong khi nó đang chạy. Thuật toán di
truyền làm cho các nhiễm sắc thể sau vòng lặp cuối cùng trở nên khá giống nhau
trừ những gen bị biến dị. Tùy theo từng tình huống cụ thể thì sự hội tụ có thể
nhanh hay chậm. Khi bài toán gặp khó khăn trong hội tụ thì có một số cách giải
quyết sau: tăng số lượng nhiễm sắc thể, tăng số lần biến dị hoặc thêm vào một số
tham số mà bạn biết từ bản chất vật lý của bài toán…
3.5. Đánh giá kết quả và hướng phát triển.
Phương pháp sử dung thuật toán di truyền là một phương pháp khó, và có
thể áp dụng cho nhiều bài toán tối ưu các loại anten khác nhau. Nhưng trong thời
gian ngắn tôi đã đạt được một số kết quả sau: Nêu ra cơ sở lý thuyết cho việc chế
tạo anten lôga – chu kỳ; lý thuyết về thuật toán di truyền và phương pháp áp dụng
cho việc thiết kế anten; và nghiên cứu việc áp dụng thuật toán di truyền cho bài
toán tối ưu thiết kế anten lôga- chu kỳ.
Trong thời gian nghiên cứu tiếp theo tôi sẽ có thể hoàn thiện tốt hơn phần
mềm tối ưu cho các loại anten thường được sử dụng trong thực tế nhằm chế tạo ra
những anten chất lượng tốt. Đồng thời chế tạo ra mẫu anten đã tối ưu theo thuật
toán, kiểm nghiệm và đo đạc mẫu để chứng minh thuật toán.
Lưu Thị Hoa Linh 38
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Tài liệu tham khảo
[1]. Phan Anh, Lý thuyết và kỹ thuật anten, Nhà xuất bản Khoa Học Kỹ Thuật
[2]. Nguyễn Hoàng Hải, Nguyễn Khắc Kiểm, Lập trình MatLab, NXB Khoa học
kỹ thuật, 2003.
[3]. A.Brindle, Genetic algorithms for function optimization, Doctoral
dissestation, Univ.Alberta, Edmonton, Canada, 1981.
[4]. Eric A.Jones, William T.joines, Design of Yagi-Uda antennas using genetic
algorithms, IEEE transaction on antennas and propagation, Vol.45, No.9,
September 1997.
[5]. Randy L.Haupt, An introduction to Genetic Algorithms for Electromagnetics,
IEEE antennas and propagation, Vol.37, No.2, April 1995.
Lưu Thị Hoa Linh 39
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
PHỤ LỤC
MÃ NGUỒN CÁC CHƯƠNG TRÌNH
%-----------------------------------------------------------------------------------------------
% Chương trình chính tính toán TTDT.
%-----------------------------------------------------------------------------------------------
global k0 so_chan_tu last lmax chuky alpha hsdh trovao
clc;
disp(‘ THIET KE ANTEN LOGA – CHU KY ’)
% doi tan so ra dang Hz
fmax=600*1e+6; %Hoặc nhập tay fmax=input('Tan so lon nhat (MHz):')
fmin=400*1e+6; %Hoặc nhập tay fmin = input('Tan so nho nhat (MHz):')
F0 = 500*1e+6; %input('Nhap tan so hoat dong giua fmax va fmin(MHz):');
last = 10; %input('Nhap so vong chay: ');
lamda = 3e+8/F0; %Bước sóng hoạt động
k0 = 2*pi/lamda;
so_chan_tu = 6;
lmax = 3e+8/(2*fmin);
No_bit = 10; %input(‘so bit ma hoa: ’);
row = 8; %số anten thử nghiệm
% Những giới hạn của các tham số anten
a0 = 0.85;
b0 = 0.95;
a1 = 10*pi/180; %alpha min = 10 o
b1 = 20*pi/180; %alpha max = 20o
gen = 2; %số gen trong 1 nhiễm sắc thể
cell = gen*No_bit; %Kích thước nhiễm sắc thể
muctieu = [];
Lưu Thị Hoa Linh 40
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Gene = round(rand(row,cell));
for ib = 1:last
disp(‘ vong lap thu ’)
ib
for ic = 1:row
cost = Gene(ic,:);
for jb = 1:No_bit
ckkc(1,jb) = cost(1,jb);
goc_d(1,jb) = cost(1,jb+No_bit);
end
%giải mã
chuky = decode(ckkc,No_bit);
alpha = decode(goc_d,No_bit);
%giới hạn các tham số
chuky = a0+chuky*(b0-a0)/(2^No_bit);
alpha = a1+alpha*(b1-a1)/(2^No_bit);
%tính hàm mục tiêu
muctieu(ic,1) = loga(chuky,alpha);
end
muctieu = -muctieu;
[muctieu,ind] = sort(muctieu);
Gene = Gene(ind(1:round(row/2)),:);
%lai ghép
cross = ceil((cell-1)*rand(round(row/2),1));
for ic = 1:2:(round(round(row/2)))
Gene(round(row/2)+ic,1:cross) = Gene(ic,1:cross);
Gene(round(row/2)+ic,cross+1:cell) = Gene(ic+1,cross+1:cell);
Lưu Thị Hoa Linh 41
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
Gene(round(row/2)+ic+1,1:cross) = Gene(ic+1,1:cross);
Gene(round(row/2)+ic+1,cross+1:cell) = Gene(ic,cross+1:cell);
end
if ib ~ = last
%biến dị
ix = ceil(row*rand);
iy = ceil(cell*rand);
Gene(ix,iy) = 1-Gene(ix,iy);
else
% vòng lặp cuối thì không biến dị, cho in các giá trị
cuoicung = Gene(1,:);
for jb = 1:No_bit
ck_cuoicung(1,jb) = cuoicung(1,jb);
goc_cuoicung(1,jb) = cuoicung(1,jb+No_bit);
end
chuky = decode(ck_cuoicung,No_bit);
alpha = decode(goc_cuoicung,No_bit);
chuky = a0+chuky*(b0-a0)/(2^No_bit);
alpha = a1+alpha*(b1-a1)/(2^No_bit);
muctieu_cuoicung = loga(chuky,alpha)
chuky
alpha = alpha*180/pi % alpha tính theo độ
trovao
hsdh_dBi = 10*log10(hsdh)/1.6 % tính ra dBi
end
end
Lưu Thị Hoa Linh 42
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
%-----------------------------------------------------------------------------------------------
% Phần chương trình tính hàm object (loga.m).
%-----------------------------------------------------------------------------------------------
function out = loga(cky,anpha)
global chky alfa k_cach lmax so_chan_tu x dong L hsdh trovao
chky = cky;
alfa = anpha;
L = [];
L(so_chan_tu) = lmax; % tính chiều dài các chấn tử.
for n = so_chan_tu-1:-1:1
L(n) = L(n+1)*chky;
end
x(1) = (L(1)/2)/tan(alfa); % tính tọa độ
for n = 2:so_chan_tu
x(n) = x(n-1)/chky;
end
for n = 1:so_chan_tu
for m = 1:so_chan_tu
if n ~ = m
k_cach(n,m) = abs(x(n)-x(m));
else k_cach(n,m) = 0;
end
end
end
bk = 4e-3;
Lưu Thị Hoa Linh 43
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
% Dẫn nạp các đoạn fiđe giữa các chấn tử liên tiếp
y = zeros(2,2,so_chan_tu+1);
for k = 2:so_chan_tu+1
for n = 1:2
for m = 1:2
y(n,m,k) = DanNap(n,m,k);
end
end
end
% Trở kháng riêng và trở kháng tương hỗ.
z = zeros(so_chan_tu,so_chan_tu);
for n = 1:so_chan_tu
for m = 1:so_chan_tu+1
if m = = so_chan_tu+1
z(n,m) = 0;
elseif n = = m
z(n,m) = tk_rieng(bk,L(n));
else
z(n,m) = trk_th(L(n),L(m),k_cach(n,m));
end
end
end
% tính các hệ số của ma trận
V0 = 1;
a = zeros(so_chan_tu);
for n = 1:so_chan_tu
for m = 1:so_chan_tu
Lưu Thị Hoa Linh 44
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
thstg = z(n,m)*(y(1,1,m+1)+y(2,2,m)+y(2,1,m+1)*z(n,m+1));
if n = = m
thstg = thstg+1;
end
if m = = 1
a(n,m) = thstg;
else
a(n,m) = thstg+y(1,2,m)*z(n,m-1);
end
end
b(n) = z(n,1)*y(2,1,2)*V0;
end
%Giải hệ điện áp trên các chấn tử.
v = a\b';
% Dong dien:
v(so_chan_tu+1) = 0;
for n = 1:so_chan_tu
v_thstg = y(1,1,n+1)*v(n)+y(1,2,n+1)*v(n+1)+y(2,2,n)*v(n);
if n = = 1
dong(n) = v_thstg+y(2,1,n);
else
dong(n) = v_thstg+y(2,1,n)*v(n-1);
end
end
% Trở kháng vào của anten
trovao = v(1)/dong(1);
t1 = dblquad('bf_hfh',0,pi,0,2*pi);
Lưu Thị Hoa Linh 45
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
hsdh = 4*pi/t1;
object = 30*hsdh;%-abs(75-real(trovao));
out = object;
%-----------------------------------------------------------------------------------------------
% Tổng bình phương hai hàm phương hướng (bf_hfh.m).
%-----------------------------------------------------------------------------------------------
function out = bf_hfh(goc1,goc2)
global so_chan_tu k0 L dong x
goc1 = 0:0.15:2*pi;
tong1 = 0;
goc2 = 0:0.15:2*pi;
tong2 = 0;
for m = 1:so_chan_tu
dau = cos(m*pi);
t_g1 = (1-cos(k0*L(m)/2))./sin(k0*L(m)/2);
e_mu1 = exp(i*k0*x(m)*cos(goc1));
tong1 = tong1+dau.*dong(m).*t_g1.*e_mu1;
t_g2 = (cos((k0*L(m)/2)*sin(goc2))-cos(k0*L(m)/2))/sin(k0*L(m)/2);
e_mu2 = exp(i*k0*x(m)*cos(goc2));
tong2 = tong2+dau.*dong(m).*t_g2.*e_mu2;
end
hfh_h = abs(tong1)/abs(tong1(1));
hfh_e = abs(tong2)/abs(tong2(1));
out = (sqr(hfh_h)+sqr(hfh_e)).*sin(goc1);
Lưu Thị Hoa Linh 46
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
%-----------------------------------------------------------------------------------------------
% Trở kháng tương hỗ (trk_th.m).
%-----------------------------------------------------------------------------------------------
function out = trk_th(l01,l02,kc)
global kcach l1 l2 k0
l1 = l01;
l2 = l02;
kcach = kc;
t1 = -30*quad('ham_tk1',-l2/2,l2/2);
t2 = -30*quad('ham_tk2',-l2/2,l2/2);
out = t1+j*t2;
%-----------------------------------------------------------------------------------------------
% hai file ham_tk1.m và ham_tk2.m:
%-----------------------------------------------------------------------------------------------
function out = ham_tk1(t)
global kcach k0 l1 l2
R = sqrt(kcach^2 + t.^2);
r1 = sqrt(kcach^2 + (t+l1/2).^2);
r2 = sqrt(kcach^2 + (t-l1/2).^2);
out = (2*sin(k0*R)*cos(k0*l1/2)/R-sin(k0*r1)/r1-sin(k0*r2)/r2)*
*sin(k0*(l2/2-abs(t)))/(sin(k0*l1/2)*sin(k0*l2/2));
%-----------------------------------------------------------------------------------------------
function out = ham_tk2(t)
global kcach k0 l1 l2
R = sqrt(kcach^2 + t.^2);
r1 = sqrt(kcach^2 + (t+l1/2).^2);
r2 = sqrt(kcach^2 + (t-l1/2).^2);
Lưu Thị Hoa Linh 47
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
out = (2*cos(k0*R)*cos(k0*l1/2)/R-cos(k0*r1)/r1-cos(k0*r2)/r2)*
*sin(k0*(l2/2-abs(t)))/(sin(k0*l1/2)*sin(k0*l2/2));
%-----------------------------------------------------------------------------------------------
% Trở kháng riêng (trk_rieng.m).
%-----------------------------------------------------------------------------------------------
function out = tk_rieng(bk,l)
format long e;
global M h k0 te_ta an_pha
M = 8; %Số vòng lặp
h = 1/M;
te_ta = k0*l/2;
an_pha = 2*bk/l;
% Gán các hệ số của ma trận R (n x m)
for i = 1:M
for j = 1:M
a(i,j) = complex(R_thuc(i,j-1),R_ao(i,j-1));
end
b(i) = complex(sin(i*te_ta*h),0);
end
% Giải hệ phức
x = a\b';
he_so = complex(0,20*h);
out = he_so/x(1);
%-----------------------------------------------------------------------------------------------
% Hàm hệ số ma trận
function R_thuc = R_thuc(m,n);
global M h k0 te_ta an_pha
Lưu Thị Hoa Linh 48
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
tam1 = 0;
tam2 = 0;
for i = 0:M
tam1 = tam1+s(i)*(1/Can(i,n)+1/Can(i,-n));
tam2 = tam2+2*s(i)/Can(i,0);
end;
r1 = cos(te_ta*Can(n,-m))/Can(n,-m);
r2 = cos(te_ta*Can(n,m))/Can(n,m);
r3 = 2*cos(te_ta*m*h)*cos(te_ta*Can(n,0))/Can(n,0);
r4 = s(n)*(r1+r2-r3);
r5 = (3/h)*log(((1-n*h+sqrt(sqr(an_pha)+sqr(1-n*h)))*
*(1+n*h+sqrt(sqr(an_pha)+sqr(1+n*h))))/sqr(an_pha));
r6 = cos(te_ta*an_pha)*(r5-tam1);
r7 = cos(te_ta*an_pha)*cos(te_ta*m*h)*
*((6/h)*log((1+sqrt(1+sqr(an_pha)))/an_pha)-tam2);
if m = = n
tam = r4+r6;
else tam = r4;
end
if n = = 0
R_thuc = tam-r7;
else R_thuc = tam;
end
% phần ảo
function R_ao = R_ao(m,n);
global te_ta an_pha h
r1 = sin(te_ta*Can(n,-m))/Can(n,-m);
Lưu Thị Hoa Linh 49
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
r2 = sin(te_ta*Can(n,m))/Can(n,m);
r3 = 2*cos(te_ta*m*h)*sin(te_ta*Can(n,0))/Can(n,0);
r = s(n)*(r1+r2-r3);
R_ao = -r;
%------------------------------------------------------------------------------
function can = Can(n,m);
global h an_pha
can = sqrt(sqr(an_pha)+h*h*sqr(n+m));
%------------------------------------------------------------------------------
% Hàm trọng số Simpson
function s = s(n);
if n = = 0|n = = 8 s = 1;
elseif rem(n,2) = = 0 s = 2;
else s = 4;
end
%-----------------------------------------------------------------------------------------------
% Dẫn nạp (DanNap.m).
%-----------------------------------------------------------------------------------------------
function out = DanNap(n1,m1,kc)
global k0 so_chan_tu x k_cach
format short e;
TKS = 120*pi;
k1 = kc;
if k1 = = so_chan_tu+1
out = 0;
elseif k1 = = 1
if n1 = = m1
Lưu Thị Hoa Linh 50
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
out = coth(i*k0*x(1))/TKS;
else
out = -csch(i*k0*x(1))/TKS;
end
elseif n1 = = m1
out = coth(i*k0*k_cach(k1,k1-1))/TKS;
else
out = -csch(i*k0*k_cach(k1,k1-1))/TKS;
end
%-----------------------------------------------------------------------------------------------
Lưu Thị Hoa Linh 51
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
MỤC LỤC
CHƯƠNG 1: ANTEN LOGA – CHU KỲ ........................................................5
1.1. Dải thông tần và dải tần công tác của anten. ...........................................5
1.1.1. Dải thông tần...................................................................................................5
1.1.2. Dải tần công tác. .............................................................................................5
1.2. Phương pháp mở rộng dải tần số của anten chấn tử. ..............................6
1.2.1. Giảm nhỏ trở kháng sóng của chấn tử. ......................................................6
1.2.2. Biến đổi từ từ thiết diện của chấn tử. .........................................................7
1.2.3. Hiệu chỉnh trở kháng vào của chấn tử. ......................................................9
1.3. Mở rộng dải tần công tác theo nguyên lý tương tự. Áp dụng cho anten
lôga - chu kỳ. .......................................................................................................9
1.3.1. Nguyên lý cấu tạo của anten lôga - chu kỳ. .............................................10
1.4. Đặc điểm kết cấu anten loga chu kỳ và phương pháp tính toán ...........13
CHƯƠNG 2: THUẬT TOÁN DI TRUYỀN...................................................19
2.1. Giới thiệu về thuật toán di truyền..........................................................19
2.2. Khái quát về thuật toán di truyền và các bước tiến hành. ....................20
2.2.1. Khái quát về thuật toán di truyền .............................................................20
2.2.2. Các bước tiến hành ......................................................................................20
2.2.3. Thuật toán di truyền cho tối ưu và thiết kế anten .................................23
2.2.4. Hàm mục tiêu và tỉ lệ phù hợp...................................................................25
2.2.5. Sự chọn lọc.....................................................................................................25
2.2.6. Mô tả anten như nhiễm sắc thể..................................................................26
2.2.7. Các toán tử di truyền...................................................................................28
CHƯƠNG 3: ỨNG DỤNG THUẬT TOÁN DI TRUYỀN...........................30
3.1. Các tham số tối ưu cho anten lôga - chu kỳ. ..........................................30
3.2. Các bước thực hiện tính toán. ................................................................30
3.2.1. Khởi tạo ngẫu nhiên quần thể gồm một số nhiễm sắc thể. ...................30
3.2.2. Thực hiện mã hóa các tham số của từng anten.......................................31
3.2.3. Giải mã các gen sang dạng thập phân và tính theo giới hạn................31
3.2.4. Thực hiện tính toán mỗi thiết kế anten........................................................32
3.2.5. Loại bỏ những nhiễm sắc thể không thỏa mãn theo tiêu chí đặt ra....34
3.2.6. Sử dụng các toán tử di truyền....................................................................34
3.2.7. Kiểm tra điều kiện dừng. ............................................................................34
3.3. Một số kết quả. ........................................................................................35
3.3.1. Trường hợp 1: chỉ tối ưu hệ số định hướng, không tối ưu trở kháng. .....35
3.3.1.1 Kết quả 1.....................................................................................35
3.3.1.2 Kết quả 2.....................................................................................35
3.3.1.3 Kết quả 3.....................................................................................36
Lưu Thị Hoa Linh 52
Đại học Công Nghệ - ĐHQG Hà Nội Khóa luận tốt nghiệp
3.3.2. Trường hợp 2: Tối ưu hệ số định hướng với hệ số là 30, đồng thời tối
ưu phần thực của trở kháng. ....................................................................................36
3.3.2.1 Kết quả 4.....................................................................................36
3.3.2.2 Kết quả 5.....................................................................................37
3.3.2.3 Kết quả 6.....................................................................................37
3.4. Nhận xét tổng quát các kết quả đạt được. .............................................37
3.5. Đánh giá kết quả và hướng phát triển. ..................................................38
Tài liệu tham khảo39
PHỤ LỤC 40
% Chương trình chính tính toán TTDT...........................................................40
% Phần chương trình tính hàm object (loga.m). .............................................43
% Tổng bình phương hai hàm phương hướng (bf_hfh.m)..............................46
% Trở kháng tương hỗ (trk_th.m). ..................................................................47
% Trở kháng riêng (trk_rieng.m). ...................................................................48
% Dẫn nạp (DanNap.m). ..................................................................................50
Lưu Thị Hoa Linh 53
Các file đính kèm theo tài liệu này:
- CPL 40.pdf