6.10 Applications

A. Introduction

Dans cette section, nous revenons sur l’application des équations différentielles dans la modélisation des systèmes d’ingénierie. En particulier, nous nous intéressons aux vibrations mécaniques et aux circuits électriques, deux grands domaines donnant lieu à l’application de systèmes d’équations différentielles.

Les équations différentielles sont d’une grande utilité dans divers domaines d’ingénierie. En génie chimique, elles sont essentielles pour modéliser la cinétique des réactions et la dynamique des processus, notamment dans des scénarios tels que les problèmes de mélanges impliquant plusieurs réservoirs et substances, des facteurs déterminants pour la conception des réacteurs et l’optimisation des processus. En génie civil, les modèles d’équations différentielles sont indispensables pour évaluer la sécurité et la longévité des structures soumises à diverses conditions de charge, comme dans l’analyse de la résistance sismique des bâtiments à plusieurs étages. L’ingénierie aérospatiale se sert de ces équations pour simuler le mouvement des avions et des engins spatiaux, en y intégrant à la fois la dynamique de translation et la dynamique de rotation. Ces connaissances sont primordiales pour élaborer des systèmes de contrôle qui améliorent la stabilité et la manœuvrabilité. L’ingénierie environnementale utilise également des modèles d’équations différentielles pour suivre la propagation des polluants, ce qui permet de concevoir des mesures efficaces de protection de l’environnement.

B. Circuits électriques

Les lois de Kirchhoff, que nous avons vues à la section 2.5, servent de base à la dérivation des équations directrices. Ces lois facilitent l’analyse des circuits en fournissant une approche systématique du calcul des courants et des tensions en divers points du circuit. Dans les circuits plus complexes, par exemple les circuits série-parallèle,

 

Exemple 6.10.1 : Circuit série RL – Système d’équations linéaires

a) Pour le schéma de circuit électrique donné, déterminer le système d’équations différentielles qui décrit les courants dans les différentes branches du circuit. On suppose que tous les courants initiaux sont nuls. b) Une fois que le système d’équations différentielles et les conditions initiales sont établis, résoudre le système pour les courants dans chaque branche du circuit.

Schéma de circuit avec une batterie de 12 V reliée à un réseau série-parallèle.

Description du schéma

 

Prenons un circuit alimenté par une batterie de 12 V CC. Une résistance de 4 Ω est branchée en série sur la borne positive de l’alimentation. Après cette résistance, le circuit se divise en deux branches parallèles. La première branche parallèle comporte une résistance de 2 Ω et la seconde un inducteur de 0,1 H. Ces deux branches convergent ensuite et le circuit continue à travers un inducteur de 0,2 H avant de revenir à la borne négative de l’alimentation. Compte tenu de cette configuration, calculer les courants I1 (aux bornes de la résistance de 4 Ω), I2 (aux bornes de la résistance de 2 Ω) et I3 (aux bornes de l’inducteur de 0,1 H). On suppose que les inducteurs sont en régime permanent.
Afficher/Masquer la solution

 

