MATLAB/simulink MATLABで磁場を計算する その１
磁場に関するクーロンの法則は、磁荷\(m_1,m_2\)、距離\(r\)とすると
\begin{align}F=k\frac{m}{r^2} \mathrm{}\end{align}
ここで\(k\)は単位に関する定数...
MATLAB/simulink
MATLAB/simulink MATLAB/simulink MATLABで磁場に関するクーロンの法則を計算する その１
磁場に関するクーロンの法則は、磁荷\(m_1,m_2\)、距離\(r\)とすると
\begin{align}F=k\frac{m_1m_2}{r^2} \mathrm{}\end{align}
ここで\(k\)は単位に...
MATLAB/simulink MATLABでフィボナッチ数列を計算する
詳しくは→【C言語】フィボナッチ数列を計算する
MATLABでフィボナッチ数列を計算する。以下コード。
N=10;
F=zeros(1,N+2);
n=0;
F(1,n+1)=0;
n=1;
F(1,n+1)=1;
fo...
MATLAB/simulink Matlabで立方体表面のメッシュを作る
Matlabで立方体表面のメッシュを作る。２枚ずつ作って組み合わせればいい。
N=50;
x1=linspace(-1,1,N);
y1=linspace(-1,1,N);
z1=linspace(-1,1,2);
x2=...
MATLAB/simulink MATLABでLpノルムを計算する
\(L_p\)ノルムは次のように定義される。
\begin{align}|| {\bf x} ||_p = \left ( |x_1|^p + |x_2|^p + \cdots + |x_n|^p \right )^{-p}...
MATLAB/simulink MATLABでヒルベルト行列を作る
MATLABでヒルベルト行列を作る。ヒルベルト行列は各要素が
\begin{align}H_{ij}=\frac{1}{i+j+1}\end{align}
の正方行列である。
clc
clear
N=5;
H=...
MATLAB/simulink MATLABで双一次変換を使った離散PID制御と連続PID制御の応答を確認する
双一次変換を使った離散PID制御と連続PID制御の応答を確認する。双一次変換は連続時間の伝達関数に対して\(s\)を
\begin{align}s=\frac{2(z-1)}{T(z+1)}\end{align}
に置...
MATLAB/simulink MatlabでPID制御のシミュレーションをする
MatlabでPID制御のシミュレーションをする。システムとPID制御器の伝達関数は
\begin{align}P&=\frac{1}{s+1} \\C&=K_P + \frac{K_I}{s} + K_D ...
MATLAB/simulink 【解数】Pythonで指数と対数のグラフを描く
指数と対数のグラフを描く。以下ソースコード
import numpy as np
import matplotlib.pyplot as plt
x = np.arange(-10.0, 10.0, 0.02)
plt.p...
MATLAB/simulink 【制御】加法的不確かさを持つモデル集合のボード線図を書く
加法的不確かさは
\begin{align}\tilde{P}=\{P+\Delta W_a: \|P\|_\infty \leq 1\}\end{align}
で与えられる。この定義に従い、次のような場合のモデル集合...
MATLAB/simulink 【制御】乗法的不確かさを持つモデル集合のボード線図を書く
乗法的不確かさは
\begin{align}\tilde{P}=\{(1+\Delta W_m)P: \|P\|_\infty \leq 1\}\end{align}
で与えられる。この定義に従い、次のような場合のモデ...
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=...
MATLAB/simulink 【MATLAB】ヘロンの公式を実装する
ヘロンの公式は三角形の三辺と面積とを結ぶ関係式で
\begin{align}S=\sqrt{s(s-a)(s-b)(s-c)} \hspace{10mm} s=\frac{a+b+c}{2}\end{align}
MA...
MATLAB/simulink 【MATLAB】２つの座標間の角速度ベクトルの関係をクォータニオンを使って定義する
２つの座標間の角速度ベクトルをクォータニオンを使って表すと
\begin{align}\dot{\boldsymbol{q}} = \boldsymbol{T}_q (\boldsymbol{q}) \boldsymbol{...