利用matlab繪制電基本振子E面方向圖和空間立體方向圖 (2).doc

微波技術(shù)與天線實(shí)驗(yàn)報(bào)告 學(xué)院： 信電學(xué)院班級(jí)： 通信2班 姓名: 學(xué)號(hào)： 2014年4月29日利用matlab繪制電基本振子E面方向圖和空間立體方向圖一、 實(shí)驗(yàn)?zāi)康?、 復(fù)習(xí)Matlab的使用。2、 利用Matlab繪制電基本振子E面方向圖。3、 利用Matlab繪制電基本振子空間立體方向圖。二、 實(shí)驗(yàn)原理電基本振子（Electric Short Dipole）又稱電流元，它是指一段理想的高頻電流直導(dǎo)線，其長度l遠(yuǎn)小于波長，其半徑a遠(yuǎn)小于l，同時(shí)振子沿線的電流I處處等幅同相。用這樣的電流元可以構(gòu) 成實(shí)際的更復(fù)雜的天線，因而電基本振子的輻射特性是研究更復(fù)雜天線輻射特性的基礎(chǔ)。三、 實(shí)驗(yàn)設(shè)備計(jì)算機(jī)、matlab7.0仿真軟件四、 實(shí)驗(yàn)內(nèi)容及步驟 1、根據(jù)電基本振子方向函數(shù)利用Matlab編程，并畫出其E方向圖。源程序如下：sita=meshgrid(eps:pi/180:pi);fai=meshgrid(eps:2*pi/180:2*pi);f=abs(sin(sita);fmax=max(max(f);a=linspace(0,2*pi);f=sin(a);subplot(1,1,1),polar(a,abs(f); title(電基本振子E平面);運(yùn)行結(jié)果如圖1（電基本振子E方向圖）。2、利用Matlab繪制電基本振子空間立體圖。代碼如下：sita=meshgrid(eps:pi/180:pi);fai=meshgrid(eps:2*pi/180:2*pi);f=abs(sin(sita);fmax=max(max(f);x,y,z=sph2cart(fai,pi/2-sita,f/fmax);subplot(1,1,1),mesh(x,y,z);axis(-1 1 -1 1 -1 1);title(電基本振子空間立體方向圖);運(yùn)行結(jié)果如圖2（電基本振子空間立體方向圖）。圖1 電基本振子E方向圖圖2 電基本振子空間立體方向圖五、 實(shí)驗(yàn)小結(jié)通過MATLAB編程，熟悉了MATLAB的基本使用方法，對(duì)電基本振子的E方向圖和空間立體圖有了深入的了解，在實(shí)驗(yàn)中遇到的編程錯(cuò)誤讓我認(rèn)識(shí)到了自己的不足，在今后的學(xué)習(xí)中一定要養(yǎng)成細(xì)心的好習(xí)慣，為以后的學(xué)習(xí)打下了堅(jiān)實(shí)的基礎(chǔ)。利用matlab繪制對(duì)稱振子天線E方向圖及空間立體方向圖一、 實(shí)驗(yàn)?zāi)康?、 復(fù)習(xí)Matlab的使用。2、 利用Matlab繪制對(duì)稱振子天線E面方向圖。3、 利用Matlab繪制對(duì)稱振子空間立體方向圖。二、 實(shí)驗(yàn)原理1、 對(duì)稱振子是中間饋電，其兩臂由兩段等長導(dǎo)線構(gòu)成的振子天線。2、 對(duì)稱振子以波腹電流歸算的方向函數(shù)為 ： 上式實(shí)際上也就是對(duì)稱振子E面的方向函數(shù)。三、 實(shí)驗(yàn)設(shè)備計(jì)算機(jī)、matlab7.0仿真軟件四、 實(shí)驗(yàn)內(nèi)容及步驟1、 根據(jù)對(duì)稱振子方向函數(shù)利用Matlab編程，并畫出其E方向圖。源程序如下：clear all;close all;sita=eps:pi/90:2*pi;fi= eps:pi/90:2*pi;f1=cos(pi/2*cos(sita)./sin(sita);f2=cos(pi/4*(3*sin(sita)+1);f=abs(f1).*abs(f2);f=f/max(f);subplot(121);polar(sita,f);title(對(duì)稱振子E平面);運(yùn)行結(jié)果如圖1所示（對(duì)稱振子E平面圖）。2、 利用Matlab繪制電基本振子空間立體圖。代碼如下：clear;m=moviein(20);for i=1:20;sita=meshgrid(eps:pi/180:pi);fai=meshgrid(eps:2*pi/180:2*pi);l=i*0.1;r=abs(cos(2.*pi.*l.*cos(sita)-cos(2*pi*l)./(sin(sita)+eps);rmax=max(max(r);x,y,z=sph2cart(fai,pi/2-sita,r/rmax);mesh(x,y,z);axis(-1 1 -1 1 -1 1);title(對(duì)稱振子空間立體方向圖);m(:,i)=getframe;end;運(yùn)行結(jié)果如圖2所示（對(duì)稱振子空間立體方向圖）。圖1 對(duì)稱振子E平面圖圖2 對(duì)稱振子空間立體方向圖五、 實(shí)驗(yàn)小結(jié)通過MATLAB編程