しろねこらぼ

RLCバンドパスフィルタの伝達関数

回路方程式が \begin{align}E _{I} =L\frac{di}{dt} + \frac{1}{C} \int i dt + Ri\end{align} のようなバンドパスフィルタを考える。 いま出力が \b […]

2021/10/29 07:19

承認と許可の違い

電気工事士法第43条には自家用電気工作物を設置するものが「許可」を受けて主任技術者の免状を受けてないものを主任技術者に選任することができることが規定されている。 電気法規にはたびたび「許可」と「承認」が出てくるので、法律 […]

2021/10/27 11:54

C言語で漸化式を解く

漸化式とは前回計算した値を使い今の値を計算するような式である。今回は簡単な漸化式 \begin{align}a_n=a_{n-1}+b\end{align} の形をした漸化式をC言語で計算する。 ソース \(N\)を変更 […]

2021/10/26 16:06

零時ホールドのボード線図

零時ホールドの伝達関数は \begin{align}H(s)=\frac{1-e^{-sT}}{s}\end{align} ボード線図を書くと次のようになる。 \(T\)を変化させるほどにゲインが下がる。これはサンプル＆ […]

2021/10/24 09:55

「制御基礎理論―古典から現代まで」を買いました

中野先生の本を買いました。 目次 自動制御 信号の伝達と伝達関数 ブロック線図の構成要素 ブロック線図の等価変換 微分・積分要素のブロック線図 等価変換の応用 シグナルフロー線図 ラプラス変換と自動制御 ラプラス変換とラ […]

2021/10/24 01:37

MATLABを使って零次ホールドを試す

matlabを使って零次ホールドを試す。対象のシステムは \begin{align}G=\frac{s}{1+s}\end{align} 零次ホールドは \begin{align}H=\dfrac{1-e^{-sT}}{ […]

2021/10/23 22:41

システムの離散化と零次ホールド

D/Aなどにより現在の出力が次の出力に移るまでの出力は一定値に保持されるのがふつうである。これを零次ホールドという。 零次ホールドの伝達関数は\(u(t)-u(t-1)\)に対応するようにすればよいので begin{al […]

2021/10/22 10:13

合格発表日

今日は合格発表日 例年だと9時半頃発表されます。

2021/10/22 08:10

MATLABで楕円を書く

楕円の座標は \begin{align}x&=a \cos \theta\y&=b \sin \theta\end{align} で計算できる。\(a=b\)の時、円になる。

2021/10/22 07:13

単層パーセプトロンをclassを使って定義しORを実装する

これの続き。 単層パーセプトロンをclass定義して、動作確認にORを計算する。 mylib これで単層パーセプトロンをいくらでも複製できるようになった。

2021/10/20 03:27

単層パーセプトロンを使ってORを作る

単層パーセプトロンを使ってORを作る。ORは A B X 0 0 0 0 1 1 1 0 1 1 1 1

2021/10/18 04:07

双一次変換を使って一次遅れ系の伝達関数を離散化する

一次遅れ系の伝達関数 \begin{align}G(s)=\frac{K}{Ts+1}\end{align} を双一次変換で離散化する。\(s\)に \begin{align}s=\frac{2(1-z^{-1})}{ […]

2021/10/12 02:33

MATLABで忘却係数付き逐次最小二乗法を実装する

忘却係数付き逐次最小二乗法の更新則は \begin{align}\hat{\theta}_{N} &= \hat{\theta}_{N-1} + \dfrac{P_{N-1} z_{N} }{\rho + z_{ […]

2021/10/11 01:05

忘却係数付き逐次最小二乗法

忘却係数付き逐次最小二乗法とは次のようなものである。 ・評価関数 \begin{align}J_{N}= \sum_{i=1}^{N} \rho^{N-i} \left ( y_{i} – z_{i}^{T} […]

2021/10/10 03:12

インボリュート曲線を描く

歯車の設計にはインボリュート曲線が使われる。今回は基礎円とそこから延びるインボリュート曲線を描く。 インボリュート曲線とは円柱に巻き付けられた糸を撓ませることなく引き出した時、糸の先端がたどる軌跡を示す曲線を言う。 媒介 […]

2021/10/09 04:05

簡単な場合の最小二乗法

回帰式を \begin{align}\hat{y}=Ax+B\end{align} とする。残差の二乗和は \begin{align}J = \sum e_i^2 = \sum \{ y_i – (Ax_i […]

2021/10/06 22:44

単層パーセプトロンを使ってANDを作る

単層パーセプトロンを使ってANDを作る。ANDは A B X 0 0 0 0 1 0 1 0 0 1 1 1 となるような演算である。一方で単層パーセプトロンとはそれぞれの入力\(x\)に重み \( w \) を乗じて和 […]

2021/10/06 02:16

オペアンプを使ってPID制御器を作る #1

以前MATLABでPID制御をシミュレーションした。今回はアナログ電子回路でPID制御を実装する。 簡単に作る場合はオペアンプを用いればよく、かなり単純である。 次回以降理論を追いつつ設計する。

2021/10/06 00:13

C言語でFizz Buzz

#include int main() { int i, N = 36; for (i […]

2021/10/04 00:07

MATLAB2021aにおけるfor文の計算速度

MATLABを使うときよく言われるのは「for文は遅い」だと思う。今回は2021aでその速度を調べてみた。 使用したソースコード 結果は for文：0.330623 秒ベクトル：0.000204 秒 最近早くなったとは聞 […]

2021/10/03 00:05