Skip to content
Snippets Groups Projects
Commit ff758e48 authored by Fredrik Bagge Carlson's avatar Fredrik Bagge Carlson
Browse files

small tweaks

parent 9f416c0e
Branches
Tags
1 merge request!3Julia1
Pipeline #668 failed
......@@ -73,10 +73,10 @@ end
y,u = prbs_experiment(duration = 10, amplitude = 1)
plot([u y], lab=["u" "y"])
bopl(G123, lab="FRA", legend=:left)
na = 4 # Order of A polynomial
nb = 2 # Order of B polynomial
arxtf, Σ = arx(h, y, u, na, nb) # Estimate transfer function with ARX method
bopl(G123, lab="Measured transfer function")
bodeconfidence!(arxtf, Σ, ω = exp10.(LinRange(0,3,200)), linecolor=:red) # The exclamation mark (!)
bodeconfidence!(arxtf, Σ, ω = exp10.(LinRange(0,3,200)), lab="ARX na=$(na)") # The exclamation mark (!)
# adds to the previous plot. If there is no plot open, remove the !
......@@ -160,7 +160,7 @@ bodeconfidence
ribbon := (lowermag, uppermag)
yscale --> :log10
xscale --> :log10
alpha --> 0.3
fillalpha --> 0.2
ω, mag
end
@series begin
......@@ -169,7 +169,7 @@ bodeconfidence
ylabel --> "Phase [deg]"
xlabel --> "Frequency [rad/s]"
xscale --> :log10
alpha --> 0.3
fillalpha --> 0.2
ω, phase
end
nothing
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment