Find the exact eigen value graph of dipole antenna using characteristic mode theory

Discussion in 'Wireless & RF Design' started by shaswatee23, Mar 12, 2013.

  1. shaswatee23

    Thread Starter New Member

    Mar 12, 2013
    1
    0
    I am doing M.Tech project and my topic is Antenna design using characteristic mode theory. For this project I have taken the impedance matrix of the antenna which is depend on only the shape and size of the antenna not the excitation level.

    So for finding the impedance matrix the halfwave dipole antenna has chossen having length 0.5lamda, 0.001lamda diameter.

    The eigen value and the eigen vector can be find out by using the formula i.e. [J,D]=eig(X,R). Here the J is eigen current and D is the eigen value of the impedace matrix where X and R are imaginary and real parts of the impedance matrix.
    The characteristic current curve i.e. eigen vector curve has come properly but in the curve of eigen value has some problem.

    The value of eigen vector are coming correctly as previously found. but the eigen values are coming very small value and the resonat frequency of corresponding eigen value is not matched.For first eigenvalue its coming high frequency and for forth eigen value its coming lower frequency.
    For the above problem the code I used posted here.

    Code ( (Unknown Language)):
    1.  
    2. pi=3.14159265;
    3. rad=pi/180;
    4. beta=2.0*pi;
    5. eta=120*pi;
    6. %mu=(4*pi)*(10^-7);
    7. %******************************************************
    8. %***  INPUT DATA
    9. %******************************************************
    10. nm=(input('NUMBER OF SUBDIVISIONS <ODD NUMBER> ='));
    11. tl1=(input('TOTAL DIPOLE LENGTH <WAVELENGTHS> ='));
    12. ra1=(input('RADIUS OF DIPOLE <WAVELENGTHS> ='));
    13. %*********************************************************
    14. %***  EXITATION OPTION # : 1 ---> MAGNETIC-FRILL
    15. %                          2 ---> DELTA-GAP
    16. %**********************************************************
    17. disp('EXCITATION :');
    18. disp('        OPTION (1):MAGNETIC-FRILL');
    19. disp('        OPTION (2): DELTA-GAP');
    20. iex=input('OPTION NUMBER :');
    21. na=1;
    22. for f=100:100:1000
    23.     c=3*(10^8);
    24.     lamda=c / (f * (10^6) );
    25.     %beta=(2*pi)/lamda;
    26. tl=tl1*lamda;
    27. ra=ra1*lamda;
    28. hl=tl/2;
    29. nmh=0.5*(nm+1);
    30. dz=2*hl/nm;
    31. zm=hl-0.5*dz;
    32. b=0.5*dz;
    33. a=-0.5*dz;
    34. n=79;
    35. hd=(b-a)/(n+1);
    36. rb=2.3*ra;
    37. tlab=2.0*log(2.3);
    38. for N=1:nm
    39.     zm=hl-(N-0.5)*dz;
    40. for I = 1: nm                
    41.            zn=hl-(I-0.5)*dz;
    42.             za1=zn-zm+a;
    43.           %******************************************************************
    44.           %     FAST ALGORITHM FORM OF THE SIMPSON'S INTEGRAL ROUTINE
    45.           %******************************************************************
    46.           recgp=sqrt(ra*ra+za1*za1);
    47.           cgp1=exp(-j*beta*recgp)*((1.0+j*beta*recgp)*(2.0*recgp*recgp-3.0*ra*ra)+(beta*ra*recgp)^2)/(4.0*pi*beta^2*recgp^5);                              
    48.           zb1=zn-zm+b;
    49.           roc=sqrt(ra*ra+zb1*zb1);
    50.           cgp2=exp(-j*beta*roc)*((1.0+j*beta*roc)*(2.0*roc*roc-3.0*ra*ra)+(beta*ra*roc)^2)/(4.0*pi*beta^2*roc^5);      
    51.           crt=cgp1+cgp2;
    52.           for k = 1: n
    53.                 xk=a+k*hd;
    54.                 zx1=zn-zm+xk;
    55.                 r=sqrt(ra*ra+zx1*zx1);
    56.                 cgp3=exp(-j*beta*r)*((1.0+j*beta*r)*(2.0*r*r-3.0*ra*ra)+(beta*ra*r)^2)/(4.0*pi*beta^2*r^5);            
    57.                 if mod(k,2)~=0
    58.                       crt=crt+4*cgp3;
    59.                 else
    60.                       crt=crt+2*cgp3;
    61.                 end
    62.            end  
    63.         crt=crt*hd*.333333;
    64.         zmn1(N,I)=crt;
    65. end
    66. end
    67. zmn1;
    68. X=imag(zmn1);
    69. R=real(zmn1);
    70. %[J,D]=eig(X,R);
    71. [J,D]=eig(imag(zmn1),real(zmn1));
    72. J;
    73. %i=[1:nm];
    74. %xi=i;
    75. %yi=J(i,1);
    76. %li=J(i,2);
    77. %ki=J(i,3);
    78. %di=J(i,4);
    79. %wi=J(i,8);
    80. %plot(xi,yi,'k',xi,li,'r',xi,ki,'g',xi,di,'y')
    81. %grid on;
    82. %xlabel('segmentation')
    83. %ylabel('EigenVector')
    84. %title('CHARACTERISTIC CURRENT DISTRIBUTION ALONG THE DIPOLE')
    85. %legend('J1','J2','J3','J4')
    86. l1(na,1)=D(1,1);
    87. l1(na,2)=D(2,2);
    88. l1(na,3)=D(3,3);
    89. l1(na,4)=D(4,4);
    90. ms1(na,1)=abs(1 / (1 + (j*l1(na,1)) ) );  %modal significans of the first eigen value
    91. na=na+1;
    92. f;
    93. end
    94. t=length(l1);
    95. k=1:t;
    96. xi=[100:100:1000];
    97. yi1=l1(k,1);
    98. yi2=l1(k,2);
    99. yi3=l1(k,3);
    100. yi4=l1(k,4);
    101. plot(xi,yi1,xi,yi2,xi,yi3,xi,yi4)
    102. legend('D1','D2','D3','D4')
    103. xlabel('Freq in MHz')
    104. ylabel('EigenValue')
    105. title('find the resonant frequency')
    106. figure;
    107. zi1=ms1(k,1);
    108. plot(xi,zi1)
    109.  
     
    Last edited by a moderator: Mar 12, 2013
Loading...