python シグモイド関数を描画する
シグモイド関数の一つをPythonを使って描画する 関数は \begin{align}y=\frac{1}{1+e^{-x}}\end{align} import numpy as np from matplotlib import pyp...
プログラミング
python 
python ルンゲクッタ法でRC回路の出力を求める
電流の関係式 \begin{align}i=\frac{dq}{dt}= C \frac{dv}{dt} \end{align} より出力電圧は \begin{align}\frac{dv}{dt} = \frac{E-v}{CR} \en...
python Pythonでテント写像を描画する
テント写像は \begin{align}x_{n+1} = 1-\left | 1 - 2 x_{n} \right |\end{align} import numpy as np from matplotlib import pyplot...
python Pythonで指定した極、零点、ゲインを実現する伝達関数を求める
これの続き。zpk2tfを使うと指定した極、零点、ゲインを持つ伝達関数を簡単に設計できるようになる。極、零点、ゲインは次のように指定する。 zero = np.array([]) pole = np.array() G = tf(*zpk2...
python Pythonで古典制御と現代制御の双方の視点からばねマスダンパ系を解析する
これの続き。Pythonで同じ解析をした。 from control.matlab import * import numpy as np from matplotlib import pyplot as plt N = 1000 t = ...
python Pythonでドップラー効果を計算する
ドップラー効果は観測者の速度を \( v_{0}\)、 音源の速度を\( v_{s}\) 、音源の周波数を\( f_{0}\) とすると \begin{align}f=f_{0} \frac{v-v_{0}}{v-v_{s}}\end{al...
python Pythonでナイキスト線図を書く
これの続き。Pythonでナイキスト線図を書いてみた。 from control.matlab import * import numpy as np from matplotlib import pyplot as plt N=100 t...
python Pythonで忘却係数付き逐次最小二乗法を実装する
これの続き。忘却係数付き逐次最小二乗法の更新則は \begin{align}\hat{\theta}_{N} &= \hat{\theta}_{N-1} + \dfrac{P_{N-1} z_{N} }{\rho + z_{N}^{T} P...
python Pythonでボード線図を書く
Pythonでボード線図を書くにはbode関数を使えばいい。 from control.matlab import * from matplotlib import pyplot as plt s = tf('s') num= den= s...
python Pythonでシステムの極を調べる
Pythonでシステムの極を調べるには syspole = pole(sys) を実行すればいい。以下コード from control.matlab import * s = tf('s') zeta = 1 omega = 2 P = o...
python PythonでPID制御をシミュレーションする
Pythonでフィードバック結合のシミュレーションをする。 \begin{align}P(s)=\frac{\omega_{n}^2 }{s^2 + 2 \zeta \omega_{n} s + \omega_{n}^2}\end{alig...
python Pythonでフィードバック結合をシミュレーションする
Pythonでフィードバック結合のシミュレーションをする。 \begin{align}P(s)=\frac{\omega_{n}^2 }{s^2 + 2 \zeta \omega_{n} s + \omega_{n}^2}\end{alig...
python Pythonで二次遅れ系のシミュレーションをする
Pythonで二次遅れ系をシミュレーションをする。 \begin{align}G(s)=\frac{\omega_{n}^2 }{s^2 + 2 \zeta \omega_{n} s + \omega_{n}^2}\end{align} コ...
python Pythonで描いたグラフにgridを追加する
前回の結果のままでは見にくいのでgridを追加する。gridは plt.grid(color='k', linestyle='dotted', linewidth=1) で追加できる。引数を指定すればいろいろ変更できる。 import nu...
python Pythonで三角関数のグラフを書く
pythonでグラフを書くにはmatplotlibを使えばいい。 今回のコードを試すには三角関数も使うのでnumpyも必要になる。 import numpy as np import matplotlib.pyplot as plt x =...
python PythonでFizz Buzz
これの続き。 N = 200 for i in range(1, N + 1): if i % 3 == 0 and i % 5 == 0: print("FizzBuzz") elif i % 3 == 0: print("Fizz") ...
MATLAB/simulink MATLABの互換にはGNU Octaveがいいらしい
GNU Octaveを試してみる。フリーなので今ある記事も置き換える予定。
MATLAB/simulink 逐次最小二乗法を使って二次遅れ系の伝達関数を推定する その２
前回の記事の続き 前回と同様に逐次最小二乗法を用いて二次遅れ系の伝達関数 \begin{align}G(s)=\frac{\omega_{n}^2 }{s^2 + 2 \zeta \omega_{n} s + \omega_{n}^2}\e...
MATLAB/simulink 逐次最小二乗法を使って二次遅れ系の伝達関数を推定する
逐次最小二乗法については前回の記事を参照。 逐次最小二乗法を用いて二次遅れ系の伝達関数 \begin{align}G(s)=\frac{\omega_{n}^2 }{s^2 + 2 \zeta \omega_{n} s + \omega_{...
C/C++/C# 二度漬けがバレると店のオヤジから怒られるプログラム
乱数を使って二度漬けして遊ぶプログラムを作りました。色々変えて遊んでみてください。 #include <stdio.h> #include <stdlib.h> #include <time.h> char* scanstr(void) {...