How can solve differential equation in Matrix form( dI/dt=c.*(M*I)+c.*(L*V) )
1 view (last 30 days)
Show older comments
vas=1;vbs=2;vcs=3;var=1.5;vbr=2.5;vcr=3.5; Rs=1.06;Rr=0.8;Lm=0.0664;Ls=0.2065;Lr=0.0810;p=2;Jg=23; S=1-(Lm^2)/(Ls*Lr);Os=1;Or=0;Om=Os-Or;L=Lm/Ls;tm=pi/2; vabcs=[vas;vbs;vcs]; Tm=[cos(tm) -sin(tm); sin(tm) cos(tm)]; Tc=2/3.*[1 -1/2 -1/2;0 sqrt(3)/2 -sqrt(3)/2 ]; vabs=Tc*vabcs;vur=u*vabcr;vDQr=Tc*vur;vabr=Tm*vDQr; syms ias(t) ibs(t) iar(t) ibr(t) I=[ias;ibs;iar;ibr]; c=1/(S*Ls*Lr); c11=-Rs*Lr;c12=Om*Lm^2;c13=Rr*Lm;c14=Om*Lm*Lr; c21=-c12;c22=c11;c23=-c14;c24=c13; c31=Rs*Lm;c32=-Om*Ls*Lm;c33=-Rr*Ls;c34=-Om*Lr*Ls; c41=-c32;c42=c31;c43=-c34;c44=c33; M=[c11 c12 c13 c14 ; c21 c22 c23 c24; c31 c32 c33 c34; c41 c42 c43 c44]; L=[Lr 0 -Lm 0; 0 Lr 0 -Lm ; -Lm 0 Ls 0; 0 -Lm 0 Ls ]; vab=[vabs;vabr]; odes=diff(I)==c.*(M*I)+ c.*(L*vab) [iasSol(t),ibsSol(t),iarSol(t),ibrSol(t)] = dsolve(odes); iasSol(t)= simplify(iasSol(t)) ibsSol(t)= simplify(ibsSol(t)) iarSol(t)= simplify(iarSol(t)) ibrSol(t)= simplify(ibrSol(t)) I can not find any results. And thank you in advance. email: driss.belkhiri@edu.uiz.ac.ma
0 Comments
Answers (0)
See Also
Categories
Find more on Error Detection and Correction in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!