update demo.jl

Author: ufechner7

examples/demo.jl

@@ -1,5 +1,5 @@
 using Pkg
-if ! haskey(Pkg.installed(), "ControlPlots")
+if ! ("ControlPlots" ∈ keys(Pkg.project().dependencies))
     using TestEnv; TestEnv.activate()
 end
 using KitePodModels, ControlPlots
@@ -9,26 +9,22 @@ const dt = 1.0 / se().sample_freq
 
 kcu = KCU(se())
 
-# rel_depower_offset = se().depower_offset/100.0
-# for step in 1:Int(round(t_end/dt))
-#     time = step * dt
-#     if step < 20
-#         set_depower_steering(kcu, rel_depower_offset, 0.0)
-#         push!(depower_set, rel_depower_offset)
-#         push!(steering_set, 0.0)
-#     else
-#         set_depower_steering(kcu, rel_depower_offset+0.5, 0.5)
-#         push!(depower_set, rel_depower_offset+0.5)
-#         push!(steering_set, 0.5)
-#     end
-#     push!(times, time)
-#     KitePodModels.on_timer(kcu)
-#     push!(depower, get_depower(kcu))
-#     push!(steering, get_steering(kcu))
-# end
-
 rel_depower_min = 0.22
 rel_depower_max = 0.60
 α_min = calc_alpha_depower(kcu, rel_depower_min)
 α_max = calc_alpha_depower(kcu, rel_depower_max)
-println("α_min = $(rad2deg(α_min)), α_max = $(rad2deg(α_max))")
+# Calculate the change of the angle between the kite and the last tether segment [rad] 
+# as function of the actual rel_depower value. 
+# println("α_min = $(rad2deg(α_min)), α_max = $(rad2deg(α_max))")
+
+# plot alpha_depower as function of rel_depower
+function plot_alpha_depower(kcu; rel_depower_min=0.22, rel_depower_max=0.60)
+    rel_depower = range(rel_depower_min, rel_depower_max, length=100)
+    α = zeros(length(rel_depower))
+    for i in eachindex(α)
+        α[i] = calc_alpha_depower(kcu, rel_depower[i])
+    end
+    plot(rel_depower, rad2deg.(α); ylabel="α [°]", xlabel="rel_depower", fig="α = f(rel_depower)")  
+end 
+
+plot_alpha_depower(kcu)