Obtener el valor de una función
Publicado por Vladimir Sousa (2 intervenciones) el 30/01/2018 00:04:22
Estimados compañeros, estoy calculando la integral de una función mediante matemática simbólica con la función (syms). En el momento de sustituir la expresión con la función (subs) obtengo una expresión nuevamente pero compuesta solo por números. Necesito por favor que me indiquen que debo hacer para que me de el valor numérico solamente y no la expresión.
Abajo le coloco el srcipt.
Agradezco mucho la colaboración.
Esto me da como solución para los valores sustituidos:
Tags =
-(3^(1/2)*((((3069529069953393*sin(3638376646725481/1125899906842624))/562949953421312 - (133*2^(1/2)*sin((163*pi)/450))/30)*((1366332973168105*cos(2365131654361039/562949953421312))/4398046511104 + (1267715505890751309*cos(3638376646725481/1125899906842624))/56294995342131200 - (2^(1/2)*exp((pi*31i)/225)*54929i)/3000 + (54929*2^(1/2)*exp((pi*163i)/450))/3000))/(120*pi) - (((3069529069953393*sin(3638376646725481/1125899906842624))/562949953421312 + (133*2^(1/2)*sin((163*pi)/450))/15)*(4430665693088175824955*cos(3638376646725481/1125899906842624) - 3602445386007332257792*2^(1/2)*sin((31*pi)/225) + 46063128596874337680000))/(17530261549539655680000*pi)))/3
De modo que no me da el número como resultado de esto.
Abajo le coloco el srcipt.
Agradezco mucho la colaboración.
1
2
3
4
5
6
7
8
9
10
11
12
13
l=input('Estado de carga; l=');
Rs=4.13;
r1=Rs*((234.5+115)/(234.5+25));
Pmecn=38.5;
Pad=19.8;
Polos=2;syms t
Vab=VabRMS(l).*sqrt(2).*sin(2*60.*pi.*t);
Vca=VcaRMS(l).*sqrt(2).*sin(2*60*pi.*t+(AngVcaRMS(l).*pi./180));
Ia=IaRMS(l).*sqrt(2).*sin(2.*60.*pi.*t+(AngIaRMS(l).*pi./180));
Ib=IbRMS(l).*sqrt(2).*sin(2*60.*pi.*t+(AngIbRMS(l).*pi./180));
Tag=((sqrt(3)*Polos)./6).*((Ia-Ib).*int(Vca+Rs.*(2.*Ia+Ib))+(2.*Ia+Ib).*int(Vab-r1.*(Ia-Ib)));
Tags=abs (subs(((sqrt(3)*Polos)./6).*((Ia-Ib).*int(Vca+Rs.*(2.*Ia+Ib))+(2.*Ia+Ib).*int(Vab-r1.*(Ia-Ib))),0))
Esto me da como solución para los valores sustituidos:
Tags =
-(3^(1/2)*((((3069529069953393*sin(3638376646725481/1125899906842624))/562949953421312 - (133*2^(1/2)*sin((163*pi)/450))/30)*((1366332973168105*cos(2365131654361039/562949953421312))/4398046511104 + (1267715505890751309*cos(3638376646725481/1125899906842624))/56294995342131200 - (2^(1/2)*exp((pi*31i)/225)*54929i)/3000 + (54929*2^(1/2)*exp((pi*163i)/450))/3000))/(120*pi) - (((3069529069953393*sin(3638376646725481/1125899906842624))/562949953421312 + (133*2^(1/2)*sin((163*pi)/450))/15)*(4430665693088175824955*cos(3638376646725481/1125899906842624) - 3602445386007332257792*2^(1/2)*sin((31*pi)/225) + 46063128596874337680000))/(17530261549539655680000*pi)))/3
De modo que no me da el número como resultado de esto.
Valora esta pregunta
0
