diff --git a/BE_Beam.m b/BE_Beam.m index 610fdba..f04f2d6 100644 --- a/BE_Beam.m +++ b/BE_Beam.m @@ -51,28 +51,30 @@ eig_tustin3 = eig(c2d(sys, Ts3, 'tustin')) ; eig_zoh1 = eig(c2d(sys, Ts1, 'zoh')) ; eig_zoh2 = eig(c2d(sys, Ts2, 'zoh')) ; eig_zoh3 = eig(c2d(sys, Ts3, 'zoh')) ; - -figure; +%% +Question4 = figure ; subplot(1,3,1) ; plot(eig_cont, 'x'); +axis square; grid on; xline(0, 'r--'); xlabel('Axe Réel'); ylabel('Axe Imaginaire'); -title('Pôles du système'); +title({'Valeurs propres en temps continu'}); subplot(1,3,2); hold on; -plot(eig_tustin1,'x', 'LineWidth',2); -plot(eig_tustin2,'x', 'LineWidth',2); -plot(eig_tustin3,'x', 'LineWidth',2); - +plot(eig_tustin1,'x', 'LineWidth',1); +plot(eig_tustin2,'x', 'LineWidth',1); +plot(eig_tustin3,'x', 'LineWidth',1); +axis square; theta = linspace(0, 2*pi, 100); -plot(cos(theta), sin(theta), 'k--', 'HandleVisibility','off'); +plot(cos(theta), sin(theta), 'k--', 'HandleVisibility','off','LineWidth',0.25); xlim([-1.1, 1.1]); ylim([-1.1, 1.1]); grid on; -title('Pôles dans le plan en Z (Discret)'); +title({'Valeurs propres en temps discret',' avec la méthode de Tustin'}); +legend("Ts1", "Ts2", "Ts3"); xlabel('Réel'); ylabel('Imaginaire'); hold off; @@ -82,20 +84,19 @@ hold on; plot(eig_zoh1,'x', 'LineWidth',1); plot(eig_zoh2,'x', 'LineWidth',1); plot(eig_zoh3,'x', 'LineWidth',1); - -legend("Eigen Values Ts1", "Eigen Values Ts2", "Eigen Values Ts3"); +axis square; +legend("Ts1", "Ts2", "Ts3"); theta = linspace(0, 2*pi, 100); -plot(cos(theta), sin(theta), 'k--', 'HandleVisibility','off'); +plot(cos(theta), sin(theta), 'k--', 'HandleVisibility','off','LineWidth',0.25); xlim([-1.1, 1.1]); ylim([-1.1, 1.1]); grid on; -title('Valeurs propres avec un bloqueur d''ordre 0 (ZOH)'); +title({'Valeurs propres en temps discret',' avec un bloqueur d''ordre 0'}); xlabel('Réel'); ylabel('Imaginaire'); hold off; -% Conlsuion : On en conclut qqc. - +exportgraphics(gcf, './latex/Illustrations/Question4.pdf') %% Question 5 @@ -103,12 +104,14 @@ t = 0:0.01:10; cmd = double (t >= 1) ; [y, t] = lsim(sys,cmd,t) ; -figure ; +Question5 = figure ; plot(t, cmd, t, y); xlabel("Temps (s)"); ylabel("Amplitude"); legend("u(t)", "y(t)"); title("Simulation en boucle ouverte") ; +exportgraphics(Question5, './latex/Illustrations/Question5.pdf') + %% Question 7 @@ -118,6 +121,7 @@ Cw = @(zeta) [zeta^2*(2*zeta^3-5*zeta^2+10)/20 ... zeta^4*(3*zeta-5)/60 ... 0, 0, 0, 0]; Cw_L = Cw(L); +%% TransientTime = [] ; for i = 0.01:0.01:10 @@ -131,12 +135,15 @@ k = min_index * 0.01 H = inv(Cw_L*inv(-A+B*C*k)*B) [y, t, x] = lsim(feedback(sys,k), H*cmd, t) ; u = H*cmd'-k*y ; -figure ; +%% +Question8 = figure ; plot(t, u, t, y, t, cmd, t, (Cw_L*x')'); xlabel("Temps (s)"); ylabel("Amplitude"); -legend("u(t)", "y(t)", "wc(t)", "w(t)"); +legend("$u(t)$", "$y(t)$", "$w_c(L,t)$", "$w(L,t)$", 'Interpreter', 'latex'); title("Simulation du retour de sortie avec pré-gain") ; +exportgraphics(Question8, './latex/Illustrations/Question8.pdf') + %% Question 9 diff --git a/latex/Illustrations/Question4.pdf b/latex/Illustrations/Question4.pdf new file mode 100644 index 0000000..2b6d0eb Binary files /dev/null and b/latex/Illustrations/Question4.pdf differ diff --git a/latex/Illustrations/Question5.pdf b/latex/Illustrations/Question5.pdf new file mode 100644 index 0000000..606a01c Binary files /dev/null and b/latex/Illustrations/Question5.pdf differ diff --git a/latex/Illustrations/Question8.pdf b/latex/Illustrations/Question8.pdf new file mode 100644 index 0000000..81a5001 Binary files /dev/null and b/latex/Illustrations/Question8.pdf differ diff --git a/latex/contents.tex b/latex/contents.tex index 1a37de7..8940606 100644 --- a/latex/contents.tex +++ b/latex/contents.tex @@ -115,10 +115,10 @@ Puis on a dans la deuxième ligne de l'équation 1 : \subsection{Exercice 4} - +\includegraphics[width=\textwidth]{Illustrations/Question4} \subsection{Exercice 5} - +\includegraphics[width=\textwidth]{Illustrations/Question5} \subsection{Exercice 6} @@ -130,6 +130,7 @@ Puis on a dans la deuxième ligne de l'équation 1 : \subsection{Exercice 8} +\includegraphics[width=\textwidth]{Illustrations/Question8} \subsection{Exercice 9} diff --git a/latex/main.pdf b/latex/main.pdf index 8498ee6..f3fa8d3 100644 Binary files a/latex/main.pdf and b/latex/main.pdf differ diff --git a/latex/main.tex b/latex/main.tex index 040d853..356fb4b 100644 --- a/latex/main.tex +++ b/latex/main.tex @@ -8,7 +8,7 @@ \newcommand{\firstcouverture}{ \parbox{\textwidth}{ \sffamily % arial - \textbf{Wissal GUARNI \\ Justin BOS \\ Nolan REYNIER-NOMER \\ Aleksander TABAN}\\ + \textbf{Wissal GUARNI \\ Justin BOS \\ Nolan REYNIER NOMER \\ Aleksander TABAN}\\ Elèves Ingénieurs de l'INSA Toulouse\\ Département GEI \\ Spécialité AE-SE \\