Thiết kế, mô phỏng anten tạo chùm và ứng dụng trong mạng MANET

NHIỆM VỤ ĐỒ ÁN “Thiết kế, mô anten tạo chùm ứng dụng mạng MANET” Thiết kế mô anten tạo chùm có đầy đủ thông số kỹ thuật phần mềm Matlab, trình bày bước cụ thể giai đoạn, kết thu anten tạo chùm ứng dụng anten tạo chùm node mạng MANET  Nội dung đồ án gồm: - Lý thuyết anten - Ứng dụng anten tạo chùm mạng MANET - Thiết kế mô anten tạo chùm Thái Nguyên, tháng năm 2016 Sinh viên Đỗ Quốc Trường LỜI CAM ĐOAN Tôi xin cam đoan nội dung đồ án ”Thiết kế, mô anten tạo chùm ứng dụng mạng MANET” tự tìm hiểu, nghiên cứu hướng dẫn thầy giáo KS.Trần Tuấn Việt Mọi trích dẫn tài liệu mà tham khảo ghi rõ nguồn gốc Nếu sai xin chịu hình thức kỉ luật trường Đại học Công nghệ thông tin Truyền thông Thái Nguyên, tháng năm 2016 Sinh viên Đỗ Quốc Trường LỜI CẢM ƠN Lời em xin gửi lời cảm ơn đến toàn thể thầy, cô giáo trường Đại học Công nghệ thông tin Truyền thông Thái Nguyên, người tận tình dạy dỗ, bảo em suốt năm năm học vừa qua nhà trường Tiếp theo em xin gửi lời cảm ơn sâu sắc đến KS Trần Tuấn Việt, người trực tiếp hướng dẫn em suốt trình học tập, nghiên cứu đồ án Bộ môn Công nghệ ô tô Hệ thống cảm biến – Trường Đại học Công nghệ thông tin Truyền thông Thái Nguyên Người truyền cho em cách tư có hệ thống, phương pháp nghiên cứu, tiếp cận thực tế - điều quý báu với em trường làm việc thực tế Em xin gửi lời cảm ơn tới toàn thể quý thầy cô Bộ môn Công nghệ ô tô Hệ thống cảm biến người dẫn dắt định hướng nghiên cứu cho em suốt năm qua Em xin tỏ lòng biết ơn chân thành tới cha mẹ, gia đình em người sinh thành, nuôi nấng, tin tưởng động viên em Xin gửi lời cảm ơn tới tất bạn bè đặc biệt tập thể lớp công nghệ ô tô K10, người cổ vũ, động viên, chia sẻ với em suốt năm qua Sinh viên thực Đỗ Quốc Trường MỤC LỤC NHIỆM VỤ ĐỒ ÁN LỜI CAM ĐOAN LỜI CẢM ƠN MỤC LỤC DANH MỤC HÌNH ẢNHTHUẬT NGỮ VIẾT TẮT LỜI NÓI ĐẦU 10 CHƯƠNG LÝ THUYẾT ANTEN 11 1.1 Tổng quan anten .11 1.1.1 Lịch sử phát triển anten 11 1.1.2 Hệ phương trình Maxwell .14 1.1.3 Quá trình vật lý xạ sóng điệntừ 15 1.1.4 Các thông số anten .16 1.2 Anten tạo chùm 23 1.2.1 Định nghĩa 23 1.2.2 Nguyên lý hoạt động .24 1.2.3 Thuật toán Hướng đến (DOA) 24 1.2.4 Anten tạo chùm thích ứng .31 1.2.5 Hiệu ứng ảnh hưởng tương hỗ .35 CHƯƠNG ỨNG DỤNG CỦA ANTEN TẠO CHÙM TRONG MẠNG MANET 37 2.1 Mạng tùy biến không dây .37 2.1.1 Định nghĩa 37 2.1.2 Tổng quan mạng tùy biến di động (MANET) 37 2.2 Vấn đề cần thiết thiết kế anten tạo chùm 39 2.3 MANET sử dụng hệ thống anten tạo chùm 42 CHƯƠNG THIẾT KẾ VÀ MÔ PHỎNG ANTEN TẠO CHÙM 48 3.1 Qúa trình thiết kế phương pháp thiết kế 48 3.1.1 Thiết kế thành phần đơn – miếng vi dải .49 3.1.2 Miếng hình chữ nhật 49 3.1.3 Thiết kế mảng 51 3.1.4 Thiết kế mảng tuyến tính .51 3.1.5 Thiết kế mạng lưới phân tử phẳng .53 3.1.6 Mạng lưới phân tử phẳng 4x4 so với Mạng lưới phân tử phẳng 8x8 54 3.1.7 Tạo chùm tia thích ứng 56 3.2 Mô với thông số cụ thể - kết cho mô anten tạo chùm .57 KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN 65 TÀI LIỆU THAM KHẢO .66 PHỤ LỤC 67 DANH MỤC HÌNH ẢNH Hình 1.1 Hệ thống anten thu phát [1] .13 Hình 1.2 Đồ thị phương hướng toạ độ cực [1] .20 Hình 1.3 Đồ thị phương hướng toạ độ góc [2] 20 Hình 1.4 Dàn phẳng có biểu diễn đồ thị độ trễ thời gian 25 Hình 1.5 Tín hiệu đến mảng phần tử [1] 30 Hình 1.6 Mảng phần tử tiếp nhận tín hiệu mong muốn θ0 = 00 SNOI θ0 = 300 [1] 32 Hình 1.7 Tính toán độ trễ pha sớm pha theo tín hiệu n(t) [1] 33 Hình 1.8 So sánh hệ số mảng có hiệu ứng ảnh hưởng tương hỗ.[1] 34 Hình 2.1 Ví dụ MANET đa chặng [1] .38 Hình 2.3 Búp sóng hướng thu thiết bị [2] 41 Hình 2.4 So sánh hiệu suất mạng sử dụng anten đẳng hướng mạng có anten thông minh [1] 43 Hình 2.5 Sơ đồ bố trí 55 node MANET sử dụng để thiết kế hệ thống anten thông minh thông lượng cao [1] 43 Hình 2.6 Ví dụ node liên hệ với nhau.[1] .44 Hình 2.7 Biểu đồ định thời gian Giao thức MAC [1] 45 Hình 3.1 Thiết kế anten miếng hình chữ nhật sử dụng mô hình khoang [1] 50 Hình 3.2 S (độ suy giảm thích ứng) cho anten miếng hình chữ nhật 20 GHz [1] 50 Hình 3.3 Mẫu biên độ Anten miếng hình chữ nhật sử dụng mô hình khoang [1] 52 Hình 3.4 Thiết kế mảng tuyến tính thành phần với miếng hình chữ nhật.[1] 52 Hình 3.5 Thiết kế mạng lưới phân tử phẳng với miếng hình chữ nhật.[1] .54 Hình 3.6 Thông lượng mạng trung bình so với tải mạng kích thước mảng khác kích thích [1] .55 Hình 3.7 Ảnh hưởng góc tín hiệu hướng đến độ rộng búp sóng [1] 56 Hình 3.8 Mặt phẳng E = 0° so sánh mẫu phóng xạ tạo chùm tia không thích ứng chùm tia thích ứng [1] 57 Hình 3.9 So sánh thông lượng mạng tạo chùm tia thích ứng không thích ứng [1] 57 Hình 3.10 Cường độ xạ tín hiệu thu anten tạo chùm mảng phẩn tử khác 58 Hình 3.11 Tín hiệu thu biểu diễn theo vị trí có độ lợi mặt tín hiệu 59 Hình 3.12 Những null búp sóng thu anten tạo chùm .60 Hình 3.13 Tín hiệu có cường độ xạ cao vào anten tạo chùm có búp sóng nhỏ chiều rộng hẹp khả truyền đến tín hiệu quan tâm .61 Hình 3.14 Hình dạng anten tạo chùm mô 3D Matlab 62 Hình 3.15 Tỷ lệ lỗi bit BER Beamforming .63 THUẬT NGỮ VIẾT TẮT CDMA Phân chia mã đa truy cập DSP Bộ xử lý tín hiệu số DOA Thuật toán hướng đến EIRP Công suất xạ đẳng hướng tương đương ESPRIT Kỹ thuật quay bất biến FSD Kỹ thuật phân tích không gian FDMA Phân chia tần số đa truy cập LMS Định hướng búp sóng tối ưu MSE Sai số toàn phương trung bình MUSIC Thuật toán phân tách đa tín hiệu ML Kỹ thuật hợp lý cực đại MANET Mạng tùy biến không dây SIR Tiêu chuẩn cực đại tỷ số tín hiệu nhiễu SIR SDMA Phân chia không gian đa truy cập SOI Tín hiệu quan tâm SNOI Tín hiệu không quan tâm TDMA Phân chia thời gian đa truy cập LỜI NÓI ĐẦU Thông tin xuất từ lâu từ người biết dùng lửa, tiếng động âm thanh, kí hiệu tượng hình để liên lạc trao đổi Trải qua trình phát triển, nhu cầu thông tin liên lạc người đòi hỏi phù hợp với thực tế nhanh, xác đến năm 1837 Samuel Morse phát minh ám hiệu truyền tin dựa cách thức đóng mở dòng điện gây nên tiếng Năm 1894 Maxwell đưa lý thuyết dạng vật chất lan truyền xa chân không Năm 1902 kiện mở kỷ nguyên vể thông tin liên lạc, tạo tiền đề cho nhiều ứng dụng viễn thông sau Đóng góp vào thông tin liên lạc không kể tới vai trò anten thiết bị dùng để truyền đạt thu nhận tín hiệu Anten xuất từ lâu nói có niên đại với thông tin liên lạc Để đáp ứng nhu cầu thông tin liên lạc ngày phát triển công nghệ anten phải phát triển theo điển hình ứng dụng truyền xa thông tin vệ tinh sử dụng anten parabol để thu phát với loại anten có độ lợi cao độ định hướng lớn Ngoài không nói đến xu hướng thời đại nhỏ gọn, đa ứng dụng Đây điều tất yếu anten phải nhỏ gọn để đáp ứng yêu cầu Để cải thiện nhược điểm có nhiều cách thức hệ thống anten mảng (array antenna) để tăng độ lợi cho anten, điều khiển búp sóng Ngoài có cải thiện vật chất cho anten cải thiện anten tạo chùm có đặc tính khác biệt Đồ án tốt nghiệp em : “Thiết kế, mô anten tạo chùm ứng dụng mạng MANET” Đồ án em chia thành chương : Chương Lý thuyết anten Chương Ứng dụng anten tạo chùm mạng MANET Chương Thiết kế mô anten tạo chùm Do thời gian kiến thức có hạn nên đồ án không tránh khỏi thiếu sót Rất mong đóng góp ý kiến giúp đỡ thầy cô giáo Em xin chân thành cảm ơn ! 10 %Plot steered response figure(3);clf surf(kx, ky, spectrumNormalized, 'edgecolor', 'none', 'FaceAlpha', 0.8) view(0, 90) axis square %Do some magic to make the figure look nice set(gcf,'color','k') cmap = colormap; cmap(1,:) = [1 1]*0.2; colormap(cmap); set(gca,'color',[0 0],'xcolor',[1 1],'ycolor',[1 1],'zcolor',[1 1]) set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[]) %set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Mi norGridColor',[1 1],'MinorGridLineStyle','-') xlabel('k_x = sin(\theta)cos(\phi)') ylabel('k_y = sin(\theta)sin(\phi)') %% 2D-array case, different source strengths, spectrum in linear scale %Relative amplitude difference between sources in decibel, %could also have written amplitudes = [3 0] amplitudes = [0 -2 -3]; % Create input signal inputSignal = createSignal(x_pos, y_pos, f, c, fs, thetaArrivalAngles, phiArrivalAngles, amplitudes); % Calculate delay-and-sum steered response [S, kx, ky] = steeredResponseDelayAndSum(x_pos, y_pos, w, inputSignal, f, c, 74 thetaScanningAngles, phiScanningAngles); %Normalise spectrum spectrumNormalized = abs(S)/max(max(abs(S))); % Plot the steered response figure(4);clf surf(kx, ky, spectrumNormalized, 'edgecolor', 'none', 'FaceAlpha', 0.8) view(0, 90) axis square %Do some magic to make the figure look nice set(gcf,'color','k') cmap = colormap; cmap(1,:) = [1 1]*0.2; colormap(cmap); set(gca,'color',[0 0],'xcolor',[1 1],'ycolor',[1 1],'zcolor',[1 1]) set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[]) %set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Mi norGridColor',[1 1],'MinorGridLineStyle','-') xlabel('k_x = sin(\theta)cos(\phi)') ylabel('k_y = sin(\theta)sin(\phi)') %% 2D-array case, different source strengths, spectrum in logarithmic scale %Convert the delay-and-sum steered response to decibel spectrumLog = 10*log10(spectrumNormalized); % Plot the steered response figure(5);clf surf(kx, ky, spectrumLog, 'edgecolor', 'none', 'FaceAlpha', 0.8) view(0, 90) 75 axis square %Do some magic to make the figure look nice set(gcf,'color','k') cmap = colormap; cmap(1,:) = [1 1]*0.2; colormap(cmap); set(gca,'color',[0 0],'xcolor',[1 1],'ycolor',[1 1],'zcolor',[1 1]) set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[]) %set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Mi norGridColor',[1 1],'MinorGridLineStyle','-') xlabel('k_x = sin(\theta)cos(\phi)') ylabel('k_y = sin(\theta)sin(\phi)') %% 2D-array case, different source strengths, spectrum in logarithmic scale with dynamic range %Dynamic range in decibels dynamicRange = 5; % Plot the steered response figure(6);clf surf(kx, ky, spectrumLog, 'edgecolor', 'none', 'FaceAlpha', 0.8) view(0, 90) axis square %Do some magic to make the figure look nice set(gcf,'color','k') cmap = colormap; cmap(1,:) = [1 1]*0.2; colormap(cmap); set(gca,'color',[0 0],'xcolor',[1 1],'ycolor',[1 1],'zcolor',[1 1]) set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[]) 76 %set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Mi norGridColor',[1 1],'MinorGridLineStyle','-') xlabel('k_x = sin(\theta)cos(\phi)') ylabel('k_y = sin(\theta)sin(\phi)') %Set the dynamic range caxis([-dynamicRange 0]) %% 2D-array case, different source strengths, use plotSteeredResponse function %Max range in decibels maxDynamicRange = 30; %Plot the steered response, no need to normalise the %spectrum or converting to decibel before input plotSteeredResponse(S, kx, ky, maxDynamicRange, 'lin', 'black', '3D') function [] = plotSteeredResponse(S, kx, ky, maxDynamicRange, scaleView, displayTheme, displayStyle) %plotSteeredResponse - plot the steered response with dynamic range slider % %Plots the steered response with a slider bar in either linear or %logarithmic scale and in black or white theme % %plotSteeredResponse(S, kx, ky, maxDynamicRange, scaleView, displayTheme, displayStyle) % %IN %S - NxM matrix of delay-and-sum steered response power %kx - 1xN vector of theta scanning angles in polar coordinates 77 %ky - 1xM vector of phi scanning angles in polar coordinates %maxDynamicRange - max dynamic range in decibels in the image %scaleView - slider scale, use 'lin' for linear and 'log' for log %displayTheme - color theme of image, use 'white' or 'black' %displayStyle - view in 2D or 3D % %OUT %[] - The figure plot % %Created by Jørgen Grythe, Norsonic AS %Last updated 2015-09-30 if ~exist('displayStyle','var') displayStyle = '2D'; end if ~exist('displayTheme','var') displayTheme = 'black'; end if ~exist('scaleView','var') scaleView = 'log'; end if ~exist('maxDynamicRange','var') maxDynamicRange = 15; end %The input is a power signal, normalize and convert to decibel S = abs(S)/max(max(abs(S))); S = 10*log10(S); 78 %Default display value defaultDisplayValue = 10; if defaultDisplayValue > maxDynamicRange defaultDisplayValue = maxDynamicRange/2; end %Plot the steered response figure; surf(kx, ky, S, 'edgecolor', 'none', 'FaceAlpha', 0.8) if isequal(displayStyle, '2D') view(0, 90) elseif isequal(displayStyle, '3D') view(40, 40) else error('Use 2D or 3D for displayStyle') end axis square if isequal(displayTheme,'black') set(gcf,'Color','k') cmap = colormap; cmap(1,:) = [1 1]*0.2; colormap(cmap); set(gca,'Color',[0 0],'xcolor',[1 1],'ycolor',[1 1],'zcolor',[1 1]) % set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Min orGridColor',[1 1],'MinorGridLineStyle','-') textColor = 'w'; elseif isequal(displayTheme,'white') set(gcf,'Color','w') cmap = colormap; cmap(1,:) = [1 1]*0.95; colormap(cmap); 79 set(gca,'Color',[1 1],'xcolor',[0 0],'ycolor',[0 0],'zcolor',[0 0]) % set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Min orGridColor',[0 0],'MinorGridLineStyle','-') textColor = 'k'; else error('Use black or white for displayStyle') end %Add labels and title set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[]) xlabel('k_x = sin(\theta)cos(\phi)', 'Color', textColor) ylabel('k_y = sin(\theta)sin(\phi)', 'Color', textColor) title(['Dynamic range: ' sprintf('%0.2f', defaultDisplayValue) ' dB'], 'FontWeight', 'normal', 'Color', textColor); %Just show color from defaultDisplayValue up to dBmax caxis([-defaultDisplayValue 0]) %Add slider that goes from 0.01 dB up to dBmax and changes color/title range = [0.01 maxDynamicRange]; if isequal(scaleView,'log') h = uicontrol('style', 'slider', 'Units', 'normalized', 'position', [0.9 0.1 0.03 0.6], 'value', log10(defaultDisplayValue), 'min', log10(range(1)), 'max', log10(range(2))); addlistener(h,'ContinuousValueChange',@(hObject,eventdata) caxis([-10^get(hObject, 'Value') 0])); 80 addlistener(h,'ContinuousValueChange',@(hObject,eventdata) title(['Dynamic range: ' sprintf('%0.2f', 10^get(hObject, 'Value')) ' dB'],'fontweight','normal')); elseif isequal(scaleView,'lin') h = uicontrol('style', 'slider', 'Units', 'normalized', 'position', [0.9 0.1 0.03 0.6], 'value', defaultDisplayValue, 'min', range(1), 'max', range(2)); addlistener(h,'ContinuousValueChange',@(hObject,eventdata) caxis([-get(hObject, 'Value') 0])); addlistener(h,'ContinuousValueChange',@(hObject,eventdata) title(['Dynamic range: ' sprintf('%0.2f', get(hObject, 'Value')) ' dB'],'fontweight','normal')); else error('Use log or lin for scale value') end %% Help figure to visualise vectors in spherical coordinates theta = 35; phi = 60; %Vector coordinates lx = sin(theta*pi/180)*cos(phi*pi/180); ly = sin(theta*pi/180)*sin(phi*pi/180); lz = cos(theta*pi/180); %Half sphere [sx, sy, sz] = sphere(100); sz(find(sz < 0)) = 0; 81 %Set color for lines ++ cmap = [1 1]; h = figure(1); set(gcf,'color','w') hold on %Plot sphere surf(sx,sy,sz,'edgecolor','none','Facecolor',cmap,'FaceAlpha',0.2) plot(cos(0:pi/50:2*pi),sin(0:pi/50:2*pi),'Color',cmap) %Plot the vector with help lines line([0 lx],[0 ly],[0 lz],'LineWidth',1.5,'Color',[1.0000 0.9059 0.0941]) line([lx lx],[ly ly],[0 lz],'Color',cmap,'LineStyle',':') line([0 lx],[0 ly],[0 0],'Color',cmap,'LineStyle',':') line([0 lx],[0 ly],[lz lz],'Color',cmap,'LineStyle',':') %Plot coordinate system lines %x-lines line([0 1],[0 0],[0 0],'Color',cmap) line([-1 1],[1 1],[0 0],'Color',cmap) %y-lines line([0 0],[0 1],[0 0],'Color',cmap) line([-1 -1],[-1 1],[0 0],'Color',cmap) %z-lines line([0 0],[0 0],[0 1],'Color',cmap) line([-1 -1],[1 1],[0 1],'Color',cmap) %text text(0.8,-0.1,'x','color',cmap) text(-0.2,0.9,'y','color',cmap) text(0,-0.1,0.9,'z','color',cmap) text(lx/10,ly/10,0.2,'\theta','color',cmap) text(0.1,0.2,0,'\phi','color',cmap) 82 %Set colors and view axis ++ view(30,30) set(gca,'color',[0 0],'xcolor',cmap,'ycolor',cmap,'zcolor',cmap) set(gca,'XTickLabel',[],'YTickLabel',[],'ZTickLabel',[]) set(gca,'XMinorGrid','on','YMinorGrid','on','ZMinorGrid','on','Min orGridColor',cmap,'MinorGridLineStyle','-') axis equal title(['\theta = ' num2str(theta) ', \phi = ' num2str(phi)],'FontWeight','Normal') function [S, kx, ky] = steeredResponseDelayAndSum(xPos, yPos, w, inputSignal, f, c, thetaScanningAngles, phiScanningAngles) %steeredResponseDelayAndSum - calculate delay and sum in frequency domain % %Calculates the steered response from the delay-and-sum algorithm in the %frequency domain based on sensor positions, input signal and scanning angles % %[S, kx, ky] = steeredResponseDelayAndSum(xPos, yPos, w, inputSignal, f, c, thetaScanningAngles, phiScanningAngles) % %IN %xPos - 1xP vector of x-positions [m] %yPos - 1xP vector of y-positions [m] %w - 1xP vector of element weights %inputSignal - PxL vector of inputsignals consisting of L samples %f - Wave frequency [Hz] %c - Speed of sound [m/s] %thetaScanningAngles - 1xN vector of theta scanning angles [degrees] %phiScanningAngles - 1xM vector of phi scanning angles [degrees] % 83 %OUT %S - NxM matrix of delay-and-sum steered response power %kx - 1xN vector of theta scanning angles in polar coordinates %ky - 1xM vector of phi scanning angles in polar coordinates % %Last updated 2015-10-06 if ~exist('thetaScanningAngles', 'var') thetaScanningAngles = -90:90; end if ~exist('phiScanningAngles', 'var') phiScanningAngles = 0:180; end nSensors = numel(xPos); nSamples = numel(inputSignal); nThetaAngles = numel(thetaScanningAngles); nPhiAngles = numel(phiScanningAngles); %Calculate steering vector for all scanning angles [e, kx, ky] = steeringVector(xPos, yPos, f, c, thetaScanningAngles, phiScanningAngles); %Multiply input signal by weighting vector inputSignal = diag(w)*inputSignal; %Calculate correlation matrix R = inputSignal*inputSignal'; R = R/nSamples; 84 %Calculate power as a function of steering vector/scanning angle (delay-and-sum) S = zeros(nThetaAngles,nPhiAngles); for angleTheta = 1:nThetaAngles for anglePhi = 1:nPhiAngles ee = reshape(e(angleTheta, anglePhi,:), nSensors, 1); S(angleTheta,anglePhi) = ee'*R*ee; end end function [ee, kx, ky] = steeringVector(xPos, yPos, f, c, theta, phi) %steeringVector - calculate steering vector of 2D array % %Calculates the steering vector for different scanning angles % %[ee, kx, ky] = steeringVector(xPos, yPos, f, c, theta, phi) % %IN %xPos - 1xP vector of x-positions [m] %yPos - 1xP vector of y-positions [m] %f - Wave frequency [Hz] %c - Speed of sound [m/s] %theta - 1xM vector of theta angles [degrees] %phi - 1xN vecor of phi angles [degrees] % %OUT %ee - MxNxP matrix of steering vectors %kx - theta scanning angles in polar coordinates %ky - phi scanning angles in polar coordinates % %Last updated 2015-09-30 if ~isvector(xPos) 85 error('X-positions of array elements must be a 1xP vector where P is number of elements') end if ~isvector(yPos) error('Y-positions of array elements must be a 1xP vector where P is number of elements') end %theta is the elevation and is the normal incidence angle from -90 to 90 if ~exist('theta', 'var') theta = -pi/2:pi/180:pi/2; else theta = theta*pi/180; end %phi is the azimuth, and is the angle in the XY-plane from to 360 if ~exist('phi', 'var') phi = 0:pi/180:2*pi; else phi = phi*pi/180; end %Wavenumber k = 2*pi*f/c; %Number of elements/sensors in the array P = size(xPos,2); %Changing wave vector to spherical coordinates kx = sin(theta)'*cos(phi); ky = sin(theta)'*sin(phi); 86 %Calculate steering vector/matrix kxx = bsxfun(@times,kx,reshape(xPos,1,1,P)); kyy = bsxfun(@times,ky,reshape(yPos,1,1,P)); ee = exp(1j*k*(kxx+kyy)); 87 NHẬN XÉT CỦA GIÁO VIÊN HƯỚNG DẪN ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… ………………………………………………………………………………… Thái nguyên, ngày … tháng năm 2016 Giáo viên hướng dẫn 88 ... động (MANET) 37 2.2 Vấn đề cần thiết thiết kế anten tạo chùm 39 2.3 MANET sử dụng hệ thống anten tạo chùm 42 CHƯƠNG THIẾT KẾ VÀ MÔ PHỎNG ANTEN TẠO CHÙM 48 3.1 Qúa trình thiết. .. cho anten cải thiện anten tạo chùm có đặc tính khác biệt Đồ án tốt nghiệp em : Thiết kế, mô anten tạo chùm ứng dụng mạng MANET Đồ án em chia thành chương : Chương Lý thuyết anten Chương Ứng dụng. .. Điều cần thiết cho anten tạo chùm thích ứng dựa DOA Đối với anten tạo chùm thích ứng dựa tham chiếu (hoặc quy chuẩn), Least Mean Square (LMS) sử dụng chương này, anten tạo chùm thích ứng không