@@ -25,7 +25,7 @@ winch = Winch(wcs, set)
2525# find equilibrium speed
2626function find_equilibrium_speed (winch, set_speed, force, n= 10000 )
2727 last_v_act = 0.0
28- for v_set in range (- 2.0 , 2 * set_speed, n)
28+ for v_set in range (- 5 , 2 * set_speed+ 5 , n)
2929 lim_speed = minimum ([v_set, set_speed])
3030 set_force (winch, force)
3131 set_v_set (winch, lim_speed)
9696function margins ()
9797 global sys
9898 margins = Float64[]
99- for v_wind in range (- 1 , 1 , length= 3 )
99+ for v_wind in range (0.2 , 1 , length= 5 )
100100 local v_set, dm
101101 v_set = 0.57 * v_wind
102102 sys = open_loop_system (winch, v_set, v_wind)
@@ -119,11 +119,11 @@ function margins()
119119end
120120margins ()
121121
122- for v_wind in range (- 1 , 1 , length= 3 )
122+ for v_wind in range (0.2 , 1 , length= 5 )
123123 global sys
124124 local v_set
125125 v_set = 0.57 * v_wind
126126 sys = open_loop_system (winch, v_set, v_wind)
127- bode_plot (sys; from= 0.76 , to= 2.85 , title= " System with LFC, v_wind=-1 ..1 m/s" )
127+ bode_plot (sys; from= - 0.3 , to= 2.85 , title= " System with LFC, v_wind=0.2 ..1 m/s" )
128128end
129129nothing
0 commit comments