Docsity
Docsity

Prepara tus exámenes
Prepara tus exámenes

Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity


Consigue puntos base para descargar
Consigue puntos base para descargar

Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium


Orientación Universidad
Orientación Universidad


Prova octave, Apuntes de Química

Asignatura: Recursos informatics, Profesor: Anna de juan, Carrera: Química, Universidad: UB

Tipo: Apuntes

2014/2015

Subido el 15/07/2015

kevinrib
kevinrib 🇪🇸

4

(208)

41 documentos

1 / 2

Toggle sidebar

Esta página no es visible en la vista previa

¡No te pierdas las partes importantes!

bg1
Prova OCTAVE 2015
1.
a.
temps=[0 100 200 400 600 800 1000]
temps =
0 100 200 400 600 800 1000
octave:6> I=[1.5e-5 6.7e-6 3.6e-6 1.4e-6 4.7e-7 2.6e-7 9.3e-8]
I =
Columns 1 through 6:
1.5000e-05 6.7000e-06 3.6000e-06 1.4000e-06 4.7000e-07 2.6000e-07
Column 7:
9.3000e-08
octave:7> Integral=trapz(temps,I)
Integral = 0.0023953
b.
tempsint=[0:10:1000];
Iinterp=interp1(temps,I,tempsint,'spline');
Integralinterp=trapz(tempsint,Iinterp)
Integralinterp = 0.0022545
2.
a.
T=[250:20:400];
V=[20:5:100];
R=0.082;
B=0.12
C=0.012
[VV,TT]=meshgrid(V,T)
P=(R*TT).*(1./VV+(B./(VV.^2))+(C./(VV.^3)));
mesh(V,T, P )
xlabel('T')
ylabel('V')
zlabel('P')
print(fig1.gif’,’-dgif’)
pf2

Vista previa parcial del texto

¡Descarga Prova octave y más Apuntes en PDF de Química solo en Docsity!

Prova OCTAVE 2015 1. a. temps=[0 100 200 400 600 800 1000] temps = 0 100 200 400 600 800 1000 octave:6> I=[1.5e-5 6.7e-6 3.6e-6 1.4e-6 4.7e-7 2.6e-7 9.3e-8] I = Columns 1 through 6: 1.5000e-05 6.7000e-06 3.6000e-06 1.4000e-06 4.7000e-07 2.6000e- Column 7: 9.3000e- octave:7> Integral=trapz(temps,I) Integral = 0. b. tempsint=[0:10:1000]; Iinterp=interp1(temps,I,tempsint,'spline'); Integralinterp=trapz(tempsint,Iinterp) Integralinterp = 0.

a. T=[250:20:400]; V=[20:5:100]; R=0.082; B=0. C=0. [VV,TT]=meshgrid(V,T) P=(RTT).(1./VV+(B./(VV.^2))+(C./(VV.^3))); mesh(V,T,P) xlabel('T') ylabel('V') zlabel('P') print(‘fig1.gif’,’-dgif’)

b.

R=0. C=0. B=0. T= P=0. vi=(RT)/P function fx=f(V) R=0.082; B=0.12; C=0.02; T=350; P=0.7; fx=(RT(1/V+(B/V^2)+(C/(V^3))))-P; endfunction* [V,fv,m]=fsolve("f", 41) Resultats numèrics V=41.120 L fv= m= c. R=0. T= V=[20:5:95] P=[1.49 1.19 0.99 0.85 0.75 0.66 0.59 0.54 0.49 0.46 0.42 0.38 0.36 0.34 0.32 0.30] rtv1=(RT)./V rtv2=(RT)./(V.^2) vindep=[rtv1;rtv2]; coefs=P/vindep coefs = 0.89110 0.39989 [A B]