29/07/2016, 14:15
A=dlmread('time_epsilon_strain_25N.txt') % leggo la matrice da file
time=A(:,1) % formo vettore tempo da matrice A
epsilon=A(:,2) % formo vettore deformazione da matrice A
scatter(time,epsilon,20) % plottaggio punti sperimentali
xlabel('time')
ylabel('epsilon')
title('PP 25N')
E=13000; % modulo elastico(valore da sostituire da quello ricavato in laib)
eta = 5500; %viscosità
vettore=[1:200];
j=1;
for i = 0:1:199
x = time(1+i);
y = epsilon(1+i);
y_appr = sigma_0/E + (sigma_0/eta)*x;
errore = abs(y-y_appr)
vettore(j)=y_appr;
j=j+1;
if errore>0.08
fprintf('non va bene')
end
end
scatter(time,epsilon,20) ; % plottaggio punti sperimentali
xlabel('time');
ylabel('epsilon');
title('Linoleum 25N');
hold on
plot(time,vettore)
Skuola.net News è una testata giornalistica iscritta al Registro degli Operatori della Comunicazione.
Registrazione: n° 20792 del 23/12/2010.
©2000—
Skuola Network s.r.l. Tutti i diritti riservati. — P.I. 10404470014.
Powered by phpBB © phpBB Group - Privacy policy - Cookie privacy
phpBB Mobile / SEO by Artodia.