Add .mat files

2021-08-11 00:28:55 +02:00 · 2021-08-11 00:28:55 +02:00 · 1792ed5fb7
commit 1792ed5fb7
parent 3f2d54dc7f
12 changed files with 35 additions and 0 deletions
--- a/matlab/mat/Khac_iff_struts.mat
+++ b/matlab/mat/Khac_iff_struts.mat
--- a/matlab/mat/Khac_iff_struts_jacobian_cok.mat
+++ b/matlab/mat/Khac_iff_struts_jacobian_cok.mat
--- a/matlab/mat/Khac_iff_struts_svd.mat
+++ b/matlab/mat/Khac_iff_struts_svd.mat
--- a/matlab/mat/Kiff_opt.mat
+++ b/matlab/mat/Kiff_opt.mat
--- a/matlab/mat/frf_iff_unbalanced_vib_table_m.mat
+++ b/matlab/mat/frf_iff_unbalanced_vib_table_m.mat
--- a/matlab/mat/frf_iff_vib_table_m.mat
+++ b/matlab/mat/frf_iff_vib_table_m.mat
--- a/matlab/mat/frf_vib_table_m.mat
+++ b/matlab/mat/frf_vib_table_m.mat
--- a/matlab/mat/identified_plants_enc_plates_masses.mat
+++ b/matlab/mat/identified_plants_enc_plates_masses.mat
--- a/matlab/mat/noise_meas_100s_20kHz.mat
+++ b/matlab/mat/noise_meas_100s_20kHz.mat
--- a/matlab/mat/sim_iff_vib_table_m.mat
+++ b/matlab/mat/sim_iff_vib_table_m.mat
--- a/matlab/mat/sim_vib_table_m.mat
+++ b/matlab/mat/sim_vib_table_m.mat
--- a/matlab/src/generateSpiralAngleTrajectory.m
+++ b/matlab/src/generateSpiralAngleTrajectory.m
@ -0,0 +1,35 @@
+function [ref] = generateSpiralAngleTrajectory(args)
+% generateSpiralAngleTrajectory -
+%
+% Syntax: [ref] = generateSpiralAngleTrajectory(args)
+%
+% Inputs:
+%    - args
+%
+% Outputs:
+%    - ref - Reference Signal
+
+arguments
+    args.R_tot  (1,1) double {mustBeNumeric, mustBePositive} = 10e-6 % [rad]
+    args.n_turn (1,1) double {mustBeInteger, mustBePositive} = 5 % [-]
+    args.Ts     (1,1) double {mustBeNumeric, mustBePositive} = 1e-3 % [s]
+    args.t_turn (1,1) double {mustBeNumeric, mustBePositive} = 1 % [s]
+    args.t_end  (1,1) double {mustBeNumeric, mustBePositive} = 1 % [s]
+end
+
+time_s = 0:args.Ts:args.n_turn*args.t_turn;
+time_e = 0:args.Ts:args.t_end;
+
+Rx = sin(2*pi*time_s/args.t_turn).*(args.R_tot*time_s/(args.n_turn*args.t_turn));
+Ry = cos(2*pi*time_s/args.t_turn).*(args.R_tot*time_s/(args.n_turn*args.t_turn));
+
+Rx = [Rx, 0*time_e];
+Ry = [Ry, Ry(end) - Ry(end)*time_e/args.t_end];
+
+t = 0:args.Ts:args.Ts*(length(Rx) - 1);
+
+ref = zeros(length(Rx), 7);
+
+ref(:, 1) = t;
+ref(:, 5) = Rx;
+ref(:, 6) = Ry;