a) Dans l’exemple 4.8.2, nous avons examiné ce circuit RL et l’avons analysé avec la transformée de Laplace. Dans cet exemple, nous résolvons le même circuit avec la méthode matricielle. Les équations du système pour le circuit sont les suivantes, avec pour conditions initiales la nullité de tous les courants à l’instant t=0.

 {(4I_1+0.1I'_3+0.2I'_1=12),(0.1I'_3-2I_2=0), (I_1-I_2-I_3=0) :} ;       I_1(0)=I_2(0)=I_3(0)=0  (6.10.1)

b) Étapes de résolution du système :

1. Le système de la section 6.10.1 est un mélange d’équations différentielles et d’équations algébriques. Il faut d’abord convertir ce système en un système d’équations différentielles linéaires en utilisant la deuxième équation pour exprimer 0.1I'_3 en tant que 2I_2 dans la première équation et en isolant les premières dérivées dans les deux premières équations. Ce qui donne :

 {(I'_1=-20I_1-10I_2+60),(I'_3=20I_2), (I_1-I_2-I_3=0) :} ;       I_1(0)=I_2(0)=I_3(0)=0 (6.10.1)

Pour créer un système d’équations différentielles linéaires à partir du système donné, il faut tenir compte du fait que I_2 n’a pas de dérivée présente. Pour contourner ce problème, il faut éliminer I_2 des équations. Pour ce faire, on réarrange la troisième équation de façon à exprimer I_2en termes deI_1 et I_3, puis on remplace cette expression par I_2 dans les autres équations.

 I_2=I_1-I_3 (6.10.2)

On simplifie ensuite le système en

 {(I'_1=-30I_1+10I_3+60),(I'_3=20I_1-20I_3) :} ;       I_1(0)=I_3(0)=0 (6.10.2)

2. On exprime ensuite le problème de valeur initiale (PVI) sous forme matricielle.

 bb(I)'=[(-30,10),(20,-20)] bb(I)+[(60),(0)],    bb(I)(0)=[(0),(0)]

3. Ensuite, on trouve une solution fondamentale au système complémentaire associé (homogène). Le polynôme caractéristique de A est donné par

 c_A(lambda)=det(A-lambda I)

 =|(-30-lambda,10),(20,-20-lambda) |

 =(-30-lambda)(-20-lambda)-200

 =lambda^2+50lambda+400

 =(lambda+10)(lambda+40)

Les valeurs propres et les vecteurs propres correspondants sont

 lambda_1=-10 :  bb(u)_1=[(1),(2)] 

 lambda_2=-40 : bb(u)_1=[(-1),(1)] 

Par conséquent, {e^(-10t) [(1),(2)], \ e^(-40t) [(-1),(1)] } est un ensemble fondamental de solutions au système complémentaire.

Ainsi, la matrice fondamentale I_c(t) pour le système complémentaire est

 I_c(t)=[(e^(-10t),-e^(-40t)),(2e^(-10t),e^(-40t))]

3. Il faut maintenant déterminer une solution particulière du système.

(i) Déterminer I_c^-1(t)

 I_c^-1(t)=1/3[(e^(10t),e^(10t)),(-2e^(40t),e^(40t))]

(ii) Déterminerbb"v"(t)en définissant une constante d’intégration nulle

 bb"v"(t)=int\ I_c^-1(t)bb"f"(t) dt

 =int \ [(e^(10t),e^(10t)),(-2e^(40t),e^(40t))] [(60),(0)] dt

 =[(2e^(10t)),(-e^(40t))]

Une solution particulière du système est donc donnée par

 bb"I"_p(t)=I_c(t).bb"v"(t)

 =[(e^(-10t),-e^(-40t)),(2e^(-10t),e^(-40t))] [(2e^(10t)),(-e^(40t))]

 =[(3),(3)]

 

4. Une solution générale du système est donc donnée par

 bb"I"=bb(I)_p+bb(I)_c

 =[(3),(3)] +c_1e^(-10t) [(1),(2)] +c_2e^(-40t) [(-1),(1)]

 

5. On applique maintenant les conditions initiales pour trouver les constantes dans la solution générale.

 bb(I)(0)=[(0),(0)]

 [(3),(3)] +c_1e^(-10t) [(1),(2)] +c_2e^(-40t) [(-1),(1)]=[(0),(0)]

On obtient ainsi un système de deux équations à deux inconnues.

 {(c_1-c_2=-3),(2c_1+c_2=-3 ) :}

En résolvant le système, on obtient

 c_1=-2,    c_2=1

La solution du problème de valeur initiale est donc

 bb(I)=[(3),(3)] -2e^(-10t) [(1),(2)] +e^(-40t) [(-1),(1)]

Cela donne, dans les expressions finales pour I_1 et I_3 :

 I_1=3-2e^(-10t)-e^(-40t)

 I_3=3-4e^(-10t)+e^(-40t)

6. Pour trouver I_2, on remplace à nouveau l’expression par I_1 et I_3dans l’équation (6.10.1), ce qui donne

I_2=I_1-I_3 

 I_2=2e^(-10t)-2e^(-40t) 

 

C. Vibrations mécaniques

L’analyse des vibrations mécaniques est cruciale pour concevoir des systèmes résistants aux charges dynamiques. Un modèle plus réaliste qui capture l’essence des systèmes mécaniques implique de considérer non seulement les masses et les ressorts, mais aussi les éléments d’amortissement et les forces externes. Cette section se concentre sur un système composé de deux masses reliées par des ressorts dans un arrangement horizontal, avec amortissement et des forces externes agissant sur les deux masses. Un tel modèle peut représenter un large éventail d’applications techniques, des suspensions de véhicules aux composants de machines. Un schéma du système est présenté à la figure 6.10.1.

Diagramme d’un système de vibration mécanique avec deux masses reliées par trois ressorts dans un arrangement horizontal, avec deux amortisseurs. La masse m1 est reliée à une paroi fixe par le ressort k1 et à la masse m2 par le ressort k2, tandis que la masse m2 est reliée à une autre paroi fixe par le ressort k3. Les amortisseurs associés à m1 et m2 sont respectivement représentés par les coefficients c1 et c2. Les forces externes F1(t) et F2(t) agissent respectivement sur m1 et m2, indiquant l’interaction dynamique du système avec son environnement.

Figure 6.10.1 Schéma d’un système couplé masse-ressort comprenant deux masses reliées par trois ressorts et deux amortisseurs.

 Hypothèses

Pour procéder à la dérivation des équations du mouvement du système, on fait les hypothèses suivantes :

  • Amortissement linéaire : chaque masse est couplée à un élément d’amortissement, caractérisé par les coefficients d’amortissement linéaire c_1pour m_1  et c_2 pour m_2. Ces coefficients quantifient la résistance au mouvement de chaque masse.
  • Forces externes : des forces externes dépendantes du temps F_1(t) et F_2(t) agissent respectivement sur m_1 et m_2 et sont considérées comme positives dans la direction droite.
  • Élasticité linéaire : les ressorts obéissent à la loi de Hooke, c’est-à-dire que la force exercée sur chaque ressort est directement proportionnelle à son déplacement par rapport à la longueur de repos.
  • Petits déplacements : l’analyse suppose de petits déplacements à partir du point d’équilibre, ce qui autorise la linéarisation du système. Les déplacements sont réputés être positifs lorsqu’ils sont orientés vers la droite.
  • Corps rigide : les masses sont traitées comme des points matériels, tandis que les ressorts et les amortissements sont censés être de masse nulle, se concentrant uniquement sur les forces et déplacements axiaux.

 Configuration du système

Nous considérons un cas général avec un système constitué de deux masses, m_1 ​ et m_2, reliées par trois ressorts aux constantes de rigidité k_1, k_2 et k_3 et augmentées par deux amortisseurs. Les ressorts extérieurs sont attachés à des parois fixes. Les forces externes agissent sur les masses et les amortisseurs contrecarrent leurs mouvements. Cette structure permet d’ajuster les forces externes et les effets d’amortissement, de façon à s’adapter aux scénarios dans lesquels ces forces sont absentes en fixant leurs valeurs respectives à zéro.

 

Exemple 6.10.2 : Vibrations mécaniques – Système amorti et forcé

Dériver le système d’équations différentielles pour le système couplé amorti et forcé décrit ci-dessus (figure 6.10.1).

Afficher/Masquer la solution

 

La dynamique de ce système amorti avec des forces externes est régie par deux équations différentielles linéaires couplées du second ordre, reflétant l’équilibre des forces sur chaque masse. Ici, on considère que la direction des forces externes est vers la droite et que les déplacements sont également positifs (sur la droite), en supposant que le déplacement de la masse 1, x_1, est supérieur au déplacement de la masse 2, x_2, de sorte que x_2-x_1<0/span>.

1) Les forces agissant sur la masse m_1sont

  • Force de rappel du ressort k_1F_(s1)=-k_1x_1.
  • Force de rappel du ressort k_2F_(s2)=-k_2(x_2-x_1), où x_2-x_1 est le déplacement du ressort du milieu.
  • Force d’amortissement F_(d1)=-c_1x_1', où c_1 est le coefficient d’amortissement de l’amortisseur 1. Si elle est présente, la force d’amortissement est proportionnelle à la vitesse de la masse, mais dans la direction opposée au mouvement.
  • Force externe F_1(t). Il s’agit de toute force externe agissant sur la masse m_1, qui peut être périodique ou aléatoire, et induisant des vibrations forcées.

