Школа системного моделирования

Идентификация

Загрузка данных из таблицы

using DataFrames, CSV;

Import_data = "/user/data.csv" # @param {type:"filepicker"}
table = CSV.read(Import_data, DataFrame);
print(".csv файл импортирован")
Показать_таблицу = true # @param {type:"boolean"}
if Показать_таблицу
    display(table)
end

.csv файл импортирован

Столбец_времени = 1 # @param {type:"integer"}
Столбец_входа = 2 # @param {type:"integer"} 
Столбец_выхода = 3 # @param {type:"integer"}
Начинать_со_строки = 1 # @param {type:"integer"} 
exp_time_data = table[Начинать_со_строки:end, Столбец_времени]; 
exp_in_data = table[Начинать_со_строки:end, Столбец_входа];
exp_out_data = table[Начинать_со_строки:end, Столбец_выхода];

Графики измеренных сигналов:

input = plot(exp_time_data, exp_in_data, title = "Вход", xlabel="t, c",legend=false)
output_measure = plot(exp_time_data, exp_out_data, title = "Выход", xlabel="t, c", legend=false)
plot(output_measure,input,layout = (2, 1))

Идентификация

Pkg.add(["ControlSystemIdentification", "ControlSystems"])

using ControlSystems
using ControlSystemIdentification;

