Iteracion dos pasos Matlab
Publicado por Marco (1 intervención) el 01/04/2019 20:38:42
hola miembros del foto, he escrito un fichero script en Matlab para hacer una iteración en dos pasos.
he definido variables y escrito las ecuaciones de los ajustes. como puedo iterar los ajustes los ajustes machacando los valores de A y B?
el Script es el siguiente:
lo que hay que hacer son iteraciones que calculen A y B usando los valores calculados en el paso anterior.
Muchas Gracias señores
he definido variables y escrito las ecuaciones de los ajustes. como puedo iterar los ajustes los ajustes machacando los valores de A y B?
el Script es el siguiente:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
%%%%%%%%%%%%
P=input('introduce los valores de P ');
Vm=input('introduce los valores de volumen molar ');
T=input('introduce los valores de T ');
R=input('introduce la constate de los gases R ');
B=input('introduce el termino inicial B ');
%Calculamos las cordenadas X,Y del primer paso.
Xa=1./(Vm.^2);
Ya=((R.*T)/(Vm-B))-P;
scatter(Xa,Ya)
%calculamos la pendiente y el termino independiente.primer ajuste
n=length(Xa)
C1=(sum(Ya)*sum(Xa.^2))-(sum(Xa)*sum(Xa.*Ya))/((n*sum(Xa.^2))-(sum(Xa)^2))
A=(n*sum(Xa.*Ya))-(sum(Xa)*sum(Ya))/((n*sum(Xa.^2))-(sum(Xa)^2))
%calculamos coordenadas segundo ajuste
Xb=P+(A./(Vm.^2));
Yb=(P+(A./(Vm.^2)))*Vm;
%calculamos la pendiente y el termino independiente. 2do ajuste
C2=(sum(Yb)*sum(Xb.^2))-(sum(Xb)*sum(Xb.*Yb))/((n*sum(Xb.^2))-(sum(Xb)^2))
B=(n*sum(Xb.*Yb))-(sum(Xb)*sum(Yb))/((n*sum(Xb.^2))-(sum(Xb)^2))
%%%%%%%%%%%%%%%%
lo que hay que hacer son iteraciones que calculen A y B usando los valores calculados en el paso anterior.
Muchas Gracias señores
Valora esta pregunta
0