1 - Enunciado
Nos dan un circuito R-L-C como el de la figura:
Dados v, f y w; hallar i1, i2.
2 - Obtención del modelo
Mesh method. 2ª Ley de Kirchoff.
{v = (z1*i1 + z3*i1) - (z3*i2)
{0 = (z2*i2 + z3*i2 + z4*i2) - (z3*i1)
3 - Relación con la teoría de clase
3.1. - Sistema de 1º Orden:
* Impulso de Dirac:
- syms A; syms a; syms s; syms g; syms y1; syms Y1; syms u;
- g = (A/(s+a));
- A = 10; a = 3;
- u = 1/s;
- Y1 = g*u;
- y1 = ilaplace(Y1)
- figure
- ezplot(y1,[0,3]), axis([0,2,0,12])
* Escalón Unidad:
- syms A; syms a; syms s; syms g; syms y2; syms Y2; syms us;
- g = (A/(s+a));
- A = 10; a = 3;
- u = 1;
- Y2 = g*u;
- y2 = ilaplace(Y2);
- figure
- ezplot(y2,[0,3]), axis([0,2,0,12])
3.2. - Sistema de 2º Orden:
* Impulso de Dirac:
- syms wn; syms e; syms s; syms u; syms g;
- e = 0.2;
- u = 1/s;
- g = wn^2/(s^2+2*E*wn*s+wn^2);
- Y1 = g*u;
- y1 = ilaplace(Y1);
- y13 = vpa(y1,3);
- figure(1);
- ezplot(y13,[0,3]),axis([0,3,0,4])
* Escalón Unidad:
- syms wn; syms e; syms s; syms u; syms g;
- e = 0.2;
- u = 1;
- g = wn^2/(s^2+2*E*wn*s+wn^2);
- Y2 = g*u;
- y2 = ilaplace(Y1);
- y23 = vpa(y2,3);
- figure(2);
- ezplot(y23,[0,3]),axis([0,3,0,4])
4 - Resolución con Matlab
4.1. - Sistema de 1º Orden:
* Impulso de Dirac:
ys = A/(s*(a + s))
* Escalón Unidad:
ys = A/(a + s)
4.2. - Sistema de 2º Orden:
* Impulso de Dirac:
* Escalón Unidad:
lunes, 28 de marzo de 2011
lunes, 14 de marzo de 2011
Práctica 3: Matlab
1 - Enunciado
Nos dan un circuito R-L-C como el de la figura:
Dados v, f y w; hallar i1, i2.
2 - Obtención del modelo
Mesh method. 2ª Ley de Kirchoff.
{v = (z1*i1 + z3*i1) - (z3*i2)
{0 = (z2*i2 + z3*i2 + z4*i2) - (z3*i1)
3 - Relación con la teoría de clase
- syms R; syms w; syms i1; syms i2;
- v = 10; L = 10e-3; C = 47e-6;
- z1 = R; z2 = 0+ (w*L)*1i; z3 = 0+(1/(1i*w*C)); z4 = 2;
- e1 = z1*i1+z3*i1-z3*i2-v;
- e2 = -z3*i1+z3*i2+z2*i2+z4*i2;
- [i1 i2] = solve (e1,e2,i1,i2)
4 - Resolución con Matlab
i1 =
(10*(6935975771714791*w^2*i + 1387195154342958200*w - 14757395258967641292800*i))/(147573952589676412928*w - 14757395258967641292800*R*i + 1387195154342958200*R*w + 6935975771714791*R*w^2*i - 29514790517935282585600*i)
i2 =
-(147573952589676412928000*i)/(147573952589676412928*w - 14757395258967641292800*R*i + 1387195154342958200*R*w + 6935975771714791*R*w^2*i - 29514790517935282585600*i)
Después de ejecutar el comando vpa():
vpa(i1,3)
ans =
(10.0*(6.94*10^15*w^2*i + 1.39*10^18*w - 1.48*10^22*i))/(1.48*10^20*w - 1.48*10^22*R*i + 1.39*10^18*R*w + 6.94*10^15*R*w^2*i - 2.95*10^22*i)
vpa(i2,3)
ans =
-(1.48*10^23*i)/(1.48*10^20*w - 1.48*10^22*R*i + 1.39*10^18*R*w + 6.94*10^15*R*w^2*i - 2.95*10^22*i)
Nos dan un circuito R-L-C como el de la figura:
Dados v, f y w; hallar i1, i2.
2 - Obtención del modelo
Mesh method. 2ª Ley de Kirchoff.
{v = (z1*i1 + z3*i1) - (z3*i2)
{0 = (z2*i2 + z3*i2 + z4*i2) - (z3*i1)
3 - Relación con la teoría de clase
- syms R; syms w; syms i1; syms i2;
- v = 10; L = 10e-3; C = 47e-6;
- z1 = R; z2 = 0+ (w*L)*1i; z3 = 0+(1/(1i*w*C)); z4 = 2;
- e1 = z1*i1+z3*i1-z3*i2-v;
- e2 = -z3*i1+z3*i2+z2*i2+z4*i2;
- [i1 i2] = solve (e1,e2,i1,i2)
4 - Resolución con Matlab
i1 =
(10*(6935975771714791*w^2*i + 1387195154342958200*w - 14757395258967641292800*i))/(147573952589676412928*w - 14757395258967641292800*R*i + 1387195154342958200*R*w + 6935975771714791*R*w^2*i - 29514790517935282585600*i)
i2 =
-(147573952589676412928000*i)/(147573952589676412928*w - 14757395258967641292800*R*i + 1387195154342958200*R*w + 6935975771714791*R*w^2*i - 29514790517935282585600*i)
Después de ejecutar el comando vpa():
vpa(i1,3)
ans =
(10.0*(6.94*10^15*w^2*i + 1.39*10^18*w - 1.48*10^22*i))/(1.48*10^20*w - 1.48*10^22*R*i + 1.39*10^18*R*w + 6.94*10^15*R*w^2*i - 2.95*10^22*i)
vpa(i2,3)
ans =
-(1.48*10^23*i)/(1.48*10^20*w - 1.48*10^22*R*i + 1.39*10^18*R*w + 6.94*10^15*R*w^2*i - 2.95*10^22*i)
Suscribirse a:
Entradas (Atom)