[ Info: Precompiling RecursiveArrayToolsMonteCarloMeasurementsExt [53dc5841-8f62-55e1-8f15-93d074dd5209] (cache misses: wrong julia version (2))
GKS: cannot open display - headless operation mode active
[ Info: Precompiling ArrayInterfaceChainRulesExt [8673388c-2cec-59d7-825e-32163c06c3e1] (cache misses: wrong julia version (2))
GKS: cannot open display - headless operation mode active
WARNING: using ControlSystemIdentification.input in module Main conflicts with an existing identifier.
┌ Error: handle_connection handler error. 
│ 
│ ===========================
│ HTTP Error message:
│ 
│ ERROR: IOError: write: broken pipe (EPIPE)
│ Stacktrace:
│   [1] uv_write(s::Sockets.TCPSocket, p::Ptr{UInt8}, n::UInt64)
│     @ Base ./stream.jl:1100
│   [2] unsafe_write(s::Sockets.TCPSocket, p::Ptr{UInt8}, n::UInt64)
│     @ Base ./stream.jl:1154
│   [3] unsafe_write
│     @ /usr/local/ijulia-core/packages/HTTP/ShTJs/src/Connections.jl:132 [inlined]
│   [4] write
│     @ ./strings/io.jl:248 [inlined]
│   [5] write
│     @ ./io.jl:794 [inlined]
│   [6] unsafe_write(http::HTTP.Streams.Stream{HTTP.Messages.Request, HTTP.Connections.Connection{Sockets.TCPSocket}}, p::Ptr{UInt8}, n::UInt64)
│     @ HTTP.Streams /usr/local/ijulia-core/packages/HTTP/ShTJs/src/Streams.jl:95
│   [7] write
│     @ ./strings/io.jl:248 [inlined]
│   [8] write(io::HTTP.Streams.Stream{HTTP.Messages.Request, HTTP.Connections.Connection{Sockets.TCPSocket}}, s::Base.CodeUnits{UInt8, String})
│     @ Base ./strings/basic.jl:807
│   [9] (::HTTP.Handlers.var"#1#2"{HTTP.Handlers.Router{typeof(HTTP.Handlers.default404), typeof(HTTP.Handlers.default405), Nothing}})(stream::HTTP.Streams.Stream{HTTP.Messages.Request, HTTP.Connections.Connection{Sockets.TCPSocket}})
│     @ HTTP.Handlers /usr/local/ijulia-core/packages/HTTP/ShTJs/src/Handlers.jl:61
│  [10] #invokelatest#2
│     @ ./essentials.jl:1055 [inlined]
│  [11] invokelatest
│     @ ./essentials.jl:1052 [inlined]
│  [12] handle_connection(f::Function, c::HTTP.Connections.Connection{Sockets.TCPSocket}, listener::HTTP.Servers.Listener{Nothing, Sockets.TCPServer}, readtimeout::Int64, access_log::Nothing, verbose::Bool)
│     @ HTTP.Servers /usr/local/ijulia-core/packages/HTTP/ShTJs/src/Servers.jl:469
│  [13] (::HTTP.Servers.var"#16#17"{HTTP.Handlers.var"#1#2"{HTTP.Handlers.Router{typeof(HTTP.Handlers.default404), typeof(HTTP.Handlers.default405), Nothing}}, HTTP.Servers.Listener{Nothing, Sockets.TCPServer}, Set{HTTP.Connections.Connection}, Int64, Nothing, ReentrantLock, Bool, Base.Semaphore, HTTP.Connections.Connection{Sockets.TCPSocket}})()
│     @ HTTP.Servers /usr/local/ijulia-core/packages/HTTP/ShTJs/src/Servers.jl:401
│   request =
│    HTTP.Messages.Request:
│    """
│    GET /api/v1/workspace/variables?field=name&order=ASC HTTP/1.1
│    Host: engee.com
│    X-Real-IP: 10.244.7.249
│    X-Forwarded-For: 10.129.0.16, 10.244.7.249
│    X-Forwarded-Proto: http
│    Connection: close
│    X-Request-ID: 3478b009c0d8ba2c36055ef809ff7c03
│    X-Forwarded-Host: engee.com
│    X-Forwarded-Port: 443
│    X-Forwarded-Scheme: https
│    X-Scheme: https
│    sec-ch-ua-platform: "Windows"
│    authorization: ******
│    engee-blocks-locale: true
│    sec-ch-ua: "Chromium";v="142", "YaBrowser";v="25.12", "Not_A Brand";v="99", "Yowser";v="2.5"
│    sec-ch-ua-mobile: ?0
│    user-agent: Mozilla/5.0 (Windows NT 10.0; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/142.0.0.0 YaBrowser/25.12.0.0 Safari/537.36
│    accept: application/json, application/x-msgpack
│    dnt: 1
│    engee-locale: ru
│    sec-fetch-site: same-origin
│    sec-fetch-mode: cors
│    sec-fetch-dest: empty
│    referer: https://engee.com/prod/user/demo54365638-svetlana/
│    accept-encoding: gzip, deflate, br, zstd
│    accept-language: ru,en;q=0.9,es;q=0.8,bg;q=0.7
│    priority: u=1, i
│    cookie: ******
│    
│    """
└ @ HTTP.Servers /usr/local/ijulia-core/packages/HTTP/ShTJs/src/Servers.jl:483

Ts = 0.1;
identdata  = iddata(exp_out_data[1:div(end, 2)], exp_in_data[1:div(end, 2)], Ts);   
plot(identdata)

valdata = iddata(exp_out_data[div(end, 2):end], exp_in_data[div(end, 2):end], Ts);   
plot(valdata)

sys = subspaceid(identdata, 2; zeroD=true, verbose=true)

[ Info: Estimating with r = 14, s1 = 14, s2 = 14, stability constraint = true, zeroD = true, focus = prediction, scaleU = true
[ Info: Fraction of variance explained: 0.9991335120383276

N4SIDStateSpace{StateSpace{Discrete{Float64}, Float64}, LinearAlgebra.Hermitian{Float64, Matrix{Float64}}, LinearAlgebra.Hermitian{Float64, Matrix{Float64}}, Matrix{Float64}, Matrix{Float64}, Matrix{Float64}, Vector{Float64}, LinearAlgebra.SVD{Float64, Float64, Matrix{Float64}, Vector{Float64}}, Float64}
A = 
 0.6618554901243604    0.08213944438423898
 0.008608619552158294  0.9940841776152444
B = 
 -0.0010922620723805732
  2.8219041220266243e-5
C = 
 -0.0606243526845628  0.005906283927498871
D = 
 0.0

Sample Time: 0.1 (seconds)
Discrete-time state-space model

simplot(identdata, sys)

simplot(valdata, sys)

Преобразование дискретной модели в непрерывную (discrete to continuous)

objct=d2c(sys.sys)

StateSpace{Continuous, Float64}
A = 
 -4.133182023925924     1.0060586787716823
  0.10543991961010012  -0.06398534390475205
B = 
 -0.013351075535161832
  0.00034471965774845683
C = 
 -0.0606243526845628  0.005906283927498871
D = 
 0.0

Continuous-time state-space model

Модуль 4. Лекция 12 - Школа системного моделирования

Школа системного моделирования

Идентификация

Загрузка данных из таблицы

Идентификация

Тип

Краткое описание

Связанные материалы

Теги

Категории

Языки

Форматы

Уровень

Статус

Опубликовано

Обновлено

Источник

Row	time	V	w(t)
	Float64	Float64	Float64
1	0.0	0.0	0.0
2	0.01	0.0	0.0
3	0.02	0.0	0.0
4	0.03	0.0	0.0
5	0.04	0.0	0.0
6	0.05	0.0	0.0
7	0.06	0.0	0.0
8	0.07	0.0	0.0
9	0.08	0.0	0.0
10	0.09	0.0	0.0
11	0.1	0.0	0.0
12	0.11	0.0	0.0
13	0.12	0.0	0.0
⋮	⋮	⋮	⋮
1990	19.89	0.0	1.05514e-19
1991	19.9	0.0	6.95588e-20
1992	19.91	0.0	4.58556e-20
1993	19.92	0.0	3.02296e-20
1994	19.93	0.0	1.99285e-20
1995	19.94	0.0	1.31375e-20
1996	19.95	0.0	8.66073e-21
1997	19.96	0.0	5.70946e-21
1998	19.97	0.0	3.76388e-21
1999	19.98	0.0	2.48128e-21
2000	19.99	0.0	1.63575e-21
2001	20.0	0.0	1.07834e-21