python Pythonで伝達関数を部分分数分解する
制御工学ではよく伝達関数の性質を調べるために部分分数分解をすることがある。部分分数分解とは分数の分母を因数分解し、それらをいくつかの分数の和に分解することを指す。例えば
\begin{align}\frac{1}{(x+p_...
古典制御
python 
MATLAB/simulink MatlabでPID制御のシミュレーションをする
MatlabでPID制御のシミュレーションをする。システムとPID制御器の伝達関数は
\begin{align}P&=\frac{1}{s+1} \\C&=K_P + \frac{K_I}{s} + K_D ...
python 【制御】Pythonで単位ステップ関数を描く
Pythonで単位ステップ関数を描画する。単位ステップ関数は
\begin{align}H(x)=\begin{cases}1 \hspace{10mm} (x \geq 0)\\0 \hspace{10mm} (x <...
制御工学 【制御】無限大ノルムの性質２
伝達関数\(G(s)\)について、\(H_{\infty}\)ノルムは
\begin{align}\parallel G(s) \parallel_{\infty} = \sup_\omega \left | G(s) \r...
制御工学 【制御】無限大ノルムの性質１
無限大ノルム
\begin{align}\parallel G(s) \parallel_{\infty} = \sup_\omega \left | G(s) \right |\end{align}
は伝達関数\(G...
python 【制御】互いに逆数の関係にあるシステムのボード線図と性質
2つのシステムが
\begin{align}G_1=\frac{1}{s^2+s+1} \hspace{10mm} G_2= \frac{1}{G_1}\end{align}
のような逆数の関係にある時、それぞれのボー...
python 【制御】直列に接続されたシステムのボード線図と性質
2つのシステム
\begin{align}G_1=\frac{1}{s} \hspace{10mm} G_2=\frac{1}{s^2+s+1}\end{align}
が直列に接続されているとき、全体のボード線図はそれ...
制御工学 【制御】安定多項式の定義
システムを示す次の伝達関数があるとする。
\begin{align}G(s)=\frac{N(s)}{D(s)}=\frac{b_{m} s^{m} + \cdots + b_{1} s + b_{0} }{s^{n} + ...
MATLAB/simulink 【制御】MATLABで離散時関系のH∞ノルムを求める
MATLABならば連続時間のときと同じ。
以下コード
s=tf('s');
sys=c2d(1/(s^2+s+1),1);
norm(sys,Inf)
MATLAB/simulink 【制御】MATLABでH∞ノルムを計算する
線形時不変なシステム
\begin{align}H(s)=\frac{1}{s^2+s+1}\end{align}
のH∞ノルムを求める。すでに用意されている関数を使えばすぐに実装できる。
以下コード
s=...
制御工学 二次遅れシステムの基本形
二次遅れ要素の例として、ばね-質量-ダンパ系の運動方程式は、
\begin{align}f(t)= m \frac{d^2x(t)}{dt^2} + c \frac{dx(t)}{dt} + k x(t)\end{align...
制御工学 一次遅れシステムの基本形
入力信号\(x(t)\)と出力信号\(y(t)\)の間に次の一階微分方程式が成り立つものを一次遅れ要素もしくは一次要素という。
次のような微分方程式を持つシステムは一次遅れ系である。
\begin{align}\tau...
制御工学 伝達関数の定義
次のようなシステムを示す\(n\)階斉次微分方程式
\begin{align}\dfrac{d^n}{dt^n} y(t) &+ a_{n-1} \dfrac{d^{n-1}}{dt^{n-1}}y(t) + \cd...
python Pythonで指定した極、零点、ゲインを実現する伝達関数を求める
これの続き。zpk2tfを使うと指定した極、零点、ゲインを持つ伝達関数を簡単に設計できるようになる。極、零点、ゲインは次のように指定する。
zero = np.array([])
pole = np.array()
G = tf(*...
python Pythonで古典制御と現代制御の双方の視点からばねマスダンパ系を解析する
これの続き。Pythonで同じ解析をした。
from control.matlab import *
import numpy as np
from matplotlib import pyplot as plt
N = ...