1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
| function essai_test
clc
clear all
close all
global val
sld = uicontrol('Style', 'slider',...
'Min',1,'Max',100,'Value',100,...
'Position', [100 350 120 20],...
'SliderStep',[0.1 0.1],...
'Callback', @slid);
val=2.5;
%test = xlsread('S11.xlsx');
test=-35*rand(1,10);
tg_perte = 0.025;
PI = pi;
C = 3.0 * 10.^8;
f = linspace(6,18,801)*10.^9;
t=7.5;
lambda_0 = (C ./ f)*10.^3;
AS11V = test(:,1);
plot(f,AS11V,'r');
hold on
h=plot(0,0);
grid
title('Mesure paramètre S11 vertical')
while val>1 % FIN pour eps-r=1
eps_r = val;
epsilon_etoile = eps_r*(1 - 1i*tg_perte);
lambda_g = lambda_0 ./ sqrt(eps_r);
gamma = (2*PI./lambda_g)*sqrt(1 - 1i*tg_perte);
%
R1p = (1 - sqrt(epsilon_etoile))/(1 + sqrt(epsilon_etoile));
R2p = (sqrt(epsilon_etoile) - 1)/(1 + sqrt(epsilon_etoile));
%
exponentiel = exp(-1i*2*gamma.*t);
%
REFLEXION = (R1p + R2p.*exponentiel)./(1 + R2p*R1p.*exponentiel);
%plot(f,20*log10(abs(REFLEXION)))
set(h,'Xdata',f,'Ydata',20*log10(abs(REFLEXION)))
pause(1e-3)
end
disp('FIN')
function slid(source,event)
global val
val=get(source,'Value')/20;
disp(val) |
Partager