% Système électromécanique élémentaire. ELE4801. Jesus Morales. Hiver 2018. close all; clear; clc %__________________________________________________________ % Definition des paramètres % circuit électrique et mécanique global r xo k K D M PHI dt r = 10; xo = 3e-3; k = 6.293e-5; K = 2667; D = 4; M = 0.055; dt = 2e-5; Tmax = 0.2; t = 0:dt:Tmax; % temps de simulation % problème non linéaire tolerance = 1e-10; % tolerance d'erreur max_iter = 20; % nombre maximum d'itérations %____________________________________________________________ % Initialisation du modèle d'espace d'état %______________________________________________________________ % Intégration numérique for n = 1:N-1 Xk = ; while error > tolerance iter = iter + 1; if iter == max_iter; disp('NON CONVERGENCE') break; end Xk_new = Xk - J\F; % Newton Raphson error = ; Xk = Xk_new; end end %_____________________________________________________________