Large update

matlab/src/initializeDisturbances.m

@@ -19,6 +19,10 @@
       args.Fty_x logical {mustBeNumericOrLogical} = true
       % Translation Stage - Z direction
       args.Fty_z logical {mustBeNumericOrLogical} = true
+      % Spindle - X direction
+      args.Frz_x logical {mustBeNumericOrLogical} = true
+      % Spindle - Y direction
+      args.Frz_y logical {mustBeNumericOrLogical} = true
       % Spindle - Z direction
       args.Frz_z logical {mustBeNumericOrLogical} = true
   end
@@ -106,6 +110,38 @@
     Fty_z = zeros(length(t), 1);
   end
 
+  % if args.Frz_x && args.enable
+  %   phi = dist_f.psd_rz;
+  %   C = zeros(N/2,1);
+  %   for i = 1:N/2
+  %     C(i) = sqrt(phi(i)*df);
+  %   end
+  %   rng(131);
+  %   theta = 2*pi*rand(N/2,1); % Generate random phase [rad]
+  %   Cx = [0 ; C.*complex(cos(theta),sin(theta))];
+  %   Cx = [Cx; flipud(conj(Cx(2:end)))];;
+  %   u = N/sqrt(2)*ifft(Cx); % Disturbance Force Rz z [N]
+  %   Frz_x = u;
+  % else
+    Frz_x = zeros(length(t), 1);
+  % end
+
+  % if args.Frz_y && args.enable
+  %   phi = dist_f.psd_rz;
+  %   C = zeros(N/2,1);
+  %   for i = 1:N/2
+  %     C(i) = sqrt(phi(i)*df);
+  %   end
+  %   rng(131);
+  %   theta = 2*pi*rand(N/2,1); % Generate random phase [rad]
+  %   Cx = [0 ; C.*complex(cos(theta),sin(theta))];
+  %   Cx = [Cx; flipud(conj(Cx(2:end)))];;
+  %   u = N/sqrt(2)*ifft(Cx); % Disturbance Force Rz z [N]
+  %   Frz_z = u;
+  % else
+    Frz_y = zeros(length(t), 1);
+  % end
+
   if args.Frz_z && args.enable
     phi = dist_f.psd_rz;
     C = zeros(N/2,1);
@@ -134,14 +170,14 @@
 
 if exist('./mat', 'dir')
     if exist('./mat/nass_disturbances.mat', 'file')
-        save('mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_z', 'Fd', 'Ts', 't', 'args', '-append');
+        save('mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_x', 'Frz_y', 'Frz_z', 'Fd', 'Ts', 't', 'args', '-append');
     else
-        save('mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_z', 'Fd', 'Ts', 't', 'args');
+        save('mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_x', 'Frz_y', 'Frz_z', 'Fd', 'Ts', 't', 'args');
     end
 elseif exist('./matlab', 'dir')
     if exist('./matlab/mat/nass_disturbances.mat', 'file')
-        save('matlab/mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_z', 'Fd', 'Ts', 't', 'args', '-append');
+        save('matlab/mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_x', 'Frz_y', 'Frz_z', 'Fd', 'Ts', 't', 'args', '-append');
     else
-        save('matlab/mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_z', 'Fd', 'Ts', 't', 'args');
+        save('matlab/mat/nass_disturbances.mat', 'Dwx', 'Dwy', 'Dwz', 'Fty_x', 'Fty_z', 'Frz_x', 'Frz_y', 'Frz_z', 'Fd', 'Ts', 't', 'args');
     end
 end