Diagramme montrant les forces agissant sur la masse m1

D’après la deuxième loi du mouvement de Newton,

 ma=sommeF

 m_1x''_1=-c_1x'_1-k_1x_1+k_2(x_2-x_1)+F_1(t)

Cette équation peut être simplifiée en

 m_1x''_1=-c_1x'_1-(k_1+k_2)x_1+k_2x_2+F_1(t)

On remarque que, comme x_2ltx_1, le ressort k_2 est comprimé et la force qu’il exerce sur la masse 1 est orientée vers la gauche (négative), le but étant de rétablir le ressort à sa longueur d’équilibre.

2) Les forces agissant sur la masse m_2 sont

  • Force de rappel du ressort k_3F_(s1)=-k_3x_2.
  • Force de rappel du ressort k_2F_(s2)=k_2(x_2-x_1), où x_2-x_1 est le déplacement du ressort du milieu.
  • Force d’amortissement F_(d1)=-c_2x_2', où c_2 est le coefficient d’amortissement de l’amortisseur 2. Si elle est présente, la force d’amortissement est proportionnelle à la vitesse de la masse, mais dans la direction opposée au mouvement.
  • Force externe F_2(t). Il s’agit de toute force externe agissant sur la masse m_2, qui peut être périodique ou aléatoire, et induisant des vibrations forcées.

