maple1




exercice 1 double masses
> with(inttrans):eqn1 := alpha*m*diff(x[1](t),t$2) = k*(x[2](t) - x[1](t)) + u(t);
> eqn2 := m*diff(x[2](t),t$2) = k*(x[1](t) - x[2](t));
> e1l := laplace(eqn1,t,s);
> Ics := {x[1](0) = 0, D(x[1])(0) = 0, x[2](0) = 1, D(x[2])(0) = 0};
> E1l := subs(Ics,e1l);
> e2l := laplace(eqn2,t,s);
> E2l := subs(Ics,e2l);
> ans := solve({E1l,E2l},map(laplace,{x[1](t),x[2](t)},t,s));
> Ans := subs(m=1,k=1,u(t) =Dirac(t),ans) ;
> impulse_response :=invlaplace(Ans,s,t);
> plot(subs(impulse_response,alpha=10,{x[1](t),x[2](t)}),t=0..10);
> map(x -> lhs(x) = int(rhs(x),t),impulse_response);
> step_response:=%:
> plot(subs(step_response,alpha=10,{x[1](t),x[2](t)}),t=0..10);
> e1 := diff(x1(t),t) = x2(t);
> e2 := diff(x2(t),t) = -k/(m*alpha)*x1(t) + k/(m*alpha)*x3(t) + u(t)/(m*alpha) ;
> e3 := diff(x3(t),t) = x4(t); e4 := diff(x4(t),t) = k/m*x1(t) - k/m*x3(t);
> eqn_list:= {seq(e||n,n=1..4)}; var_list := {seq(evaln(x||n)(t),n=1..4)};
> ics :={x1(0)=0,x2(0)=0,x4(0)=0,x3(0)=1} ; # initial length of sprign = 1
> dsolve(subs(u(t) =1, eqn_list union ics), var_list); # try a step - no good
> map(laplace,eqn_list ,t,s); # try laplace again
> subs(ics, %);
> solve(%%,map(laplace,var_list,t,s));
> subs(u(t)=Dirac(t),%);
> I_response := eval(%);
> invlaplace(%,s,t);
> ans:= subs(ics,%);
> plot(subs(ans,m=1,k=1,alpha=10,{x3(t),x1(t)}),t=0..10);

Commentaires (3)

1. arjaz 18/02/2006

MAD GANT TIRADGH ITYARAN DA

2. amdenar 28/05/2006

salut
je vois que tu utilise maple
tu ne programmes pas en java ?

ps: je suis de douar amdenar

3. hassani 03/12/2006

A mat3nam à ait tmazirt n oguensou ou awd win berra.
je m'excuse je suis nouveau si vouv voulez bien nous expliquer pour comprendre ce charabia ( excusez le terme )
Je suis de taltnamart, mais je bosse au nord du Maroc et cela fait plus de 4 ans que je n'ai pas vu le bled. TAGHI TGHOUFI N TMAZIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIRT AWI. AKN A3AWN RABBI .

Ajouter un commentaire

Vous utilisez un logiciel de type AdBlock, qui bloque le service de captchas publicitaires utilisé sur ce site. Pour pouvoir envoyer votre message, désactivez Adblock.

MOBILISONS NOUS POUR LE DEVELOPPEMENT DE SARGHINE EN AIDANT L'ASSOCIATION AMSIRAR.

Créer un site gratuit avec e-monsite - Signaler un contenu illicite sur ce site