I have simulated an structure on IE3D from which i got S11 and s21 ..................

then i want to use MATLAB to calculate effective values of permeability and permittivity in order to check type of metamaterial (Negative refractive index ).......I am using Method explained in physical review paper by d r smith..........

***************please check am i going right with this formula or I have to change the approach ****************

% extracted values of frequency %

f=[1.00000000,1.67346939, 1.89795918,2.34693878, 2.79591837, 3.02040816, 3.24489796, 3.69387755 ,4.36734694 ,4.81632653 ,5.26530612,5.71428571,6.38775510 , 7.06122449 , 7.51020408 , 8.85714286 , 9.30612245 , 9.75510204 ,10.20408163 ,10.65306122 , 10.87755102 , 11.10204082 , 11.77551020 ];

% s11 parameter %

s11=[-10.2096766593,-7.6772215391,-7.5990323576,-8.8890690294,-19.7071254719,-11.1259166031, -3.3145373349 , -0.7839962323 , -3.4421138099 , -3.9468612419, -5.3320037502, -7.5585127979, -14.8891887486 , -19.9105565933 , -12.7098476122 ,-14.0834418139 , -32.4189096035 , -10.1206097795, -5.5714529494 ,-2.2400924818 , -1.2289516205 , -0.8503614427 , -5.7961227275];

%extracted s21 parameter from simulation%

s21=[-0.7915865485,-1.6875595060,-1.9007912686,-2.0296233434,-1.5034860028,-1.2149884970 , -5.4510835084, -26.7766971832 , -5.9835155484 , -5.2473000329 , -4.2659963387 , -3.1587227525 , -1.3988610332 , -0.3081867986, -1.3291504276 , -1.9306769742 , -0.6652145887 , -2.6538876439, -7.1033229681, -12.0406416672 , -11.7075593833 , -8.5077506090 ,-1.7398624551];

% formula from IEEE paper %

n=(acos((1./(2.*s21)).*(1-(s11).^2-(s21).^2)))./1.6; % Index %

z=sqrt(((1+s11).^2)-(s21).^2)./(((1-s11).^2)-(s21).^2); % Impedence %

subplot(3,1,1)

plot(f, real(n), 'blue', f, imag(n), 'r');

title('index');

u=n.*z; % permeability %

e=n./(z); % permittivity %

subplot(3,1,2)

plot(f, real(u), 'blue', f, imag(u), 'r');

title('permeability');

u=n.*z;

subplot(3,1,3);

plot(f, real(e), 'blue', f, imag(e), 'r');

title('permittivity');

Similar questions and discussions