From 45fac2c8834d77bc71d222dfec70887bbbded389 Mon Sep 17 00:00:00 2001
From: JB <bogdangheorghe98@gmail.com>
Date: Thu, 27 Mar 2025 23:58:35 +0200
Subject: [PATCH] code prepared for course
---
sec3_predefined_utilities.m | 11 +++++------
1 file changed, 5 insertions(+), 6 deletions(-)
diff --git a/sec3_predefined_utilities.m b/sec3_predefined_utilities.m
index cbceec1..cc33107 100644
--- a/sec3_predefined_utilities.m
+++ b/sec3_predefined_utilities.m
@@ -58,18 +58,17 @@ F = integrator('F', 'cvodes', ode, opts);
figure; grid on; hold on;
xi = [2; 0];
-trajectory_vals = [xi];
+traj_vals = [xi];
for i = 1:floor(T / T_sample)
if (T_sample * i < T_star)
- trajectory = F('x0', trajectory_vals(:, end), 'p', 1);
+ trajectory = F('x0', traj_vals(:, end), 'p', -1);
else
- trajectory = F('x0', trajectory_vals(:, end), 'p', -1);
+ trajectory = F('x0', traj_vals(:, end), 'p', 1);
end
- trajectory_vals = [trajectory_vals, full(trajectory.xf)];
+ traj_vals = [traj_vals, full(trajectory.xf)];
end
%get the trajectory samples
-plot(trajectory_vals(1, :), trajectory_vals(2, :));
-
+plot(traj_vals(1, :), traj_vals(2, :));
%% If else instruction
% Saturate function
--
GitLab