1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
| %évolution de etat_capt en fonction de theta_x pour plusieurs valeurs de
%alpha_c donné
%declaration des constantes
eta_c=0.9;%rendement optique du capteur
k_r=3;
rho_r=0.4;
alpha_c=[40 60 80 100];%valeurs de alpha_c donné
T_0=35;%temperature ambiante
T_r=300;%temperatude du recepteur
p_x=160;
%initialisation
for i=1:1:4
theta_r=T_r/T_0;
eta_capt(i)=eta_c-(k_r/alpha_c(i))*(theta_r-1)-(rho_r/alpha_c(i))*(theta_r^4-1);
theta_x(i)=theta_r-eta_capt(i)*(alpha_c(i)/p_x);
end
plot(eta_capt(i),theta_x(i)) |
Partager