dipoleantenna_frequence.m

一个用matlab编写

clear all;close all;
a=[1e-7,1e-5,1e-3];%三个天线半径
L=1;%天线长
lamda0=2;%波长为此数时,为半波天线
n=21;%段数,必须为奇数
e=8.85e-12;%介电常数
u=4e-7*pi;%磁导率
eta=120*pi;%波阻抗
delta=L/(n+1);
for m=1:3
    for t=1:300
        lamda=200*L/t;
        k=2*pi/lamda;%波数
        omega=2*pi*3e8/lamda;%角速度
        for mm=1:n
            for nn=1:n
                R(mm,nn)=abs((mm-nn)*delta);%mm,nn之间的距离
                R1(mm,nn)=abs((mm-nn+1)*delta);%mm+,nn-之间的距离
                R2(mm,nn)=abs((mm-nn-1)*delta);%mm-,nn+之间的距离
                if mm==nn %计算psi(mm,nn)矩阵
                    psi(mm,nn)=log(delta/a(m))/(2*pi*delta)-j*k/(4*pi);
                else
                    psi(mm,nn)=exp(-j*k*R(mm,nn))/(4*pi*R(mm,nn));
                end
                if mm==nn-1 %计算psi(mm+,nn-)矩阵
                    psi1(mm,nn)=log(delta/a(m))/(2*pi*delta)-j*k/(4*pi);
                else
                    psi1(mm,nn)=exp(-j*k*R1(mm,nn))/(4*pi*R1(mm,nn));
                end
                if mm==nn+1%计算psi(mm-,nn+)矩阵
                    psi2(mm,nn)=log(delta/a(m))/(2*pi*delta)-j*k/(4*pi);
                else
                    psi2(mm,nn)=exp(-j*k*R2(mm,nn))/(4*pi*R2(mm,nn));
                end
                Z(mm,nn)=j*omega*u*delta*delta*psi(mm,nn)+(2*psi(mm,nn)-psi1(mm,nn)-psi2(mm,nn))/(j*omega*e); %计算阻抗矩阵
            end
        end
        voltage=zeros(n,1); %对电压矩阵进行赋值
        voltage((n+1)/2,1)=1;
        I_n=inv(Z)*voltage;    %计算电流矩阵
        Zin=1/I_n((n+1)/2,1) ;   %输入阻抗
        zin(t)=Zin;
    end
    mm=0.01:0.01:3;
    figure(1);
    if(m==1)
        plot(mm,real(zin),'r');
    end
    if(m==2)
        plot(mm,real(zin),'b');
    end
    if(m==3)
        plot(mm,real(zin),'g');
    end
    grid on;hold on;
    title('输入电阻随频率变化关系图')
    xlabel('f/fo');ylabel('输入电阻(\Omega)');
    legend('1e-7','1e-5','1e-3');
    figure(2);
    if(m==1)
        plot(mm,imag(zin),'r');
    end
    if(m==2)
        plot(mm,imag(zin),'b');
    end
    if(m==3)
        plot(mm,imag(zin),'g');
    end
    grid on;hold on;
    title('输入电抗随频率变化关系图')
    xlabel('f/fo');ylabel('输入电抗(\jOmega)');
    legend('1e-7','1e-5','1e-3');
end