【入門】状態空間モデル(DCモータ)Scilab【数値計算】

MATLAB、Python、Scilab、Julia比較ページはこちら
https://www.simulationroom999.com/blog/comparison-of-matlab-python-scilab/

はじめに
DCモータ状態空間モデル
Scilabコード
シミュレーション
まとめ

はじめに

の、
MATLAB,Python,Scilab,Julia比較その70【状態空間モデル㉘】

を書き直したもの。

DCモータ状態空間モデルをScilabでシミュレーションする。

DCモータ状態空間モデル

以下の状態空間モデルをScilabでシミュレーションする。

状態方程式

\(
\begin{bmatrix}
\dot{\theta}(t) \\
\dot{\omega}(t) \\
\dot{I}(t)
\end{bmatrix}=
\begin{bmatrix}
0 && 1 && 0 \\
0 && 0 && K/J \\
0 && -K/L && -R/L
\end{bmatrix}
\begin{bmatrix}
\theta(t) \\
\omega(t) \\
I(t)
\end{bmatrix}+
\begin{bmatrix}
0 \\
0 \\
1/L
\end{bmatrix}
E(t)
\)

出力方程式

\(
\boldsymbol{y}=
\begin{bmatrix}
1 && 0 && 0 \\
0 && 1 && 0 \\
0 && 0 && 1 \\
\end{bmatrix}
\begin{bmatrix}
\theta(t) \\
\omega(t) \\
I(t)
\end{bmatrix}+
\begin{bmatrix}
0 \\
0 \\
0
\end{bmatrix}
E(t)
\)

Scilabコード

シミュレーション用のコードは以下となる。

function [x,y] = statespacemodel(A, B, C, D, u, dt, x)
    % 様態方程式
    x = x + (A*x + B*u) * dt;
    
    % 出力方程式
    y = C*x + D*u;
end

K=0.016;
J=0.000000919;
R=1.34;
L=0.00012;

A=[0,1,0 ; 0,0,K/J ; 0,-K/L,-R/L];
B=[0 ; 0; 1/L];
C=[1,0,0;0,1,0;0,0,1];
D=[0;0;0];

dt = 0.0001;
t = linspace(0, 10, 100000); // 時間(横)軸
u = zeros(1,100000);         // 入力信号生成
u(1,50000:100000)=1;          // 5秒後に0から1へ
y = zeros(3,length(t));
x = zeros(3,1);

for i = 1:length(t)
    [x,y(:,i)] = statespacemodel(A,B,C,D,u(:,i),dt,x);
end

subplot(3,1,1);
plot(t,u);
p=gca();
p.tight_limits(:)="on"
p.data_bounds(:,2)=[-0.05;1.1];
p.data_bounds(:,1)=[4.95;5.25];
subplot(3,1,2);
plot(t,y(1:2,:));
p=gca();
p.tight_limits(:)="on"
p.data_bounds(:,2)=[-1;80];
p.data_bounds(:,1)=[4.95;5.25];
subplot(3,1,3);
plot(t,y(3,:));
p=gca();
p.tight_limits(:)="on"
p.data_bounds(:,1)=[4.95;5.25];

演算部分に関してはMATLABと同一。
グラフ表示の微調整に差がある。