Pulling data from a for loop

3 vues (au cours des 30 derniers jours)
Justin Hayes
Justin Hayes le 27 Mar 2020
Commenté : Justin Hayes le 28 Mar 2020
I want to subtract Tskin at t = 10 from Tskin at t = 9. How do I do this?
T = 29; % %C
Tambient = 20; % C
Tfilm = (T + Tambient)/2;
for t = 1:1:100
k_air_2 = (2.4131 + 0.0078 * T) * 10^-2; % W/m*C
rho_air = 1.2893 - 0.004 * Tfilm; % kg/m^3
mu_air = (0.0047 * Tfilm + 1.7164) * 10^-5; % kg/m*s
k_air_ambient= (2.4131 + 0.0078 * Tfilm) * 10^-2; % W/m*C
Re = Diameter_forearm * V * rho_air / mu_air;
h_air = (k_air_ambient/Diameter_forearm) .* (0.3 + ((0.62 * Re.^.5 * Pr.^.333)/((1 + (0.4 * Pr).^.667).^.25)).*(1 + ((Re/282000).^(5/8))).^.8) % J/s*K*m^2
Uskin = (k_air_2/thickness_air) + (k_cotton/thickness_shirt) + (h_air) % J/s*K*m^2
Uambient = (k_air_ambient/thickness_air) + (k_cotton/thickness_shirt) + (h_air)
Tskin = [(T * Area * Cp_air * rho_air)/(t-(t-1)) + Qgen + (Uambient * Area * Tambient)]/((Area * Cp_air * rho_air / ((t-(t-1)))) + Uskin * Area);
T = Tskin
Tfilm = (T + Tambient)/2
fprintf('Y is %.2f C after %d seconds ' ,Tskin,t)
end

Réponse acceptée

Walter Roberson
Walter Roberson le 27 Mar 2020
T = 29; % %C
Tambient = 20; % C
Tfilm = (T + Tambient)/2;
Tskin = nan;
for t = 1:1:100
k_air_2 = (2.4131 + 0.0078 * T) * 10^-2; % W/m*C
rho_air = 1.2893 - 0.004 * Tfilm; % kg/m^3
mu_air = (0.0047 * Tfilm + 1.7164) * 10^-5; % kg/m*s
k_air_ambient= (2.4131 + 0.0078 * Tfilm) * 10^-2; % W/m*C
Re = Diameter_forearm * V * rho_air / mu_air;
h_air = (k_air_ambient/Diameter_forearm) .* (0.3 + ((0.62 * Re.^.5 * Pr.^.333)/((1 + (0.4 * Pr).^.667).^.25)).*(1 + ((Re/282000).^(5/8))).^.8) % J/s*K*m^2
Uskin = (k_air_2/thickness_air) + (k_cotton/thickness_shirt) + (h_air) % J/s*K*m^2
Uambient = (k_air_ambient/thickness_air) + (k_cotton/thickness_shirt) + (h_air)
old_Tskin = Tskin;
Tskin = [(T * Area * Cp_air * rho_air)/(t-(t-1)) + Qgen + (Uambient * Area * Tambient)]/((Area * Cp_air * rho_air / ((t-(t-1)))) + Uskin * Area);
T = Tskin
Tfilm = (T + Tambient)/2
fprintf('Y is %.2f C after %d seconds ' ,Tskin,t)
if t==10
fprintf('Tskin(10) - Tskin(9) = %g\n', Tskin - old_Tskin);
end
end
  1 commentaire
Justin Hayes
Justin Hayes le 28 Mar 2020
Thank you very much, it works!

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Startup and Shutdown dans Help Center et File Exchange

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by