# 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;
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