Corsi di Laurea

[Uchiha]

Ragazzi sto svolgendo la prova del 29 maggio 2014 solo che nell'andare a risolvere il sistema per t>0 mi compare questo errore " Warning: Explicit solution could not be found." Sto provando a riscrivere le equazioni in tutti i modi possibili ma è sempre la stessa storia.

 

Help Please 

 

ecco il codice:

 
clc
close all
clear all

%Dati

Rx=1;

C=10e-6; L=5e-6;

G11=0.15;G22=0.15;G12=-0.5;G21=-0.5;

G=[G11,G12;G21,G22];

ampvs=2; ampis=10;
w=50;
phi=pi/4;
syms t

VS=ampvs*cos(w*t+phi);
IS=ampis*cos(2*w*t+phi);

%quesito 1) Vr2_t e ir3_t per tE[-inf,+inf]

%analisi del circuito per t<0: interruttore aperto e circuito a regime
%la presenza dell'interruttore nello stato aperto implica che R4=0, poichè
%per esso non passa corrente ( ir4=0 )

syms ka vL iL ic vc vs is i1 i2 v1 v2

eq1=Rx*ka+vL;
eq2=-vL+Rx*(ka-iL)-vs;
eq3=vs+vc+Rx*(ic-i1);
eq4=-i1+G11*(Rx*(ic-i1))+G12*v2;
eq5=-is+G21*(Rx*(ic-i1))+G22*v2;

sol_m=solve(eq1,eq2,eq3,eq4,eq5,'ka,vL,ic,i1,v2');

ic=sol_m.ic;
vL=sol_m.vL;

dvcdt_m=ic/C;
dildt_m=vL/L;

A11_m=subs(dvcdt_m,[vc,iL,vs,is],[1,0,0,0]);
A12_m=subs(dvcdt_m,[vc,iL,vs,is],[0,1,0,0]);
A21_m=subs(dildt_m,[vc,iL,vs,is],[1,0,0,0]);
A22_m=subs(dildt_m,[vc,iL,vs,is],[0,1,0,0]);

A_m=[A11_m,A12_m;A21_m,A22_m];

B11_m=subs(dvcdt_m,[vc,iL,vs,is],[0,0,1,0]);
B12_m=subs(dvcdt_m,[vc,iL,vs,is],[0,0,0,1]);
B21_m=subs(dildt_m,[vc,iL,vs,is],[0,0,1,0]);
B22_m=subs(dildt_m,[vc,iL,vs,is],[0,0,0,1]);

B_m=[B11_m,B12_m;B21_m,B22_m];

%sovrapposizione degli effetti per il calcolo delle variabili di stato

%calcolo delle variabili di stato per Vs

Fvs=ampvs*exp(j*phi);
Fu_vs=[Fvs;0];
Id=eye(2);
Fxvs=inv(j*w*Id-A_m)*B_m*Fu_vs;
vcvs_m=abs(Fxvs(1))*cos(w*t+angle(Fxvs(1)));
iLvs_m=abs(Fxvs(2))*cos(w*t+angle(Fxvs(2)));

%calcolo delle variabili di stato per Is

Fis=ampis*exp(j*phi);
Fu_is=[0;Fis];
Fxis=inv(j*2*w*Id-A_m)*B_m*Fu_is;
vcis_m=abs(Fxis(1))*cos(w*t+angle(Fxis(1)));
iLis_m=abs(Fxis(2))*cos(w*t+angle(Fxis(2)));

vc_m=(vcvs_m+vcis_m)
iL_m=(iLvs_m+iLis_m)


figure(1)
subplot(211),ezplot(vc_m,[0,10*pi/w]),axis auto,grid on,title('vc per t<0')
subplot(212),ezplot(iL_m,[0,10*pi/w]),axis auto,grid on,title('iL per t<0')

%analisi per t>0: interruttore chiuso e presenza del resistore R4 in
%parallelo con il generatore di corrente sinusoidale Is

syms  ka v1 v2 i1 i2 kc

eq1p=Rx*ka+vL;
eq2p=-vL+Rx*(ka-iL)-vs;
eq3p=vs+vc+Rx*(ic-i1);
eq4p=-i1+G11*(Rx*(ic-i1))+G12*(Rx*(i2+is));
eq5p=-i2+G21*(Rx*(ic-i1))+G22*(Rx*(i2+is));

sol_p=solve(eq1p,eq2p,eq3p,eq4p,eq5p,'ka,vL,ic,i1,i2');
 
Warning: Explicit solution could not be found. 
> In solve at 179
  In esamemaggio at 95 

[angela DL]

ragazzi qualcuno ha svoltao la prova del 26 maggio 2014, il monofase?