Diagramme montrant les forces agissant sur la masse m2

D’après la deuxième loi du mouvement de Newton,

 ma=sommeF

 m_2x''_2=-c_2x'_2-k_3x_2-k_2(x_2-x_1)+F_2(t)

Cette équation peut être simplifiée en

 m_2x''_2=-c_2x'_2+k_2x_1 -(k_2+k_3)x_2+F_2(t)

On remarque que, comme x_2ltx_1, le ressortk_2 est comprimé et la force qu’il exerce sur la masse 2 est orientée vers la droite et devrait être positive, ce qui est cohérent avec le signe deimage.

Par conséquent, les déplacements dépendants du temps des masses sont décrits par le système d’équations différentielles

 m_1x''_1=-c_1x'_1-(k_1+k_2)x_1+k_2x_2+F_1(t)

 m_2x''_2=-c_2x'_2+k_2x_1 -(k_2+k_3)x_2+F_2(t)

 

Exemple 6.10.3 : Vibrations mécaniques – Réécrire un système en équations du premier ordre

a) Réécrire le système dérivé d’équations différentielles de l’exemple 6.10.2 en un système d’équations du premier ordre. b) Écrire le système sous forme matricielle.

Afficher/Masquer la solution

 

a) Les équations régissant la système masse-ressort de la figure 6.10.1 sont dérivées dans l’exemple 6.10.2.

 m_1x''_1=-c_1x'_1-(k_1+k_2)x_1+k_2x_2+F_1(t)

 m_2x''_2=-c_2x'_2+k_2x_1 -(k_2+k_3)x_2+F_2(t)

La section 6.4-C a expliqué comment convertir des équations différentielles d’ordre supérieur en un système d’équations du premier ordre. On introduit de nouvelles variables comme suit :

y_1=x_1y_2=x'_1 , y_3=x_2 , et y_4=x'_2

Les équations peuvent alors être écrites comme suit :

y'_1=y_2

 m_1y'_2=-c_1y_2-(k_1+k_2)y_1+k_2y_3+F_1(t)

y'_3=y_4

 m_2y'_4=-c_2y_4+k_2y_1 -(k_2+k_3)y_3+F_2(t)

b) Le système sous forme matricielle est

 bb(y)'(t)=Abb(y)(t)+bb(f)(t)

 bb(y)'(t)=[(0,1,0,0),(-(k_1+k_2)/(m_1),(c_1)/(m_1),(k_2)/(m_1),0),(0,0,0,1),((k_2)/(m_2),0,-(k_2+k_3)/(m_2) ,-(c_2)/(m_2))] bb(y)(t)+[(0),((F_1(t))/(m_1)),(0),((F_2(t))/(m_2)) ]

 

Exemple 6.10.4 : Résoudre un problème de valeur initiale : vibration libre non amortie

Considérons un système masse-ressort couplé, tel celui de l’exemple 6.10.2, avec les paramètres suivants : les deux masses m_1 ​et m_2 pèsent 1 kg et toutes les constantes du ressortk_1, k_2 et k_3 ​sont de 1 N/m. Le système est isolé contre les forces externes et les effets d’amortissement. Au départ, les déplacements et les vitesses sont x_1(0)=0 m, x'_1(0)=0 m/s, x_2(0)=1 m et x'_2(0)=0 m/s. Résoudre le problème de valeur initiale pour déterminer les déplacements des masses en fonction du temps. Étant donné la complexité des calculs, utiliser un logiciel d’algèbre matricielle pour trouver les valeurs propres et les vecteurs propres.

