Compte rendu d'examen:
Réalisé par:QAIDI AYMANE (35)
Encadré par:[Link]
exercise1:
PID:
PI:
EX2:
function [model_morning, model_noon, model_evening] =
identify_sun_model(time, temp, humidity)
fprintf('Entrez l\'heure de début de la matinée (en heures) : ');
morning_start = input('');
fprintf('Entrez l\'heure de fin de la matinée (en heures) : ');
morning_end = input('');
fprintf('Entrez l\'heure de début de la mi-journée (en heures) : ');
noon_start = input('');
fprintf('Entrez l\'heure de fin de la mi-journée (en heures) : ');
noon_end = input('');
fprintf('Entrez l\'heure de début de la soirée (en heures) : ');
evening_start = input('');
fprintf('Entrez l\'heure de fin de la soirée (en heures) : ');
evening_end = input('');
morning_data = select_data(time, temp, humidity, morning_start,
morning_end);
noon_data = select_data(time, temp, humidity, noon_start,
noon_end);
evening_data = select_data(time, temp, humidity, evening_start,
evening_end);
model_morning = armax_model(morning_data);
model_noon = armax_model(noon_data);
model_evening = armax_model(evening_data);
disp('Modèle de la matinée :');
disp(model_morning);
disp('Modèle de la mi-journée :');
disp(model_noon);
disp('Modèle de la soirée :');
disp(model_evening);
end
function selected_data = select_data(time, temp, humidity, start_time,
end_time)
selected_data = [temp(time >= start_time & time <= end_time), ...
humidity(time >= start_time & time <= end_time)];
end
function model = armax_model(data)
p = 1; % ordre AR
q = 1; % ordre MA
n = length(data);
y = data(:, 1); % Température (ou autre sortie)
u = data(:, 2); % Humidité (ou autre entrée)
Phi = [y(1:n-1), u(1:n-1)];
theta = inv(Phi' * Phi) * Phi' * y(2:n); % Estimation des paramètres
model = theta;
end
