From 6b589b2f39cc5996b7c14006d5be8010afc67e89 Mon Sep 17 00:00:00 2001
From: Nao Pross <np@0hm.ch>
Date: Fri, 31 May 2024 01:39:04 +0200
Subject: Update DK iteration and model structure

---
 uav_performance_hinf.m | 45 ++++++++++++++++++++++++++-------------------
 1 file changed, 26 insertions(+), 19 deletions(-)

(limited to 'uav_performance_hinf.m')

diff --git a/uav_performance_hinf.m b/uav_performance_hinf.m
index ba2fd3b..4a6de8a 100644
--- a/uav_performance_hinf.m
+++ b/uav_performance_hinf.m
@@ -19,41 +19,43 @@ function [perf] = uav_performance_hinf(params, do_plots)
 s = tf('s');
 
 % Bandwitdhs
-bw_alpha = .7 * params.actuators.ServoNominalAngularVelocity;
-bw_omega = 8;
+T_alpha = params.actuators.ServoSecondsTo60Deg;
+T_omega = 0.1;
 
-bw_xy = .1;
-bw_z = .4;
+T_xy = 5;
+T_z = 10;
 
-bw_xydot = .5;
-bw_zdot = .1;
+% Inverse performance functions
+W_Palpha = make_weight(1/T_alpha, 10, 2);
+W_Pomega = make_weight(1/T_omega, 10, 2);
 
-bw_phitheta = bw_xy;
-bw_psi = .08;
+W_Rxy = 1 / ((s * T_xy)^2 + 2 * T_xy * .8 * s + 1);
+W_Rz = 1 / (s * T_z + 1);
 
-% Inverse performance functions
-W_Palpha = .25 / (s / bw_alpha + 1); 
-W_Pomega = .1 / (s / bw_omega + 1);
+W_Pxy = tf(5);
+W_Pz = tf(1);
 
-W_Pxy = 5 * bw_xy^2 / (s^2 + 2 * .9 * bw_xy * s + bw_xy^2);
-W_Pz = bw_z^2 / (s^2  + 2 * 1 * bw_z * s + bw_z^2);
+W_Pxydot = tf(.01);
+W_Pzdot = tf(.01);
 
-W_Pxydot = tf(.1); % .2 / (s / bw_xydot + 1);
-W_Pzdot = tf(.1); % .5 / (s / bw_zdot + 1);
- 
-W_Pphitheta = .01 / (s / bw_phitheta + 1);
-W_Ppsi = tf(.1); % .1 / (s / bw_psi + 1);
+W_Pphitheta = .001 / (s * T_xy + 1);
+W_Ppsi = tf(.1);
 
 % Construct performance vector by combining xy and z
+W_ref = blkdiag(W_Rxy * eye(2), W_Rz);
+
 W_PP = blkdiag(W_Pxy * eye(2), W_Pz);
 W_PPdot = blkdiag(W_Pxydot * eye(2), W_Pzdot);
 W_PTheta = blkdiag(W_Pphitheta * eye(2), W_Ppsi);
+
 perf = struct(...
   'FlapAngle', W_Palpha * eye(4), ...
   'Thrust', W_Pomega, ...
+  'ReferenceFilter', W_ref, ...
   'Position', W_PP, ...
   'Velocity', W_PPdot, ...
-  'Angles', W_PTheta);
+  'Angles', W_PTheta ...
+);
 
 if do_plots
   % Bode plots of performance requirements
@@ -61,6 +63,8 @@ if do_plots
 
   bodemag(W_Palpha);
   bodemag(W_Pomega);
+  bodemag(W_Rxy);
+  bodemag(W_Rz);
   bodemag(W_Pxy);
   bodemag(W_Pz);
   bodemag(W_Pxydot);
@@ -70,6 +74,7 @@ if do_plots
 
   grid on;
   legend('$W_{P,\alpha}$', '$W_{P,\omega}$', ...
+    '$W_{R,xy}$', '$W_{R,z}$', ...
     '$W_{P,xy}$', '$W_{P,z}$', ...
     '$W_{P,\dot{x}\dot{y}}$', '$W_{P,\dot{z}}$', ...
     '$W_{P,\phi\theta}$', '$W_{P,\psi}$', ...
@@ -80,11 +85,13 @@ if do_plots
   % Step response of position requirements
   figure; hold on;
   step(W_Pxy); step(W_Pz);
+  step(W_Rxy); step(W_Rz);
   step(W_Pxydot); step(W_Pzdot);
   step(W_Palpha);
   step(W_Pomega);
   grid on;
   legend('$W_{P,xy}$', '$W_{P,z}$', ...
+    '$W_{R,xy}$', '$W_{R,z}$', ...
     '$W_{P,\dot{x}\dot{y}}$', '$W_{P,\dot{z}}$', ...
     '$W_{P,\alpha}$', '$W_{P,\omega}$', ...
     'interpreter', 'latex', 'fontSize', 8);
-- 
cgit v1.2.1