Preguntas urgentes sobre la implementación de circuitos en matlab
G1=1/10;G2=1/20;G3=1/10;G4=1/10;
us=20;is=1; p>
Gn=; %Establecer matriz de coeficientes B
X1=A\B*[ is1;0];Uoc=X1(3)
X2 =A \B*[0;0;1];Req=X2(3)
RL=Req;P=Uoc^2*RL/(Req+RL)^2 %Encontrar su potencia máxima
RL=0:0.01:10;P=(RL*Uoc./(Req+RL)).*Uoc./(Req+RL); %Establece la secuencia RL para encontrar su potencia
figura(1),plot(RL,P), grid %Trazar la curva de cambio del consumo de energía con RL
Convolución
>> p=0.1;
t=0:p:10;
f1=exp( -2*t).*u(t);
f2=u(t); )-u (t-4);
f=conv(f1,f2);
subtrama(1,3,1);
trama( t,f1 ,'r');
title('f1(t)=e^-2*t*u(t)');
xlabel('t( seg)' );
ylabel('f1(t)');
subplot(1,3,2);
plot(t,f2 ,'g ');
title('f2(t)=u(t)-u(t-4)');
xlabel('t(sec)' );
ylabel('f2(t)');
subplot(1,3,3);
plot(f);
title('f(t)=f1(t)*f2(t)');
xlabel('t(sec )');
ylabel( 'f( t)');
Dinámica de primer orden
1, respuesta de entrada cero de secuencia RC
U0=2;R=2;C= 0.5; % Ingrese los parámetros dados
U1=3;R1=3;C1=0.5; %Ingrese los parámetros dados
t=[0:0.1 :5]; el rango de tiempo
Uc1=U0*exp(-t/(R*C)); Uc2=U1*exp(-t/(R*C)); %Calcule el valor de voltaje del capacitor
Ur1=U0*exp(-t/(R*C)); Ur2=U1*exp(-t/(R*C) % Calcular el valor de voltaje de resistencia
I1=U0/ R*exp(-t /(R*C)); I2=U1/R*exp(-t/(R*C)); %Calcular valor actual
Pc1=U0^2/R *exp(-2* t/(R*C));Pc2=U1^2/R*exp(-2*t/(R*C)); Calcular el valor de potencia del condensador
Pr1 =U0^2/R* exp(-2*t/(R*C)); Pr2=U1^2/R*exp(-2*t/(R*C)); % valor de potencia de resistencia calculado
Figura
subplot(5,1,1);plot(t,Uc1,t,Uc2)
title('Uc (t) forma de onda')
subtrama (5,1,2);plot(t,Ur1,t,Ur2)
título('Ur(t) forma de onda')
subtrama (5,1,3) ;plot(t,I1,t,I2)
title('I(t) forma de onda')
subplo
t(5,1,4);plot(t,Pc1,t,Pc2)
title('Forma de onda Pc(t)')
subplot(5,1, 5 );plot(t,Pr1,t,Pr2)
title('Forma de onda de Pr(t)')
2. Respuesta de entrada cero paralelo RL
I0=2;R=2;L=0.5; %Ingresar parámetros dados
I1=3;R1=3;L1=0.5; %Ingresar parámetros dados
t= [0:0.05:1.5]; %Determinar el rango de tiempo
IL1=I0*exp(-t*R/L);IL2=I1*exp(-t*R/L); el valor de corriente del inductor
Ir1=I0*exp(-t*R/L); Ir2=I1*exp(-t*R/L %Calcular el valor de corriente de la resistencia
<); p>U1=I0*R*exp(-t*R/L);U2=I1*R*exp(-t*R/L);% valor de voltaje calculadoPL1=I0^2 *R*exp(-2*t*R/L);PL2=I1^2*R*exp(-2*t*R/L); %Calcular el valor de potencia del inductor
Pr1= I0 ^2*R*exp(-2*t*R/L); Pr2=I1^2*R*exp(-2*t*R/L % Calcular el valor de potencia de la resistencia
subplot(5,1,1);plot(t,IL1,t,IL2)
title('Gráfico de forma de onda de IL(t)')
subplot(5,1,2);plot(t,Ir1,t,Ir2)
title('Gráfico de forma de onda de Ir(t)')
subplot(5,1,3); plot(t,U1,t,U2)
title('Gráfico de forma de onda de U(t)')
subplot(5 ,1 ,4);plot(t,PL1,t,PL2)
title('Forma de onda de PL(t)')
subplot(5,1,5);plot(t, Pr1, t,Pr2)
title('Forma de onda de Pr(t)')
3, respuesta de estado cero de la serie RC
Us=2; R=2 ;C=0.5; %Ingrese los parámetros dados
Us1=3;R1=3;C1=0.5;
t=[0:0.1:10]; /p>
Uc1=Us*(1-exp(-t/(R*C))); Uc2=Us1*(1-exp(-t/(R*C)));
Ur1=Us*exp(-t/(R*C ))); Ur2=Us1*exp(-t/(R*C));
I1=Us/R*exp); (-t /(R*C)); I2=Us1/R*exp(-t/(R*C));
Pc1=Us^2/R*(exp(-t/ (R* C)) -exp(-2*t/(R*C)));
Pc2=Us1^2/R*(exp(-t/(R*C))- exp(- 2*t/(R*C)));
Pr1=Us^2/R*exp(-2*t/(R*C));Pr2=Us1^2/ R*exp (-2*t/(R*C));
Figura
subplot(5,1,1);plot(t,Uc1,t,Uc2)
title('Gráfico de forma de onda de Uc(t)')
subplot(5,1,2);plot(t,Ur1,t,Ur2)
title ('Gráfico de forma de onda de U
r(t)')
subplot(5,1,3);plot(t,I1,t,I2)
title('Forma de onda de I(t)')
subplot(5,1,4);plot(t,Pc1,t,Pc2)
title('Diagrama de forma de onda de Pc(t)')
subplot(5,1,5);plot(t,Pr1,t,Pr2)
title('Diagrama de forma de onda de Pr(t)')
4, RL Respuesta de estado cero paralelo
Is=2;R=2;L=0.5; %Parámetros de entrada dados
Is1=3;R1=3;L1=0.5;
t=[0:0.05:1.5];
IL1=Es*(1-exp(-t*R/L));IL2=Es1*(1- exp (-t*R/L));
Ir1=Es*exp(-t*R/L);Ir2=Es1*exp(-t*R/L);
U1=Es*R*exp(-t*R/L);U2=Es1*R*exp(-t*R/L);
PL1=Es^2*R* ( exp(-t*R/L)-exp(-2*t*R/L));PL2=Is1^2*R*(exp(-t*R/L)-exp(-2*t* R /L));
Pr1=Es^2*R*exp(-2*t*R/L);Pr2=Es1^2*R*exp(-2*t*R/L);< / p>
figura
subplot(5,1,1);plot(t,IL1,t,IL2)
title('Gráfico de forma de onda de IL(t) ' )
subplot(5,1,2);plot(t,Ir1,t,Ir2)
title('Forma de onda de Ir(t)')
subplot(5,1,3);plot(t,U1,t,U2)
title('Forma de onda de U(t)')
subplot(5,1,4) ;plot(t,PL1,t,PL2)
title('Forma de onda de PL(t)')
subplot(5,1,5);plot(t, Pr1 ,t,Pr2)
title('Forma de onda de Pr(t)')
5, respuesta completa de la serie RC
U0=2;Us= 3 ;R=2;C=0.5; %Ingrese los parámetros dados
U1=2.5;Us1=3;R1=3;C1=0.5;
t=[0: 0.1 :5];
Uc1=U0*exp(-t/(R*C))+Us*(1-exp(-t/(R*C)));
Uc2=U1*exp(-t/(R*C))+Us1*(1-exp(-t/(R*C)));
Ur1=Us*exp( - t/(R*C))-U0*exp(-t/(R*C));
Ur2=Us1*exp(-t/(R*C));
Ur2=Us1*exp(-t/(R*C))-U1*exp(-t/(R*C));
I1=(Us-U0)/R * exp(-t/(R*C));I2=(Us1-U1)/R*exp(-t/(R*C));
figura(1)
subplot(3,1,1);plot (t,Uc1,t,Uc2)
title('Gráfico de forma de onda de Uc(t)')
subplot( 3 ,1,2);plot(t,Ur1,t,Ur2)
t
itle('Gráfico de forma de onda de Ur(t)')
subplot(3,1,3);plot(t,I1,t,I2)
title('Gráfico de forma de onda de I(t)')
6, RL respuesta completa paralela
I0=2;Is=3;R=2;L=0.5 %Parámetros dados de entrada<; /p> p>
I1=2.5;Is1=3;R1=3;L1=0.5;
t=[0.01:1.5];
IL1= I0*exp (-t*R/L)+Is*(1-exp(-t*R/L));
IL2=I1*exp(-t*R/L)+Is1 *(1 -exp(-t*R/L));
Ir1=Es*exp(-t*R/L)-I0*exp(-t*R/L);
Ir2=Is1*exp(-t*R/L)-I1*exp(-t*R/L)
U1=(Es-I0)*R*exp(); -t* R/L); U2=(Is1-I1)*R*exp(-t*R/L);
figura(1)
subtrama(3, 1,1 );plot(t,IL1,t,IL2)
title('Forma de onda de IL(t)')
subplot(3,1,2);plot (t, Ir1,t,Ir2)
title('Forma de onda de Ir(t)')
subplot(3,1,3);plot(t,U1,t ,U2)
title('Diagrama de forma de onda de U(t)')
7. Descomposición de respuesta completa
U0=2.5;Us=3.5;I0 =2; Is=3;R=2;L=0.5;C=1;%Ingrese los parámetros dados
t=[0:0.01:10];
Uc= U0*exp (-t/(R*C))+Us*(1-exp(-t/(R*C))))
Uc1=U0*exp(-t/(R *C) ); Uc2=Us*(1-exp(-t/(R*C)));
Uc3=Uc4=(U0-Us)*exp(-t/( R*C ));
IL=I0*exp(-t*R/L)+Is*(1-exp(-t*R/L)); IL1=I0 *exp(-t*R/L); IL2=Is*(1-exp(-t*R/L));
IL3= IL4=(I0-Is) *exp( -t*R/L);
figura(1)
subplot(4,1,1);plot(t,Uc,t,Uc1,t, Uc2)
title('Uc(t) respuesta completa = entrada cero + estado cero')
subplot(4,1,2); IL1,t ,IL2)
title('IL(t) respuesta completa = entrada cero + estado cero')
subplot(4,1,3);plot(t, Uc,t,Uc3,t,Uc4)
title('Uc( t respuesta completa = componente transitorio + componente de estado estacionario')
subplot(4,1,4) ;plot( t,IL,t,IL3,t,IL4)
title('IL(t) respuesta completa = componente transitorio + componente de estado estacionario')
8. Serie RC seno Excitación respuesta de estado cero
Usm=2;w =pi;R=2;C=0.5;h=atan(w*C*R);z=sqrt((w*R* C)^ 2+1);%Ingrese los parámetros dados
t=[0:0.01:4];
Us=Usm*cos(w*t+pi/2);
Uc=Usm/z*cos(w*t+pi/ 2-h)-Usm/ z*cos(pi/2-h)*exp(-t/(R*C));
Uc1=Usm/z*cos(pi/2-h)*exp(-t/ (R*C));
Uc2=Usm/z*cos(w*t+pi/2-h);
Ur=1/(R*C )* Usm/z*cos(pi/2-h)*exp(-t/(R*C))-Usm*sin(h)*sin(w*t+pi/2-h);
Ur1=1/(R*C)*Usm/z*cos(pi/2-h)*exp(-t/(R*C));
Ur2=Usm*sin( h)*sin(w*t+pi/2-h);
I=Ur/R;I1=Ur1/R;I2=Ur2/R;
Figura ( 1)
subplot(2,1,1);plot(t,Uc,t,Uc1,t,Uc2)
title('Seno cero de Uc(t) respuesta de estado = componente transitorio + componente de estado estable')
subplot(2,1,2);plot(t,IL,t,IL1,t,IL2)
título ( 'Respuesta sinusoidal de estado cero de IL(t) = componente transitorio + componente de estado estacionario')
11., Respuesta de impulso en serie RC
R=2;C=0,5 ; % Introduzca los parámetros dados
t=[0:0.01:4];
Uc=1/(R*C)*exp(-t/(R*C)
I=-1/(R^2*C)*exp(-t/(R*C));
figura(1)
subplot(2,1,1);plot(t,Uc)
title('Forma de onda de Uc(t)')
subplot(2,1,2) ;plot ( t,I)
title('Forma de onda de I(t)')
12, respuesta de impulso paralelo RL
R=2;L =0.5 ; %Ingrese los parámetros dados
t=[0:0.01:4];
IL=R/L*exp(-t*R/L);
U=-R^2/L *exp(-t*R/L);
figura(1)
subtrama(2,1,1) ; plot (t,IL)
title('Forma de onda de IL(t)')
subplot(2,1,2);plot(t,U)
title(' Forma de onda de U(t)')