CAPÃTULO 1 - Universidad de Sevilla

More documents

Recommendations

Info

Departamento de Física Aplicada III. Universidad de Sevilladisp (' ')disp (' ')disp ('disp (' ')disp ('disp ('disp ('disp ('ELIJA EL PERFIL DE CURVA DESEADA: ')1) POLÍGONO ')2) POLINOMIO GRADO N QUE PASE POR TODOS LOS PUNTOS ')3) INTERPOLACIÓN POLINOMIAL ')4) INTERPOLACIÓN SPLINE ')disp (' ');opcion=input('OPCIÓN = ');format long; %MAXIMA PRECISION PARA RESOLVER EL SSIT. ECS.switch (opcion)case 1 %POLIGONOSif (tiempos(length(tiempos)) ~= duracion) %ULTIMO PUNTO DISTINTO ADURACION%CALCULAMOS LA PENDIENTE DEL ÚLTIMO TRAMO DE RECTApendiente=(temps(length(temps))-temps(length(temps)-1))/(tiempos(length(tiempos))-tiempos(length(tiempos)-1));%CALCULAMOS EL VALOR DE LA TEMPERAURA FINAL POR INTERPOLACION DELULTIMO TRAMO DE RECTATfinal=pendiente * (duracion-tiempos(length(tiempos))) +temps(length(temps));%AÑADIMOS EL PUNTO FINAL DE LA CURVAtiemposlong=[tiempos, duracion];temps=[temps,Tfinal];clear Tfinal;endif (tiempos(1) ~= 0) %PRIMER PUNTO DISTINTO DE CERO%CALCULAMOS LA PENDIENTE DEL PRIMER TRAMO DE RECTApendiente=(temps(1)-temps(2))/(tiempos(1)-tiempos(2));Tinicial=pendiente * (0-tiempos(1)) + temps(1);tiemposlong=[0,tiemposlong];temps=[Tinicial,temps];clear Tinicial;end%INTERPOLAMOST=interp1(tiemposlong,temps,t);plot(t,T);clear tiemposlong;clear pendiente;case 2 %POLINOMIO GRADO N%SE PLANTEA COMO UN SISTEMA DE ECUACIONES EN EL QUE LAS INCOGNITAS%SON LOS N COEFICIENTES DEL POLINOMIO Y LOS DATOS LOS PUNTOSINTRODUCIDOSY=temps'; %PARA RESOLVER EL SIST TRABAJAMOS CON VECTOR COLUMNAx=[];for (i=1:1:length(Y)-1) % SUBMATRIZ NECESARIA PARA RESOLVER ELSISTEMAx=[x;tiempos.^i];end138
Departamento de Física Aplicada III. Universidad de Sevillaglobal X;X=[ones(1,length(Y));x]'; % MATRIZ DE VALORES DEL TIEMPOclear x;global a;a=X\Y; %RESOLVEMOS EL SISTEMA Y OBTENEMOS LOS COEFICIENTES DELPOLINOMIOclear X Y;T=[];for (i=1:length(t)-1)%UN PUNTO PARA CADA VALOR DE tT(i)=0;for (j=1:length(a)) %MULTIPLICAMOS PORCADA COEFICIENTET(i)=T(i)+a(j)*(t(i)^(j-1));endT=[T,T(i)];endplot(t,T);clear a;clear i,j;case 3 %INTERPOLACIÓN POLINOMIALdisp (' ')disp (' ')disp ('INTRODUZCA EL GRADO DEL POLINOMIO DE INTERPOLACIÓN ')Grado = ');N=input('a=polyfit(tiempos,temps,N);T=[];for (i=1:length(t)-1) %UN PUNTO PARA CADA VALOR DE tT(i)=0;for (j=length(a):-1:1)T(i)=T(i)+a(j)*(t(i)^(length(a)-j));endT=[T,T(i)];endplot(t,T);clear N;clear a;clear i,j;case 4 %SPLINET=spline(tiempos,temps,t);plot(t,T);end %end switchhold off;clear option;clear temps;tiempos=[];concents=[];disp (' ');disp (' ');disp (' ');%REPETIMOS PARA LAS CONCENTRACIONES139
Page 3 and 4:
Departamento de Física Aplicada II
Page 5 and 6:
Page 7 and 8:
Page 11 and 12:
Page 13 and 14:
Page 15 and 16:
Page 17 and 18:
Page 19 and 20:
Page 21 and 22:
Page 23 and 24:
Page 25 and 26:
Page 27 and 28:
Page 29 and 30:
Page 31 and 32:
Page 33 and 34:
Page 35 and 36:
Page 37 and 38:
Page 39 and 40:
Page 41 and 42:
Page 43 and 44:
Page 45 and 46:
Page 47 and 48:
Page 49 and 50:
Page 51 and 52:
Page 53 and 54:
Page 55 and 56:
Page 57 and 58:
Page 59 and 60:
Page 61 and 62:
Page 63 and 64:
Page 65 and 66:
Page 67 and 68:
Page 69 and 70:
Page 71 and 72:
Page 73 and 74:
Page 75 and 76:
Page 77 and 78:
Page 79 and 80:
Page 81 and 82:
Page 83 and 84:
Page 85 and 86:
Page 87 and 88: Departamento de Física Aplicada II
Page 137: Departamento de Física Aplicada II
Page 189 and 190:
Page 191 and 192:
Page 193 and 194:
Page 195 and 196:
Page 197 and 198:
show all

CAPÃTULO 1 - Universidad de Sevilla

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?

CAPÃTULO 1 - Universidad de Sevilla