1. Estimation of Sensitivity Function
+1. Estimation of Sensitivity Function
1.1. Load Data
+1.1. Load Data
Two scans are performed: @@ -99,7 +99,7 @@ t = linspace(Ts, Ts*length(ol_drx), length(ol_d
-By comparison the frequency content of the crystal orientation errors between mode B and mode C, it is possible to estimate the Sensitivity transfer function (Figure 1). +By comparison the frequency content of the crystal orientation errors between mode B and mode C, it is possible to estimate the Sensitivity transfer function (Figure 1).
win = hanning(ceil(1/Ts));
@@ -113,7 +113,7 @@ By comparison the frequency content of the crystal orientation errors between mo
Figure 1: Estimation of the sensitivity transfer function magnitude
@@ -122,8 +122,8 @@ By comparison the frequency content of the crystal orientation errors between mo1.2. Controller
+1.2. Controller
load('X_tal_cage_PID.mat', 'K'); @@ -131,7 +131,7 @@ By comparison the frequency content of the crystal orientation errors between mo
Figure 2: Bode Plot of the Controller
@@ -139,8 +139,8 @@ By comparison the frequency content of the crystal orientation errors between mo1.3. Test
+1.3. Test
Ts = 5e-3; @@ -170,12 +170,12 @@ t = linspace(Ts, Ts*length(ol_drx), length(ol_d
2. System Identification
+2. System Identification
2.1. Identification
+2.1. Identification
ur = load('FJPUR_step.mat'); @@ -321,7 +321,7 @@ win = hanning(ceil(5*Fs));
Figure 3: Coherence
@@ -329,7 +329,7 @@ win = hanning(ceil(5*Fs)); -
Figure 4: Bode Plot of the DCM dynamics in the frame of the fast jack.
@@ -359,7 +359,7 @@ Kb = eye(3)*(2* -
Figure 5: Loop gain
@@ -367,7 +367,7 @@ Kb = eye(3)*(2* +
Figure 6: Loop gain
@@ -377,7 +377,7 @@ Kb = eye(3)*(2* -
Figure 7: Comparison of sensitivity functions
@@ -396,7 +396,7 @@ Lb = zeros(3, 3, length(f));
Figure 8: Root Locus
@@ -405,8 +405,8 @@ Lb = zeros(3, 3, length(f));2.2. Identification - New
+2.2. Identification - New
ur = load('FJPUR_step_new.mat'); @@ -577,7 +577,7 @@ Kb = eye(3)*(2* -+-
Figure 9: Comparison of sensitivity functions
@@ -596,7 +596,7 @@ Lb = zeros(3, 3, length(f));+-
Figure 10: Root Locus
@@ -605,8 +605,8 @@ Lb = zeros(3, 3, length(f));-2.3. Identification - White noise
++2.3. Identification - White noise
-ur = load('fjpur_white_noise.mat'); @@ -725,14 +725,14 @@ win = hanning(ceil(0.5*Fs));+-
Figure 11: description
+
Figure 12: Bode Plot of the DCM dynamics in the frame of the fast jack.
@@ -762,14 +762,14 @@ Kb = eye(3)*(2* -+-
Figure 13: Loop gain
+
Figure 14: Loop gain
@@ -779,7 +779,7 @@ Kb = eye(3)*(2* -+-
Figure 15: Comparison of sensitivity functions
@@ -798,7 +798,7 @@ Lb = zeros(3, 3, length(f));+-
Figure 16: Root Locus
@@ -807,8 +807,8 @@ Lb = zeros(3, 3, length(f));-2.4. Test
++2.4. Test
-%% Notch @@ -851,8 +851,8 @@ Kb = 0.8*eye(3)*-2.5. New controller - Higher bandwidth
++2.5. New controller - Higher bandwidth
%% Previously used controller @@ -916,7 +916,7 @@ Kb = 0.9*eye(3)* -+-
Figure 17: description
@@ -937,14 +937,14 @@ Lb_new = zeros(3, 3, length(f));+-
Figure 18: n
+
Figure 19: description
@@ -954,8 +954,8 @@ Lb_new = zeros(3, 3, length(f)); --2.6. Added gain
++2.6. Added gain
-%% Notch @@ -1017,21 +1017,21 @@ Kb_gain = 0.9*eye(3)+-
Figure 20: description
+-
Figure 21: description
+-
Figure 22: nyquist plot
@@ -1041,12 +1041,12 @@ Kb_gain = 0.9*eye(3)-3. Noise Budgeting
++3. Noise Budgeting
--3.1. No Displacement
++3.1. No Displacement
@@ -1275,8 +1275,8 @@ data_70_deg.drx = data_70_deg.allValues(:,6)
-3.2. Scans
++3.2. Scans
@@ -1437,8 +1437,8 @@ CPS_70_10_d = flip(-cumtrapz(flip(f), flip(pxx -
-3.3. Noise budgeting - No rotation
++3.3. Noise budgeting - No rotation
-First, we look at the position errors when the bragg axis is not moving @@ -1488,7 +1488,7 @@ win = hanning(ceil(1/Ts));
+-
Figure 23: Amplitude Spectral Density
@@ -1511,7 +1511,7 @@ CPS_dz = cumtrapz(f, pxx_ol_dz);+-
Figure 24: Cumulative Amplitude Spectrum
@@ -1520,17 +1520,17 @@ CPS_dz = cumtrapz(f, pxx_ol_dz);--3.4. Noise budgeting - Bragg rotation
++3.4. Noise budgeting - Bragg rotation
-4. Test Mode C
++4. Test Mode C
--4.1. Mode B and Mode C
++4.1. Mode B and Mode C
-data_B = extractDatData(sprintf("/home/thomas/mnt/data_id21/22Jan/blc13491/id21/test_regul_220119/%s","lut_const_fj_vel_19012022_1450.dat"), ... @@ -1620,12 +1620,12 @@ win = hanning(ceil(1/Ts));-5. Export numerator and denominator
++5. Export numerator and denominator
--5.1. Export
++5.1. Export
-K_order = 10; @@ -1673,8 +1673,8 @@ fclose(fileID);-5.2. Verify
++5.2. Verify
diff --git a/dcm_feedback.org b/dcm-feedback-control.org similarity index 99% rename from dcm_feedback.org rename to dcm-feedback-control.org index bfb5ea7..629e5f6 100644 --- a/dcm_feedback.org +++ b/dcm-feedback-control.org @@ -44,7 +44,7 @@ #+begin_export htmlK_data = importdata('X_tal_cage_PID_20Hz.dat'); @@ -1687,7 +1687,7 @@ K = tf(K_data(1,:), K_data(2, -Created: 2022-02-15 mar. 14:15
+Created: 2022-02-15 mar. 14:18
-This report is also available as a pdf.
+This report is also available as a pdf.
#+end_export diff --git a/index.html b/index.html index 60e50a8..d936273 120000 --- a/index.html +++ b/index.html @@ -1 +1 @@ -dcm_feedback.html \ No newline at end of file +dcm-feedback-control.html \ No newline at end of file