Afficher/Masquer la solution

 

Informations données :

  •  m_1=m_2=1\ "kg"
  • k_1=k_2=k_3=1\ "N"//"m"
  • Pas d’amortissement : c_1=c_2=0
  • Pas de forces externes : F_1=F_2=0
  • Conditions initiales :x_1(0)=0, x_2(0)=1\ "m"x'_1(0)=x_2'(0)=0

Dans l’exemple 6.10.3, on a converti un système du second ordre régissant le système masse-ressort couplé en un système du premier ordre et on l’a exprimé en notation matricielle.

 bb(y)'(t)=Abb(y)(t)+bb(f)(t)

 bb(y)'(t)=[(0,1,0,0),(-(k_1+k_2)/(m_1),(c_1)/(m_1),(k_2)/(m_1),0),(0,0,0,1),((k_2)/(m_2),0,-(k_2+k_3)/(m_2) ,-(c_2)/(m_2))] bb(y)(t)+[(0),((F_1(t))/(m_1)),(0),((F_2(t))/(m_2)) ]

 bb(y)=[(y_1),(y_2),(y_3),(y_4)]=[(x_1),(x'_1),(x_2),(x'_2)]

En remplaçant les valeurs données, le problème de valeur initial devient

 bb(y)'(t)=[(0,1,0,0),(-2,0,1,0),(0,0,0,1),(1,0,-2 ,0)] bb(y)(t), \ \ bb(y)(0)=[(0),(0),(1),(0) ]

En utilisant un logiciel d’algèbre matricielle, on trouve les valeurs propres et les vecteurs propres de la matrice coefficient. Les valeurs propres et les vecteurs propres se présentent sous la forme des paires complexes conjuguées suivantes :

Vecteur propre pour lambda_(1,2)=+-ibb(u)_(1,2)=[(0),(1),(0),(1)]+-i [(-1),(0),(-1),(0)]

Vecteur propre pour lambda_(3,4)=+-sqrt(3)ibb(u)_(3,4)=[(0),(-3),(0),(3)]+-i [(sqrt(3)),(0),(-sqrt(3)),(0)]

La solution pour chaque paire est donnée par l’équation présentée à la section 6.7.1.

Pour la première paire, lambda_1,2=+-i, les deux solutions linéairement indépendantes sont

 bb"y"_1=[(0),(1),(0),(1)] cos(t)-[(-1),(0),(-1),(0)] sin( t)= [(sin(t)),(cos(t)),(sin(t)),(cos(t))]

 bb"y"_2=[(0),(1),(0),(1)] sin(t)+[(-1),(0),(-1),(0)]cos( t)= [(-cos(t)),(sin(t)),(-cos(t)),(sin(t))]

Pour la seconde paire, lambda_3,4=+-sqrt(3)i, les deux solutions linéairement indépendantes sont

 bb"y"_3=[(0),(-3),(0),(3)] cos(sqrt(3) t)-[(sqrt(3)),(0),(-sqrt(3)),(0)] sin( sqrt(3) t)= [(-sqrt(3) sin(sqrt(3) t)),(-3cos(sqrt(3) t)),(sqrt(3) sin(sqrt(3) t)),(3cos(sqrt(3) t))]

 bb"y"_4=[(0),(-3),(0),(3)] sin(sqrt(3) t)+[(sqrt(3)),(0),(-sqrt(3)),(0)] cos( sqrt(3) t)= [(sqrt(3) cos(sqrt(3) t)),(-3sin(sqrt(3) t)),(-sqrt(3) cos(sqrt(3) t)),(3sin(sqrt(3) t))]

La matrice fondamentale du système est donc

 Y(t)=  [(bb(y)_1,bb(y)_2, bb(y)_3, bb(y)_4)]

 = [(sin(t),-cos(t),-sqrt(3) sin(sqrt(3) t),sqrt(3) cos(sqrt(3) t)),(cos(t),sin(t),-3 cos(sqrt(3) t),-3sin(sqrt(3) t)),(sin(t),-cos(t),sqrt(3) sin(sqrt(3) t),-sqrt(3) cos(sqrt(3) t)),(cos(t),sin(t),3 cos(sqrt(3) t),3 sin(sqrt(3) t))]

