MATLAB/simulink simulinkで状態方程式をシミュレーションする 前回の記事の続き 前回のパラメータを使ってsimulinkでシミュレーションした。モデルはこれ。 結果は同じになる。 2021.07.26 MATLAB/simulinkプログラミング制御工学古典制御現代制御
古典制御 微分要素と積分要素のボード線図 微分要素と積分要素の伝達関数は \begin{align}G_1(s)=\frac{1}{s} \hspace{5mm} G_2 (s) = s\end{align} で表される。\(s=j \omega\)を代入すれば \begin{al... 2021.07.24 古典制御
MATLAB/simulink 古典制御と現代制御の双方の視点からばねマスダンパ系を解析する ばねマスダンパ系について、古典制御と現代制御の各変数の相互関係とそれらの解析結果が一致することを示す。 ばねマスダンパ系の運動方程式は \begin{align}M \ddot{y}(t) + C \dot{y} (t) + ky (t) ... 2021.07.22 MATLAB/simulinkプログラミング制御工学古典制御現代制御
MATLAB/simulink クラーク変換とパーク変換をとりあえず試す 細かいことは気にせずクラーク変換を実装する。クラーク変換は \begin{align}\begin{bmatrix}\alpha \\ \beta \\ \gamma\end{bmatrix}=k \begin{bmatrix}1 & -\... 2021.07.20 MATLAB/simulinkプログラミング電気機器
MATLAB/simulink 線形時不変なシステムのH∞ノルムを計算する 線形時不変なシステムのH∞ノルムを計算する。適当な伝達関数 \begin{align}G(s)=\frac{1}{(s+1)(s+2)}\end{align} を定義する。H∞ノルムの定義は \begin{align}\| G(j \ome... 2021.07.20 MATLAB/simulinkプログラミング制御工学古典制御現代制御
MATLAB/simulink 線形時不変なシステムの可制御行列と可観測行列を調べる 線形時不変なシステム(A,B,C,D)の可制御性と可観測性を調べる。以前証明した結果から 可制御性行列 \begin{align}M_c=\begin{bmatrix}B & AB & \cdots & A^{n-1}B\end{bmatr... 2021.07.20 MATLAB/simulinkプログラミング制御工学現代制御
MATLAB/simulink MALTABを使ってボード線図を書く MATLABでボード線図を描くとき bode(sys) とすれば即座に得られるが、今回は別の方法で描画する。 ボード線図はゲイン線図と位相線図からなり、その定義は \begin{align}gain &: 20 \log_{10} |G(j... 2021.07.19 MATLAB/simulinkプログラミング制御工学古典制御
制御工学 関数とノルム 信号などの大きさを図る手法としてノルムが定義されている。普段よく使う絶対値 \begin{align}|x|= \begin{cases}x \hspace{5mm} &( x \geq 0) \\-x &( x < 0)\end{case... 2021.07.19 制御工学
MATLAB/simulink MATLABでもっとも単純な数値積分と数値微分を実装する 伝達関数 \begin{align}P=\frac{1}{s^2+1}\end{align} の単純な数値積分と数値微分を考える。 始めに確認のため、連続時間での結果を確認する。連続時間での微分と積分はラプラス演算子を用いて \begin{... 2021.07.18 MATLAB/simulinkプログラミング制御工学古典制御
MATLAB/simulink ボード線図の示す意味と応答への影響 伝達関数とボード線図の関係を考察する。今伝達関数が \begin{align}G(s)=\frac{s^2+1}{s^2+6s+8}\end{align} で与えられているとする。この時の分子多項式に\(s=-j \omega\)を代入した... 2021.07.16 MATLAB/simulinkプログラミング制御工学古典制御
MATLAB/simulink 離散時間モデルのボード線図とstep応答 MATLABを使って離散時間モデルのボード線図の描画とstep応答を調べる。MATLABのbode関数は離散時間モデルにも対応しているので簡単。 ボード線図 step応答 ソース T = 0.1; s=tf("s"); Gs = 2/(3*... 2021.07.15 MATLAB/simulinkプログラミング制御工学古典制御
MATLAB/simulink MATLABで離散モデルに対しPID制御器を作りシミュレーションする 前回に続いてMATLABで離散モデルに対しPID制御器を構築する。与えられたモデルを離散化するには G=c2d(omega^2/(s^2+2*zeta*omega*s+omega^2),Ts); を使えばいい。例では二次遅れ系をサンプル時間... 2021.07.13 MATLAB/simulinkプログラミング制御工学古典制御
MATLAB/simulink MATLABでPID制御器を作りシミュレーションする MATLABでPID制御器を作って実装する。MATLABにはモデルを与えればそれっぽくPID制御器を構築してくれる便利な関数pidtuneがあるので今回はこれを使う。pidtuneの構文は下記の通り。 = pidtune(G,'PID');... 2021.07.13 MATLAB/simulinkプログラミング制御工学古典制御
制御工学 現代制御の視点から見た線形時不変な多入力多出力なシステムのPID制御 線形時不変な多入力多出力なシステム \begin{align}\dot{x}(t)&=Ax(t)+Bu(t) \\y(t)&=Cx(t)\end{align} をPID制御則は \begin{align}u(t)=K_{P} (r-y(t)... 2021.07.12 制御工学現代制御
制御工学 線形システムに対するリアプノフの定理の必要十分条件 線形システムの安定性を考える。このとき \begin{align}V(x)=x^{T} Px\end{align} と考えればその微分は \begin{align}\dot{V} (x)=\dot{x}^{T} P x + x^{T} P ... 2021.07.11 制御工学現代制御
制御工学 状態方程式の同値変換 次の状態方程式で表現されるシステム \begin{align}\dot{x}(t) &=Ax(t)+Bu(t)\\y(t)&=Cx(t)\end{align} において任意の正則な行列\(T¥)を使って \begin{align}\over... 2021.07.10 制御工学現代制御
C/C++/C# C言語での多次元配列の渡し方 C言語で配列を宣言する場合 int mat; などとする。このように定義された配列を関数に渡すときは先頭のアドレスを渡せばよく、実際の関数の宣言は void func(int *mat) となる。 一方で多次元配列のこれが多次元になると話が... 2021.07.09 C/C++/C#プログラミング
制御工学 線形時不変システムの可制御性の十分条件 線形時不変な状態方程式 \begin{align}\dot{x}(t) &= Ax(t)+Bu(t) \\y(t)&=Cx(t)+Du(t)\end{align} について \begin{align}\mathrm{rank} M_c = ... 2021.07.06 制御工学現代制御
制御工学 線形時不変システムの可観測性 線形時不変な状態方程式 \begin{align}\dot{x}(t) &= Ax(t)+Bu(t) \\y(t)&=Cx(t)+Du(t)\end{align} が不可観測であるとする。不可観測性から求められる異なった二つの初期値を\(x... 2021.07.05 制御工学現代制御
制御工学 直達行列を含む線形時不変な状態方程式の出力を求める 線形時不変な状態方程式 \begin{align}\dot{x}(t) &= Ax(t)+Bu(t) \\y(t)&=Cx(t)+Du(t)\end{align} の初期値を\(x(0)\)とすればその解は \begin{align}x(t... 2021.07.05 制御工学現代制御