La solution générale du système sous forme matricielle est

 bb(y)(t)= [(sin(t),-cos(t),-sqrt(3) sin(sqrt(3) t),sqrt(3) cos(sqrt(3) t)),(cos(t),sin(t),-3 cos(sqrt(3) t),-3sin(sqrt(3) t)),(sin(t),-cos(t),sqrt(3) sin(sqrt(3) t),-sqrt(3) cos(sqrt(3) t)),(cos(t),sin(t),3 cos(sqrt(3) t),3 sin(sqrt(3) t))] [(c_1),(c_2),(c_3),(c_4)]

En appliquant les conditions initiales, on obtient

 bb(y)(0)=[(0),(0),(1),(0) ]

[(0,-1,0,sqrt(3)),(1,0,-3,0),(0,-1,0,-sqrt(3)),(1,0,3 ,0)] [(c_1),(c_2),(c_3),(c_4)]  =[(0),(0),(1),(0) ]

En résolvant les coefficients, on trouve

 c_1=0,\ c_2=-1/2,\ c_3=0,\ c_4=-sqrt3/6

La solution du problème de valeur initiale est donc

 bb(y)(t)= [(sin(t),-cos(t),-sqrt(3) sin(sqrt(3) t),sqrt(3) cos(sqrt(3) t)),(cos(t),sin(t),-3 cos(sqrt(3) t),-3sin(sqrt(3) t)),(sin(t),-cos(t),sqrt(3) sin(sqrt(3) t),-sqrt(3) cos(sqrt(3) t)),(cos(t),sin(t),3 cos(sqrt(3) t),3 sin(sqrt(3) t))] [(0),(-1/2),(0),(-sqrt(3)/6)]

Elle peut s’écrire comme suit

 bb(y)=[(y_1),(y_2),(y_3),(y_4)]  = [(1/2cos(t)-1/2cos(sqrt(3)t)),(-1/2sin(t)+(sqrt(3))/2sin(sqrt(3)t)),(1/2cos(t)+1/2cos(sqrt(3)t)),(-1/2sin(t)-(sqrt(3))/2sin(sqrt(3)t))]

Rappelons que, dans l’exemple 10.6.3, on a introduit des variables bb(y)_n pour représenter les déplacements et les vitesses des masses.

 [(y_1),(y_2),(y_3),(y_4)]=[(x_1),(x'_1),(x_2),(x'_2)]

Étant donné cette conversion, les déplacements de la masse 1 (x_1) et de la masse 2 (x_2) en fonction du temps sont déterminés par

 x_1=y_1=1/2cos(t)-1/2cos(sqrt(3)t)

 x_2=y_3=1/2cos(t)+1/2cos(sqrt(3)t)

La représentation graphique ci-dessous montre les déplacements des deux masses au fil du temps dans un système masse-ressort couplé dans cet exemple, sur un graphique avec le temps sur l’axe horizontal (allant de 0 à 10 secondes) et le déplacement en mètres sur l’axe vertical. La ligne pour la masse 1 oscille, soit un modèle de mouvement qui varie entre des déplacements positifs et des déplacements négatifs, suggérant un mouvement harmonique complexe. La ligne de la masse 2 suit un schéma oscillatoire similaire, mais avec des différences de phase et d’amplitude par rapport à la masse 1, reflétant l’interaction entre les deux masses au moyen du système de ressorts.

 

Section 6.10 Exercices

  1. Considère un système masse-ressort couplé, tel celui de l’exemple 6.10.2, avec les paramètres suivants : les masses m_1 ​et m_2 pèsent 1 kg et toutes les constantes du ressortk_1, k_2 et k_3 ​sont de 1 N/m. Le système est isolé contre les forces externes et les effets d’amortissement. Au départ, les déplacements et les vitesses sont x_1(0)=0 m, x'_1(0)=2 m/s, x_2(0)=0 m et x'_2(0)=0 m/s. Résoudre le problème de valeur initiale pour déterminer les déplacements des masses en fonction du temps. Utilise un logiciel d’algèbre matricielle pour trouver les valeurs propres et les vecteurs propres.
    Afficher/Masquer la réponse

     x_1=sin(t)+sqrt(3)/3sin(sqrt(3)t)

     x_2=sin(t)-sqrt(3)/3sin(sqrt(3)t)

Licence

Équations différentielles© par AMIR TAVANGAR. Tous droits réservés.

Partagez ce livre